The Radio Signal From Extensive Air Showers


The field of ultra-high energy cosmic rays made a lot of progresses last years with large area experiments such as the Pierre Auger Observatory, HiRes and the Telescope Array. A suppression of the cosmic ray flux at energies above  eV is observed at a very high level of significance but the origin of this cut-off is not established: it can be due to the Greisen-Zatsepin-Kuzmin suppression but it can also reflect the upper limit of particle acceleration in astrophysical objects. The key characteristics to be measured on cosmic rays is their composition. Upper limits are set above  eV on primary photons and neutrinos and primary cosmic rays are expected to be hadrons. Identifying the precise composition (light or heavy nuclei) will permit to solve the puzzle. It has been proven that the radio signal emitted by the extensive air showers initiated by ultra-high energy cosmic rays reflects their longitudinal profile and can help in constraining the primary particle. We review in this paper the emission mechanisms as a function of the frequency of the electric field.

1 Observing extensive air showers

In the current large area ultra-high energy cosmic rays experiments (Auger [1, 2] and TA [3, 4]), extensive air showers are studied through the secondary particles reaching the ground level and by collecting the fluorescence light emitted by the atmospheric diazote molecules excited by the passage of the charged particles. Surface detectors (SD) arrays are composed of particle detectors (Cerenkov water tanks, scintillators…) having a high duty cycle (around  %) and that allow to estimate the arrival direction of the shower and its core position. During dark nights (duty cycle of  % in Auger), the fluorescence detectors (FD) detect the longitudinal profile of the showers through a calorimetric measurement during their development in the air. A FD directly measures the shower energy and estimates the nature of the primary cosmic ray by the determination of the atmospheric depth corresponding to the maximum number of secondary particles in the shower. Hybrid events, detected by both SD and FD, allow to calibrate the SD using the FD reconstruction of the energy. The need for a better event-by-event identification of the nature of the primary particle have speed up the research and development of additionnal detection techniques such as the measurement of the radio signal emitted by the electrons and positrons of the showers.

2 Electric field from air showers

The radio signal associated to cosmic rays was detected for the first time in 1965[5] but the detailed characteristics of the electric field have been understood in the last years using the data from several experiments (CODALEMA [6], LOPES [7], LOFAR [8], AERA [9]…). It has been shown that the emission is coherent for wavelengths larger than or of the order of the typical dimensions of air showers: the scale of their longitudinal development ( km,  kHz), their lateral spread ( km,  MHz) and the thickness of their front ( m,  MHz). The individual electric fields emitted by all high energy ( MeV) electrons and positrons of air shower (we neglect the emission from muons and ions, being much heavier) interfere with each other: if the observation wavelength is larger than the size of the emitting region, the fields add-up coherently so that the total electric field is proportionnal to the number of electrons and positrons, and then, to the primary energy. On the contrary, if the wavelength is smaller than the size of the emitting region, the fields interfere destructively. We should therefore expect a cut-off frequency in the full frequency spectrum. The radio signal has a clear coherent contribution over a large frequency range: from tens of kHz up to some hundreds of MHz. It has been shown that the signal in the frequency range 30-80 MHz allows to reconstruct precisely the incoming direction of the shower (angular resolution below , see [10]). The deep understanding of the radio signal we reached today allows to use simulations to estimate the and LOFAR [11] obtained a resolution on this parameter of , similar to the resolution achieved by the FD. The energy estimation is underway, for instance in the AERA experiment.

2.1 Emission mechanisms

In 1967, Kahn and Lerche [12] studied the dominant mechanism: the emission from the time varying transverse currents induced by the interaction between the geomagnetic field and the charged particles (hereafter the geomagnetic mechanism). The systematic opposite drift of electrons and positrons when they propagate in the geomagnetic field associated to the interactions with the medium results into an average transverse current (with respect to the shower axis). This current induces a macroscopic electric field linearly polarized along the direction, where is the direction of the shower axis and the direction of the geomagnetic field. The orientation of this electric field is independent on the observer’s location.

A secondary emission process is due to the excess of electrons with respect to positrons in air showers (Askaryan [13]). This net negative charge excess is due to the fast in-flight positrons annihilation and to direct enrichment in electrons extracted from the medium by the Compton, Bhabha and Moeller diffusions. The charge-excess mechanism leads to a radial polarization pattern of the electric field in the plane transverse to the shower axis.

An observer will detect the superposition of both electric fields. These fields depend on time through the varying number of electrons and positrons in the showers during its development in the air. This superposition leads to an asymetry in the total electric field strength according to the observer’s location with respect to the shower core and explains the departure from the radially symetric lateral distribution function proposed by Allan [14] in 1971. The electric field value cannot be simply modelized by a 1D lateral function. Figure 1 presents the superposition of both mechanisms.

Figure 1: Left: in green, electric field generated by the geomagnetic mechanism. The polarization is linear and the same for all observers. Center: in red, electric field generated by the charge-excess mechanism. The polarization is linear with a radial pattern with respect to the shower core in the plane transverse to the shower axis. Right: in blue, total electric field measured by an observer. Its value and direction depends on the observer’s location.

Note that other mechanisms can also lead to the emission of electric fields: the presence of electric dipoles in the shower, the emission from the ions that are left behind the shower, the individual geosynchrotron emission but all of them are largely subdominant.

2.2 Simulation approaches

The simulation of the electric field emitted by air showers recently reached a mature and robust state. Different approaches were chosen in the past, using microscopic or macroscopic descriptions, in the time or frequency domains. The codes MGMR [15] and EVA [16] use the macroscopic transverse current as the main ingredient to the Maxwell’s equations. The code SELFAS [17, 18] uses the energy, impulsion, angular, position distributions of electrons and positrons from the shower universality description and computes the resulting electric field as the summation of all individual contributions (see Eq. 1). The code ZHAireS [19] uses the ZHS formalism and computes the electric field in both time and frequency domains. CoREAS [20] uses the end-point formalism and computes all individual electric fields directly inside a full CORSIKA simulated shower. As an example, the following formula used in SELFAS gives the electric field emitted by a single particle of charge and lifetime (its charge is taken as before and after ):


where is the observer’s position, is the unit vector between the particle and the observer, is the air refractive index, the particle velocity and the distance between the particle and the observer. The field is evalutated at time . We obtain the complete electric field after summation over all electrons and positrons of the shower. The first term corresponds to the Coulombian contribution; it is negligible with respect to the two others. The global contribution of the second term for a system with equal numbers of electrons and positrons vanishes. Due to the net excess of electrons in air showers, this term is not negligible and constitutes the charge excess contribution. The third term contains the particle velocity and its time derivative is proportionnal to the Lorentz force: this is the geomagnetic term. The particle velocity direction is close, at first order, to the shower axis direction. We note the presence of the air refractive index  leading to a boost of the electric field when the angle between the particle velocity and the line of sight is close to the Cherenkov angle. This is not actual Cherenkov radiation but this effect is important for observers located close to the shower axis. For these observers, in the time domain, the electric field pulse is sharpened making it dominant at higher frequencies (above  MHz) and still important below  MHz and this explains the flat electric profiles observed by LOPES [21] and CODALEMA in the past. The electric field is amplified along a ring of radius  m,  km being is the typical distance from the shower’s maximum for  eV primary energies and is the Cherenkov angle in air. Observation of this ring would be of great interest because it probes the distance to the source of the electric field and therefore could help in differentiating proton or iron initiated showers. The ring radius being small, it would be necessary to detect shower with a dense array of radio stations which is not a good solution (not cheap).

3 Radio signal in the range  MHz

The frequency band  MHz is intensively studied since the beginning of years 2000. It has been established that the main emission mechanism in this band is the geomagnetic contribution. Beside this contribution, the CODALEMA and AERA experiments observed the secondary electric field emitted by the charge excess mechanism using different observables. The CODALEMA experiment compared the shower core positions obtained using the data from the radio array to the core positions obtained using the particle array. A significant statistical discrepancy between the two core positions was found: on average, the radio core positions is shifted by  m to the east with respect to the particle core position. The radio core position is estimated using a radio lateral distribution function (LDF) of the electric field amplitude depending only on the distance to the shower axis; ie a radially symmetric LDF. Using SELFAS simulations including only the geomagnetic contribution, it has been shown that the core positions estimated from the radio and particle arrays were in good agreement. Using SELFAS with both mechanisms permitted to reproduced the same shift toward the east. This was the first strong indication of the presence of the charge excess mechanism in the data [22].

The AERA experiment used the polarization information to exhibit the charge excess contribution. Following Figure 1, an observer measures the superposition of two electric fields with different polarization patterns. If we consider only the geomagnetic mechanism, we can compute the polarization angle (with respect to the geographic east) in the horizontal plane using the north-south (NS) and east-west (EW) components of the predicted electric field:

Using the data, we can compute for each radio station in a given event the measured polarization angle in the horizontal plane (the antennas used in AERA do measure the electric field in both EW and NS directions). This angle is given by:

where is the electric field amplitude in both EW and NS directions. The correlation between and in the case of a pure geomagnetic contribution is presented in Figure 2(left). For a given event with a known incoming direction (ie we know the axis direction ) involving radio stations, we have a single value for but we have values of the measured polarization angles . This explains the bunches of different values of at a given value of . The correlation is very clear: the Pearson correlation coefficient is at 95% CL (with ) which demonstrates that the geomagnetic is indeed dominant. As we know that the charge excess mechanism is also there, we can refine our predicted value of the polarization angles . As discussed before, the electric field from this mechanism has a radial linar polarization in the plane transverse to the shower axis. The predicted polarization angle at each radio station can be deduced using the following formula:

where is the azimuth of the geomagnetic contribution, the angle between the shower axis and the geomagnetic field and is the azimuth of the charge excess contribution at the location of the radio station . The parameter is the ratio of the amplitude of the electric field from the charge excess contribution to the amplitude of the electric field from the geomagnetic contribution , modulated by (in order to cancel the geomagnetic contribution when the shower incoming direction is parallel to the geomagnetic field):

The average value of was computed using showers detected by AERA and we obtained: (see [23]). This value should not be considered as the fraction of charge excess in the shower, as it depends also on the angular distance between the shower axis and the geomagnetic field through . Using this value for , we can finally compute the predicted polarization angle for each radio station involved in a given event, taking into account both geomagnetic and charge excess contributions. The correlation between and is presented in Figure 2(right). The Pearson correlation coefficient rises up to at 95% CL (with ). The agreement between the data and the model is much better when considering both mechanisms.

Figure 2: Left: correlation between the measured polarization angle and the predicted polarization angle in the case of a pure geomagnetic contribution. Right: same than the left figure but adding the contribution of the charge excess mechanism.

In conclusion, the radio signal in the range  MHz is very well described by the superposition of two mechanisms of electric field emission: the geomagnetic and the charge excess contributions.

4 Radio signal in the GHz domain

The first signal we can think about in this frequency region is the GHz counterpart of the usual signal made of the geomagnetic and charge excess mechanisms, enhanced by the Cherenkov-like effect due to the air index effect as discussed in section 2.2; this emission has a steeply falling spectrum due the the incoherence of the fields at these frequencies. The Cherenkov-like effect results in a highly forward-beamed emission. As seen previously, this electric field is polarized accordingly to the superposition of both mechanisms. Figure 3(left) presents the Fourier spectra of the simulated electric field for a vertical shower at  eV and antennas located at different distances from the shower axis. For the antennas located at and  m from the axis, we see the coherent part of the spectrum; this corresponds to the coherent geomagnetic and charge excess electric field between  MHz and some hundreds of MHz. The signal is incoherent at higher frequencies because these frequencies corrresponds to wavelength much smaller than the emittive zone.

Figure 3: Contribution at high frequencies of the electric field from the geomagnetic and charge excess mechanisms. Left: electric field power as a function of frequency. The coherent part lies at smaller frequencies as can be seen for antennas close to the shower axis ( m and  m). Right: amplitude of the same electric field between  MHz and  GHz a a function of the observer location at the groud level. Both plots were obtained for a simulated vertical shower (using SELFAS) initiated by an iron nucleus at  eV.

Figure 3(right) shows the electric field amplitude between  MHz and  GHz as a function of the observer position at the ground level. The shower core is at the origin. At these frequencies, the Cherenkov-like ring clearly appears. This ring is due to the enhancement of the electric field when the particles velocity and the line of sight is close to the Cherenkov angle in air (see section 2.2).

The frequency domain was roughly separated in this paper between ”below GHz” and ”above GHz”. This because another emission mechanism was studied and measured by Gorham et al[24] in 2004 using a test beam at Argonne and SLAC: it is the molecular Bremsstrahlung radiation (MBR). This emission is due to the deceleration of the low-energy electrons ( eV) in the plasma created by the high energy electrons and positrons of the shower. We remind that the electrons and positrons of the shower mainly loose energy by exciting or ionizing the medium. The excited N then emits UV light detected by fluoresence telescopes but the ionized molecules provide low energy electrons forming a weakly ionized and stationary plasma. These electrons emit a Bremsstrahlung radiation that is expected to be non-polarized, isotropic which is potentially an excellent feature for a detection at large distance from the shower axis (the radiation is not focalized). The MBR implies a microwave continuum emission at the GHz level and is expected to have a direct relationship with the shower energy. Finally, this radiation is detectable with a duty cycle close to %. The experiment of Gorham et al. measured a non-polarized flux density of for a setup corresponding to a shower at  eV at  m from the axis. Extrapolation to a typical observation distance of  km leads to a threshold shower energy for a detection at of  eV and  eV assuming a quadratic or linear dependence of the power with respect to the shower energy respectively. The MAYBE [25] (electron beam) and AMY [26] (electron or positron beam) experiments also aimed at detecting and characterizing the MBR. The results disagree with those of Gorham. Recently, Conti et al. [27], using a low energy electron beam (in order to be below the Cherenkov threshold), obtained a linear scaling of the signal with the number of electrons and an asymmetric emission pattern, in opposition to the isotropic feature reported by Gorham et al. In conclusion on the beam experiments, the situation is very unclear but the high MBR signal reported by Gorham et al. is not confirmed by any of the other experiments.

AMBER [24], MIDAS [28] and EASIER [29] are three prototypes installed at the Pierre Auger Observatory aiming at detecting the GHz emission from air showers. MIDAS excludes at quadratic (resp. linear) scaling of the power flux with the primary energy as reported in at a confidence level of (resp. ). EASIER detected five (up to summer 2014) GHz events, in 3 years of data acquisition. All of those were detected by a single antenna and at a small distance from the shower axis (below  m) and for high-energy showers. The simulation with SELFAS of the first detected EASIER shower was in good agreement with the data: it shows that the GHz signal can be explained with the usual geomagnetic and charge excess mechanisms in the GHz range, for this event. The others have not been simulated.

The CROME [30] experiment, installed at the KASCADE experiment detected showers in the GHz range. The 35 detected showers in coincidence with the KASCADE setup after  months of data taking presented an east-west asymmetry in their electric field strength. High values of the electric field strength were measured following a ring pattern at the ground level. The detected GHz emission is strongly polarized with a pattern in agreement with the geomagnetic and charge excess expectations. The conclusion of CROME is that the GHz signal associated to air showers is fully consistent with geomagnetic and charge excess mechanisms as sole emission processes at these frequencies and the data rule out the MBR mechanism.

Name location year freq. (GHz) scaling emission pattern
Gorham SLAC 2004 - quadratic isotropic
MAYBE Argonne 2012 - linear isotropic, Gorham flux
AMY Frascati 2012 - MBR much smaller than Cherenkov
Conti et al. Padova 2014 linear peaked forward
AMBER Auger 2011 no CR detection
MIDAS Auger 2012 no CR detection
EASIER Auger 2011-2012 -

5 CRs detected at high energy, very close to the shower axis

CROME KASCADE 2011-2012 -

35 CRs detected in agreement with geomagnetic and charge excess mechanisms, no MBR required

Table 1: Summary of the results obtained in the GHz domain, seeking for the characteristics of the molecular Bremsstrahlung radiation (MBR).

In conclusion, the GHz signal strength and characteristics as reported in [24] are not confirmed by any other experiments. The GHz emission from showers is compatible with the GHz counterpart of the usual geomagnetic and charge excess mechanisms. No MBR emission associated with EAS is clearly observed up to now and clearly, detecting EAS using the GHz signal is not efficient.

5 Radio signal below  MHz

On the other side of the frequency spectrum, below  MHz, we can of course think to the contributions of the geomagnetic and charge excess mechanisms. Nevertheless, these contributions have not been carefully measured at the MHz level because of the atmospheric noise which increases with decreasing frequency, this explains why people favoured higher frequencies. In the past, several experiments reported the observation of radio pulses in coincidence with EAS. A complete review of the experimental status can be found in Revenu & Marin [31]. The frequencies probed at that time were between some hundreds of kHz up to  MHz. The conclusions were that a very strong signal (1 or 2 orders of magnitude higher than in the usual band - MHz) associated to air showers was undoubtly detected. People thought about several new mechanisms to explain this signal: the interaction between the ionization electrons in the air with the atmospheric electric field, the transition radiation front the electrons when the shower front hits the ground or the longitudinal and transverse emissions assuming full coherence but none of them could explain the reported huge values of the electric field. In [31], we propose a new emission mechanism: the coherent deceleration of the electrons and positrons of the shower front when it hits the ground level. The sudden deceleration of these particles when hitting the ground (which we call the sudden death of the shower) produces a coherent Bremsstrahlung emission for wavelengths larger than the typical size of the shower front, or equivalently for frequencies smaller than  MHz. We can also understand this emission using a macroscopic point of view. When the shower front disappears below the ground level, the macroscopic charge density and current vary very quickly and their space and time derivatives, providing electric fields according to Maxwell’s equations, reach strong values. We therefore expect a significant electric field emission and we suspect that the previous experiments did measure this signal. We used the simulation code SELFAS to characterize the electric field emitted by the sudden death mechanism. The basic formula, derived from Eq. 1 and using the Coulomb gauge reads:


where, as in Eq. 1, is the air refractive index, and are the line of sight and the distance between the observer and the particle , the velocity of this particle and its electric charge. The summation is done over the total number of particles that emitted an electric field detected by the observer at time . All these quantities are evaluated at the retarded time , related to the observer’s time through . Figure 4 presents the vertical component of the electric received by two observers located at  m and  m of the shower axis for vertical showers with a primary energy of  eV and  eV. The origin of time corresponds to the time when the shower front hits the ground.

Figure 4: Vertical component of the electric field as a function of time for two observers located at  m and  m of the shower core for vertical simulated showers with primary energy of  eV and  eV. The electric field emitted during the shower development in the air is clearly visible around  ns. The sudden death pulse, generated when the shower front hits the ground, appears at  ns and  ns for an observer at  m and  m, respectively, from the shower core.

In this figure, the electric field pulse created by the disappearance of the shower front in the ground has a very specific shape. Its maximum is located at a time corresponding to the time needed for the signal to reach the observer from the shower core. The reception time of this signal is therefore delayed by with respect to the time at the shower core, ie  ns at  m and  ns at  m. The shape of the pulse can be understood as the shape of the ground distribution of the electrons and positrons. The simulation shows that the sudden death pulse is polarized according to Eq. 2 in the direction of because at first order, the of the particles at the ground level is parallel to the shower axis . The vertical component should therefore be important. The sudden death pulse amplitude scales linearly with the primary energy and, a very important feature, is that this amplitude decreases as where is the distance to the shower core. It means that it should be possible to detect showers up to some kilometers from the core, contrarily to the signal in - MHz. Figure 5 summarizes our current understanding of the emission mechanisms as a function of the frequency.

Figure 5: Summary of the current status of the experiments and of the understood mechanisms of electric field emission in showers, as a function of the frequency.

6 Conclusion

At the time of this conference, the radio signal emitted by extensive air showers is understood at a very high level of precision. This signal, between 30 MHz and 80 MHz, allows to reconstruct precisely the arrival direction of the primary cosmic ray. It also allows to estimate its nature following the results from the LOFAR experiment which claims a resolution of  g cm. This value is comparable to the resolution of a fluorescence detector but with an uptime close to 100% which is very interesting for improving our knowledge on ultra-high energy cosmic rays. The drawback of this signal is the relatively small range of the order of some hundreds of meters. The contribution from the geomagnetic and charge excess mechanisms can explain the data in a large range of frequencies, between 30 MHz up to some GHz. At low frequency, below 20 MHz, showers have already been observed in the past. A very large electric field amplitude was observed with no satisfactory underlying mechanism. We are currently installing the EXTASIS experiment in Nançay, dedicated to the low frequency signal of the showers when the front hits the ground. The range of this sudden death signal is much higher than the range in 30-80 MHz; this could be a very important feature for the detection of cosmic ray with a large efficiency.



  1. I. Allekotte et al. Nucl. Instrum. Meth. A, 586:409–420, 2008.
  2. J. Abraham et al. Nucl. Instrum. Meth. A, 620:227–251, 2010.
  3. T. Abu-Zayyad et al. Nucl. Instrum. Meth. A, 689(0):87– 97, 2012.
  4. P. Tinyakov. Nucl. Instrum. Meth. A, 742(0):29—- 34, 2014. 4th RICAP
  5. J.V. Jelley et al. Nature, 205:327–328, 1965.
  6. D. Ardouin et al. Nucl. Instrum. Meth. A, (555):148–163, 2005.
  7. H. Falcke et al. (LOPES Collaboration). Nature, (435):313, 2005.
  8. A. Nelles and the LOFAR collaboration. 33rd ICRC, number 0558, July 2013.
  9. J. Neuser and The Pierre Auger Collaboration. 6 ARENA, Annapolis, USA, June 2014.
  10. B. Revenu. 32th ICRC, 2011.
  11. S. Buitink and the LOFAR collaboration. 33rd ICRC, number 0579, July 2013.
  12. F.D. Kahn and I. Lerche. Proc. Roy. Soc. A, 289:206, 1966.
  13. G.A. Askaryan. J. Exp. Theor. Phy., 21:658, 1962.
  14. H. R. Allan. Progress in elementary part. and cosmic ray phys.. North Holland, 1971.
  15. K.D. de Vries et al. Astropart. Phys., 34:267, 2010.
  16. K. Werner, K. D. de Vries, and O. Scholten. Astropart. Phys., 37:5–16, 2012
  17. V. Marin and B. Revenu. Astropart. Phys., 35:733–741, 2012.
  18. V. Marin. 5 ARENA, Erlangen, Germany, June 2012.
  19. J. Alvarez Muiz, W. R. Carvalho Jr, and E. Zas. Astropart. Phys., 35:325–341, 2012.
  20. T. Huege and C. James. 33rd ICRC, July 2013.
  21. W.D. Apel et al. Astropart.Phys., 32:294–303, 2010.
  22. V. Marin and the CODALEMA Collaboration. 32nd ICRC, 2011.
  23. A. Aab et al. (The Pierre Auger Collaboration). Phys. Rev. D, 89:052002, March 2014.
  24. P. W. Gorham et al. Phys. Rev. D, 78:032007, Aug 2008.
  25. C. Williams et al. European Physical Journal Web of Conferences, 53:8008, June 2013.
  26. M. Bohacova and the AMY collaboration. 33rd ICRC, July 2013.
  27. E. Conti, G. Collazuol, and G. Sartori. ArXiv e-prints, August 2014.
  28. J. Alvarez-Muñiz et al. Nucl. Instrum. Meth. A, 719:70–80, August 2013.
  29. R. Gaïor and The Pierre Auger collaboration. 33rd ICRC, July 2013.
  30. F. Werner. 6 ARENA, Annapolis, USA, July 2014.
  31. B. Revenu and V. Marin. 33rd ICRC, number 398, July 2013.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description