-
We measure the energy emitted by extensive air showers in the form of radio
emission in the frequency range from 30 to 80 MHz. Exploiting the accurate
energy scale of the Pierre Auger Observatory, we obtain a radiation energy of
15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV
arriving perpendicularly to a geomagnetic field of 0.24 G, scaling
quadratically with the cosmic-ray energy. A comparison with predictions from
state-of-the-art first-principle calculations shows agreement with our
measurement. The radiation energy provides direct access to the calorimetric
energy in the electromagnetic cascade of extensive air showers. Comparison with
our result thus allows the direct calibration of any cosmic-ray radio detector
against the well-established energy scale of the Pierre Auger Observatory.
-
The Auger Engineering Radio Array (AERA) is part of the Pierre Auger
Observatory and is used to detect the radio emission of cosmic-ray air showers.
These observations are compared to the data of the surface detector stations of
the Observatory, which provide well-calibrated information on the cosmic-ray
energies and arrival directions. The response of the radio stations in the 30
to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of
the incoming electric field. For the latter, the energy deposit per area is
determined from the radio pulses at each observer position and is interpolated
using a two-dimensional function that takes into account signal asymmetries due
to interference between the geomagnetic and charge-excess emission components.
The spatial integral over the signal distribution gives a direct measurement of
the energy transferred from the primary cosmic ray into radio emission in the
AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air
shower arriving perpendicularly to the geomagnetic field. This radiation energy
-- corrected for geometrical effects -- is used as a cosmic-ray energy
estimator. Performing an absolute energy calibration against the
surface-detector information, we observe that this radio-energy estimator
scales quadratically with the cosmic-ray energy as expected for coherent
emission. We find an energy resolution of the radio reconstruction of 22% for
the data set and 17% for a high-quality subset containing only events with at
least five radio stations with signal.
-
Neutrinos in the cosmic ray flux with energies near 1 EeV and above are
detectable with the Surface Detector array of the Pierre Auger Observatory. We
report here on searches through Auger data from 1 January 2004 until 20 June
2013. No neutrino candidates were found, yielding a limit to the diffuse flux
of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound
predictions. Neutrino identification is attempted using the broad
time-structure of the signals expected in the SD stations, and is efficiently
done for neutrinos of all flavors interacting in the atmosphere at large zenith
angles, as well as for "Earth-skimming" neutrino interactions in the case of
tau neutrinos. In this paper the searches for downward-going neutrinos in the
zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for
upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L.
single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an
$E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times
10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~
s^{-1}~ sr^{-1}}$.
-
A measurement of the cosmic-ray spectrum for energies exceeding
$4{\times}10^{18}$ eV is presented, which is based on the analysis of showers
with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger
Observatory between 1 January 2004 and 31 December 2013. The measured spectrum
confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$
eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with
index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by
a smooth suppression region. For the energy ($E_\text{s}$) at which the
spectral flux has fallen to one-half of its extrapolated value in the absence
of suppression, we find
$E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$
eV.
-
Interplanetary Coronal Mass Ejections are the manifestation of solar
transient eruptions, which can significantly modify the plasma and magnetic
conditions in the heliosphere. They are often preceded by a shock, and a
magnetic flux rope is detected in situ in a third to half of them. The main aim
of this study is to obtain the best quantitative shape for the flux rope axis
and for the shock surface from in situ data obtained during spacecraft
crossings of these structures. We first compare the orientation of the flux
ropes axes and shock normals obtained from independent data analyses of the
same events, observed in situ at 1AU from the Sun. Then, we carry out an
original statistical analysis of axes/shock normals by deriving the statistical
distributions of their orientations. We fit the observed distributions using
the distributions derived from several synthetic models describing these
shapes. We show that the distributions of axis/shock orientations are very
sensitive to their respective shape. One classical model, used to analyze
interplanetary imager data, is incompatible with the in situ data. Two other
models are introduced, for which the results for axis and shock normals lead to
very similar shapes; the fact that the data for MCs and shocks are independent
strengthen this result. The model which best fit all the data sets has an
ellipsoidal shape with similar aspect ratio values for all the data sets. These
derived shapes for the flux rope axis and shock surface have several potential
applications. First, these shapes can be used to construct a consistent ICME
model. Second, these generic shapes can be used to develop a quantitative model
to analyze imager data, as well as constraining the output of numerical
simulations of ICMEs. Finally, they will have implications for space weather
forecasting, in particular for forecasting the time arrival of ICMEs at the
Earth.
-
Flux ropes are twisted magnetic structures, which can be detected by in situ
measurements in the solar wind. However, different properties of detected flux
ropes suggest different types of flux-rope population. As such, are there
different populations of flux ropes? The answer is positive, and is the result
of the analysis of four lists of flux ropes, including magnetic clouds (MCs),
observed at 1 AU. The in situ data for the four lists have been fitted with the
same cylindrical force-free field model, which provides an estimation of the
local flux-rope parameters such as its radius and orientation. Since the
flux-rope distributions have a large dynamic range, we go beyond a simple
histogram analysis by developing a partition technique that uniformly
distributes the statistical fluctuations over the radius range. By doing so, we
find that small flux ropes with radius R<0.1 AU have a steep power-law
distribution in contrast to the larger flux ropes (identified as MCs), which
have a Gaussian-like distribution. Next, from four CME catalogs, we estimate
the expected flux-rope frequency per year at 1 AU. We find that the predicted
numbers are similar to the frequencies of MCs observed in situ. However, we
also find that small flux ropes are at least ten times too abundant to
correspond to CMEs, even to narrow ones. Investigating the different possible
scenarios for the origin of those small flux ropes, we conclude that these
twisted structures can be formed by blowout jets in the low corona or in
coronal streamers.
-
Contributions of the Pierre Auger Collaboration to the 33rd International
Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
-
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.
-
A large amount of magnetized plasma is frequently ejected from the Sun as
coronal mass ejections (CMEs). Some of these ejections are detected in the
solar wind as magnetic clouds (MCs) that have flux rope signatures. Magnetic
clouds are structures that typically expand in the inner heliosphere. We derive
the expansion properties of MCs in the outer heliosphere from one to five
astronomical units to compare them with those in the inner heliosphere. We
analyze MCs observed by the Ulysses spacecraft using insitu magnetic field and
plasma measurements. The MC boundaries are defined in the MC frame after
defining the MC axis with a minimum variance method applied only to the flux
rope structure. As in the inner heliosphere, a large fraction of the velocity
profile within MCs is close to a linear function of time. This is indicative
of} a self-similar expansion and a MC size that locally follows a power-law of
the solar distance with an exponent called zeta. We derive the value of zeta
from the insitu velocity data. We analyze separately the non-perturbed MCs
(cases showing a linear velocity profile almost for the full event), and
perturbed MCs (cases showing a strongly distorted velocity profile). We find
that non-perturbed MCs expand with a similar non-dimensional expansion rate
(zeta=1.05+-0.34), i.e. slightly faster than at the solar distance and in the
inner heliosphere (zeta=0.91+-0.23). The subset of perturbed MCs expands, as in
the inner heliosphere, at a significantly lower rate and with a larger
dispersion (zeta=0.28+-0.52) as expected from the temporal evolution found in
numerical simulations. This local measure of the expansion also agrees with the
distribution with distance of MC size,mean magnetic field, and plasma
parameters. The MCs interacting with a strong field region, e.g. another MC,
have the most variable expansion rate (ranging from compression to
over-expansion).
-
Observations of magnetic clouds (MCs) are consistent with the presence of
flux ropes detected in the solar wind (SW) a few days after their expulsion
from the Sun as coronal mass ejections (CMEs). Both the \textit{in situ}
observations of plasma velocity profiles and the increase of their size with
solar distance show that MCs are typically expanding structures. The aim of
this work is to derive the expansion properties of MCs in the inner heliosphere
from 0.3 to 1 AU.We analyze MCs observed by the two Helios spacecraft using
\textit{in situ} magnetic field and velocity measurements. We split the sample
in two subsets: those MCs with a velocity profile that is significantly
perturbed from the expected linear profile and those that are not. From the
slope of the \textit{in situ} measured bulk velocity along the Sun-Earth
direction, we compute an expansion speed with respect to the cloud center for
each of the analyzed MCs. We analyze how the expansion speed depends on the MC
size, the translation velocity, and the heliocentric distance, finding that all
MCs in the subset of non-perturbed MCs expand with almost the same
non-dimensional expansion rate ($\zeta$). We find departures from this general
rule for $\zeta$ only for perturbed MCs, and we interpret the departures as the
consequence of a local and strong SW perturbation by SW fast streams, affecting
the MC even inside its interior, in addition to the direct interaction region
between the SW and the MC. We also compute the dependence of the mean total SW
pressure on the solar distance and we confirm that the decrease of the total SW
pressure with distance is the main origin of the observed MC expansion rate. We
found that $\zeta$ was $0.91\pm 0.23$ for non-perturbed MCs while $\zeta$ was
$0.48\pm 0.79$ for perturbed MCs, the larger spread in the last ones being due
to the influence of the environment conditions on the expansion.
-
Relating in situ measurements of relativistic solar particles to their parent
activity in the corona requires understanding the magnetic structures that
guide them from their acceleration site to the Earth. Relativistic particle
events are observed at times of high solar activity, when transient magnetic
structures such as Interplanetary Coronal Mass Ejections (ICMEs) often shape
the interplanetary magnetic field (IMF). They may introduce interplanetary
paths that are longer than nominal, and magnetic connections rooted far from
the nominal Parker spiral. We present a detailed study of the IMF
configurations during ten relativistic solar particle events of the 23rd
activity cycle to elucidate the actual IMF configuration guiding the particles
to Earth, where they are measured by neutron monitors. We use magnetic field
(MAG) and plasma parameter measurements (SWEPAM) from ACE, and determine
interplanetary path lengths of energetic particles through a modified version
of the velocity dispersion analysis based on energetic particle measurements
with SoHO/ERNE. We find that the majority (7/10) of the events is detected in
the vicinity of an ICME. Their interplanetary path lengths are found to be
longer (1.5-2.6 AU) than those of the two events propagating in the slow solar
wind (1.3 AU). The largest apparent path length is found in an event within the
fast solar wind, probably due to enhanced pitch angle scattering. The derived
path lengths imply that the first energetic and relativistic protons are
released at the Sun at the same time as electron beams emitting type III radio
bursts. The timing of the first high-energy particle arrival at Earth is
dominantly determined by the type of IMF in which the particles propagate.
Initial arrival times are as expected from Parker's model in the slow solar
wind, and significantly larger in or near transient structures such as ICMEs.
We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.
The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.
Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.
A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.
Interplanetary Coronal Mass Ejections are the manifestation of solar transient eruptions, which can significantly modify the plasma and magnetic conditions in the heliosphere. They are often preceded by a shock, and a magnetic flux rope is detected in situ in a third to half of them. The main aim of this study is to obtain the best quantitative shape for the flux rope axis and for the shock surface from in situ data obtained during spacecraft crossings of these structures. We first compare the orientation of the flux ropes axes and shock normals obtained from independent data analyses of the same events, observed in situ at 1AU from the Sun. Then, we carry out an original statistical analysis of axes/shock normals by deriving the statistical distributions of their orientations. We fit the observed distributions using the distributions derived from several synthetic models describing these shapes. We show that the distributions of axis/shock orientations are very sensitive to their respective shape. One classical model, used to analyze interplanetary imager data, is incompatible with the in situ data. Two other models are introduced, for which the results for axis and shock normals lead to very similar shapes; the fact that the data for MCs and shocks are independent strengthen this result. The model which best fit all the data sets has an ellipsoidal shape with similar aspect ratio values for all the data sets. These derived shapes for the flux rope axis and shock surface have several potential applications. First, these shapes can be used to construct a consistent ICME model. Second, these generic shapes can be used to develop a quantitative model to analyze imager data, as well as constraining the output of numerical simulations of ICMEs. Finally, they will have implications for space weather forecasting, in particular for forecasting the time arrival of ICMEs at the Earth.
Flux ropes are twisted magnetic structures, which can be detected by in situ measurements in the solar wind. However, different properties of detected flux ropes suggest different types of flux-rope population. As such, are there different populations of flux ropes? The answer is positive, and is the result of the analysis of four lists of flux ropes, including magnetic clouds (MCs), observed at 1 AU. The in situ data for the four lists have been fitted with the same cylindrical force-free field model, which provides an estimation of the local flux-rope parameters such as its radius and orientation. Since the flux-rope distributions have a large dynamic range, we go beyond a simple histogram analysis by developing a partition technique that uniformly distributes the statistical fluctuations over the radius range. By doing so, we find that small flux ropes with radius R<0.1 AU have a steep power-law distribution in contrast to the larger flux ropes (identified as MCs), which have a Gaussian-like distribution. Next, from four CME catalogs, we estimate the expected flux-rope frequency per year at 1 AU. We find that the predicted numbers are similar to the frequencies of MCs observed in situ. However, we also find that small flux ropes are at least ten times too abundant to correspond to CMEs, even to narrow ones. Investigating the different possible scenarios for the origin of those small flux ropes, we conclude that these twisted structures can be formed by blowout jets in the low corona or in coronal streamers.
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.
A large amount of magnetized plasma is frequently ejected from the Sun as coronal mass ejections (CMEs). Some of these ejections are detected in the solar wind as magnetic clouds (MCs) that have flux rope signatures. Magnetic clouds are structures that typically expand in the inner heliosphere. We derive the expansion properties of MCs in the outer heliosphere from one to five astronomical units to compare them with those in the inner heliosphere. We analyze MCs observed by the Ulysses spacecraft using insitu magnetic field and plasma measurements. The MC boundaries are defined in the MC frame after defining the MC axis with a minimum variance method applied only to the flux rope structure. As in the inner heliosphere, a large fraction of the velocity profile within MCs is close to a linear function of time. This is indicative of} a self-similar expansion and a MC size that locally follows a power-law of the solar distance with an exponent called zeta. We derive the value of zeta from the insitu velocity data. We analyze separately the non-perturbed MCs (cases showing a linear velocity profile almost for the full event), and perturbed MCs (cases showing a strongly distorted velocity profile). We find that non-perturbed MCs expand with a similar non-dimensional expansion rate (zeta=1.05+-0.34), i.e. slightly faster than at the solar distance and in the inner heliosphere (zeta=0.91+-0.23). The subset of perturbed MCs expands, as in the inner heliosphere, at a significantly lower rate and with a larger dispersion (zeta=0.28+-0.52) as expected from the temporal evolution found in numerical simulations. This local measure of the expansion also agrees with the distribution with distance of MC size,mean magnetic field, and plasma parameters. The MCs interacting with a strong field region, e.g. another MC, have the most variable expansion rate (ranging from compression to over-expansion).
Observations of magnetic clouds (MCs) are consistent with the presence of flux ropes detected in the solar wind (SW) a few days after their expulsion from the Sun as coronal mass ejections (CMEs). Both the \textit{in situ} observations of plasma velocity profiles and the increase of their size with solar distance show that MCs are typically expanding structures. The aim of this work is to derive the expansion properties of MCs in the inner heliosphere from 0.3 to 1 AU.We analyze MCs observed by the two Helios spacecraft using \textit{in situ} magnetic field and velocity measurements. We split the sample in two subsets: those MCs with a velocity profile that is significantly perturbed from the expected linear profile and those that are not. From the slope of the \textit{in situ} measured bulk velocity along the Sun-Earth direction, we compute an expansion speed with respect to the cloud center for each of the analyzed MCs. We analyze how the expansion speed depends on the MC size, the translation velocity, and the heliocentric distance, finding that all MCs in the subset of non-perturbed MCs expand with almost the same non-dimensional expansion rate ($\zeta$). We find departures from this general rule for $\zeta$ only for perturbed MCs, and we interpret the departures as the consequence of a local and strong SW perturbation by SW fast streams, affecting the MC even inside its interior, in addition to the direct interaction region between the SW and the MC. We also compute the dependence of the mean total SW pressure on the solar distance and we confirm that the decrease of the total SW pressure with distance is the main origin of the observed MC expansion rate. We found that $\zeta$ was $0.91\pm 0.23$ for non-perturbed MCs while $\zeta$ was $0.48\pm 0.79$ for perturbed MCs, the larger spread in the last ones being due to the influence of the environment conditions on the expansion.
Relating in situ measurements of relativistic solar particles to their parent activity in the corona requires understanding the magnetic structures that guide them from their acceleration site to the Earth. Relativistic particle events are observed at times of high solar activity, when transient magnetic structures such as Interplanetary Coronal Mass Ejections (ICMEs) often shape the interplanetary magnetic field (IMF). They may introduce interplanetary paths that are longer than nominal, and magnetic connections rooted far from the nominal Parker spiral. We present a detailed study of the IMF configurations during ten relativistic solar particle events of the 23rd activity cycle to elucidate the actual IMF configuration guiding the particles to Earth, where they are measured by neutron monitors. We use magnetic field (MAG) and plasma parameter measurements (SWEPAM) from ACE, and determine interplanetary path lengths of energetic particles through a modified version of the velocity dispersion analysis based on energetic particle measurements with SoHO/ERNE. We find that the majority (7/10) of the events is detected in the vicinity of an ICME. Their interplanetary path lengths are found to be longer (1.5-2.6 AU) than those of the two events propagating in the slow solar wind (1.3 AU). The largest apparent path length is found in an event within the fast solar wind, probably due to enhanced pitch angle scattering. The derived path lengths imply that the first energetic and relativistic protons are released at the Sun at the same time as electron beams emitting type III radio bursts. The timing of the first high-energy particle arrival at Earth is dominantly determined by the type of IMF in which the particles propagate. Initial arrival times are as expected from Parker's model in the slow solar wind, and significantly larger in or near transient structures such as ICMEs.