Evaluation of cosmic ray rejection algorithms on single-shot exposures

To maximise data output from single-shot astronomical images, the rejection of cosmic rays is important. We present the results of a benchmark trial comparing various cosmic ray rejection algorithms. The procedures assess relative performances and characteristics of the processes in cosmic ray detection, rates of false detections of true objects, and the quality of image cleaning and reconstruction. The cosmic ray rejection algorithms developed by Rhoads (2000, PASP, 112, 703), van Dokkum (2001, PASP, 113, 1420), Pych (2004, PASP, 116, 148), and the IRAF task xzap by Dickinson are tested using both simulated and real data. It is found that detection efficiency is independent of the density of cosmic rays in an image, being more strongly affected by the density of real objects in the field. As expected, spurious detections and alterations to real data in the cleaning process are also significantly increased by high object densities. We find the Rhoads' linear filtering method to produce the best performance in the detection of cosmic ray events; however, the popular van Dokkum algorithm exhibits the highest overall performance in terms of detection and cleaning.

Measurement of the radiation energy in the radio signal of extensive air showers as a universal estimator of cosmic-ray energy.

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 +/- 0.7 (stat) +/- 6.7 (syst) 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-principles 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.

Investigating the cosmic-ray ionization rate in the galactic interstellar medium through observations of H3+

Observations of H3+ in the Galactic diffuse interstellar medium (ISM) have led to various surprising results, including the conclusion that the cosmic-ray ionization rate (zeta_2) is about 1 order of magnitude larger than previously thought. The present survey expands the sample of diffuse cloud sight lines with H3+ observations to 50, with detections in 21 of those. Ionization rates inferred from these detections are in the range (1.7+-1.0)x10^-16 s^-1 < zeta_2 < (10.6+-6.8)x10^-16 s^-1 with a mean value of zeta_2=(3.3+-0.4)x10^-16 s^-1. Upper limits (3sigma) derived from non-detections of H3+ are as low as zeta_2 < 0.4x10^-16 s^-1. These low upper-limits, in combination with the wide range of inferred cosmic-ray ionization rates, indicate variations in zeta_2 between different diffuse cloud sight lines. Calculations of the cosmic-ray ionization rate from theoretical cosmic-ray spectra require a large flux of low-energy (MeV) particles to reproduce values inferred from observations. Given the relatively short range of low-energy cosmic rays --- those most efficient at ionization --- the proximity of a cloud to a site of particle acceleration may set its ionization rate. Variations in zeta_2 are thus likely due to variations in the cosmic-ray spectrum at low energies resulting from the effects of particle propagation. To test this theory, H3+ was observed in sight lines passing through diffuse molecular clouds known to be interacting with the supernova remnant IC 443, a probable site of particle acceleration. Where H3+ is detected, ionization rates of zeta_2=(20+-10)x10^-16 s^-1 are inferred, higher than for any other diffuse cloud. These results support both the concept that supernova remnants act as particle accelerators, and the hypothesis that propagation effects are responsible for causing spatial variations in the cosmic-ray spectrum and ionization rate. Future observations of H3+ near other supernova remnants and in sight lines where complementary ionization tracers (OH+, H2O+, H3O+) have been observed will further our understanding of the subject.

Atmospheric effects on extensive air showers observed with the surface detector of the Pierre Auger observatory

Atmospheric parameters, Such as pressure (P), temperature (T) and density (rho proportional to P/T), affect the development of extensive air showers initiated by energetic cosmic rays. We have Studied the impact of atmospheric variations on extensive air showers by means of the surface detector of the Pierre Auger Observatory. The rate of events shows a similar to 10% seasonal modulation and similar to 2% diurnal one. We find that the observed behaviour is explained by a model including the effects associated with the variations of P and rho. The former affects the longitudinal development of air showers while the latter influences the Moliere radius and hence the lateral distribution of the shower particles. The model is validated with full simulations of extensive air showers using atmospheric profiles measured at the site of the Pierre Auger Observatory. (C) 2009 Elsevier B.V. All rights reserved.

Contribution from normal and starburst galaxies to the extragalactic gamma-ray background (EGRB)

The extragalactic diffuse emission at gamma-ray energies has interesting cosmological implications since these photons suffer little or no attenuation during their propagation from the site of origin. The emission could originate from either truly diffuse processes or from unresolved point sources such as AGNs, normal galaxies and starburst galaxies. Here, we examine the unresolved point source origin of the extragalactic gamma-ray background emission from normal galaxies and starburst galaxies. gamma-ray emission from normal galaxies is primarily coming from cosmic-ray interactions with interstellar matter and radiation (similar to 90%) along with a small contribution from discrete point sources (similar to 10%). Starburst galaxies are expected to have enhanced supernovae activity which leads to higher cosmic-ray densities, making starburst galaxies sufficiently luminous at gamma-ray energies to be detected by the current gamma-ray mission(Fermi Gamma-ray Space Telescope).

The relative abundances of Sn, Te, Xe, Ba, and Ce in the cosmic radiation

Elements with even atomic number (Z) in the interval 50 ≤ Z ≤ 58 have been resolved in the cosmic radiation using the Heavy Nuclei Experiment on the HEAO-3 satellite. Their relative abundances have been compared with the results expected from pure r-process material, pure s-process material, and solar system material, both with and without a modification due to possible first ionization potential effects. Such effects may be the result of the preferential acceleration, and hence enhancement in the cosmic rays, or those elements having low first ionization potentials. We find that our measurements are inconsistent with pure r-process material at the greater than 98% confidence level whether or not the first ionization potential adjustments are made.

In addition, we have compared our results with mixtures having varying ratios of pure r-process material to pure s-process material. We find that, if no first ionization potential effects are included,

(r/s)CRS/(r/s)SS = 0.20^+0.18_-0.14

where CRS refers to the cosmic ray source and SS refers to the solar system, consistent with having an almost pure s-process source. If the first ionization potential adjustments are applied

(r/s)CRS/(r/s)SS = 1.5^+1.1_-0.7

consistent with a solar system mixture.

The abundances of ultraheavy elements in the cosmic radiation

Analysis of the data from the Heavy Nuclei Experiment on the HEAO-3 spacecraft has yielded the cosmic ray abundances of odd-even element pairs with atomic number, Z, in the range 33 ≤ Z ≤60, and the abundances of broad element groups in the range 62 ≤ Z ≤83, relative to iron. These data show that the cosmic ray source composition in this charge range is quite similar to that of the solar system provided an allowance is made for a source fractionation based on first ionization potential. The observations are inconsistent with a source composition which is dominated by either r-process or s-process material, whether or not an allowance is made for first ionization potential. Although the observations do not exclude a source containing the same mixture of r- and s-process material as in the solar system. the data are best fit by a source having an r- to s-process ratio of 1.22^(+0.25)_(0.21), relative to the solar system The abundances of secondary elements are consistent with the leaky box model of galactic propagation, implying a pathlength distribution similar to that which explains the abundances of nuclei with Z<29.

The energy spectra of the even elements in the range 38 ≤ Z ≤ 60 are found to have a deficiency of particles in the range ~1.5 to 3 GeV/amu, compared to iron. This deficiency may result from ionization energy loss in the interstellar medium, and is not predicted by propagation models which ignore such losses. ln addition, the energy spectra of secondary elements are found to be different to those of the primary elements. Such effects are consistent with observations of lighter nuclei, and are in qualitative agreement with galactic propagation models using a rigidity dependent escape length. The energy spectra of secondaries arising from the platinum group are found to be much steeper than those of lower Z. This effect may result from energy dependent fragmentation cross sections.

A filamentation instability for streaming cosmic rays

We demonstrate that cosmic rays form filamentary structures in the precursors of supernova remnant shocks due to their self-generated magnetic fields. The cosmic ray filamentation results in the growth of a long-wavelength instability, and naturally couples the rapid non-linear amplification on small scales to larger length-scales. Hybrid magnetohydrodynamics-particle simulations are performed to confirm the effect. The resulting large-scale magnetic field may facilitate the scattering of high-energy cosmic rays as required to accelerate protons beyond the knee in the cosmic ray spectrum at supernova remnant shocks. Filamentation far upstream of the shock may also assist in the escape of cosmic rays from the accelerator.

The transport of cosmic rays in self-excited magnetic turbulence

The process of diffusive shock acceleration relies on the efficacy with which hydromagnetic waves can scatter charged particles in the precursor of a shock. The growth of self-generated waves is driven by both resonant and non-resonant processes. We perform high-resolution magnetohydrodynamic simulations of the non-resonant cosmic ray driven instability, in which the unstable waves are excited beyond the linear regime. In a snapshot of the resultant field, particle transport simulations are carried out. The use of a static snapshot of the field is reasonable given that the Larmor period for particles is typically very short relative to the instability growth time. The diffusion rate is found to be close to, or below, the Bohm limit for a range of energies. This provides the first explicit demonstration that self-excited turbulence reduces the diffusion coefficient and has important implications for cosmic-ray transport and acceleration in supernova remnants.

On the maximum energy of shock-accelerated cosmic rays at ultra-relativistic shocks

The maximum energy to which cosmic rays can be accelerated at weakly magnetised ultra-relativistic shocks is investigated. We demonstrate that for such shocks, in which the scattering of energetic particles is mediated exclusively by ion skin-depth scale structures, as might be expected for a Weibel-mediated shock, there is an intrinsic limit on the maximum energy to which particles can be accelerated. This maximum energy is determined from the requirement that particles must be isotropized in the downstream plasma frame before the mean field transports them far downstream, and falls considerably short of what is required to produce ultra-high-energy cosmic rays. To circumvent this limit, a highly disorganized field is required on larger scales. The growth of cosmic ray-induced instabilities on wavelengths much longer than the ion-plasma skin depth, both upstream and downstream of the shock, is considered. While these instabilities may play an important role in magnetic field amplification at relativistic shocks, on scales comparable to the gyroradius of the most energetic particles, the calculated growth rates have insufficient time to modify the scattering. Since strong modification is a necessary condition for particles in the downstream region to re-cross the shock, in the absence of an alternative scattering mechanism, these results imply that acceleration to higher energies is ruled out. If weakly magnetized ultra-relativistic shocks are disfavoured as high-energy particle accelerators in general, the search for potential sources of ultra-high-energy cosmic rays can be narrowed.

Empirical evidence for a nonlinear effect of galactic cosmic rays on clouds

Galactic cosmic ray (GCR) changes have been suggested to affect weather and climate, and new evidence is presented here directly linking GCRs with clouds. Clouds increase the diffuse solar radiation, measured continuously at UK surface meteorological sites since 1947. The ratio of diffuse to total solar radiation-the diffuse fraction, (DF)-is used to infer cloud, and is compared with the daily mean neutron count rate measured at Climax; Colorado from 1951-2000, which provides a globally representative indicator of cosmic rays. Across the UK, oil days of high cosmic ray flux (above 3600 X 10(2) neutron counts h(-1), which occur 87% of the time on average) compared with low cosmic ray flux, (i) the chance of an overcast day increases by (19 +/- 4)%; and (ii) the diffuse fraction increases by (2 +/- 0.3)%. During sudden transient reductions in cosmic rays (e.g. Forbush events), simultaneous decreases occur in the diffuse fraction. The diffuse radiation changes are; therefore; unambiguously due to cosmic rays. Although the statistically significant nonlinear cosmic ray effect is small, it will have a considerably larger aggregate effect on longer timescale (e.g. centennial) climate variations when day-to-day variability averages out.

Cosmic rays, clouds, and climate

It has been proposed that Earth's climate could be affected by changes in cloudiness caused by variations in the intensity of galactic cosmic rays in the atmosphere. This proposal stems from an observed correlation between cosmic ray intensity and Earth's average cloud cover over the course of one solar cycle. Some scientists question the reliability of the observations, whereas others, who accept them as reliable, suggest that the correlation may be caused by other physical phenomena with decadal periods or by a response to volcanic activity or El Niño. Nevertheless, the observation has raised the intriguing possibility that a cosmic ray–cloud interaction may help explain how a relatively small change in solar output can produce much larger changes in Earth's climate. Physical mechanisms have been proposed to explain how cosmic rays could affect clouds, but they need to be investigated further if the observation is to become more than just another correlation among geophysical variables.

Cloud base height and cosmic rays

Cosmic rays modify current flow in the global atmospheric electrical circuit. Charging at horizontal layer cloud edges has been observed to be consistent with global circuit vertical current flow through the cloud, which can modify the properties of small and pure water droplets. Studies have been hampered by the absence of cloud edge observations, hence cloud base height information is investigated here. Cloud base height measured at the Lerwick Observatory, Shetland, UK, is analysed using threshold tests and spectral analysis. The cloud base height distributions for low cloud (cloud base less than 800 m) are found to vary with cosmic ray conditions. Further, 27 day and 1.68 year periodicities characteristic of cosmic ray variations are present, weakly, in the cloud base height data of stratiform clouds, when such periodicities are present in neutron monitor cosmic ray data. These features support the idea of propagation of heliospheric variability into layer clouds, through the global atmospheric electric circuit.

Heliospheric modulation of galactic cosmic rays during grand solar minima: past and future variations

Galactic cosmic ray flux at Earth is modulated by the heliospheric magnetic field. Heliospheric modulation potential, Φ, during grand solar minima is investigated using an open solar flux (OSF) model with OSF source based on sunspot number, R, and OSF loss on heliospheric current sheet inclination. Changing dominance between source and loss means Φ varies in- (anti-) phase with R during strong (weak) cycles, in agreement with Φ estimates from ice core records of 10Be concentration, which are in-phase during most of the last 300 years, but anti-phase during the Maunder Minimum. Model results suggest “flat” OSF cycles, such as solar cycle 20 result from OSF source and loss terms temporarily balancing throughout the cycle. Thus even if solar activity continues to decline steadily, the long-term drop in OSF through SC21 to SC23 may plateau during SC24, though reemerge in SC25 with the inverted phase relation.