Search

In this work we analyze the ground-level enhancement of the cosmic-ray intensity due to solar energetic particles as observed on 29 September 1989, by using two pairs of standard and lead-free neutron monitors. This enables one to separate spectral and anisotropy effects. This has been done previously by several authors for other events, but in this paper we make use of the large size and long duration of this event, as well as the fact that it is perhaps the best-observed one in the whole data base since 1942. It is shown that the method is more sensitive than the standard method that uses neutron monitors at different locations. The analysis provides a prototype for what can potentially be achieved by a new generation of mini neutron monitors.

Neutron monitors are the premier instruments for precise measurements of time variations (e.g., of solar origin) in the galactic cosmic ray (GCR) flux in the range of ∼1–100 GeV. However, it has proven challenging to accurately determine the yield function (effective area) versus rigidity in order to relate a neutron monitor's count rate with those of other monitors worldwide and the underlying GCR spectrum. Monte Carlo simulations of the yield function have been developed, but there have been few opportunities to validate these observationally, especially regarding the particular environment surrounding each monitor. Here we have precisely measured the count rate of a calibration neutron monitor near the Princess Sirindhorn Neutron Monitor (PSNM) at Doi Inthanon, Thailand (18.59∘N, 98.49∘E, 2560 m altitude), which provides a basis for comparison with count rates of other neutron monitors worldwide that are similarly calibrated. We directly measured the effect of surrounding structure by operating the calibrator outside and inside the building. Using Monte Carlo simulations, we clarify differences in response of the calibrator and PSNM, as well as the calibrator outside and inside the building. The dependence of the calibrator count rate on surrounding structure can be attributed to its sensitivity to neutrons of 0.5–10 MeV and a shift of sensitivity to nucleons of higher energy when placed inside the building. Simulated calibrator to PSNM count rate ratios inside and outside agree with observations within a few percent, providing useful validation and improving confidence in our ability to model the yield function for a neutron monitor station.

We study the temporal intensity profile, or pulse shape, of cosmic ray ground level enhancements (GLEs) by calculating the rise $(\tau_\mathrm{r})$ and decay $(\tau_\mathrm{d})$ times for a small subset of all available events. Although these quantities show very large inter-event variability, a linear dependence of $\tau_\mathrm{d} \approx 3.5 \tau_\mathrm{r}$ is found. We interpret these observational findings in terms of an interplanetary transport model, thereby including the effects of scattering (in pitch-angle) as these particles propagate from (near) the Sun to Earth. It is shown that such a model can account for the observed trends in the pulse shape, illustrating that interplanetary transport must be taken into account when studying GLE events, especially their temporal profiles. Furthermore, depending on the model parameters, the pulse shape of GLEs may be determined entirely by interplanetary scattering, obscuring all information regarding the initial acceleration process, and hence making a classification between impulsive and gradual events, as is traditionally done, superfluous.

Transverse-integrated nodal diffusion methods currently represent the standard in full core neutronic simulation. The primary shortcoming of this approach is the utilization of the quadratic transverse leakage approximation. This approach, although proven to work well for typical LWR problems, is not consistent with the formulation of nodal methods and can cause accuracy and convergence problems. In this work, an improved, consistent quadratic leakage approximation is formulated, which derives from the class of higher-order nodal methods developed some years ago. Further, a number of iteration schemes are developed around this consistent quadratic leakage approximation which yields accurate node average results in much improved calculational times. The most promising of these iteration schemes results from utilizing the consistent leakage approximation as a correction method to the standard quadratic leakage approximation. Numerical results are demonstrated on a set of benchmark problems and further applied to a realistic reactor problem, particularly the SAFARI-1 reactor, operating at Necsa, South Africa. The final optimal solution strategy is packaged into a standalone module which may simply be coupled to existing nodal diffusion codes.

This contribution to the 100th commemoration of the discovery of cosmic rays (6–8 August, 2012 in Bad Saarow, Germany) is about observations of those cosmic rays that are sensitive to the structure and the dynamics of the heliosphere. This places them in the energy range of 107–1010 eV. For higher energies the heliosphere becomes transparent; below this energy range the particles become strictly locked into the solar wind. Rather than give a strict chronological development, the paper is divided into distinct topics. It starts with the Pioneer/Voyager missions to the outer edges of the heliosphere, because the most recent observations indicate that a distinct boundary of the heliosphere might have been reached at the time of the meeting. Thereafter, the Ulysses mission is described as a unique one because it is still the only spacecraft that has explored the heliosphere at very high latitudes. Next, anomalous cosmic rays, discovered in 1972–1974, constitute a separate component that is ideally suited to study the acceleration and transport of energetic particles in the heliosphere. At this point the history and development of ground-based observations is discussed, with its unique contribution to supply a stable, long-term record. The last topic is about solar energetic particles with energies up to ∼1010 eV http://www.journals.elsevier.com/astroparticle-physics/

In a recent paper McCracken et al. (J. Geophys. Res. 113:A12101, 2008) proposed that the Ground Level Enhancement (GLE) of 20 January 2005 may have been produced by more than one acceleration mechanism, with the first acceleration due to the solar flare and the second one due to the CME associated with that event. They also noted several other GLEs with similar multiple pulse structures. This paper systematically investigates all the GLEs of solar cycle 23, from GLE 55 on 6 November 1997 to GLE 70 on 13 December 2006, to study their morphology and pulse structure, and to determine whether the multiple structures that may be found in these events are qualitatively similar to that of the GLE of 20 January 2005. We use all the data of all NMs that saw each event, to have as much directional and spectral information as possible. It is shown that three of these 16 events do contain such double-pulse structures, and the properties of these three are discussed in some detail.

This paper is a follow-on of that of Kruger et al. [Kruger, H., Moraal, H., Bieber, J.W., Clem, J.M., Evenson, P.A., Pyle, K.R., Duldig, M.L., Humble, J.E. A calibration neutron monitor: energy response and instrumental temperature sensitivity. J. Geophys. Res. 113, A08101, doi:10.1029/2008JA013229, 2008], that describes the characteristics of a pair of calibration neutron monitors that were developed to intercalibrate the count rates of the world’s neutron monitors against each other. Such an intercalibration will allow the calculation of energy (rigidity) spectra, which will enhance the quality of the neutron monitor data. Kruger et al. (2008) investigated the energy and temperature response of the calibrators. This paper studies the statistical accuracy of the calibration procedure, its repeatability, and the sensitivity to its environment. The paper concludes with a calibration procedure that can minimise the uncertainties caused by these five effects, or at least correct for them.

Voyager 1 has observed strong increases in the intensities of 2-160 MeV electrons since crossing the termination shock of the heliosphere in 2004 December. Before this time these intensities were submerged below the detector background, except for occasional transient events. These increases are large compared to the concurrent increases of positive ions such as H, He, and O. A significant part is probably due to temporal effects as the heliosphere was recovering to solar minimum conditions from 2005 to early 2010. The intensity observed by Voyager 2 since its crossing of the shock in 2007 September is 5-10 times lower than that observed by Voyager 1, which is so low that the electron intensity may still be below the background produced by high-energy protons in the detector. This points to a large north-south asymmetry in the properties of the heliosheath. It is shown that the observations suggest that these electrons are not freshly accelerated on the termination shock, but rather that they are of galactic origin-while they may be re-accelerated by that shock. In this paper, these intensities are modeled with numerical solutions of the cosmic-ray transport equation. It is shown that because they are relativistic, the electrons are muchmore » more sensitive to the form of the diffusion coefficient at low rigidities than ions, and that this can explain the asymmetry.« less

On January 20, 2005 there was an X 7.1 solar flare at 0636 UT with an accompanied halo coronal mass ejection (CME). The resultant interplanetary shock impacted earth ∼36 h later. Near earth, the Advanced Composition Explorer (ACE) spacecraft observed two impulses with a staircase structure in density and pressure. The estimated earth-arrival times of these impulses were 1713 UT and 1845 UT on January 21, 2005. Three MINIature Spectrometer (MINIS) balloons were aloft on January 21st; one in the northern polar stratosphere and two in the southern polar stratosphere. MeV relativistic electron precipitation (REP) observed by all three balloons is coincident (

Recent accurate measurements of the cosmic-ray intensity ratio ^(51)V/^(51)Cr below 1 GeV nucleon^(-1) provide a powerful new tool to study cosmic-ray modulation in the heliosphere. This paper describes how energy changes during this modulation process influence this ratio. In particular, our model includes acceleration at the solar wind termination shock, and we find that this mechanism significantly enhances the ^(51)V/^(51)Cr ratio at 1 AU. It is also shown that this acceleration makes the ratio more sensitive to the form of local low-energy interstellar spectra, below ~100 MeV nucleon^(-1), than without it. Specifically, this acceleration provides an independent confirmation of the consensus that low-energy spectra should be flatter than their high-energy power-law forms.

During the solar minimum period of 1954 the cosmic-ray diurnal variation as observed by neutron monitors and muon telescopes underwent a dramatic swing in its direction of maximum intensity, from the normal value between 16:00 and 18:00 local time to as early as 08:00. It is shown that this swing can be explained as being due to a negative radial density gradient of cosmic rays in the inner heliosphere and that this negative gradient is caused by large radial and latitudinal diffusion mean free paths that bring in particles from high latitudes. In principle, such large diffusion mean free paths should simultaneously cause high intensities, as were observed in 1954.

New insights regarding the form of the accelerated anomalous cosmic ray spectrum at the solar wind termination shock, obtained from solutions of the cosmic ray transport equation, are presented. A simple analytical expression for the spectrum on the shock is derived, and its dependence on the acceleration parameters is shown. Particular attention is paid to the high-energy cutoff of this spectrum. This expression and its parameterization should be applicable for the acceleration of charged particles in any spherical shock.

We review the observed properties of anomalous cosmic rays and the present status of our knowledge of the processes by which they originate. We compiled a comprehensive set of ACR energy spectral data from various spacecraft throughout the heliosphere during the passes of Ulysses over the poles of the Sun and present first results of a detailed modeling effort. In several contributions, we discuss the questions of injection and possible pre-acceleration of pickup ions, summarize new observations on the ionic charge composition, and present new results on the composition of minor ions in ACRs.

[1] Since the middle 1950s, neutron monitors have provided a continuous record of the intensity of secondary atmospheric particles produced by the primary cosmic radiation above the atmosphere. The number of counts due to these secondary particles is related to the primary spectrum above the atmosphere through the so-called atmospheric yield function. This yield function includes the secondary particles produced in the atmosphere, as well as inside the particular detector. In this paper the primary focus is to recalculate the yield function for neutron monitors. The motivation for this study is that the quality of the experimental observations has increased to such an extent that it has become possible to reduce uncertainties in the yield function down to approximately 10%. It thus becomes possible to refine our knowledge of the atmospheric cascade process considerably. We parameterize the yield functions in a simple way; we also make reference to the yield and response functions of muon, Cherenkov and stratospheric balloon detectors, and we compare their response to that of the neutron monitors.

The temporal variation of the cosmic-ray intensity in the heliosphere is called cosmic-ray modulation. The main periodicity is the response to the 11-year solar activity cycle. Other variations include a 27-day solar rotation variation, a diurnal variation, and irregular variations such as Forbush decreases. General awareness of the importance of this cosmic-ray modulation has greatly increased in the last two decades, mainly in communities studying cosmogenic nuclides, upper atmospheric physics and climate, helio-climatology, and space weather, where corrections need to be made for these modulation effects. Parameterized descriptions of the modulation are even used in archeology and in planning the flight paths of commercial passenger jets. The qualitative, physical part of the modulation is generally well-understood in these communities. The mathematical formalism that is most often used to quantify it is the so-called Force-Field approach, but the origins of this approach are somewhat obscure and it is not always used correct. This is mainly because the theory was developed over more than 40 years, and all its aspects are not collated in a single document. This paper contains a formal mathematical description intended for these wider communities. It consists of four parts: (1) a description of the relations between four indicators of “energy”, namely energy, speed, momentum and rigidity, (2) the various ways of how to count particles, (3) the description of particle motion with transport equations, and (4) the solution of such equations, and what these solutions mean. Part (4) was previously described in Caballero-Lopez and Moraal (J. Geophys. Res, 109: A05105, doi:10.1029/2003JA010358 2004). Therefore, the details are not all repeated here. The style of this paper is not to be rigorous. It rather tries to capture the relevant tools to do modulation studies, to show how seemingly unrelated results are, in fact, related to one another, and to point out the historical context of some of the results. The paper adds no new knowledge. The summary contains advice on how to use the theory most effectively. http://dx.doi.org/10.1007/s11214-011-9819-3

[1] The solar minimum of 2009 was characterized by a prolonged increase toward the maximum cosmic ray intensity, which was higher than it was during the maxima of 22 and 44 years ago. In the previous two so-called qA

[1] Worldwide observations of the cosmic ray ground level enhancement (GLE) of 20 January 2005 are used to investigate a commonly observed but poorly understood feature of this class of event. It is argued that the GLE comprised two distinctly different cosmic ray populations. The first resulted in an impulsive, highly anisotropic, field-aligned pulse with a relatively hard rigidity spectrum and significant velocity dispersion. The characteristics of the anisotropy were almost identical to those for similar impulsive increases observed during GLEs in 1960, 1978, and 1989. The π0γ ray observations from the RHESSI and CORONAS-F spacecraft and Type III radio emissions yield a path length of 1.76 ± 0.1 AU to Earth for the first pulse. After the highest energies in the initial anisotropic pulse had passed Earth, another field-aligned but mildly anisotropic cosmic ray pulse developed slowly worldwide, exhibiting the characteristics of the conventional GLE. The risetime and anisotropy of this second population indicate substantial scattering, apparently at variance to the essentially scatter-free nature of the initial pulse. We show that the coexisting scatter-free initial impulsive increase and the diffusive character of the second pulse are consistent with the standard quasi-linear theory of pitch angle diffusion. Throughout the GLE, the anisotropy remained field-aligned, and a third maximum, seen by some stations, is shown to be due to changes in the direction of the heliospheric magnetic field (HMF). Examination of 22 large (>20%) GLEs in the historical record shows that the impulsive pulse never occurs after the commencement of the P2 pulse, indicating that the impulsive-gradual combination is not due to a chance sampling of differing scattering regions of the HMF. It is further shown that impulsive pulses, or their equivalents, have been observed in 13 out of the 15 GLEs associated with solar activity in the solar longitude range 24°–98°W, leading us to propose that the event of 20 January 2005 should be regarded as the defining example of the GLE. The observations lead us to propose two separate acceleration episodes in the typical GLE: (1) acceleration directly associated with the flare itself and located in the lower corona and (2) acceleration by a supercritical shock driven by the associated coronal mass ejection, located at ∼3–5 solar radii and farther in the upper corona. A one-to-one association with so-called impulsive and gradual solar energetic particle events at lower energies is proposed. On the basis of these observations, a generic model for the GLE is proposed.

[1] Neutron monitors are integral detectors of secondary cosmic rays. Since each of them has its own detection efficiency, energy spectra cannot readily be derived from their observations. To circumvent this problem, latitudinal surveys have been conducted for many years with mobile neutron monitors to derive such spectra. Another way, however, is to use the worldwide stationary neutron monitor network, but then the counting rates of these monitors must be normalized sufficiently accurately against one another. For this reason, two portable calibration neutron monitors were built at the Potchefstroom campus of the North-West University, South Africa. When calibrations of an adequate number of the worldwide neutron monitors have been done, this paper demonstrates that the differential response functions derived from them will provide experimental data for modulation studies in the rigidity range of 1 to 15 GV. Several properties of these calibrators were investigated, in order to achieve sufficient calibration accuracy. The energy response over the cutoff rigidity range from the poles to the equator, as well as the instrumental temperature sensitivity, are described in this paper. The main conclusion is that the calibrator has a difference of almost 4% in its energy response over the cutoff rigidity range 0–16 GV in comparison with a standard 3NM64 neutron monitor. Furthermore, it is shown that not only the calibrator, but also the NM64 and IGY monitors, have fairly large instrumental temperature sensitivities. Correction coefficients for these effects are given.

[1] During January, 2005, there were several large X-class solar flares and associated solar energetic particle (SEP) events. Coincidentally, the MINIS balloon campaign had multiple payloads aloft in the stratosphere above Antarctica measuring dc electric fields, conductivity and x-ray flux. One-to-one increases in the electrical conductivity and decreases to near zero of both the vertical and horizontal electric field components were observed in conjunction with an increase in particle flux at SEP onset. Combined with an atmospheric electric field mapping model, these data are consistent with a shorting out of the global electric circuit and point toward substantial ionospheric convection modifications. Additionally, two subsequent, rapid changes were detected in the vertical electric field component several hours after SEP onset. These changes result in similar fluctuations in the calculated vertical current density. We will describe how rigidity cut-off dynamics may be crucial in understanding these sudden jumps in the vertical electric field.

[1] Adiabatic energy losses form a well-understood part of cosmic ray modulation in the heliosphere. However, shock acceleration at the solar wind termination shock and the modulation in the heliosheath beyond this shock alter the net cosmic ray energy change from that in standard noshock modulation models. This paper compares cosmic ray energy changes in several versions of the modulation process, from the simplest Force-Field model to one with a shock and heliosheath included. Implications for observations and their interpretation are noted. In particular, the shock acceleration mechanism seems to provide a sensitive probe of interstellar spectra at low energies (

The ground-level enhancement (GLE) of 1989 September 29 is one of the largest of 71 solar energetic particle events observed by neutron monitors on Earth. It was smaller than the record-breaking GLE 5 of 1956 February 23, but by some measures it was larger than GLE 69 of 2005 January 20. It is also the most extensively studied of the 71 GLEs, and it was observed by more than 50 ground-based detectors in the worldwide network. This paper contains another study of the event, with the main difference from previous studies that all the existing observations are employed, instead of the usual selection of stations. An effort is made to represent all the information graphically. This reveals new insight in the event, mainly about its time profile. The main conclusion is that the event is the best example available of a "classical" GLE that has a gradual increase toward peak intensity and does not contain two or more distinct peaks as inferred previously. It does, however, suggest that there were two acceleration or release mechanisms: a prompt, rapid one and a delayed, slower one. This conclusion is based on a detailed comparison with GLE 69 of 2005 January 20, which is the best-known example of a double-peaked event with a "prompt" component. It is also found that the rigidity spectrum was probably softer than derived in several previous studies, and that the decay phase of the event reveals that the cosmic-ray diffusion coefficient in the neutron monitor range is proportional to rigidity http://iopscience.iop.org/0004-637X http://dx.doi.org/10.1088/0004-637X/790/2/154 South African NRF grant SNA2011110300007 and Mexican PAPIIT-UNAM grant IN110413.

Charged particles propagating through a structured magnetic field are subject to drift motion. The primary aim of the present paper is therefore to investigate the effects of gradient, curvature, and neutral sheet drift on the evolution of the electron and positron spectra in a pulsar wind nebula, where the drift motion is a direct result of the magnetic field having an Archimedean spiral structure. In order to investigate the evolution of the spectra, the steady-state, axisymmetric Fokker-Planck transport equation is solved numerically using a finite-difference scheme. Apart from drift motion, the transport processes of convection and diffusion, along with the energy loss processes of adiabatic cooling and synchrotron radiation, are also included in the model. It is found that drift, particularly neutral sheet drift, can lead to a quantitative difference in the evolution of the electron and positron spectra. This difference may be of importance when interpreting the positron excess observed by PAMELA and AMS-02 near Earth.

The world’s neutron monitor network was initiated about 50 years ago. It grew to a peak of almost 100 stations towards the end of the 1960s and at present about half of these are still active. Many of the original questions about the production of atmospheric secondaries, geomagnetic effects, and neutron monitor response in general, have been settled satisfactorily. Due to their long-term reliability and automated data acquisition, the remaining neutron monitors in the network are well suited for important future contributions in several areas. Amongst these are (a) spectral measurements, which require an optimal distribution along cutoff rigidities; (b) anisotropy studies, which require a set of neutron monitors with well-defined, narrow cones of acceptance for charged particles, covering all directions as evenly as possible; and (c) solar neutron measurements, which primarily require an even distribution in longitude, at high altitude and near the equator. Steps have already been taken to improve the network with these goals in mind, and to standardize methods. This contribution describes some of these steps and suggests further strategies to achieve the most optimal network.

The global processes that determine cosmic ray modulation are reviewed. The essential elements of the theory which describes cosmic ray behavior in the heliosphere are summarized, and a series of discussions is presented which compare the expectations of this theory with observations of the spatial and temporal behavior of both galactic cosmic rays and the anomalous component; the behavior of cosmic ray electrons and ions; and the 26-day variations in cosmic rays as a function of heliographic latitude. The general conclusion is that the current theory is essentially correct. There is clear evidence, in solar minimum conditions, that the cosmic rays and the anomalous component behave as is expected from theory, with strong effects of gradient and curvature drifts. There is strong evidence of considerable latitude transport of the cosmic rays, at all energies, but the mechanism by which this occurs is unclear. Despite the apparent success of the theory, there is no single choice for the parameters which describe cosmic ray behavior, which can account for all of the observed temporal and spatial variations, spectra, and electron vs. ion behavior.

A possible way to calculate particle spectra as a function of position in pulsar wind nebulae is to solve a Fokker–Planck transport equation. This paper presents numerical solutions to the transport equation with the processes of convection, diffusion, adiabatic losses, and synchrotron radiation included. In the first part of the paper, the steady-state version of the transport equation is solved as a function of position and energy. This is done to distinguish the various effects of the aforementioned processes on the solutions to the transport equation. The second part of the paper deals with a time-dependent solution to the transport equation, specifically taking into account the effect of a moving outer boundary. The paper highlights the fact that diffusion can play a significant role in reducing the amount of synchrotron losses, leading to a modification in the expected particle spectra. These modified spectra can explain the change in the photon index of the synchrotron emission as a function of position. The solutions presented in this paper are not limited to pulsar wind nebulae, but can be applied to any similar central source system, e.g., globular clusters. http://iopscience.iop.org/0004-637X/765/1/30/article http://iopscience.iop.org/0004-637X http://dx.doi.org/10.1088/0004-637X/765/1/30

We have studied short-lived (21 minute average duration), highly anisotropic pulses of cosmic rays that constitute the first phase of 10 large ground-level enhancements (GLEs), and which extend to rigidities in the range 5‐20 GV. We provide a set of constraints that must be met by any putative acceleration mechanism for this type of solar-energetic-particle (SEP) event. The pulses usually have very short rise-times (three to five minutes) at all rigidities, and exhibit the remarkable feature that the intensity drops precipitously by 50% to 70% from the maximum within another three to five minutes. Both the rising and falling phases exhibit velocity dispersion, which indicates that there are particles with rigidities in the range 1 90 MeV gamma-ray bursts, indicating that freshly accelerated SEPs had impinged on higher-density matter in the chromosphere prior to the departure of the SEP pulse for Earth. This study was based on an updated archive of the 71 GLEs in the historic record, which is now available for public use.

A latitude survey of the cosmic ray intensity at sea level was conducted during the 1986/1987 solar minimum period on commercial vessels of the South African Marine Corporation (SAFMARINE). The results show that the differential response function for the 1986/1987 solar minimum agrees well with that measured in 1965. Both these response functions are significantly lower than those for 1976 and 1954. This result supports the 22-year modulation cycle as predicted, for example, by models including drift effects of the charged cosmic ray particles in the large-scale interplanetary magnetic field. A cross-over of the spectra at rigidities of about 7 GV was also observed. Such a crossover is necessary to explain both the stationary neutron monitor counting rates and the lower-energy balloon and space observations in consecutive solar cycles.