Application of CORSIKA Simulation Code to Study Lateral and Longitudinal Distribution of Fluorescence Light in Cosmic Ray Extensive Air Showers

In this paper, we used CORSIKA code to understand the characteristics of cosmic ray induced showers at extremely high energy as a function of energy, detector distance to shower axis, number, and density of secondary charged particles and the nature particle producing the shower. Based on the standard properties of the atmosphere, lateral and longitudinal development of the shower for photons and electrons has been investigated. Fluorescent light has been collected by the detector for protons, helium, oxygen, silicon, calcium and iron primary cosmic rays in different energies. So we have obtained a number of electrons per unit area, distance to the shower axis, shape function of particles density, percentage of fluorescent light, lateral distribution of energy dissipated in the atmosphere and visual field angle of detector as well as size of the shower image. We have also shown that location of highest percentage of fluorescence light is directly proportional to atomic number of elements. Also we have shown when the distance from shower axis increases and the shape function of particles density decreases severely. At the first stages of development, shower axis distance from detector is high and visual field angle is small; then with shower moving toward the Earth, angle increases. Overall, in higher energies, the fluorescent light method has more efficiency. The paper provides standard calibration lines for high energy showers which can be used to determine the nature of the particles.     

Coherent Cherenkov radio pulses from hadronic showers up to EeV energies

The Cherenkov radio pulse emitted by hadronic showers of energies in the EeV range in ice is calculated for the first time using full three dimensional simulations of both shower development and the coherent radio pulse emitted as the excess charge develops in the shower. A Monte Carlo, ZHAireS, has been developed for this purpose combining the high energy hadronic interaction capabilities of AIRES, and the dense media propagation capabilities of TIERRAS, with the precise low energy tracking and specific algorithms developed to calculate the radio emission in ZHS. A thinning technique is implemented to allow the simulation of radio pulses induced by showers up to 10 EeV in ice. The code is validated comparing the results for electromagnetic and hadronic showers to those obtained with GEANT4 and ZHS codes. The contribution to the pulse of other shower particles in addition to electrons and positrons, mainly protons, pions and muons, is found to be below 3% for 10 PeV and above proton induced showers. The characteristics of hadronic showers and the corresponding Cherenkov frequency spectra are compared with those from purely electromagnetic showers. The dependence of the spectra on shower energy and high-energy hadronic model is addressed and parameterizations for the radio emission in hadronic showers in ice are given for practical applications.     

Cherenkov τ shower earth-skimming method for PeV–EeV ντ observation with Ashra

We describe a method of observation for PeV–EeV τ neutrinos using Cherenkov light from the air showers of decayed τs produced by τ neutrino interactions in the Earth. Aiming for the realization of neutrino astronomy utilizing the Earth-skimming τ neutrino detection technique, highly precise determination of arrival direction is key due to the following issues: (1) clear identification of neutrinos by identifying those vertices originating within the Earth’s surface and (2) identification of very high energy neutrino sources. The Ashra detector uses newly developed light collectors which realize both a 42°-diameter field-of-view and arcminute resolution. Therefore, it has superior angular resolution for imaging Cherenkov air showers. In this paper, we estimate the sensitivity of and cosmic-ray background resulting from application of the Ashra-1 Cherenkov τ shower observation method. Both data from a commissioning run and a long-term observation (with fully equipped trigger system and one light collector) are presented. Our estimates are based on a detailed Monte Carlo simulation which describes all relevant shower processes from neutrino interaction to Cherenkov photon detection produced by τ air showers. In addition, the potential to determine the arrival direction of Cherenkov showers is evaluated by using the maximum likelihood method. We conclude that the Ashra-1 detector is a unique probe into detection of very high energy neutrinos and their accelerators.     

The optimum distance at which to determine the size of a giant air shower

To determine the size of an extensive air shower it is not necessary to have knowledge of the function that describes the fall-off of signal size from the shower core (the lateral distribution function). In this paper an analysis with a simple Monte Carlo model is used to show that an optimum ground parameter can be identified for each individual shower. At this optimal core distance, r_opt, the fluctuations in the expected signal, S(r_opt), due to a lack of knowledge of the lateral distribution function are minimised. Furthermore it is shown that the optimum ground parameter is determined primarily by the array geometry, with little dependence on the energy or zenith angle of the shower or choice of lateral distribution function. For an array such as the Pierre Auger Southern Observatory, with detectors separated by 1500 m in a triangular configuration, the optimum distance at which to measure this characteristic signal is close to 1000 m.     

Correlated features of arrival time and angle-of-incidence distributions of EAS muons

Various features and correlations of the arrival time and angle-of- incidence distributions of muons of extensive air showers (EAS) are studied by analyses based on Monte-Carlo simulations of the EAS development by using the air shower simulation code CORSIKA. Trends and dependencies of the temporal dispersion of the EAS muon component on shower size and distance from the shower core are displayed by the distribution of the arrival time and angle of incidence of the first muon, of the mean and median of the single shower distribution. Special attention is called to multi-correlations in observations at different radial distances from the shower core. These ‘radial’ correlations provide additional information for the discrimination of different EAS primaries, while the correlation of muon arrival time and angle-of-incidence is shown to improve the mass separation only insignificantly. This feature does not basically change, when arrival times and angles of incidence are displayed by transformed quantities like ‘muon production heights’.     

Asteroid Apophis: Evaluating the impact hazards of such bodies

Soon after the discovery of asteroid 99942 Apophis, it was classified as a potentially hazardous object with a high probability of an impact on the Earth in 2029. Although subsequent observations have substantially reduced the probability of a collision, it has not been ruled out; moreover, similar-sized asteroids in orbits intersecting the Earth’s orbit may well be discovered in the near future. We conduct a numerical simulation of an atmospheric passage and an impact on the Earth’s surface of a stony cosmic body with a diameter of 300 m and kinetic energy of about 1000 Mt, which roughly corresponds to the parameters of the asteroid Apophis, at atmospheric entry angles of 90° (vertical stroke), 45°, and 30°. The simulation is performed by solving three-dimensional equations of hydrodynamics and radiative transfer equations in the approximations of radiative heat conduction and volume emission. The following hazards are considered: an air shock wave, ejecta from the crater, thermal radiation, and ionospheric disturbances. Our calculations of the overpressure and wind speed on the Earth’s surface show that the zone of destruction of the weakest structures can be as large as 700–1000 km in diameter; a decrease in the flight path angle to the surface leads to a marked increase in the area affected by the shock wave. The ionospheric disturbances are global in nature and continue for hours: at distances of several thousand kilometers at altitudes of more than 100 km, air density disturbances are tens of percent and the vertical and horizontal velocity components reach hundreds of meters per second. The impact of radiation on objects on the Earth’s surface is estimated by solving the equation of radiative transfer along rays passing through a luminous area. In clear weather, the size of the zone where thermal heating may ignite wood can be as large as 200 km, and the zone of individual fire outbreaks associated with the ignition of flammable materials can be twice as large. In the 100-km central area, which is characterized by very strong thermal damage, there is ignition of structures, roofs, clothes, etc. The human hazardous area increases with the decrease in the trajectory angle, and people may experience thermal effects at distances of up to 250–400 km from the crater.     

One-dimensional hybrid approach to extensive air shower simulation

An efficient scheme for one-dimensional extensive air shower simulation and its implementation in the program conex are presented. Explicit Monte Carlo simulation of the high-energy part of hadronic and electro-magnetic cascades in the atmosphere is combined with a numeric solution of cascade equations for smaller energy sub-showers to obtain accurate shower predictions. The developed scheme allows us to calculate not only observables related to the number of particles (shower size) but also ionization energy deposit profiles which are needed for the interpretation of data of experiments employing the fluorescence light technique. We discuss in detail the basic algorithms developed and illustrate the power of the method. It is shown that Monte Carlo, numerical, and hybrid air shower calculations give consistent results which agree very well with those obtained within the corsika program.     

Interpretation of the time structure of the EAS disc measured by the GREX/COVER_PLASTEX experiment

Simulation results for the time structure of the extensive air shower disc are presented and compared with data from the GREX/COVER_PLASTEX experiment. The distribution of the arrival times at various distances from the shower core and the contributions from the secondary particles to the shape of the distribution are described. The main parameters of the distribution, the mean time of arrival τ and the standard deviation σ, reflect the shower disc profile and thickness. The dependence of the shower profile and thickness on the energy and mass of the primary particle initiating the shower as well as on its inclination angle is discussed. The influence of the experimental conditions on the disc profile and thickness measured by the GREX/COVER_PLASTEX experiment is analysed and a parametrization of the average profile and thickness is given.     

Coherent radio emission from cosmic ray air-showers in the turbulent atmosphere

Coherent electromagnetic erenkov radiation is produced by cosmic ray air showers passing through the atmosphere. This radiation is detected by radio telescopes. We demonstrate here that the effect of random spatial fluctuations in the refractive index of air, about a mean exceeding unity, causes the airshower to emit not only the spontaneous coherent radio emission described elsewhere by Kahn and Lerche, but also an induced radiation field which can exceed the spontaneous field in certain frequency bands. Further the conditions for emission of the coherent radio erenkov radiation are altered by the presence of the refractive index fluctuations. And the Earth's magnetic field gives rise to the dominant term in the far-field radiation, be it spontaneous or induced, since it causes a systematic separation of electrons and positrons in the shower which, for parameters currently acceptable for air showers, is the major factor in determining the far-field radiation pattern. Also we suggest that the coherent 500 Mc/sec radiation seen from occasional showers is probably a reflection of an atmospheric correlation length of order 15 cm at the time the shower passes through the atmosphere.     

The air shower maximum probed by Cherenkov effects from radio emission

Radio detection of cosmic-ray-induced air showers has come to a flight the last decade. Along with the experimental efforts, several theoretical models were developed. The main radio-emission mechanisms are established to be the geomagnetic emission due to deflection of electrons and positrons in Earth’s magnetic field and the charge-excess emission due to a net electron excess in the air shower front. It was only recently shown that Cherenkov effects play an important role in the radio emission from air showers. In this article we show the importance of these effects to extract quantitatively the position of the shower maximum from the radio signal, which is a sensitive measure for the mass of the initial cosmic ray. We also show that the relative magnitude of the charge-excess and geomagnetic emission changes considerably at small observer distances where Cherenkov effects apply.     

The shape of the radio wavefront of extensive air showers as measured with LOFAR

Extensive air showers, induced by high energy cosmic rays impinging on the Earth’s atmosphere, produce radio emission that is measured with the LOFAR radio telescope. As the emission comes from a finite distance of a few kilometers, the incident wavefront is non-planar. A spherical, conical or hyperbolic shape of the wavefront has been proposed, but measurements of individual air showers have been inconclusive so far. For a selected high-quality sample of 161 measured extensive air showers, we have reconstructed the wavefront by measuring pulse arrival times to sub-nanosecond precision in 200 to 350 individual antennas. For each measured air shower, we have fitted a conical, spherical, and hyperboloid shape to the arrival times. The fit quality and a likelihood analysis show that a hyperboloid is the best parameterization. Using a non-planar wavefront shape gives an improved angular resolution, when reconstructing the shower arrival direction. Furthermore, a dependence of the wavefront shape on the shower geometry can be seen. This suggests that it will be possible to use a wavefront shape analysis to get an additional handle on the atmospheric depth of the shower maximum, which is sensitive to the mass of the primary particle.     

Muon density measurements with the KASCADE central detector

Frequency distributions of local muon densities in high-energy extensive air showers (EAS) are presented as signature of the primary cosmic ray energy spectrum in the knee region. Together with the gross shower variables like shower core position, angle of incidence, and the shower sizes, the KASCADE experiment is able to measure local muon densities for two different muon energy thresholds. The spectra have been reconstructed for various core distances, as well as for particular subsamples, classified on the basis of the shower size ratio N_μ/N_e. The measured density spectra of the total sample exhibit clear kinks reflecting the knee of the primary energy spectrum. While relatively sharp changes of the slopes are observed in the spectrum of EAS with small values of the shower size ratio, no such feature is detected at EAS of large N_μ/N_e ratio in the energy range of 1–10 PeV. Comparing the spectra for various thresholds and core distances with detailed Monte Carlo simulations the validity of EAS simulations is discussed.     

On the transfer of radiation at asteroidal surfaces in relation to their orbit deflection — II

The efficiency of absorption of X-rays generated by a nuclear explosion at the surface of an asteroid, estimated earlier, is used to calculate the explosion yield needed to deflect the orbit of an asteroid. Following the work of Ahrens &38; Harris, it is shown that a recoil velocity of 1 cm s⁻¹ is required to deflect an asteroid from a collision course with the Earth, and the necessary yield of explosion energy is estimated. If it is assumed that the scaling law between the energy and the diameter of the resulting crater, obtained from experiments carried out on the Earth, remains valid on the asteroid surface, where gravity is much weaker, an explosion energy of 8 and 800 megaton (Mton) equivalent of TNT would be required for asteroids of diameter 1 and 10 km respectively. If, on the other hand, the crater diameter is proportional to a certain power of the gravity g , the power being determined from a dimension analysis, 130 kton and 12 Mton would be required to endow asteroids of diameters 1 and 10 km with the required velocity, respectively. The result indicates that in order to estimate the required explosion energy, a better understanding of cratering under gravity much weaker than on the Earth would be required.     

Radar absorption due to a corotating interaction region encounter with Mars detected by MARSIS

Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) is a subsurface and topside ionosphere radar sounder aboard the European Space Agency spacecraft Mars Express, in orbit at Mars since 25 December 2003, and in operation since 17 June 2005. The ionospheric sounding mode of MARSIS is capable of detecting the reflection of the sounding wave from the martian surface. This ability has been used in previous work to show that the surface reflection is absorbed and disappears during periods when high fluxes of energetic particles are incident on the ionosphere of Mars. These absorption events are believed to be the result of increased collisional damping of the sounding wave, caused by increased electron density below the spacecraft, in turn caused by impact ionization from the impinging particles. In this work we identify two absorption events that were isolated during periods when the surface reflection is consistently visible and when Mars is nearly at opposition. The visibility of the surface reflection is viewed in conjunction with particle and photon measurements taken at both Mars and Earth. Both absorption events are found to coincide with Earth passing through solar wind speed and ion flux signatures indicative of a corotating interaction region (CIR). The two events are separated by an interval of approximately 27 days, corresponding to one solar rotation. The first of the two events coincides with abruptly enhanced particle fluxes seen in situ at Mars. Simultaneous with the particle enhancement there are an abrupt decrease in the intensity of electron oscillations, typically seen by the Mars Express particle instrument ASPERA-3 between the magnetic pileup boundary and the martian bow shock, and a sharp drop in the solar wind pressure, seen in the proxy quantity based on MGS magnetometer observations. The decrease in oscillation intensity is therefore the probable effect of a relaxation of the martian bow shock. The second absorption event does not show a particle enhancement and complete ASPERA-3 data during that time are unavailable. Other absorption events are the apparent result of solar X-ray and XUV enhancements. We conclude that surface reflection absorption events are sometimes caused by enhanced ionospheric ionization from high energy particles accelerated by the shocks associated with a CIR. A full statistical analysis of CIRs in relation to observed absorption events in conjunction with a quantitative analysis of the deposition of ionization during space weather events is needed for a complete understanding of this phenomenon. If such analyses can be carried out, radar sensing of the martian ionosphere might be useful as a space weather probe.     

Searching for tau neutrinos with Cherenkov telescopes

Cherenkov telescopes have the capability of detecting high energy tau neutrinos in the energy range of 1–1000 PeV by searching for very inclined showers. If a tau lepton, produced by a tau neutrino, escapes from the Earth or a mountain, it will decay and initiate a shower in the air which can be detected by an air shower fluorescence or Cherenkov telescope. In this paper, we present detailed Monte Carlo simulations of corresponding event rates for the VERITAS and two proposed Cherenkov Telescope Array sites: Meteor Crater and Yavapai Ranch, which use representative AGN neutrino flux models and take into account topographic conditions of the detector sites. The calculated neutrino sensitivities depend on the observation time and the shape of the energy spectrum, but in some cases are comparable or even better than corresponding neutrino sensitivities of the IceCube detector. For VERITAS and the considered Cherenkov Telescope Array sites the expected neutrino sensitivities are up to factor 3 higher than for the MAGIC site because of the presence of surrounding mountains.     

Predictions for the Aurigid Outburst of 2007 September 1

The September 2007 encounter of Earth with the 1-revolution dust trail of comet C/1911 N1 (Kiess) is the most highly anticipated dust trail crossing of a known long period comet in the next 50 years. The encounter was modeled to predict the expected peak time, duration, and peak rate of the resulting outburst of Aurigid shower meteors. The Aurigids will radiate with a speed of 67 km/s from a radiant at R.A. = 92°, Decl. = +39° (J2000) in the constellation Auriga. The expected peak time is 11:36 ± 20 min UT, 2007 September 1, and the shower is expected to peak at Zenith Hourly Rate = 200/h during a 10-min interval, being above half this value during 25 min. The meteor outburst will be visible by the naked eye from locations in Mexico, the Western provinces of Canada, and the Western United States, including Hawaii and Alaska. A concerted observing campaign is being organized. Added in proof: first impression of the shower. Prepared as a contribution to the conference proceedings of “Meteoroids 2007”, to be published in the journal “Earth, Moon, and Planets”.     

The center of lateral iso-density contours for inclined cosmic air showers

The horizontal lateral density of a cosmic air shower with a non-zero zenith angle is asymmetric. The asymmetry consist of a stretching of the iso-density contours to ellipses and to a shift of the center of the elliptic contours with respect to the core of the shower. The shift is caused by atmospheric attenuation. The modeling of the attenuation results in an equation for the shift as a function of zenith angle and the size of the iso-density contours. A more accurate equation is obtained by investigating the shift in lateral densities of simulated showers. It is shown how the shift can be incorporated in an elliptic lateral density function. A linear approximation for the shift allows for an analytical solution for the shifted elliptic density. Its predictions for the polar variations of the density are compared with data of simulated showers.     

Measurements of attenuation and absorption lengths with the KASCADE experiment

The attenuation of the electron shower size beyond the shower maximum is studied with the KASCADE extensive air shower (EAS) experiment in the primary energy range of about 10¹⁴–10¹⁶ eV. Attenuation and absorption lengths are determined by applying different approaches, including the method of constant intensity, the decrease of the flux of EASs with increasing zenith angle, and its variation with ground pressure. We observe a significant dependence of the results on the applied method. The determined values of the attenuation length ranges from 175 to 196 g/cm² and of the absorption length from 100 to 120 g/cm². The origin of these differences is discussed emphasizing the influence of intrinsic shower fluctuations.     

Observations of the 2001 Leonid meteor shower using VHF meteor radar

Every year the Earth crosses or passes near one of the dust trails left by Comet 55P/Tempel-Tuttle in its pass through the Solar System every 33.2 years. This produces a meteor shower Commonly called the Leonid. The 2001 Leonid meteor shower is one of the strongest in recent years. We present observations made by the 50 MHz all-sky meteor radar located at the Platteville Atmospheric Observatory in Colorado (40° N, 105° W). The spatial and temporal distributions of the meteor activity detected by the radar during the 2001 Leonid shower differs from the observed sporadic activity detected by VHF radars. Estimation of the radiant flux of the meteor shower of the shower by a well-known methodology is presented, and the intensity of the phenomena is discussed.