Solar energetic particle access to distant longitudes through turbulent fieldline meandering
^{1}
Jeremiah Horrocks Institute, University of Central Lancashire,
Preston,
UK
email:
tlmlaitinen@uclan.ac.uk
^{2}
Université Libre de Bruxelles, Service de Physique Statistique et des Plasmas, CP
231, 1050
Brussels,
Belgium
^{3}
Department of Mathematics, University of Waikato,
Hamilton, New
Zealand
^{4}
Department of Physics and KIPAC, Stanford
University, Stanford,
CA
94305,
USA
Received: 21 November 2015
Accepted: 23 March 2016
Context. Current solar energetic particle (SEP) propagation models describe the effects of interplanetary plasma turbulence on SEPs as diffusion, using a FokkerPlanck (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 crossfield diffusion.
Aims. We study whether the recently suggested early nondiffusive 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 nondiffusive propagation mode has been shown to dominate the SEP crossfield 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 fieldline random walk in the interplanetary medium and requires an SEP transport model that properly describes the nondiffusive early phase of SEP crossfield propagation.
Key words: Sun: particle emission / turbulence / Sun: heliosphere / magnetic fields / diffusion
© ESO, 2016
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0),
which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Solar energetic particles (SEPs) are accelerated up to relativistic energies during solar eruptions. Their propagation in interplanetary space is controlled by the largescale Parker spiral magnetic field. SEPs can also propagate across the mean field as a result of largescale drifts (Marsh et al. 2013) and the turbulent magnetic field fluctuations superposed on the mean field.
SEP propagation in interplanetary space is typically modelled by solving a FokkerPlanck equation that describes the effect of the plasma turbulence on the particles as diffusion (Parker 1965; Jokipii 1966). The fieldaligned propagation is affected by smallscale inhomogeneities (Jokipii 1966), whereas the crossfield propagation is described as diffusion caused by randomwalking field lines (Jokipii 1966; Matthaeus et al. 2003; Shalchi 2010). These theories give the ratio of the crossfield and fieldaligned diffusion coefficients κ_{⊥}/κ_{∥}~0.01 near Earth, and these values are supported by cosmic ray observations (Burger et al. 2000; Potgieter et al. 2014) and particle simulations (Giacalone & Jokipii 1999). Recently, several works have studied the effect of crossfield diffusion on SEP evolution (Zhang et al. 2009; Dröge et al. 2010; He et al. 2011; Giacalone & Jokipii 2012; Qin et al. 2013) and showed up to 360° heliolongitudinal extents for wide source regions with onsets within days of the SEP injection.
Recent multispacecraft SEP observations, made with the nearEarth spacecraft SOHO and ACE and the two STEREO spacecraft, show that SEPs from a single solar eruption have fast access to a very wide range of longitudes, even 180° from the eruption location (e.g. Dresing et al. 2012; Richardson et al. 2014). Similar event extents as measured by longitudinal peak intensity distribution width have been observed both in gradual and impulsive events (Wiedenbeck et al. 2013; Cohen et al. 2014; Richardson et al. 2014), which in traditional view have very different acceleration origins (e.g. Reames 1999). Detailed analysis of two SEP multispacecraftobserved events, with SEP transport fitted with a 3D FP transport model, favour a small source region of about 20° for the SEPs and a high diffusion coefficient ratio κ_{⊥}/κ_{∥} ~ 0.1 (Dresing et al. 2012; Dröge et al. 2014). High ratios in the range of 0.1−1 have also been inferred in other studies of heliospheric particles (e.g. Zhang et al. 2003; Dwyer et al. 1997).
High values of κ_{⊥}/κ_{∥} are not supported by our current theoretical understanding of SEP transport in plasma turbulence, however. At 1 AU, the amplitude of transverse magnetic fluctuations, normalised with the background magnetic field, δB_{⊥}/B, is observed to be between 0.1 and 0.5 (e.g. Burlaga & Turner 1976), for which fullorbit particle simulations suggest κ_{⊥}/κ_{∥}< 0.01Giacalone & Jokipii (1999). Current theoretical understanding also supports much lower values for the ratio (e.g. Pei et al. 2010).
In this paper, we introduce a new model for heliospheric SEP propagation that is capable of reproducing the fast transport in longitude that is observed at the beginning of SEP events. Our work is based on the notion that SEPs spread initially much faster than allowed by a diffusion description, and that this is a result of fieldline meandering (Laitinen et al. 2013). The crossfield propagation of SEPs along meandering fields has been analysed previously to study SEP intensity dropouts (Giacalone et al. 2000) and to determine diffusion coefficients (Kelly et al. 2012). Laitinen et al. (2013) showed for a uniform magnetic field that particle propagation across the mean field is initially a fast nondiffusive process along the meandering fieldlines that only later relaxes to diffusion. Our new approach combines fieldline meandering and timeasymptotic diffusion in a new description of earlytime SEP propagation in the Parker field, the FP+FLRW model. We compare the new model to the traditional diffusion approach (the FP model).
2. Models
2.1. FP model
In the traditional approach, SEP propagation is solved using the FokkerPlanck equation for the distribution function f (Roelof 1969; Skilling 1971; Isenberg 1997; Zhang et al. 2009)
(1)where v, p and μ are the particle speed, momentum, and pitchangle cosine, V_{sw} and b the solar wind velocity and a unit vector along the local mean magnetic field, respectively, and L = −B/ (∂B/∂s) the focusing length of the particles, with s the arc length along the field line, and Q is the injection function. The terms of the order (dV_{sw}/ dt) /v (Isenberg 1997) are typically neglected in SEP transport studies as small compared to the other terms (e.g. Zhang et al. 2009; Dröge et al. 2010), and they vanish for the constant radial solar wind considered in our study (see Sect. 2.3).
The pitchangle diffusion coefficient D_{μμ} and the crossfield diffusion coefficient κ_{⊥} are given as (Jokipii 1966; Matthaeus et al. 2003)
where S(k) is the turbulence power spectrum, S_{∥} the spectrum of slab turbulence (see Sect. 2.3), Ω and r_{L} are the particle Larmor frequency and radius, respectively, the free parameter (Matthaeus et al. 2003), and the parallel mean free path (4)The fieldaligned spatial diffusion tensor is given as (5)that is, without parallel diffusion, because this is accounted for with D_{μμ}.
The changes in SEP energy during propagation in the interplanetary space are given by the terms in the third line of Eq. (1). They include momentum changes that are due to the diverging solar wind, and in an approximate way, the effect of drifts on the particle momentum (Dalla et al. 2015). However, the equation does not contain the drifts in its convective term (second term in the first line), which means that the transport picture presented by Eq. (1) is incomplete. The effect of latitudinal drifts on the particle energy was recently discussed by Dalla et al. (2015), who found that in 100 hours a 10 MeV proton can lose almost half of its energy through drifting along the V_{sw} × B electric field and adiabatic deceleration. An improved approach, the drift kinetic theory, introduces drifts within the convective term (e.g. le Roux & Webb 2009; Webb et al. 2009).
The energy changes that are due to drifts and the diverging solar wind may be significant when performing a detailed fitting of the SEP event onset and the decay phase (e.g. Ruffolo 1995; Dalla et al. 2015). However, we here consider neither, but only the first ten hours of the event. For simplicity, the energy change term is therefore neglected in this study. Because the spatial transport associated with drift is small for a 10 MeV proton within the early phase of an SEP event (Dalla et al. 2015), drifts are neglected within the convective term.
Particle diffusion parallel and across the mean field direction, as understood by current theories, are not independent of each other. The crossfield diffusion coefficient in Eq. (3) depends on the parallel diffusion coefficient, resulting in compound diffusion (e.g. Kóta & Jokipii 2000; Qin et al. 2002; Matthaeus et al. 2003). The pitch angle dependence of the crossfield diffusion coefficient has also been discussed recently, with suggestions for using κ_{⊥} ∝  μ  (e.g. Qin & Shalchi 2014) or κ_{⊥} ∝ r_{L}, the particle Larmor radius (Dröge et al. 2010). While the suggested different forms of pitch angle dependence of κ_{⊥} influence particle crossfield propagation (e.g. Strauss & Fichtner 2015), it has not been clearly established which type of dependence best represents particle behaviour in turbulent fields. We therefore use the conventional pitchangle independent form of κ_{⊥} here.
We solved Eq. (1) using stochastic differential equations (SDEs, e.g. Gardiner 2009), which solve the equations using pseudoparticles. The method is briefly described in Appendix A, further details are given in Kopp et al. (2012). We used enhanced pitch angle scattering across μ = 0 by using a pitch angle diffusion coefficient of the form (6)suggested by Beeck & Wibberenz (1986), where κ_{0} contains the μindependent terms in Eq. (2), q is the spectral index of the slab turbulence, and H = 0.1 is a parameter that enhances scattering between pitch angle hemispheres. The FP equation is solved in a Parker spiral magnetic field with the magnitude B given as (7)where B_{0} = 5 nT is the magnetic field at heliocentric distance r_{0} = 1 AU, and a = V_{sw}/ (Ω_{⊙}sinθ), where V_{sw} = 400 km s^{1}, Ω_{⊙} = 2.8631 × 10^{6} rad/s is the solar rotation rate and θ the colatitude.
2.2. FP+FLRW model
Laitinen et al. (2013) demonstrated that the earlytime particle transport across the mean magnetic field is dominated by the particles propagating along meandering field lines. With respect to the mean magnetic field direction, the earlytime crossfield propagation of the particles is therefore only weakly stochastic. At later times, Laitinen et al. (2013) showed that the crossfield extent of the particle population begins to widen, resulting in an asymptotic diffusive behaviour.
In this study, we used an FP+FLRW model that explicitly introduces the random walk of field lines across the Parker spiral magnetic field and particle propagation along these randomwalking field lines. A similar model was successfully used by Laitinen et al. (2013) for Cartesian geometry. In addition to propagating particles along meandering field lines, the particles are diffused across them with the timeasymptotic crossfield particle diffusion coefficient κ_{⊥} to facilitate the latetime widening of the particle population. Laitinen et al. (2013) showed that such a model reproduces both the initial wide extent of the SEP event (where the FLRW dominates) and the asymptotic diffusive behaviour at later times (where the particle diffusion dominates) well and that is significantly better than the FLRW or FP on their own.
In our model, the field line random walk is facilitated by the turbulent interplanetary magnetic field, which is described by a power spectrum. A mix of 2D and slab modes in turbulent magnetic field gives rise to diffusive spreading of field lines (Matthaeus et al. 1995). This allows us to describe the path of a particle following a meandering field line using SDEs, with (8)where dr_{⊥} is the displacement across the Parker field direction for advance dr_{∥} along the field, D_{FL} is the fieldline diffusion coefficient, and W_{⊥} a Gaussian random number with zero mean and unit variance. The randomwalking path is calculated for each simulated particle before the pseudoparticle trajectory is integrated. The pseudoparticles are propagated along and diffused from the meandering field lines according to Eq. (1), with crossfield diffusion given by Eq. (3), pitch angle diffusion as given in Eqs. (2) and (6), and focusing, where we used the Parker field geometry when calculating the focusing length.
The fieldline diffusion coefficient was obtained as presented in Matthaeus et al. (1995), for composite slab and 2D turbulence. The contribution from the slab component is proportional to S_{∥}(k_{∥} = 0), which in our turbulence model (see below) vanishes. We therefore only considered the 2D contribution to the fieldline diffusion, giving (Matthaeus et al. 1995) (9)It should be noted that particle propagation as a beam along field lines that randomwalk according to the fieldline diffusion coefficient in Eq. (9) causes transport across the mean field much faster than diffusion with the coefficient given by Eq. (3). The latter aims to describe the crossfield spreading of particles at late times, when the pitchangle distribution of the particles is isotropic. The two descriptions converge for a particle beam in turbulence with vanishing slab component, a_{NLGC} = 1 and a scatterfree velocity correlation in Eq. (3) of Matthaeus et al. (2003). The two descriptions of particle crossfield propagation are therefore consistent with each other in the appropriate limit.
Consistent with the SDE approach in the traditional model, the effect of the crossfield step length is not considered in the particle propagation time. This affects the accurate timing of particle propagation (e.g. Strauss & Fichtner 2014).
Fig. 1 Left: parallel and perpendicular scattering mean free paths for 10 MeV protons obtained with Eqs. (2)–(4) and λ_{⊥} = 3κ_{⊥}/v. Right: diffusion coefficient ratios. For the red curve, the crossfield diffusion coefficient is obtained as κ_{⊥} = vD_{FL}, representing a crossfield diffusion coefficient that is due to fieldline meandering of unscattered particles. 

Open with DEXTER 
2.3. Turbulence model
To compare the FP and FP+FLRW models, it is essential that D_{μμ}, κ_{⊥} and D_{FL}, as defined in Eqs. (2), (3) and (9), are consistent with each other. We ensured this by deriving them for a simple interplanetary turbulence model. We considered only wave modes parallel (slab) and perpendicular (2D) to the mean magnetic field (Gray et al. 1996) and defined the power spectrum as(10)where k_{∥} and k_{⊥} = k_{⊥} are magnitudes of the wave number. The spectra are given as broken power laws
where L_{∥} and L_{⊥} are the spectral breakpoint scales with a value of 0.007 AU used for both scales, consistent with solar wind observations (e.g. Tu & Marsch 1995, and references therein). The largest scale in the turbulence model, L_{0}(r), is taken to be r, a natural choice for scaling in spherically expanding solar wind. The turbulence amplitude δB^{2} is related to the normalisation factors C_{∥} and C_{⊥} to give . The turbulence amplitudes are normalised to give the parallel scattering mean free path λ_{∥} = 0.3 AU for a 10 MeV proton at 1 AU heliocentric distance (e.g. Palmer 1982), and (Bieber et al. 1996), giving δB^{2}/B^{2} = 0.04 at 1 AU, in line with interplanetary turbulence observations (e.g. Burlaga & Turner 1976).
The plasma turbulence evolves as the solar wind propagates from the Sun to the interplanetary medium (e.g. Tu & Marsch 1995). The amplitude of the fluctuations changes as a function of radial distance as the properties of the background plasma change. The turbulence also evolves nonlinearly, with energy cascading towards smaller scales. The radial and spectral evolution of turbulence in the heliosphere have been studied in several works over several decades. Many of the works, however, considered isotropic or slab turbulence (e.g. Tu et al. 1984; Tu 1987; Zhou & Matthaeus 1990; Vainio et al. 2003) instead of the more realistic anisotropic geometry, or the composite mode used in this and many other recent SEP transport studies. Some formulations exist for spectral transport separately to small parallel and perpendicular scales (e.g. Cranmer & van Ballegooijen 2005; Laitinen 2005). These works, however, have not been considered with respect to particle transport coefficients. The recent detailed heliospheric turbulence transport model by Zank et al. (2012) provides the evolution of the total turbulence energy and correlation lengths instead of the turbulence spectrum, and the application of their model to SEP transport modelling is a challenging task that has so far not been attempted.
Developing an SEP transport model that considers the spectral evolution of composite turbulence in the heliosphere is therefore an undertaking that has not been considered before, and it warrants a separate study. In this work, we take a simpler approach and only consider the radial evolution of the turbulence. We used the WKB description (Richter & Olbers 1974; Tu et al. 1984) for the slab and 2D components. For simplicity, we did not consider the evolution of the spectral shape nor the wave refraction, and we kept the wave components in their pure slab and 2D geometries. Furthermore, we did not consider the change in k that is due to the varying wave phase speed because assuming a constant inertialframe wave phase speed agrees well with observations (Vainio et al. 2003). With these simplifications, the turbulence evolves as (13)where v_{Ar} is the radial component of Alfvén velocity, V_{r} = V_{sw,r} + v_{Ar} and n_{e} the solar wind electron number density, with the subscript 0 denoting the values at reference distance r_{0}. We used a constant solar wind velocity V_{sw,r0} and electron density . Using these, we find (14)which is consistent with the observed approximately ∝r^{3} trend in interplanetary space (e.g. Bavassano et al. 1982). We here used V_{sw,r0} = 400 km s^{1} and v_{a,r0} = 30 km s^{1} to represent the values at r_{0} = 1 AU.
It should be noted that our model for the fieldline random walk differs from that of Giacalone (2001), where the fieldline random walk is caused by the motion of the magnetic field footpoints at the Sun that are due to supergranulation. The Giacalone (2001) model provides an explanation for the source of the fieldline random walk at the Sun, but it does not allow for turbulence evolution in interplanetary space, which further modifies the turbulence as the magnetic fields are convected with the solar wind (e.g. Bruno & Carbone 2005, and references therein). In our model, the fieldline meandering is derived from the turbulence spectrum and the corresponding field line diffusion coefficient, with no limiting assumptions on the interplanetary evolution of the turbulence.
2.4. Transport parameters
After specifying a turbulence model, we now calculate the particle transport parameters that are fully consistent with the properties of the turbulence for both the FP and FP+FLRW models. The transport parameters are shown in Fig. 1 for 10 MeV protons. In the left panel we show the particle transport parameters as parallel and crossfield mean free paths, which are the same for both the FP and the FP+FLRW model, and the right panel shows the ratio of the crossfield and parallel diffusion coefficients, κ_{⊥}/κ_{∥}. The values of the particle diffusion coefficients are broadly similar to those used in other studies (e.g. Zhang et al. 2009; Pei et al. 2010).
Fig. 2 Sample path (red curve) of a 10 MeV proton injected at the Sun (located at origin) with the FP (left panel) and FP+FLRW model (right panel), respectively. The black curve shows the Parker spiral field, while the dashed blue curve shows a meandering field line in the FP+FLRW model (right panel). 

Open with DEXTER 
At small heliocentric distances, our model gives high values for the parallel mean free path (left panel of Fig. 1), allowing the particles to propagate along the field lines essentially without scattering from the Sun to 0.4 AU. A similar transport regime is also present in models that use a constant radial diffusion coefficient to parametrise the parallel diffusion coefficient (e.g. Zhang et al. 2009; Dröge et al. 2010; He et al. 2011; Giacalone & Jokipii 2012; Qin et al. 2013).
In addition to κ_{⊥}/κ_{∥}, we show also the ratio for particles that propagate as a beam along the meandering field line, v·D_{FL}, in the right panel of Fig. 1. As discussed in Sect. 2.2, v·D_{FL} is always larger than κ_{⊥}. In addition, as shown by Fraschetti & Jokipii (2011) and Fraschetti & Giacalone (2012), for typical interplanetary turbulence parameters the rate of decoupling of particles from their field lines is always slower than the spreading of particles across the mean field. Thus, at early times the SEP propagation along randomwalking field lines dominates over both the asymptotic perpendicular diffusion and the particle decoupling from the field.
The crossfield mean free path for both particles and field lines is very small close to the Sun, of the order of 10^{5} AU. The value is deceptive, however: a low value close to the Sun can result in large longitudinal spread. In radial geometry, the fieldline diffusion coefficient can be written in terms of a longitudinal diffusion coefficient, D_{FL} = ^{⟨}Δx^{2}^{⟩}/ (2r) = r^{2}^{⟨}Δφ^{2}^{⟩}/ (2r) ≡ r^{2}D_{FLφ}. In a radial geometry, a constant longitudinal diffusion coefficient is therefore represented by a diffusion coefficient that is proportional to r^{2}. For our turbulence model, we have D_{FLφ} = (40°)^{2}/AU close to the Sun, and D_{FLφ} = (10°)^{2}/AU at 1 AU.
It should be noted that a simple parametrisation with κ_{⊥}/κ_{∥} a constant or ∝ r does not reproduce the behaviour of the diffusion coefficients obtained with a full model of the turbulence, such as the one used in this study and in Pei et al. (2010). The radial evolution of turbulence and its effect on the diffusion coefficients will be discussed in more detail in a separate study.
3. Results and discussion
We here study a simple injection profile, for E_{0} = 10 MeV protons and trace the propagation of pseudoparticles through interplanetary space according to the model equations. The resulting profile can also be convolved with different injection profiles to produce response for time and longitudeextended injections, which will be a subject of a future full parameter study.
The propagation of the particles in FP and FP+FLRW models is shown in Fig. 2 as a 2D projection of the 3D orbits. The red curves depict the path of a simulated pseudoparticle, and the black curves show the Parker spiral. In the right panel, the blue dashed curve depicts the diffusively meandering field line. In both models, the particles scatter along the field line and diffuse away from it in the crossfield direction. In the FP+FLRW model, however, the field line itself meanders away from the Parker Spiral field. As discussed above, the diffusion of the field line is faster than the diffusion of a particle across the field line, thus the particle orbit tends to follow the meandering field line (blue dashed curve in Fig. 2), and as a result particles spread faster across the mean Parker spiral field than in the FP case.
In Fig. 3 we compare the particle distributions predicted by the diffusion (FP) (left) and our new FP+FLRW (right) model, one (top row) and three (bottom row) hours after injection. The distribution contours are calculated for particles near the ecliptic plane, between latitudes ± 10°. In the FP model, particles diffusively spread in longitude as a function of time, and as they are focused rapidly away from the Sun in the diverging magnetic field, the most significant spreading of the particles takes place at distances >0.5 AU. In the FP+FLRW model, however, the particles have a much wider distribution in longitude and show significant spreading close to the Sun.
To compare the models with observed SEP events, we show in Fig. 4 the time evolution of the simulated intensities at three different longitudes for the two models, compared to an event that took place during November 3–4, 2011. GómezHerrero et al. (2015) have analysed this event using a 1D transport model, which, based on anisotropy observations, suggests a source region with 270° longitudinal extent. However, they found no supporting evidence for such a wide source and suggested a disturbed interplanetary magnetic field as one of the possible causes for the event extent.
Fig. 3 SEP particle density in arbitrary units for particles between latitudes −10° and 10°, one hour (top row) and three hours (bottom row) after an impulsive injection at 1 solar radius and longitude φ = 0°. Left panels show the result of the FP model, while the right panels show the FP+FLRW model. The black curve shows the Parker fieldline starting at φ = 0°. The field line reaches a radial distance of 1 AU (red circle) at φ = −62°. 

Open with DEXTER 
Fig. 4 SEP proton intensities as a function of time at 1 AU. The green dashed curve shows 610 MeV SEP proton intensities observed during a solar eruption on Nov. 3–4, 2011, at SOHO (top left panel), STEREO B (top right) and STEREO A (bottom left) spacecraft. The blue and red curves show the 10 MeV proton intensities simulated using the FP and FP+FLRW models, respectively, using an impulsive injection at flare maximum, t_{0} = 22:41 – 500 s UT on Nov. 3 (dashed vertical line) at the flare longitude. The spacecraft locations are shown in the bottom right panel by filled green circle, yellow triangle, and blue square for SOHO, STEREO B, and STEREO A, respectively. The blue arrow denotes the SEP injection longitude, the thick solid line and red circle are the same as in Fig. 3, and the dashed curves show the Parker spiral fieldlines connecting the spacecraft to the solar surface. 

Open with DEXTER 
The SEP event is related to a solar eruption at longitude 50° east of the STEREO B (STB) spacecraft, with an M4.7 flare maximum at t_{max} = 22:41 UT on November 3 at N20E62 from Earth’s point of view (Richardson et al. 2014; GómezHerrero et al. 2015). Using the longitude angle φ as defined in Fig. 3, where the particle source is at φ = 0°, STB was at longitude φ_{STB} = 50°, (see lower right panel of Fig. 4). In this coordinate system, the SEP source is connected along the Parker field line (solid black spiral curve in the lower right panel of Fig. 4) to longitude φ = −62° at radial distance of 1 AU (red circle). During the event, the STEREO A (STA) and SOHO spacecraft were located at longitudes φ_{STA} = −103° and φ_{SOHO} = 152°, respectively.
Observed intensities of 6–10 MeV protons are shown in Fig. 4 with dashed green curves. They were obtained with ERNE onboard SOHO (Torsti et al. 1995) (top left panel), LET onboard STB (Mewaldt et al. 2008) (top right), and the LET instrument onboard STA (bottom left). The blue and red curves in Fig. 4 show the intensities from the FP and FP+FLRW models, respectively, for the injection function Q with t_{0} = t_{max}−500 s, (subtracting 500 s of lighttravelling time from Sun to 1 AU). The injection intensity is fitted so that the highest intensities of the models match the observed SEP highest intensity at STA.
The SEP event was best observed with LET/STA, located 41° clockwise from the wellconnected longitude of φ = −62°. The 6–10 MeV proton intensity increased by 5 orders of magnitude within a few hours (bottom left panel of Fig. 4). The event was also observed by LET/STB and ERNE/SOHO at 112° and 146°, respectively, from the wellconnected fieldline. Even with such poor connection, SEP intensities began to rise within an hour of the observed SEP onset at STA, to rise by 3 orders of magnitude above the preevent background. It is therefore clear that the SEPs had rapid access to a wide range of longitudes during this solar event.
While we did not perform a full transport fitting in this report, our study implies that the rapid access of SEPs during the November 3–4, 2011 event to the whole inner heliosphere can be achieved with our FP+FLRW model even with a narrow source region. The FP+FLRWmodelled particles are observed within an hour at the longitudes corresponding to the November 3, 2011 spacecraft locations. The FP model can reproduce the rapid onset of the SEP event at the STA spacecraft, with a delay of half an hour compared to the FP+FLRW model. It is not capable of reproducing the SOHO or STB observations, however, because it fails to spread particles to these longitudes to intensity levels above the preevent background during the first ten hours of the event.
The inability of the FP model to spread particles to the longitudes of the SOHO and STB spacecraft within the first ten hours is shown in Fig. 5, where we show the evolution of the simulated SEP event as it would have been measured at 1 AU at different heliographic longitudes, with t_{0} = 0. For the FP model the intensity falls very sharply as a function of longitude, explaining the failure of the FP model in Fig. 4 to describe the November 2011 event without resorting to extremely wide source regions (GómezHerrero et al. 2015). As shown in Fig. 5, the SEPs spread quickly for the FP+FLRW model, and the SEP onset within a wide range of longitudes, between φ = −150° and 50°, is delayed only by an hour compared to the onset at the wellconnected longitude, φ = −62°. Such a rapid spreading of first particles is a key feature of observations: the earliest SEP onsets at 100° from the flare longitude are delayed only an hour compared to onsets at bestconnected longitude (Richardson et al. 2014).
Fig. 5 Contour plot of particle density in arbitrary units as a function of time and longitude at 1 AU radial distance from the Sun (red circle in Figs. 1 and 2). The top panel shows the FP+FLRW model result, while the bottom panel shows results for the FP model. 

Open with DEXTER 
The longitudinal extent of an SEP event is frequently described by means of a Gaussian fit, f(φ) = f_{0}exp {−(φ − φ_{0})^{2}/ (2σ^{2})}, to the peak SEP intensities observed by spacecraft at different longitudes. Using STEREO and nearEarth spacecraft, several studies have reported values of σ = 30°−50° to best represent the range of SEP event widths (Lario et al. 2006, 2013; Wiedenbeck et al. 2013; Dresing et al. 2014; Richardson et al. 2014). To compare our models to observations, we performed a Gaussian fit to the peak intensities of the simulated SEP events at 1 AU for the two models presented in this work. The longitudinal distributions and the corresponding fits are shown in Fig. 6 with solid and dashed curves, respectively. The grey area depicts the observational range of σ = 30°−50°. The FP+FLRW model (red curve) can reproduce the observed SEP event widths, with width σ = 34°, while the FP model (blue curve) can not; it has a very narrow SEP width, with only σ = 11°. This strong difference between the FP+FLRW and FP models is consistent with Laitinen et al. (2013), who found for Cartesian geometry that at early times the mean square width in fullorbit simulations exceeded the FP model result by an order of magnitude and that an FP+FLRW model reproduces the behaviour seen in fullorbit simulations well.
Fig. 6 Peak intensity at 1 AU as a function of longitude for the FP (blue solid curve) and FP+FLRW (red solid curve) models. The dashed lines show the fits of Gaussian profiles to the peak intensities, with σ = 11° for the FP model, and σ = 34° for the FP+FLRW model. The grey area depicts the range of Gaussian widths of 30−50° obtained from observations. 

Open with DEXTER 
Our results show the importance of modelling the early phase of an SEP event using a proper physical description. Both the FP and FP+FLRW models aim to describe the effect of the turbulence on the particles, including the crossfield spreading that is due to fieldline meandering. In the FP model, the fieldline meandering is included only timeasymptotically using a diffusion description. However, as shown by Laitinen et al. (2013), the initial spreading of particles is fast and nondiffusive, relaxing to diffusive behaviour only at later times. This is accounted for in our new FP+FLRW model and results in particles reaching much broader longitudinal extents than predicted by an FP model using the same turbulence parameters.
Propagation along meandering field lines implies that similar event extents should also be reached with other particle species. This is indeed found in observations: Lario et al. (2013) and Richardson et al. (2014) have found a similar range of peak intensity widths for electrons and protons at different energies, and Wiedenbeck et al. (2013) and Cohen et al. (2014) for heavy ions during impulsive events. While these observations do not exclude other SEPspreading mechanisms, such as transport in the corona, symphatetic flaring and extended particle sources (see, e.g., Wiedenbeck et al. 2013; GómezHerrero et al. 2015), the effect of interplanetary fieldline meandering offers a direct explanation for the SEP event widths that are observed to be independent of energy and particle species and the gradual or impulsive classification. The influence of the source height, size, and different interplanetary conditions on the resulting SEP event extent will be explored in a future full parametric study. We will also address the varying degrees of observed SEP anisotropies, which typically indicate crossfield transport in the interplanetary space (Dresing et al. 2014), but in some cases indicate rapid access of particles throughout the heliosphere (GómezHerrero et al. 2015).
4. Conclusions
We developed a new SEP transport model that takes the nondiffusive propagation of SEPs early in the event history into account for a Parker spiral geometry. We showed that the early onset of SEPs over a wide range of longitudes can be explained by fieldline random walk and requires an SEP transport model that properly describes the nondiffusive early phase of SEP crossfield propagation. Our FP+FLRW model is the first model that is capable of reproducing the observed fast access of SEPs to distant longitudes, when the particle and fieldline diffusion coefficients are consistently derived from an interplanetary turbulence model. When the FLRW is not included (in the FP model), a much narrower crossfield extent of the SEP event is produced. We conclude that introducing fieldline wandering into SEP modelling has the potential of resolving the problem of fast access of SEPs to a wide range of longitudes.
Acknowledgments
We acknowledge the STEREO/IMPACT team for providing the data used in this study. The SOHO/ERNE data were obtained from the Space Research Laboratory, University of Turku (http://www.srl.utu.fi/erne_data/). T.L. and S.D. acknowledge support from the UK Science and Technology Facilities Council (STFC) (grants ST/J001341/1 and ST/M00760X/1), M.M. from the European Commission FP7 Project COMESEP (263252), and F.E. from NASA grant NNX14AG03G. The contribution of A.K. benefitted from financial support through project He 3279/151, funded by the Deutsche Forschungsgmeinschaft (DFG), at the CAU Kiel, where large parts of this work where carried out. T.L., F.E. and S.D. acknowledge support from the International Space Science Institute as part of international team 297. We thank Horst Fichtner for helpful discussions. Access to the University of Central Lancashire’s High Performance Computing Facility is gratefully acknowledged.
References
 Bavassano, B., Dobrowolny, M., Mariani, F., & Ness, N. F. 1982, J. Geophys. Res., 87, 3617 [NASA ADS] [CrossRef] [Google Scholar]
 Beeck, J., & Wibberenz, G. 1986, ApJ, 311, 437 [NASA ADS] [CrossRef] [Google Scholar]
 Bieber, J. W., Wanner, W., & Matthaeus, W. H. 1996, J. Geophys. Res., 101, 2511 [NASA ADS] [CrossRef] [Google Scholar]
 Bruno, R., & Carbone, V. 2005, Liv. Rev. Sol. Phys., 2, 4 [Google Scholar]
 Burger, R. A., Potgieter, M. S., & Heber, B. 2000, J. Geophys. Res., 105, 27447 [NASA ADS] [CrossRef] [Google Scholar]
 Burlaga, L. F., & Turner, J. M. 1976, J. Geophys. Res., 81, 73 [NASA ADS] [CrossRef] [Google Scholar]
 Cohen, C. M. S., Mason, G. M., Mewaldt, R. A., & Wiedenbeck, M. E. 2014, ApJ, 793, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Cranmer, S. R., & van Ballegooijen, A. A. 2005, ApJS, 156, 265 [NASA ADS] [CrossRef] [Google Scholar]
 Dalla, S., Marsh, M. S., & Laitinen, T. 2015, ApJ, 808, 62 [NASA ADS] [CrossRef] [Google Scholar]
 Dresing, N., GómezHerrero, R., Klassen, A., et al. 2012, Sol. Phys., 281, 281 [NASA ADS] [Google Scholar]
 Dresing, N., GómezHerrero, R., Heber, B., et al. 2014, A&A, 567, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dröge, W., Kartavykh, Y. Y., Klecker, B., & Kovaltsov, G. A. 2010, ApJ, 709, 912 [NASA ADS] [CrossRef] [Google Scholar]
 Dröge, W., Kartavykh, Y. Y., Dresing, N., Heber, B., & Klassen, A. 2014, J. Geophys. Res., 119, 6074 [NASA ADS] [CrossRef] [Google Scholar]
 Dwyer, J. R., Mason, G. M., Mazur, J. E., et al. 1997, ApJ, 490, L115 [NASA ADS] [CrossRef] [Google Scholar]
 Effenberger, F., Fichtner, H., Scherer, K., et al. 2012, ApJ, 750, 108 [NASA ADS] [CrossRef] [Google Scholar]
 Fraschetti, F., & Giacalone, J. 2012, ApJ, 755, 114 [NASA ADS] [CrossRef] [Google Scholar]
 Fraschetti, F., & Jokipii, J. R. 2011, ApJ, 734, 83 [NASA ADS] [CrossRef] [Google Scholar]
 Gardiner, C. W. 2009, in Stochastic Methods, 4th edn. (Berlin, Heidelberg: SpringerVerlag), 13 [Google Scholar]
 Giacalone, J. 2001, J. Geophys. Res., 106, 15881 [NASA ADS] [CrossRef] [Google Scholar]
 Giacalone, J., & Jokipii, J. R. 1999, ApJ, 520, 204 [NASA ADS] [CrossRef] [Google Scholar]
 Giacalone, J., & Jokipii, J. R. 2012, ApJ, 751, L33 [NASA ADS] [CrossRef] [Google Scholar]
 Giacalone, J., Jokipii, J. R., & Mazur, J. E. 2000, ApJ, 532, L75 [NASA ADS] [CrossRef] [Google Scholar]
 GómezHerrero, R., Dresing, N., Klassen, A., et al. 2015, ApJ, 799, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Gray, P. C., Pontius, D. H., & Matthaeus, W. H. 1996, Geophys. Res. Lett., 23, 965 [NASA ADS] [CrossRef] [Google Scholar]
 He, H.Q., Qin, G., & Zhang, M. 2011, ApJ, 734, 74 [NASA ADS] [CrossRef] [Google Scholar]
 Isenberg, P. A. 1997, J. Geophys. Res., 102, 4719 [NASA ADS] [CrossRef] [Google Scholar]
 Jokipii, J. R. 1966, ApJ, 146, 480 [NASA ADS] [CrossRef] [Google Scholar]
 Kelly, J., Dalla, S., & Laitinen, T. 2012, ApJ, 750, 47 [NASA ADS] [CrossRef] [Google Scholar]
 Kopp, A., Büsching, I., Strauss, R. D., & Potgieter, M. S. 2012, Comput. Phys. Comm., 183, 530 [NASA ADS] [CrossRef] [Google Scholar]
 Kóta, J., & Jokipii, J. R. 2000, ApJ, 531, 1067 [NASA ADS] [CrossRef] [Google Scholar]
 Laitinen, T. 2005, Astrophys. Space Sci. Trans., 1, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Laitinen, T., Dalla, S., & Marsh, M. S. 2013, ApJ, 773, L29 [NASA ADS] [CrossRef] [Google Scholar]
 Lario, D., Kallenrode, M.B., Decker, R. B., et al. 2006, ApJ, 653, 1531 [NASA ADS] [CrossRef] [Google Scholar]
 Lario, D., Aran, A., GómezHerrero, R., et al. 2013, ApJ, 767, 41 [NASA ADS] [CrossRef] [Google Scholar]
 le Roux, J. A., & Webb, G. M. 2009, ApJ, 693, 534 [NASA ADS] [CrossRef] [Google Scholar]
 Marsh, M. S., Dalla, S., Kelly, J., & Laitinen, T. 2013, ApJ, 774, 4 [NASA ADS] [CrossRef] [Google Scholar]
 Matthaeus, W. H., Gray, P. C., Pontius, Jr., D. H., & Bieber, J. W. 1995, Phys. Rev. Lett., 75, 2136 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Matthaeus, W. H., Qin, G., Bieber, J. W., & Zank, G. P. 2003, ApJ, 590, L53 [NASA ADS] [CrossRef] [Google Scholar]
 Mewaldt, R. A., Cohen, C. M. S., Cook, W. R., et al. 2008, Space Sci. Rev., 136, 285 [NASA ADS] [CrossRef] [Google Scholar]
 Palmer, I. D. 1982, Rev. Geophys. Space Phys., 20, 335 [NASA ADS] [CrossRef] [Google Scholar]
 Parker, E. N. 1965, Planet. Space Sci., 13, 9 [NASA ADS] [CrossRef] [Google Scholar]
 Pei, C., Bieber, J. W., Burger, R. A., & Clem, J. 2010, J. Geophys. Res., 115, 12107 [Google Scholar]
 Potgieter, M. S., Vos, E. E., Boezio, M., et al. 2014, Sol. Phys., 289, 391 [NASA ADS] [CrossRef] [Google Scholar]
 Qin, G., & Shalchi, A. 2014, Appl. Phys. Res., 6, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Qin, G., Matthaeus, W. H., & Bieber, J. W. 2002, ApJ, 578, L117 [NASA ADS] [CrossRef] [Google Scholar]
 Qin, G., Wang, Y., Zhang, M., & Dalla, S. 2013, ApJ, 766, 74 [NASA ADS] [CrossRef] [Google Scholar]
 Reames, D. V. 1999, Space Sci. Rev., 90, 413 [NASA ADS] [CrossRef] [Google Scholar]
 Richardson, I. G., von Rosenvinge, T. T., Cane, H. V., et al. 2014, Sol. Phys., 289, 3059 [NASA ADS] [CrossRef] [Google Scholar]
 Richter, A. K., & Olbers, D. J. 1974, Astrophys. Space Sci., 26, 95 [NASA ADS] [CrossRef] [Google Scholar]
 Roelof, E. C. 1969, in Lectures in HighEnergy Astrophysics, eds. H. Ögelman, & J. R. Wayland, 111 [Google Scholar]
 Ruffolo, D. 1995, ApJ, 442, 861 [NASA ADS] [CrossRef] [Google Scholar]
 Shalchi, A. 2010, ApJ, 720, L127 [NASA ADS] [CrossRef] [Google Scholar]
 Skilling, J. 1971, ApJ, 170, 265 [NASA ADS] [CrossRef] [Google Scholar]
 Strauss, R. D., & Fichtner, H. 2014, A&A, 572, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Strauss, R. D., & Fichtner, H. 2015, ApJ, 801, 29 [NASA ADS] [CrossRef] [Google Scholar]
 Torsti, J., Valtonen, E., Lumme, M., et al. 1995, Sol. Phys., 162, 505 [NASA ADS] [CrossRef] [Google Scholar]
 Tu, C.Y. 1987, Sol. Phys., 109, 149 [NASA ADS] [CrossRef] [Google Scholar]
 Tu, C.Y., & Marsch, E. 1995, Space Sci. Rev., 73, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Tu, C.Y., Pu, Z.Y., & Wei, F.S. 1984, J. Geophys. Res., 89, 9695 [NASA ADS] [CrossRef] [Google Scholar]
 Vainio, R., Laitinen, T., & Fichtner, H. 2003, A&A, 407, 713 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Webb, G. M., Le Roux, J. A., & Zank, G. P. 2009, in AIP Conf. Ser. 1183, eds. X. Ao, & G. Z. R. Burrows, 85 [Google Scholar]
 Wiedenbeck, M. E., Mason, G. M., Cohen, C. M. S., et al. 2013, ApJ, 762, 54 [NASA ADS] [CrossRef] [Google Scholar]
 Zank, G. P., Dosch, A., Hunana, P., et al. 2012, ApJ, 745, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, M., Jokipii, J. R., & McKibben, R. B. 2003, ApJ, 595, 493 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, M., Qin, G., & Rassoul, H. 2009, ApJ, 692, 109 [NASA ADS] [CrossRef] [Google Scholar]
 Zhou, Y., & Matthaeus, W. H. 1990, J. Geophys. Res., 95, 10291 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Stochastic differential equations in the FP and FP+FLRW models
The method of stochastic differential equations (SDEs) is a tool for solving FokkerPlancktype transport equations. Instead of solving this equation directly, as is the case when applying finitedifference methods, for example, the SDEs trace socalled pseudoparticles, that is, phasespace elements, through the phasespace. These pseudoparticles obey the corresponding Langevin equation. The simplest example of a transport equation is (A.1)with the Langevin equation (A.2)where f(x,t) is the distribution function that depends on space x and time t. v and D denote velocity and diffusion coefficient, dt is the time step and and W_{t} a Gaussian random number with zero mean and unit variance. The solution of the entire equation is obtained by binning the results for an appropriate number of pseudoparticles. Although these phasespace elements represent the behaviour of real particles in most cases quite well, they
should not be confused with the latter. In contrast to the traditional approach, the SDE method also allows solving the transport equation, with a different Langevin equation, backwards in time. For details on these points and on the method in general, we refer to Kopp et al. (2012) and the literature cited there.
In the present case we solved the transport equation, Eq. (1), under the assumption of a pure radial and constant solar wind (V_{sw} = V_{sw}e_{r}) and a vanishing latitudinal component of the magnetic field, which in particular significantly reduces the Skilling term R (the term in curved brackets in the second line of Eq. (1)). We first solved the equation in a local Cartesian system with z being the direction along the magnetic field, and transformed the transport quantities back to the global system before the actual integration was carried out (Effenberger et al. 2012). In this local system, the Langevin equations read (A.3)Here, u_{∥} = V_{sw}b_{r}, where b_{r} is the radial component of the unit vector along the local magnetic field line and . The two local unit vectors perpendicular to the magnetic field are e_{x} = e_{θ} × e_{z} and e_{y} = e_{z} × e_{x}, and the corresponding particle crossfield diffusion coefficients κ_{⊥ ,x} and κ_{⊥ ,y}, respectively. For the FP model, is the unit vector along the local Parker field direction, whereas for the FP+FLRW model the pseudoparticle is first propagated (virtually) along the meandering field line, and itself is computed through the difference between the old and new points.
The transport coefficients in Eq. (A.3) are transformed into the global system by means of the tensor T_{αβ} = e_{α,β}, so that the (spatial) diffusion tensor becomes nondiagonal. In this case, the square root of a tensor has to be computed, see Kopp et al. (2012).
It is important to note that calculating the meandering field line using Eq. (8) is performed before the particle is injected into the meandering field line. The particle crossfield propagation, while diffusive relative to the meandering field line, is therefore not diffusive relative to the mean magnetic field direction until the crossfield distance that is due to integration of Eq. (A.3) is comparable to that due to the integration of Eq. (8). At early times, the perpendicular pseudoparticle propagation is therefore only weakly stochastic with respect to the mean magnetic field, making our model consistent with the results presented by Laitinen et al. (2013).
All Figures
Fig. 1 Left: parallel and perpendicular scattering mean free paths for 10 MeV protons obtained with Eqs. (2)–(4) and λ_{⊥} = 3κ_{⊥}/v. Right: diffusion coefficient ratios. For the red curve, the crossfield diffusion coefficient is obtained as κ_{⊥} = vD_{FL}, representing a crossfield diffusion coefficient that is due to fieldline meandering of unscattered particles. 

Open with DEXTER  
In the text 
Fig. 2 Sample path (red curve) of a 10 MeV proton injected at the Sun (located at origin) with the FP (left panel) and FP+FLRW model (right panel), respectively. The black curve shows the Parker spiral field, while the dashed blue curve shows a meandering field line in the FP+FLRW model (right panel). 

Open with DEXTER  
In the text 
Fig. 3 SEP particle density in arbitrary units for particles between latitudes −10° and 10°, one hour (top row) and three hours (bottom row) after an impulsive injection at 1 solar radius and longitude φ = 0°. Left panels show the result of the FP model, while the right panels show the FP+FLRW model. The black curve shows the Parker fieldline starting at φ = 0°. The field line reaches a radial distance of 1 AU (red circle) at φ = −62°. 

Open with DEXTER  
In the text 
Fig. 4 SEP proton intensities as a function of time at 1 AU. The green dashed curve shows 610 MeV SEP proton intensities observed during a solar eruption on Nov. 3–4, 2011, at SOHO (top left panel), STEREO B (top right) and STEREO A (bottom left) spacecraft. The blue and red curves show the 10 MeV proton intensities simulated using the FP and FP+FLRW models, respectively, using an impulsive injection at flare maximum, t_{0} = 22:41 – 500 s UT on Nov. 3 (dashed vertical line) at the flare longitude. The spacecraft locations are shown in the bottom right panel by filled green circle, yellow triangle, and blue square for SOHO, STEREO B, and STEREO A, respectively. The blue arrow denotes the SEP injection longitude, the thick solid line and red circle are the same as in Fig. 3, and the dashed curves show the Parker spiral fieldlines connecting the spacecraft to the solar surface. 

Open with DEXTER  
In the text 
Fig. 5 Contour plot of particle density in arbitrary units as a function of time and longitude at 1 AU radial distance from the Sun (red circle in Figs. 1 and 2). The top panel shows the FP+FLRW model result, while the bottom panel shows results for the FP model. 

Open with DEXTER  
In the text 
Fig. 6 Peak intensity at 1 AU as a function of longitude for the FP (blue solid curve) and FP+FLRW (red solid curve) models. The dashed lines show the fits of Gaussian profiles to the peak intensities, with σ = 11° for the FP model, and σ = 34° for the FP+FLRW model. The grey area depicts the range of Gaussian widths of 30−50° obtained from observations. 

Open with DEXTER  
In the text 