Search

Aims. We model the energetic storm particle (ESP) event of 14 July 2012 using the energetic particle acceleration and transport model named PArticle Radiation Asset Directed at Interplanetary Space Exploration'(PARADISE), together with the solar wind and coronal mass ejection (CME) model named EUropean Heliospheric FORcasting Information Asset'(EUHFORIA). The simulation results illustrate both the capabilities and limitations of the utilised models. We show that the models capture some essential structural features of the ESP event; however, for some aspects the simulations and observations diverge. We describe and, to some extent, assess the sources of errors in the modelling chain of EUHFORIA and PARADISE and discuss how they may be mitigated in the future. Methods. The PARADISE model computes energetic particle distributions in the heliosphere by solving the focused transport equation in a stochastic manner. This is done using a background solar wind configuration generated by the ideal magnetohydrodynamic module of EUHFORIA. The CME generating the ESP event is simulated by using the spheromak model of EUHFORIA, which approximates the CME's flux rope as a linear force-free spheroidal magnetic field. In addition, a tool was developed to trace CME-driven shock waves in the EUHFORIA simulation domain. This tool is used in PARADISE to (i) inject 50 keV protons continuously at the CME-driven shock and (ii) include a foreshock and a sheath region, in which the energetic particle parallel mean free path, λ¥, decreases towards the shock wave. The value of λ¥ at the shock wave is estimated from in situ observations of the ESP event. Results. For energies below ∼ 1 MeV, the simulation results agree well with both the upstream and downstream components of the ESP event observed by the Advanced Composition Explorer. This suggests that these low-energy protons are mainly the result of interplanetary particle acceleration. In the downstream region, the sharp drop in the energetic particle intensities is reproduced at the entry into the following magnetic cloud, illustrating the importance of a magnetised CME model., We acknowledge the use of ACE data from the ACE Science Center and of SOHO/ERNE data from sepem.eu and of the SEPEM Reference Data Set version 2.00, European Space Agency (2016). N.W. acknowledges funding from the Research Foundation – Flanders (FWO – Vlaanderen, fellowship no. 1184319N). A.A. and B.S. acknowledge the support by the Spanish Ministerio de Ciencia e Innovación (MICINN) under grant PID2019-105510GB-C31 and through the “Centre of Excellence María de Maeztu 2020–2023” award to the ICCUB (CEX2019-000918-M). C.S. acknowledges the NASA Living With a Star Jack Eddy Postdoctoral Fellowship Program, administered by UCAR’s Cooperative Programs for the Advancement of Earth System Science (CPAESS) under award no. NNX16AK22G. D.L. acknowledges support from NASA Living With a Star (LWS) programs NNH17ZDA001N-LWS and NNH19ZDA001N-LWS, and the Goddard Space Flight Center Heliophysics Innovation Fund (HIF) program. The work in the University of Turku was performed under the umbrella of Finnish Centre of Excellence in Research of Sustainable Space. This project has received funding from the European Union’s Horizon 2020 research and innovation programs under grant agreement No. 870405 (EUHFORIA 2.0). These results were also obtained in the framework of the ESA project “Heliospheric modelling techniques” (Contract No. 4000133080/20/NL/CRS). Computational resources and services used in this work were provided by the VSC (Flemish Supercomputer Centre), funded by the FWO and the Flemish Government-Department EWI.

Aims. The first relativistic solar proton event of solar cycle 25 (SC25) was detected on 28 October 2021 by neutron monitors (NMs) on the ground and particle detectors onboard spacecraft in the near-Earth space. This is the first ground level enhancement (GLE) of the current cycle. A detailed reconstruction of the NM response together with the identification of the solar eruption that generated these particles is investigated based on in-situ and remote-sensing measurements. Methods. In-situ proton observations from a few MeV to $\sim$500 MeV were combined with the detection of a solar flare in soft X-rays (SXRs), a coronal mass ejection (CME), radio bursts and extreme ultraviolet (EUV) observations to identify the solar origin of the GLE. Timing analysis was performed and a relation to the solar sources was outlined. Results. GLE73 reached a maximum particle rigidity of $\sim$2.4 GV and is associated with type III, type II, type IV radio bursts and an EUV wave. A diversity of time profiles recorded by NMs was observed. This points to an anisotropic nature of the event. The peak flux at E$>$10 MeV was only $\sim$30 pfu and remained at this level for several days. The release time of $\geq$1 GV particles was found to be $\sim$15:40 UT. GLE73 had a moderately hard rigidity spectrum at very high energies ($\gamma$ $\sim$5.5). Comparison of GLE73 to previous GLEs with similar solar drivers is performed, Comment: accepted for publication in A&A Let

Kollhoff, A., et al., Context. On 2020 November 29, the first widespread solar energetic particle (SEP) event of solar cycle 25 was observed at four widely separated locations in the inner (≲ 1 AU) heliosphere. Relativistic electrons as well as protons with energies > 50 MeV were observed by Solar Orbiter (SolO), Parker Solar Probe, the Solar Terrestrial Relations Observatory (STEREO)-A and multiple near-Earth spacecraft. The SEP event was associated with an M4.4 class X-ray flare and accompanied by a coronal mass ejection and an extreme ultraviolet (EUV) wave as well as a type II radio burst and multiple type III radio bursts. Aims. We present multi-spacecraft particle observations and place them in context with source observations from remote sensing instruments and discuss how such observations may further our understanding of particle acceleration and transport in this widespread event. Methods. Velocity dispersion analysis (VDA) and time shift analysis (TSA) were used to infer the particle release times at the Sun. Solar wind plasma and magnetic field measurements were examined to identify structures that influence the properties of the energetic particles such as their intensity. Pitch angle distributions and first-order anisotropies were analyzed in order to characterize the particle propagation in the interplanetary medium. Results. We find that during the 2020 November 29 SEP event, particles spread over more than 230° in longitude close to 1 AU. The particle onset delays observed at the different spacecraft are larger as the flare-footpoint angle increases and are consistent with those from previous STEREO observations. Comparing the timing when the EUV wave intersects the estimated magnetic footpoints of each spacecraft with particle release times from TSA and VDA, we conclude that a simple scenario where the particle release is only determined by the EUV wave propagation is unlikely for this event. Observations of anisotropic particle distributions at SolO, Wind, and STEREO-A do not rule out that particles are injected over a wide longitudinal range close to the Sun. However, the low values of the first-order anisotropy observed by near-Earth spacecraft suggest that diffusive propagation processes are likely involved., The CAU Kiel team acknowledges support from the German Federal Ministry for Economic Affairs and Energy, the German Space Agency (Deutsches Zentrum für Luft- und Raumfahrt e.V., DLR) under grants 50OT0901, 50OT1202, 50OT1702, and 50OT2002. A.K. acknowledges financial support from the ANR COROSHOCK project (ANR-17-CE31-0006-01). D.L. acknowledges support from NASA-HGI grant NNX16AF73G and the NASA Program NNH17ZDA001N-LWS. This study has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 101004159 (SERPENTINE). The U. Turku team acknowledges funding by the Academy of Finland (Grant No. 336809). The UAH team acknowledges financial support by the Spanish Ministerio de Ciencia, Innovación y Universidades FEDER/MCIU/AEI Projects ESP2017-88436-R and PID2019-104863RB-I00/AEI/10.13039/501100011033. [...] The Swedish National Space Agency, ESA-PRODEX and all the participating institutes. I.G.R. acknowledges support from NASA programs NNH19ZDA001N-HSR and NNH19ZDA001N-LWS, and the STEREO mission. IRFU team acknowledges support from the Swedish National Space Agency grant 20/136. A.A. acknowledges the support of the Spanish Ministerio de Ciencia e Innovación (MICINN) under grant PID2019- 105510GB-C31 and through the ‘Center of Excellence María de Maeztu 2020-2023’ award to the ICCUB (CEX2019-000918- M). F.C. acknowledges the financial support by the Spanish MINECO-FPI-2016 predoctoral grant with FSE. E.A. would like to acknowledge the financial support by the Academy of Finland (Postdoctoral Grant No 322455). V.K. acknowledges the support by NASA under grants 18-2HSWO218_2-0010 and 19-HSR-19_2-0143. Solar Orbiter magnetometer operations are funded by the UK Space Agency (grant ST/T001062/1). T.S.H. is supported by STFC grant ST/S000364/1. T.L. and S.D. acknowledge support from the UK Science and Technology Facilities Council (STFC; grant ST/R000425/1).

Solar activity in all its varied manifestations is driven by the magnetic field. Two global quantities are particularly important for many purposes, the Sun’s total and open magnetic flux, which can be computed from sunspot number records using models. Such sunspot-driven models, however, do not take into account the presence of magnetic flux during grand minima, such as the Maunder minimum. Here we present a major update of a widely used simple model, which now takes into account the observation that the distribution of all magnetic features on the Sun follows a single power law. The exponent of the power law changes over the solar cycle. This allows for the emergence of small-scale magnetic flux even when no sunspots have been present for multiple decades and leads to non-zero total and open magnetic flux also in the deepest grand minima, such as the Maunder minimum, thus overcoming a major shortcoming of the earlier models. The results of the updated model compare well with the available observations and reconstructions of the solar total and open magnetic flux. This opens up the possibility of improved reconstructions of the sunspot number from time series of the cosmogenic isotope production rate.

Wimmer-Schweingruber, R. F., et al., [Context] Solar Orbiter strives to unveil how the Sun controls and shapes the heliosphere and fills it with energetic particle radiation. To this end, its Energetic Particle Detector (EPD) has now been in operation, providing excellent data, for just over a year. [Aims] EPD measures suprathermal and energetic particles in the energy range from a few keV up to (near-) relativistic energies (few MeV for electrons and about 500 MeV nuc-1 for ions). We present an overview of the initial results from the first year of operations and we provide a first assessment of issues and limitations. In addition, we present areas where EPD excels and provides opportunities for significant scientific progress in understanding how our Sun shapes the heliosphere. [Methods] We used the solar particle events observed by Solar Orbiter on 21 July and between 10 and 11 December 2020 to discuss the capabilities, along with updates and open issues related to EPD on Solar Orbiter. We also give some words of caution and caveats related to the use of EPD-derived data. [Results] During this first year of operations of the Solar Orbiter mission, EPD has recorded several particle events at distances between 0.5 and 1 au from the Sun. We present dynamic and time-averaged energy spectra for ions that were measured with a combination of all four EPD sensors, namely: the SupraThermal Electron and Proton sensor (STEP), the Electron Proton Telescope (EPT), the Suprathermal Ion Spectrograph (SIS), and the High-Energy Telescope (HET) as well as the associated energy spectra for electrons measured with STEP and EPT. We illustrate the capabilities of the EPD suite using the 10 and 11 December 2020 solar particle event. This event showed an enrichment of heavy ions as well as 3He, for which we also present dynamic spectra measured with SIS. The high anisotropy of electrons at the onset of the event and its temporal evolution is also shown using data from these sensors. We discuss the ongoing in-flight calibration and a few open instrumental issues using data from the 21 July and the 10 and 11 December 2020 events and give guidelines and examples for the usage of the EPD data. We explain how spacecraft operations may affect EPD data and we present a list of such time periods in the appendix. A list of the most significant particle enhancements as observed by EPT during this first year is also provided., We acknowledge the Spanish Ministerio de Ciencia, Innovación y Universidades for their invaluable support on the development of the EPD suite and the ICU under grants FEDER/MCIU – Agencia Estatal de Investigación/Projects ESP2015-68266-R, ESP2017-88436-R and PID2019-104863RB-I00/AEI/10.13039/501100011033. We also thank acknowledge the financial support by the Spanish MINECO-FPI-2016 predoctoral grant with FSE. We thank the German Federal Ministry for Economic Affairs and Energy and the German Space Agency (Deutsches Zentrum für Luft- und Raumfahrt, e.V., (DLR)) for their unwavering support of STEP, EPT, and HET under grants numbers 50OT0901, 50OT1202, 50OT1702, and 50OT2002; and ESA for supporting the build of SIS under contract number SOL.ASTR.CON.00004, as well as the University of Kiel and the Land Schleswig-Holstein for their support of SIS. We thank NASA headquarters and the GSFC Solar Orbiter project office for support of SIS at APL under contract NNN06AA01C. A.A. acknowledges the support by the Spanish Ministerio de Ciencia e Innovación (MICINN) under grant PID2019- 105510GB-C31 and through the ‘Center of Excellence María de Maeztu 2020-2023’ award to the ICCUB (CEX2019-000918- M). Solar Orbiter magnetometer operations are funded by the UK Space Agency (grant ST/T001062/1). Tim Horbury is supported by STFC grant ST/S000364/1. We acknowledge funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 101004159 (SERPENTINE).

Aran, A., et al., [Context] In mid-June 2020, the Solar Orbiter (SolO) mission reached its first perihelion at 0.51 au and started its cruise phase, with most of the in situ instruments operating continuously. [Aims] We present the in situ particle measurements of the first proton event observed after the first perihelion obtained by the Energetic Particle Detector (EPD) suite on board SolO. The potential solar and interplanetary (IP) sources of these particles are investigated. Methods. Ion observations from ∼20 keV to ∼1 MeV are combined with available solar wind data from the Radio and Plasma Waves (RPW) instrument and magnetic field data from the magnetometer on board SolO to evaluate the energetic particle transport conditions and infer the possible acceleration mechanisms through which particles gain energy. We compare > 17-20 MeV ion count rate measurements for two solar rotations, along with the solar wind plasma data available from the Solar Wind Analyser (SWA) and RPW instruments, in order to infer the origin of the observed galactic cosmic ray (GCR) depressions. [Results] The lack of an observed electron event and of velocity dispersion at various low-energy ion channels and the observed IP structure indicate a local IP source for the low-energy particles. From the analysis of the anisotropy of particle intensities, we conclude that the low-energy ions were most likely accelerated via a local second-order Fermi process. The observed GCR decrease on 19 June, together with the 51.8-1034.0 keV nuc-1 ion enhancement, was due to a solar wind stream interaction region (SIR). The observation of a similar GCR decrease in the next solar rotation favours this interpretation and constitutes the first observation of a recurrent GCR decrease by SolO. The analysis of the recurrence times of this SIR suggests that it is the same SIR responsible for the He events previously measured in April and May. Finally, we point out that an IP structure more complex than a common SIR cannot be discarded, mainly due to the lack of solar wind temperature measurements and the lack of a higher cadence of solar wind velocity observations., EPD was built with funding from Ministerio de Ciencia e Innovación (Spain), Deutsches Zentrum für Luft- und Raumfahrt (Germany), and the European Space Agency (ESA); operations are funded by FEDER/MCI/AEI Projects ESP2017-88436-R and PID2019-104863RBI00/AEI/10.13039/501100011033 (Spain), 50OT 2002 (DLR, Germany), and NASA contract 80MSFC19F0002. The UB team acknowledges the support by the Spanish Ministerio de Ciencia e Innovación (MICINN) under grant PID2019-105510GB-C31 and through the “Center of Excellence María de Maeztu 2020-2023” award to the ICCUB (CEX2019-000918-M). The CAU Kiel team acknowledges support by the German Federal Ministry for Economic Affairs and Energy, the German Space Agency (Deutsches Zentrum für Luft- und Raumfahrt e.V., DLR) under grants 50OT0901, 50OT1202, 50OT1702, 50OT2002, and 50OC1702. The Solar Orbiter magnetometer was funded by the UK Space Agency (grant ST/T001062/1). M.L. and S.B. acknowledge financial support by the Italian MIUR-PRIN grant 2017APKP7T on Circumterrestrial Environment: Impact of Sun-Earth Interaction. N.W. acknowledges funding from the Research Foundation – Flanders (FWO – Vlaanderen, fellowship no. 1184319N). D.L. acknowledges the support from the NASA-HGI grant NNX16AF73G and the NASA Programs NNH17ZDA001N-LWS. I.G.R. and D.L. acknowledge support from NASA program NNH19ZDA001N-LWS. E.S. was supported by a PhD grant awarded by the Royal Observatory of Belgium. V.K. acknowledges the support by NASA under grants 18-2HSWO218_2-0010 and 19-HSR-19_2-0143. S.P. is supported by the projects C14/19/089 (C1 project Internal Funds KU Leuven), G.0D07.19N (FWO-Vlaanderen), SIDEX (ESA Prodex-12), and EUHFORIA 2.0 (funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 870405). IRFU team acknowledges support from the Swedish National Space Agency grant 20/136. F.C. acknowledges the financial support by the Spanish MINECO-FPI-2016 predoctoral grant with FSE. The JHU/APL team is supported under NASA contract NNN06AA01C and thanks NASA headquarters and the NASA/GSFC Solar Orbiter project office for their continuing support. Solar Orbiter Solar Wind Analyser (SWA) data are derived from scientific sensors which have been designed and created, and are operated under funding provided in numerous contracts from the UK Space Agency (UKSA, most recently grant ST/T001356/1), the UK Science and Technology Facilities Council (STFC), the Agenzia Spaziale Italiana (ASI), the Centre National d’Etudes Spatiales (CNES, France), the Centre National de la Recherche Scientifique (CNRS, France), the Czech contribution to the ESA PRODEX programme and NASA.

In the literature different so-called $\kappa$-distribution functions are discussed to fit and model the velocity (or energy) distributions of solar wind species, pickup ions or magnetospheric particles. Here we introduce a generalized (isotropic) $\kappa$-distribution as a "cookbook", which admits as special cases, or "recipes", all the other known versions of $\kappa$-models. A detailed analysis of the generalized distribution function is performed, providing general analytical expressions for the velocity moments, Debye length, and entropy, and pointing out a series of general requirements that plasma distribution functions should satisfy. From a contrasting analysis of the recipes found in the literature, we show that all of them lead to almost the same macroscopic parameters with a small standard deviation between them. However, one of these recipes called the regularized $\kappa$-distribution provides a functional alternative for macroscopic parameterization without any constraint for the power-law exponent $\kappa$., Comment: 20 pager, 9 Figures

Large magnetic structures are launched away from the Sun during solar eruptions. They are observed as (interplanetary) coronal mass ejections (ICMEs or CMEs) with coronal and heliospheric imagers. A fraction of them are observed insitu as magnetic clouds (MCs). Fitting these structures properly with a model requires a better understanding of their evolution. In situ measurements are done locally when the spacecraft trajectory crosses the magnetic configuration. These observations are taken for different elements of plasma and at different times, and are therefore biased by the expansion of the magnetic configuration. This aging effect leads to stronger magnetic fields measured at the front than at the rear of MCs, an asymmetry often present in MC data. However, can the observed asymmetry be explained quantitatively only from the expansion? Based on self-similar expansion, we derive a method to estimate the expansion rate from observed plasma velocity. We next correct for the aging effect both the observed magnetic field and the spatial coordinate along the spacecraft trajectory. This provides corrected data as if the MC internal structure was observed at the same time. We apply the method to 90 best observed MCs near Earth (1995-2012). The aging effect is the main source of the observed magnetic asymmetry only for 28\% of MCs. After correcting the aging effect, the asymmetry is almost symmetrically distributed between MCs with a stronger magnetic field at the front and those at the rear of MCs. The proposed method can efficiently remove the aging bias within insitu data of MCs, and more generally of ICMEs. This allows one to analyse the data with a spatial coordinate, such as in models or remote sensing observations., Comment: 13 pages, 8 figures

Magnetic clouds (MCs) are large-scale interplanetary transient structures in the heliosphere that travel from the Sun into the interplanetary medium. The internal magnetic field lines inside the MCs are twisted, forming a flux rope (FR). This magnetic field structuring is determined by its initial solar configuration, by the processes involved during its eruption from the Sun, and by the dynamical evolution during its interaction with the ambient solar wind. One of the most important properties of the magnetic structure inside MCs is the twist of the field lines forming the FR (the number of turns per unit length). The detailed internal distribution of twist is under debate mainly because the magnetic field (B) in MCs is observed only along the spacecraft trajectory, and thus it is necessary to complete observations with theoretical assumptions. Estimating the twist from the study of a single event is difficult because the field fluctuations significantly increase the noise of the observed B time series and thus the bias of the deduced twist. The superposed epoch applied to MCs has proven to be a powerful technique, permitting the extraction of their common features, and removing the peculiarity of individual cases. We apply a superposed epoch technique to analyse the magnetic components in the local FR frame of a significant sample of moderately asymmetric MCs observed at 1 au. From the superposed profile of B components in the FR frame, we determine the typical twist distribution in MCs. The twist is nearly uniform in the FR core (central half part), and it increases moderately, up to a factor two, towards the MC boundaries. This profile is close to the Lundquist field model limited to the FR core where the axial field component is above about one-third of its central value., 13 pages, 11 figures

International audience; The magnetometer instrument on the Solar Orbiter mission is designed to measure the magnetic field local to the spacecraft continuously for the entire mission duration. The need to characterise not only the background magnetic field but also its variations on scales from far above to well below the proton gyroscale result in challenging requirements on stability, precision, and noise, as well as magnetic and operational limitations on both the spacecraft and other instruments. The challenging vibration and thermal environment has led to significant development of the mechanical sensor design. The overall instrument design, performance, data products, and operational strategy are described.

In this paper we present a statistical study of the kinematics of 28894 coronal mass ejections (CMEs) recorded by the Large Angle and Spectrometric Coronagraph (LASCO) on board the Solar and Heliospheric Observatory spacecraft from 1996 until mid-2017. The initial acceleration phase is characterized by a rapid increase in CME velocity just after eruption in the inner corona. This phase is followed by a non-significant residual acceleration (deceleration) characterized by an almost constant speed of CMEs. We demonstrate that the initial acceleration is in the range 0.24-2616 ms-2 with median (average) value of 57 ms-2 (34 ms-2 ) and it takes place up to a distance of about 28 solar radius with median (average) value of 7.8 solar radius (6 solar radius). Additionally, the initial acceleration is significant in the case of fast CMEs (V > 900 kms-1 ), where the median (average) values are about 295 ms-2 (251 ms-2 ), respectively, and much weaker in the case of slow CMEs (V < 250 kms-1 ), where the median (average) values are about 18 ms-2 (17 ms-2 ), respectively. We note that the significant driving force (Lorentz force) can operate up to a distance of 6 solar radius from the Sun during the first 2 hours of propagation. We found a significant anti-correlation between the initial acceleration magnitude and the acceleration duration, whereas the residual acceleration covers a range from -1224 to 0 ms-2 with a median (average) value of -34 ms-2 (-17 ms-2 ). One intriguing finding is that the residual acceleration is much smaller during the 24th cycle in comparison to the 23rd cycle of solar activity. Our study has also revealed that the considered parameters, initial acceleration (ACC INI ), residual acceleration (ACC RES ), maximum velocity (V MAX ), and time at maximum velocity (Time MAX ) mostly follow solar cycles and the intensities of the individual cycle., 12 pages, 14 figures

International audience; Non-recurrent short-term variations of the galactic cosmic-ray (GCR) flux above 70 MeV n−1 were observed between 2016 February 18 and 2017 July 3 on board the European Space Agency LISA Pathfinder (LPF) mission orbiting around the Lagrange point L1 at 1.5 × 106 km from Earth. The energy dependence of three Forbush decreases is studied and reported here. A comparison of these observations with others carried out in space down to the energy of a few tens of MeV n−1 shows that the same GCR flux parameterization applies to events of different intensity during the main phase. FD observations in L1 with LPF and geomagnetic storm occurrence are also presented. Finally, the characteristics of GCR flux non-recurrent variations (peaks and depressions) of duration

The behavior of the 2D turbulence power spectrum at the lowest wavenumbers has a significant effect on the perpendicular diffusion coefficients of charged particles in the heliosphere derived from various scattering theories, and subsequently used to model the transport of cosmic rays (CRs) and solar energetic particles. In this regard, the lengthscale at which the energy-containing range begins, as opposed to that at which the inertial range commences, is of particular interest. This 2D outerscale has, however, never before been directly observed. Recently, direct measurements of magnetic islands in the solar wind have been reported by various authors. Assuming that these may provide an estimate of the 2D ultrascale, the direct calculation of the 2D outerscale becomes possible, should an observationally motivated form for the 2D turbulence power spectrum be employed. This study presents the results of such a calculation and provides comparisons of these with previous estimates of the 2D outerscale. Furthermore, the sensitivity of galactic CR electron intensities, calculated using a 3D ab initio CR modulation model, is demonstrated, and conclusions are drawn therefrom

A complete theory for the complex interaction between solar energetic particles and the turbulent interplanetary magnetic field remains elusive. In this work we aim to contribute towards such a theory by modelling the propagation of solar energetic particle electrons in plasma wave turbulence. We specify a background turbulence spectrum, as constrained through observations, calculate the transport coefficients from first principles, and simulate the propagation of these electrons in the inner heliosphere. We have also, for the first time, included dynamical effects into the perpendicular diffusion coefficient. We show that such a "physics-first" approach can lead to reasonable results, when compared qualitatively to observations. In addition, we include the effect of wave growth/damping due to streaming electrons and show that these particles can significantly alter the turbulence levels close to the Sun for the largest events., Comment: Submitted to ApJ

Heliospheric modulation conditions were unusually quite during the last solar minimum activity between Solar Cycles 23/24. Fortunately, the PAMELA space-experiment measured six-month averaged Galactic positron spectra for the period July 2006 to December 2009, over an energy range of 80 MeV to 30 GeV, which is important for solar modulation. The highest level of Galactic positrons was observed at Earth during the July-December 2009 period. A well-established, comprehensive three-dimensional (3D) numerical model is applied to study the modulation of the observed positron spectra. This model had been used previously to understand the modulation of Galactic protons and electrons also measured by PAMELA for the same period. First, a new very local interstellar spectrum for positrons is constructed, using the well-known GALPROP code together with the mentioned PAMELA observations. The 3D model is used to distinguish between the dominant mechanisms responsible for the heliospheric modulation of Galactic positrons, and to understand the effect of particle drift during this unusual minimum in particular, which is considered diffusion dominant, even though particle drift still had a significant role in modulating positrons. Lastly, the expected intensity of Galactic positrons during an A>0 polarity minimum, with similar heliospheric conditions than for 2006-2009, is predicted to be higher than what was observed by PAMELA for the 2006-2009 unusual minimum., 13 Pages, 12 figures, Submitted to Astrophysical Journal

The local interstellar spectra (LIS's) for galactic cosmic rays (CR's) cannot be directly observed at the Earth below certain energies, because of solar modulation in the heliosphere. With Voyager 1 crossing the heliopause in 2012, in situ experimental LIS data below 100 MeV/nuc can now constrain computed galactic CR spectra. Using galactic propagation models, galactic electron, proton and light nuclei spectra can be computed, now more reliably as very LIS's. Using the Voyager 1 observations made beyond the heliopause, and the observations made by the PAMELA experiment in Earth orbit for the 2009 solar minimum, as experimental constraints, we simultaneously reproduced the CR electron, proton, Helium and Carbon observations by implementing the GALPROP code. Below about 30 GeV/nuc solar modulation has a significant effect and a comprehensive three-dimensional (3D) numerical modulation model is used to compare the computed spectra with the observed PAMELA spectra possible at these energies. Subsequently the computed LIS's can be compared over as wide a range of energies as possible. The simultaneous calculation CR spectra with a single propagation model allows the LIS's for positrons, Boron and Oxygen to also be inferred. This implementation of the most comprehensive galactic propagation model (GALPROP), alongside a sophisticated solar modulation model to compute CR spectra for comparison with both Voyager 1 and PAMELA observations over a wide energy range, allows us to present new self-consistent very LIS's (and expressions) for electrons, positrons, protons, Helium, Carbon, Boron and Oxygen for the energy range of 3 MeV/nuc to 100 GeV/nuc.

Insights into the processes of Solar Energetic Particle (SEP) propagation are essential for understanding how solar eruptions affect the radiation environment of near-Earth space. SEP propagation is influenced by turbulent magnetic fields in the solar wind, resulting in stochastic transport of the particles from their acceleration site to Earth. While the conventional approach for SEP modelling focuses mainly on the transport of particles along the mean Parker spiral magnetic field, multi-spacecraft observations suggest that the cross-field propagation shapes the SEP fluxes at Earth strongly. However, adding cross-field transport of SEPs as spatial diffusion has been shown to be insufficient in modelling the SEP events without use of unrealistically large cross-field diffusion coefficients. Recently, Laitinen et al. [ApJL 773 (2013b); A&A 591 (2016)] demonstrated that the early-time propagation of energetic particles across the mean field direction in turbulent fields is not diffusive, with the particles propagating along meandering field lines. This early-time transport mode results in fast access of the particles across the mean field direction, in agreement with the SEP observations. In this work, we study the propagation of SEPs within the new transport paradigm, and demonstrate the significance of turbulence strength on the evolution of the SEP radiation environment near Earth. We calculate the transport parameters consistently using a turbulence transport model, parametrised by the SEP parallel scattering mean free path at 1 AU, l k ,and show that the parallel and cross-field transport are connected, with conditions resulting in slow parallel transport corresponding to wider events. We find a scaling smax1=l k +1=4 for the Gaussian fitting of the longitudinal distribution of maximum intensities. The longitudes with highest intensities are shifted towards the west for strong scattering conditions. Our results emphasise the importance of understanding both the SEP transport and the interplanetary turbulence conditions for modelling and predicting the SEP radiation environment at Earth., SCOPUS: ar.j, info:eu-repo/semantics/published

Aims: The heliospheric current sheet (HCS) has been observed to be southward shifted in the late declining to minimum phase of the solar cycle. Here we study the existence of a simultaneous shift in the heliosphere and in the corona using a robust new method. Methods: We use the synoptic maps of the photospheric field of the Wilcox Solar Observatory (WSO) and the Mount Wilson Observatory (MWO) together with the potential field source surface (PFSS) model to calculate the coronal magnetic field and compare it with the simultaneous heliospheric magnetic field of the NASA/NSSDC OMNI 2 dataset. We divide the magnetic field into the two sectors, towards (T) and away (A) from the Sun, and calculate how often the sector polarities at 1 AU and in the corona match each other. We divide the sectors both at 1 AU and in the corona. We also calculate the annual (T − A)/(T + A) ratios of sector occurrence both at 1 AU and in the corona. Results: We verify that the HCS/neutral line is southward shifted both in the corona and heliosphere. We find that the coronal shift is systematically larger than the simultaneous heliospheric shift. Conclusions: The fact that the southward shift of the coronal neutral line is larger than the simultaneous shift of the heliospheric current sheet at 1 AU implies that the radial evolution of the magnetic field between the two sites is different between the northern and southern hemispheres.

The final publication is available at IOS Press through http://dx.doi.org/10.3847/1538-4357/aaa774 Galactic cosmic-ray (GCR) energy spectra observed in the inner heliosphere are modulated by the solar activity, the solar polarity and structures of solar and interplanetary origin. A high counting rate particle detector (PD) aboard LISA Pathfinder, meant for subsystems diagnostics, was devoted to the measurement of GCR and solar energetic particle integral fluxes above 70 MeV n-1 up to 6500 counts s-1. PD data were gathered with a sampling time of 15 s. Characteristics and energy dependence of GCR flux recurrent depressions and of a Forbush decrease dated 2016 August 2 are reported here. The capability of interplanetary missions, carrying PDs for instrument performance purposes, in monitoring the passage of interplanetary coronal mass ejections is also discussed.

The multitude of recent multi-point spacecraft observations of solar energetic particle (SEP) events have made it possible to study the longitudinal distribution of SEPs in great detail. SEPs, even those accelerated during impulsive events, show a much wider than expected longitudinal extent, bringing into question the processes responsible for their transport perpendicular to the local magnetic field. In this paper we examine some aspects of perpendicular transport by including perpendicular diffusion into a numerical SEP transport model that simulates the propagation of impulsively accelerated SEP electrons in the ecliptic plane. We find that: (i) The pitch-angle dependence of the perpendicular diffusion coefficient is an important, and currently mainly overlooked, transport parameter. (ii) SEP intensities are generally asymmetric in longitude, being enhanced towards the west of optimal magnetic connection to the acceleration region. (iii) The maximum SEP intensity may also be shifted (parameter dependently) away from the longitude of best magnetic connectivity at 1 AU. We also calculate the maximum intensity, the time of maximum intensity, the onset time and the maximum anisotropy as a function of longitude at Earth's orbit and compare the results, in a qualitative fashion, to recent spacecraft observations.

Magnetic clouds (MCs) are twisted magnetic structures ejected from the Sun and probed by in situ instruments. They are typically modeled as flux ropes (FRs). The determination of the FR global characteristics requires the estimation of the FR axis orientation. Among the developed methods, the minimum variance (MV) is the most flexible, and features only a few assumptions. However, as other methods, MV has biases. We aim to investigate the limits of the method and extend it to a less biased method. We first identified the origin of the biases by testing the MV method on cylindrical and elliptical models with a temporal expansion comparable to the one observed in MCs. Then, we developed an improved MV method to reduce these biases. In contrast with many previous publications we find that the ratio of the MV eigenvalues is not a reliable indicator of the precision of the derived FR axis direction. Next, we emphasize the importance of the FR boundaries selected since they strongly affect the deduced axis orientation. We have improved the MV method by imposing that the same amount of azimuthal flux should be present before and after the time of closest approach to the FR axis. We emphasize the importance of finding simultaneously the FR axis direction and the location of the boundaries corresponding to a balanced magnetic flux, so as to minimize the bias on the deduced FR axis orientation. This method can also define an inner flux-balanced sub-FR. We show that the MV results are much less biased when a compromise in size of this sub-FR is achieved. For weakly asymmetric field temporal profiles, the improved MV provides a very good determination of the FR axis orientation. The main remaining bias is moderate (lower than 6 degrees) and is present mostly on the angle between the flux rope axis and the plane perpendicular to the Sun-Earth direction., 15 pages, 8 figures

The photospheric magnetic field is the source of the coronal and heliospheric magnetic fields (HMF), but their mutual correspondence is non-trivial and depends on the phase of the solar cycle. The photospheric field during the HMF sector crossings observed at 1 AU has been found to contain enhanced field intensities and definite polarity ordering, forming regions called Hale boundaries. Here we separately study the structure of the photospheric field during the HMF sector crossings during Solar Cycles 21 – 24 for the four phases of each solar cycle. We use a refined version of Svalgaard’s list of major HMF sector crossings, mapped to the Sun using the solar wind speed observed at Earth, and the daily level-3 magnetograms of the photospheric field measured at the Wilcox Solar Observatory in 1976 – 2016. We find that the structure of the photospheric field corresponding to the HMF sector crossings and the existence and properties of the corresponding Hale bipolar regions varies significantly with solar cycle, solar cycle phase, and hemisphere. The Hale boundaries in more than half of the ascending, maximum, and declining phases are clear and statistically significant. The clearest Hale boundaries are found during the (+,−) HMF crossings in the northern hemisphere of odd Cycles 21 and 23, but less systematical during the (+,−) crossings in the southern hemisphere of even Cycles 22 and 24. No similar difference between odd and even cycles is found for the (−,+) crossings. This shows that the northern hemisphere has a more organized Hale pattern overall. The photospheric field distribution also depicts a larger area for the field of the northern hemisphere during the declining and minimum phases, in a good agreement with the bashful ballerina phenomenon.

The transport of cosmic rays (CRs) in the heliosphere is determined by the properties of the solar wind plasma. The heliospheric plasma environment has been probed by spacecraft for decades and provides a unique opportunity for testing transport theories. Of particular interest for the three-dimensional (3D) heliospheric CR transport are structures such as corotating interaction regions (CIRs), which, due to the enhancement of the magnetic field strength and magnetic fluctuations within and due to the associated shocks as well as stream interfaces, do influence the CR diffusion and drift. In a three-fold series of papers, we investigate these effects by modeling inner-heliospheric solar wind conditions with the numerical magnetohydrodynamic (MHD) framework Cronos (Wiengarten et al., referred as Paper I), and the results serve as input to a transport code employing a stochastic differential equation approach (this paper). While, in Paper I, we presented results from 3D simulations with Cronos, the MHD output is now taken as an input to the CR transport modeling. We discuss the diffusion and drift behavior of Galactic cosmic rays using the example of different theories, and study the effects of CIRs on these transport processes. In particular, we point out the wide range of possible particle fluxes at a given point in space resulting from these different theories. The restriction of this variety by fitting the numerical results to spacecraft data will be the subject of the third paper of this series.

Context. Galactic cosmic ray (GCR) intensities around solar minimum times are modulated by magnetic drifts that depend on the overall solar polarity. GCR intensities reach a higher but more narrow peak during negative minima than during positive minima. However, despite these higher intensities, the variation of GCRs over timescales of solar rotation is smaller during negative minima than during positive minima. Aims. We study the variation of GCR intensity over the 27-day synodic solar rotation and over the 14-day half-rotation, in particular the long-term trend and cyclic pattern of this variation, and propose a unifying explanation for the observations. Methods. We used two high-latitude neutron monitors, Oulu and Apatity, which are most sensitive to the low-energy part of the GCR spectrum and thereby more strongly affected by the changes in the conditions of the local heliosphere. We calculated the yearly mean amplitudes of the GCR intensity variation during the full solar rotation ( A 27 ) and half-rotation ( A 14 ) in 1964–2016. Results. We verify that the A 27 and A 14 amplitudes exhibit a clear 22-yr Hale cycle during solar minima at both stations, with larger amplitudes in positive minima. We find that the mean amplitude of the Hale cycle is about 30–45% of the mean amplitude for A 14 , while is only about 15–30% for A 27 . We also find that all amplitudes depict a declining long-term trend, which we suggest is due to the weakening of solar polar magnetic fields during the last four solar cycles and the ensuing latitudinal widening of the heliospheric current sheet (HCS) region. An exceptionally wide HCS region during the last solar minimum, when A 14 reached its all-time minimum, is demonstrated by Ulysses probe observations. Conclusions. Our results emphasize the effect of polarity-dependent drift and the properties of the HCS in modulating the variation of GCR intensity during solar rotation in solar minimum times. The second rotation harmonic yields a larger Hale amplitude than the first because it is more probable for the Earth to be outside the HCS only once during the rotation than twice or more, which more strongly reduces A 14 during negative polarity times than A 27 . With the HCS region widening from minimum to minimum, the decrease in A 14 is relatively faster than in A 27 .

In this work, we analyze the long-term cosmic-ray modulation observed by the Hermanus neutron monitor, which is the detector with the longest cosmic-ray record, from 1957 July. For our study we use the force-field approximation to the cosmic-ray transport equation, and the newest results on the mean free paths from the scattering theory. We compare the modulation parameter (phgr) with different rigidity (P) dependences: P, P 2, and P 2/3. We correlate them with solar and interplanetary parameters. We found that (1) these rigidity dependences properly describe the modulation, (2) long-term cosmic-ray variations are better correlated with the magnitude of the heliospheric magnetic field (HMF) than the sunspot number, solar wind speed, and tilt angle of the HMF, and (3) the theoretical dependence of the parallel mean free path on the magnetic field variance is in agreement with the modulation parameter and therefore with the neutron monitor record. We also found that the force-field approximation is not able to take into account the effects of three-dimensional particle transport, showing a poor correlation with the perpendicular mean free path

Recent studies of particle acceleration in the heliosphere have revealed a new mechanism that can locally energize particles up to several MeV nucleon–1. Stream–stream interactions, as well as the heliospheric current sheet (CS)—stream interactions, lead to formation of large magnetic cavities, bordered by strong CSs, which in turn produce secondary CSs and dynamical small-scale magnetic islands (SMIs) of ~0.01 au or less owing to magnetic reconnection. It has been shown that particle acceleration or reacceleration occurs via stochastic magnetic reconnection in dynamical SMIs confined inside magnetic cavities observed at 1 au. The study links the occurrence of CSs and SMIs with characteristics of intermittent turbulence and observations of energetic particles of keV–MeV nucleon–1 energies at ~5.3 au. We analyze selected samples of different plasmas observed by Ulysses during a widely discussed event, which was characterized by a series of high-speed streams of various origins that interacted beyond Earth's orbit in 2005 January. The interactions formed complex conglomerates of merged interplanetary coronal mass ejections, stream/corotating interaction regions, and magnetic cavities. We study properties of turbulence and associated structures of various scales. We confirm the importance of intermittent turbulence and magnetic reconnection in modulating solar energetic particle flux and even local particle acceleration. Coherent structures, including CSs and SMIs, play a significant role in the development of secondary stochastic particle acceleration, which changes the observed energetic particle flux time–intensity profiles and increases the final energy level to which energetic particles can be accelerated in the solar wind

Since 2011 May, the Alpha Magnetic Spectrometer (AMS-02) on board the International Space Station has provided monthly cosmic proton fluxes for various low-rigidity levels (P

Interplanetary imagers provide 2D projected views of the densest plasma parts of interplanetary coronal mass ejections (ICMEs) while in situ measurements provide magnetic field and plasma parameter measurements along the spacecraft trajectory, so along a 1D cut. As such, the data only give a partial view of their 3D structures. By studying a large number of ICMEs, crossed at different distances from their apex, we develop statistical methods to obtain a quantitative generic 3D shape of ICMEs. In a first approach we theoretically obtain the expected statistical distribution of the shock-normal orientation from assuming simple models of 3D shock shapes, including distorted profiles, and compare their compatibility with observed distributions. In a second approach we use the shock normal and the flux rope axis orientations, as well as the impact parameter, to provide statistical information across the spacecraft trajectory. The study of different 3D shock models shows that the observations are compatible with a shock symmetric around the Sun-apex line as well as with an asymmetry up to an aspect ratio around 3. Moreover, flat or dipped shock surfaces near their apex can only be rare cases. Next, the sheath thickness and the ICME velocity have no global trend along the ICME front. Finally, regrouping all these new results and the ones of our previous articles, we provide a quantitative ICME generic 3D shape, including the global shape of the shock, the sheath and the flux rope. The obtained quantitative generic ICME shape will have implications for several aims. For example, it constrains the output of typical ICME numerical simulations. It is also a base to develop deeper studies of the transport of high energy solar and cosmic particles during an ICME propagation, as well as for modeling and forecasting space weather conditions near Earth., 14 pages, 13 figures

Interplanetary coronal mass ejections (ICMEs) are the interplanetary manifestations of solar eruptions. The overtaken solar wind forms a sheath of compressed plasma at the front of ICMEs. Magnetic clouds (MCs) are a subset of ICMEs with specific properties (e.g. the presence of a flux rope). When ICMEs pass near Earth, ground observations indicate that the flux of galactic cosmic rays (GCRs) decreases. The main aims of this paper are to find: common plasma and magnetic properties of different ICME sub-structures, and which ICME properties affect the flux of GCRs near Earth. We use a superposed epoch method applied to a large set of ICMEs observed \insitu\ by the spacecraft ACE, between 1998 and 2006. We also apply a superposed epoch analysis on GCRs time series observed with the McMurdo neutron monitors. We find that slow MCs at 1 AU have on average more massive sheaths. We conclude that it is because they are more effectively slowed down by drag during their travel from the Sun. Slow MCs also have a more symmetric magnetic field and sheaths expanding similarly as their following MC, while in contrast, fast MCs have an asymmetric magnetic profile and a compressing sheath in compression. In all types of MCs, we find that the proton density and the temperature, as well as the magnetic fluctuations can diffuse within the front of the MC due to 3D reconnection. Finally, we derive a quantitative model which describes the decrease of cosmic rays as a function of the amount of magnetic fluctuations and field strength. The obtained typical profiles of sheath/MC/GCR properties corresponding to slow, mid, and fast ICMEs, can be used for forecasting/modelling these events, and to better understand the transport of energetic particles in ICMEs. They are also useful for improving future operative space weather activities., Comment: 13 pages, 6 figures, paper accepted in A&A

We investigate the evolution of a well-observed, long-lived, low-latitude coronal hole (CH) over 10 solar rotations in the year 2012. By combining extreme ultraviolet (EUV) imagery from the Solar TErrestrial RElations Observatories (STEREO-A/B) and the Solar Dynamics Observatory (SDO), we are able to track and study the entire evolution of the CH having a continuous 360° coverage of the Sun. The remote sensing data are investigated together with in situ solar wind plasma and magnetic field measurements from STEREO-A/B, the Advanced Composition Explorer (ACE), and WIND. From this, we obtain how different evolutionary states of the CH as observed in the solar atmosphere (changes in EUV intensity and area) affect the properties of the associated high-speed stream measured at 1 au. Most distinctly pronounced for the CH area, three development phases are derived: (a) growing, (b) maximum, and (c) decaying phase. During these phases the CH area (a) increases over a duration of around three months from about 1 · 1010 km2 to 6 · 1010 km2, (b) keeps a rather constant area for about one month of >9 · 1010 km2, and (c) finally decreases in the following three months below 1 · 1010 km2 until the CH cannot be identified anymore. The three phases manifest themselves also in the EUV intensity and in in situ measured solar wind proton bulk velocity. Interestingly, the three phases are related to a different range in solar wind speed variations, and we find for the growing phase a range of 460–600 km s−1, for the maximum phase 600–720 km s−1, and for the decaying phase a more irregular behavior connected to slow and fast solar wind speeds of 350–550 km s−1.

Context. Current solar energetic particle (SEP) propagation models describe the effects of interplanetary plasma turbulence on SEPs as diffusion, using a Fokker-Planck (FP) equation. However, FP models cannot explain the observed fast access of SEPs across the average magnetic field to regions that are widely separated in longitude within the heliosphere without using unrealistically strong cross-field diffusion. Aims. We study whether the recently suggested early non-diffusive phase of SEP propagation can explain the wide SEP events with realistic particle transport parameters. Methods. We used a novel model that accounts for the SEP propagation along field lines that meander as a result of plasma turbulence. Such a non-diffusive propagation mode has been shown to dominate the SEP cross-field propagation early in the SEP event history. We compare the new model to the traditional approach, and to SEP observations. Results. Using the new model, we reproduce the observed longitudinal extent of SEP peak fluxes that are characterised by a Gaussian profile with σ = 30-50°, while current diffusion theory can only explain extents of 11° with realistic diffusion coefficients. Our model also reproduces the timing of SEP arrival at distant longitudes, which cannot be explained using the diffusion model. Conclusions. The early onset of SEPs over a wide range of longitudes can be understood as a result of the effects of magnetic field-line random walk in the interplanetary medium and requires an SEP transport model that properly describes the non-diffusive early phase of SEP cross-field propagation., SCOPUS: ar.j, info:eu-repo/semantics/published

The higher charge states found in slow ($, Comment: Submitted to the Astrophysical Journal

Precision measurements of the electron component of cosmic radiation provide important information about the origin and propagation of cosmic rays in the Galaxy not accessible from the study of cosmic-ray nuclear components due to their differing diffusion and energy-loss processes. However, when measured near Earth, the effects of propagation and modulation of Galactic cosmic rays in the heliosphere, particularly significant for energies up to at least 30 GeV, must be properly taken into account. In this paper the electron (e(-)) spectra measured by the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics down to 70 MeV from 2006 July to 2009 December over six-month time intervals are presented. Fluxes are compared with a state-of-the-art three-dimensional model of solar modulation that reproduces the observations remarkably well.

The last solar minimum activity period, and the consequent minimum modulation conditions for cosmic rays, was unusual. The highest levels of Galactic protons were recorded at Earth in late 2009 in contrast to expectations. A comprehensive model was used to study the proton modulation for the period from 2006 to 2009 in order to determine what basic processes were responsible for solar modulation during this period and why it differs from proton modulation during previous solar minimum modulation periods. This established model is now applied to studying the solar modulation of electron spectra as observed for 80 MeV–30 GeV by the PAMELA space detector from mid-2006 to the end of 2009. Over this period the heliospheric magnetic field had decreased significantly until the end of 2009 while the waviness of the heliospheric current sheet decreased moderately and the observed electron spectra increased by a factor of ∼1.5 at 1.0 GeV to ∼3.5 at 100 MeV. In order to reproduce the modulation evident from seven consecutive semesters, the diffusion coefficients had to increase moderately while maintaining the basic rigidity dependence. It is confirmed that the main diffusion coefficients are independent of rigidity below ∼0.5 GV, while the drift coefficient had to be reduced below this value. The 2006–2009 solar minimum epoch indeed was different than previously observed minima, at least since the beginning of the space exploration era. This period could be called “diffusion-dominated” as was also found for the modulation of protons.

Based on a hybrid galactic cosmic-ray transport model, which incorporated MHD global heliospheric data into Parker's cosmic-ray transport equation, we studied the behavior of the transport of galactic cosmic rays and the corresponding gradients in their flux near the heliopause (HP). We found that, (1) by increasing the ratio of the parallel diffusion coefficient to the perpendicular diffusion coefficient in the interstellar magnetic field of the outer heliosheath, the simulated radial flux near the HP increases as well. As the ratio multiplying factor reached 1010, the radial flux experienced a sudden jump near the HP, similar to what Voyager 1 observed in 2012. (2) The effect of changing the diffusion coefficients' ratio on the radial flux variation depends on the energy of the cosmic rays, the lower the energy, the more pronounced the effect is. (3) The magnitude of the diffusion coefficients also affect the radial flux near the HP, the modulation beyond the HP varies by adjusting the magnitude multiplying factor.

During the recent prolonged solar minimum of cycle 23/24, the PAMELA detector measured 27-day averaged Galactic proton energy spectra over the energy range that is important for solar modulation. By comparing these spectra to computed spectra from a three-dimensional model that contains all of the important heliospheric modulation processes, the recent minimum can be studied in detail from a modulation perspective. This was done by setting up a realistic heliosphere in the model, and reproducing a representative selection of seven intermittent PAMELA spectra, separated by approximately six months, from 2006 July to 2009 December. Additionally, a new very local interstellar proton spectrum was constructed using measurements below 600 MeV from Voyager 1, taken beyond the heliopause, combined with PAMELA and AMS-02 measurements above 30 GeV at the Earth. As a result of the extreme minimum modulation conditions that governed the recent solar minimum, the highest ever Galactic cosmic ray spectrum at Earth was observed by PAMELA at the end of 2009. It was found that, apart from the self-consistent changes in the heliospheric current sheet and the heliospheric magnetic field over time, additional increases in the mean free paths during this period were required below ~4 GV in order to reproduce the intensities observed by PAMELA.

Aims. The phase scintillation of the European Space Agency's (ESA) Venus Express (VEX) spacecraft telemetry signal was observed at X-band (\lambda = 3.6 cm) with a number of radio telescopes of the European VLBI Network (EVN) in the period 2009-2013. Methods. We found a phase fluctuation spectrum along the Venus orbit with a nearly constant spectral index of -2.42 +/-0.25 over the full range of solar elongation angles from 0{\deg} to 45{\deg}, which is consistent with Kolmogorov turbulence. Radio astronomical observations of spacecraft signals within the solar system give a unique opportunity to study the temporal behaviour of the signal's phase fluctuations caused by its propagation through the interplanetary plasma and the Earth's ionosphere. This gives complementary data to the classical interplanetary scintillation (IPS) study based on observations of the flux variability of distant natural radio sources. Results. We present here our technique and the results on IPS. We compare these with the total electron content (TEC) for the line of sight through the solar wind. Finally, we evaluate the applicability of the presented technique to phase-referencing Very Long Baseline Interferometry (VLBI) and Doppler observations of currently operational and prospective space missions., Comment: 7 pages, 6 figures, A&A in press

Context. Magnetic clouds (MCs) are a subset of interplanetary coronal mass ejections (ICMEs). One property of MCs is the presence of a magnetic flux rope. Is the difference between ICMEs with and without MCs intrinsic or rather due to an observational bias? Aims. As the spacecraft has no relationship with the MC trajectory, the frequency distribution of MCs versus the spacecraft distance to the MCs’ axis is expected to be approximately flat. However, Lepping & Wu (2010, Ann. Geophys., 28, 1539) confirmed that it is a strongly decreasing function of the estimated impact parameter. Is a flux rope more frequently undetected for larger impact parameter? Methods. In order to answer the questions above, we explore the parameter space of flux rope models, especially the aspect ratio, boundary shape, and current distribution. The proposed models are analyzed as MCs by fitting a circular linear force-free field to the magnetic field computed along simulated crossings. Results. We find that the distribution of the twist within the flux rope and the non-detection due to too low field rotation angle or magnitude only weakly affect the expected frequency distribution of MCs versus impact parameter. However, the estimated impact parameter is increasingly biased to lower values as the flux rope cross section is more elongated orthogonally to the crossing trajectory. The observed distribution of MCs is a natural consequence of a flux rope cross section flattened on average by a factor 2 to 3 depending on the magnetic twist profile. However, the faster MCs at 1 AU, with V > 550 km s-1, present an almost uniform distribution of MCs vs. impact parameter, which is consistent with round-shaped flux ropes, in contrast with the slower ones. Conclusions. We conclude that the sampling of MCs at various distances from the axis does not significantly affect their detection. The large part of ICMEs without MCs could be due to a too strict criteria for MCs or to the fact that these ICMEs are encountered outside their flux rope or near the leg region, or they do not contain a flux rope. Fil: Démoulin, Pascal. Centre National de la Recherche Scientifique. Observatoire de Paris; Francia Fil: Dasso, Sergio Ricardo. Consejo Nacional de Investigaciónes Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Astronomía y Física del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Astronomía y Física del Espacio; Argentina Fil: Janvier, M.. Centre National de la Recherche Scientifique. Observatoire de Paris; Francia

The energy spectra of galactic cosmic rays carry fundamental information regarding their origin and propagation. These spectra, when measured near Earth, are significantly affected by the solar magnetic field. A comprehensive description of the cosmic radiation must therefore include the transport and modulation of cosmic rays inside the heliosphere. During the end of the last decade, the Sun underwent a peculiarly long quiet phase well suited to study modulation processes. In this paper we present proton spectra measured from 2006 July to 2009 December by PAMELA. The large collected statistics of protons allowed the time variation to be followed on a nearly monthly basis down to 400 MV. Data are compared with a state-of-the-art three-dimensional model of solar modulation. © 2013. The American Astronomical Society. All rights reserved..

The last solar minimum activity period, and the consequent minimum modulation conditions for cosmic rays, was unusual. The highest levels of galactic protons were recorded at Earth in late 2009 in contrast to expectations. Proton spectra observed for 2006 to 2009 from the PAMELA cosmic ray detector on-board the Resurs-DK1 satellite are presented together with the solutions of a comprehensive numerical model for the solar modulation of cosmic rays. The model is used to determine what mechanisms were mainly responsible for the modulation of protons during this period, and why the observed spectrum for 2009 was the highest ever recorded. From mid-2006 until December 2009 we find that the spectra became significantly softer because increasingly more low energy protons had reached Earth. To simulate this effect, the rigidity dependence of the diffusion coefficients had to decrease significantly below ˜ 3 GeV. The modulation minimum period of 2009 can thus be described as relatively more `diffusion dominated' than previous solar minima. However, we illustrate that drifts still had played a significant role but that the observable modulation effects were not as well correlated with the waviness of the heliospheric current sheet as before. Protons still experienced global gradient and curvature drifts as the heliospheric magnetic field had decreased significantly until the end of 2009, in contrast to the moderate decreases observed during previous minimum periods. We conclude that all modulation processes contributed to the observed increases in the proton spectra for this period, exhibiting an intriguing interplay of these major mechanisms.

Coronal mass ejections (CMEs) are routinely tracked with imagers in the interplanetary space while magnetic clouds (MCs) properties are measured locally by spacecraft. However, both imager and insitu data do not provide direct estimation on the global flux rope properties. The main aim of this study is to constrain the global shape of the flux rope axis from local measurements, and to compare the results from in-situ data with imager observations. We perform a statistical analysis of the set of MCs observed by WIND spacecraft over 15 years in the vicinity of Earth. With the hypothesis of having a sample of MCs with a uniform distribution of spacecraft crossing along their axis, we show that a mean axis shape can be derived from the distribution of the axis orientation. In complement, while heliospheric imagers do not typically observe MCs but only their sheath region, we analyze one event where the flux-rope axis can be estimated from the STEREO imagers. From the analysis of a set of theoretical models, we show that the distribution of the local axis orientation is strongly affected by the global axis shape. Next, we derive the mean axis shape from the integration of the observed orientation distribution. This shape is robust as it is mostly determined from the global shape of the distribution. Moreover, we find no dependence on the flux-rope inclination on the ecliptic. Finally, the derived shape is fully consistent with the one derived from heliospheric imager observations of the June 2008 event. We have derived a mean shape of MC axis which only depends on one free parameter, the angular separation of the legs (as viewed from the Sun). This mean shape can be used in various contexts such as the study of high energy particles or space weather forecast., Comment: 13 pages, 12 figures

Drift along the wavy heliospheric neutral sheet is believed to play an important role in cosmic-ray modulation and can explain the peaked versus flat intensity profiles during consecutive solar magnetic epochs. Modulation models are becoming more and more realistic and in order to determine the role of the wavy neutral sheet more accurately, we revisit a previous calculation for drift along it. While mathematically correct, we argue that the previous expression for neutral sheet drift, which follows naturally from the standard expression for gradient and curvature drift, must be adapted in order for the drift speed to be less than particle speed. We compare the effect of both the previous and the current more accurate version of neutral sheet drift on cosmic-ray modulation with results obtained by other methods. http://dx.doi.org/10.1088/0004-637X/760/1/60 http://iopscience.iop.org/0004-637X/760/1/60/pdf/0004-637X_760_1_60.pdf

The spatial diffusion of cosmic rays in turbulent magnetic fields can, in the most general case, be fully anisotropic, i.e. one has to distinguish three diffusion axes in a local, field-aligned frame. We reexamine the transformation for the diffusion tensor from this local to a global frame, in which the Parker transport equation for energetic particles is usually formulated and solved. Particularly, we generalize the transformation formulas to allow for an explicit choice of two principal local perpendicular diffusion axes. This generalization includes the 'traditional' diffusion tensor in the special case of isotropic perpendicular diffusion. For the local frame, we motivate the choice of the Frenet-Serret trihedron which is related to the intrinsic magnetic field geometry. We directly compare the old and the new tensor elements for two heliospheric magnetic field configurations, namely the hybrid Fisk and the Parker field. Subsequently, we examine the significance of the different formulations for the diffusion tensor in a standard 3D model for the modulation of galactic protons. For this we utilize a numerical code to evaluate a system of stochastic differential equations equivalent to the Parker transport equation and present the resulting modulated spectra. The computed differential fluxes based on the new tensor formulation deviate from those obtained with the 'traditional' one (only valid for isotropic perpendicular diffusion) by up to 60% for energies below a few hundred MeV depending on heliocentric distance., Comment: 8 pages, 6 figures, accepted in ApJ

Magnetic Clouds (MCs) are a particular subset of Interplanetary Coronal Mass Ejections (ICMEs), forming large scale magnetic flux ropes. In this work we analyze the evolution of a particular MC (observed on March 1998) using {\it in situ} observations made by two spacecraft approximately aligned with the Sun, the first one at 1 AU from the Sun and the second one at 5.4 AU. We study the MC expansion, its consequent decrease of magnetic field intensity and mass density, and the possible evolution of the so-called global ideal-MHD nvariants. We describe the magnetic configuration of the MC at both spacecraft using different models and compute relevant global quantities (magnetic fluxes, helicity and energy) at both helio-distances. We also track back this structure to the Sun, in order to find out its solar source. We find that the flux rope is significantly distorted at 5.4 AU. However, we are able to analyze the data before the flux rope center is over-passed and compare it with observations at 1 AU. From the observed decay of magnetic field and mass density, we quantify how anisotropic is the expansion, and the consequent deformation of the flux rope in favor of a cross section with an aspect ratio at 5.4 AU of $\approx 1.6$ (larger in the direction perpendicular to the radial direction from the Sun). We quantify the ideal-MHD invariants and magnetic energy at both locations, and find that invariants are almost conserved, while the magnetic energy decays as expected with the expansion rate found. The use of MHD invariants to link structures at the Sun and the interplanetary medium is supported by the results of this multispacecraft study. We also conclude that the local dimensionless expansion rate, that is computed from the velocity profile observed by a single spacecraft, is very accurate for predicting the evolution of flux ropes in the solar wind., 16 two-column pages, 8 figures. Accepted for publication in A&A

著者人数: 150名, Accepted: 2011-04-08, 資料番号: SA1002992000

The propagation of energetic charged particles in the heliospheric magnetic field is one of the fundamental problems in heliophysics. In particular, the structure of the heliospheric magnetic field remains an unsolved problem and is discussed as a controversial topic. The first successful analytic approach to the structure of the heliospheric magnetic field was the Parker field. However, the measurements of the Ulysses spacecraft at high latitudes revealed the possible need for refinements of the existing magnetic field model during solar minimum. Among other reasons, this led to the development of the Fisk field. This approach is highly debated and could not be ruled out with magnetic field measurements so far. A promising method to trace this magnetic field structure is to model the propagation of electrons in the energy range of a few MeV. Employing three-dimensional and time-dependent simulations of the propagation of energetic electrons, this work shows that the influence of a Fisk-type field on the particle transport in the heliosphere leads to characteristic variations of the electron intensities on the timescale of a solar rotation. For the first time it is shown that the Ulysses count rates of 2.5-7 MeV electrons contain the imprint of a Fisk-type heliospheric magnetic field structure. From a comparison of simulation results and the Ulysses count rates, realistic parameters for the Fisk theory are derived. Furthermore, these parameters are used to investigate the modeled relative amplitudes of protons and electrons, including the effects of drifts. http://dx.doi.org/10.1088/0004-637X/741/1/23 http://iopscience.iop.org/0004-637X/741/1/23/pdf/0004-637X_741_1_23.pdf

Two of the paradigms in modeling the transport of galactic cosmic rays are that the modulation boundary is the heliopause and that the local interstellar spectra are identical to the galactic cosmic ray spectra. Here we demonstrate that the proton spectrum is already modulated due to an altered interstellar diffusion in the outer heliosheath as a consequence of the heliospheric “obstacle” in the interstellar flow. The main modulation effect however is adiabatic energy losses during a “confinement time” of cosmic rays inside the heliosphere. http://dx.doi.org/10.1088/0004-637X/735/2/128 http://iopscience.iop.org/0004-637X/735/2/128/pdf/0004-637X_735_2_128.pdf

We present a newly developed numerical modulation model to study the transport of galactic and Jovian electrons in the heliosphere. The model employs stochastic differential equations (SDEs) to solve the corresponding transport equation in five dimensions (time, energy, and three spatial dimensions) which is difficult to accomplish with the numerical schemes used in finite difference models. Modeled energy spectra for galactic electrons are compared for the two drift cycles to observations at Earth. Energy spectra and radial intensity profiles of galactic and Jovian electrons are compared successfully to results from previous studies. In line with general drift considerations, it is found that most 100 MeV electrons observed at Earth enter the heliosphere near the equatorial regions in the A > 0 cycle, while they enter mainly over the polar regions in theA < 0 cycle. Our results indicate that 100 MeV electrons observed at Earth originate at the heliopause with ∼600 MeV undergoing adiabatic cooling during their transport to Earth. The mean propagation time of these particles varies between ∼180 and 300 days, depending on the drift cycle. For 10 MeV Jovian electrons observed at Earth, a mean propagation time of ∼40 days is obtained. During this time, the azimuthal position of the Jovian magnetosphere varies by ∼1◦. At a 50 AU observational point, the mean propagation time of these electrons increases to ∼370 days with an azimuthal position change of Jupiter of ∼20◦. The SDE approach is very effective in calculating these propagation times. http://dx.doi.org/10.1088/0004-637X/735/2/83 http://iopscience.iop.org/0004-637X/735/2/83/pdf/0004-637X_735_2_83.pdf

It is generally accepted that the solar wind is significantly heated beyond ∼10 AU by the turbulent decay of pickup ion generated Alfvenic fluctuations. Here, we present a detailed and general calculation of the pickup ion ionization frequencies, and we evaluate these quantities within the solar wind termination shock along the stagnation line. For this supersonic solar wind region, inside of the solar wind termination shock, our results compare well with earlier estimates of these frequencies. When, in the future, turbulence transport models are extended into the heliosheath, the methodology outlined in this paper can be used to calculate ionization frequencies in this hot and dense plasma region where the resulting calculations become more complex.

The time-dependent modulation of galactic cosmic rays in the heliosphere is studied over different polarity cycles by computing 2.5 GV proton intensities using a two-dimensional, time-dependent modulation model. By incorporating recent theoretical advances in the relevant transport parameters in the model, we showed in previous work that this approach gave realistic computed intensities over a solar cycle. New in this work is that a time dependence of the solar wind termination shock (TS) position is implemented in our model to study the effect of a dynamic inner heliosheath thickness (the region between the TS and heliopause) on the solar modulation of galactic cosmic rays. The study reveals that changes in the inner heliosheath thickness, arising from a time-dependent shock position, does affect cosmic-ray intensities everywhere in the heliosphere over a solar cycle, with the smallest effect in the innermost heliosphere. A time-dependent TS position causes a phase difference between the solar activity periods and the corresponding intensity periods. The maximum intensities in response to a solar minimum activity period are found to be dependent on the time-dependent TS profile. It is found that changing the width of the inner heliosheath with time over a solar cycle can shift the time of when the maximum or minimum cosmic-ray intensities occur at various distances throughout the heliosphere, but more significantly in the outer heliosphere. The time-dependent extent of the inner heliosheath, as affected by solar activity conditions, is thus an additional time-dependent factor to be considered in the long-term modulation of cosmic rays.