Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

With LOFAR we have been able to image the development of lightning flashes with meter-scale accuracy and unprecedented detail. We discuss the primary steps behind our most recent lightning mapping method. To demonstrate the capabilities of our technique we show and interpret images of the first few milliseconds of two intra-cloud flashes. In all our flashes the negative leaders propagate in the charge layer below the main negative charge. Among several interesting features we show that in about 2~ms after initiation the Primary Initial Leader triggers the formation of a multitude (more than ten) negative leaders in a rather confined area of the atmosphere. From these only one or two continue to propagate after about 30~ms to extend over kilometers horizontally while another may propagate back to the initiation point. We also show that normal negative leaders can transition into an initial-leader like state, potentially in the presence of strong electric fields. In addition, we show some initial breakdown pulses that occurred during the primary initial leader, and even during two "secondary" initial leaders that developed out of stepped leaders.

We use the Low Frequency ARray (LOFAR) to probe the dynamics of the stepping process of negatively-charged plasma channels (negative leaders) in a lightning discharge. We observe that at each step of a leader, multiple pulses of VHF (30~--~80 MHz) radiation are emitted in short-duration bursts ($<10\ \mu$s). This is evidence for streamer formation during corona flashes that occur with each leader step, which has not been observed before in natural lightning and it could help explain X-ray emission from lightning leaders, as X-rays from laboratory leaders tend to be associated with corona flashes. Surprisingly we find that the stepping length is very similar to what was observed near the ground, however with a stepping time that is considerably larger, which as yet is not understood. These results will help to improve lightning propagation models, and eventually lightning protection models.

The limited knowledge of atmospheric parameters like humidity, pressure, temperature, and the index of refraction has been one of the important systematic uncertainties in reconstructing the depth of the shower maximum from the radio emission of air showers. Current air shower Monte Carlo simulation codes like CORSIKA and the radio plug-in CoREAS use various averaged parameterized atmospheres. However, time-dependent and location-specific atmospheric models are needed for the cosmic ray analysis method used for LOFAR data. There, dedicated simulation sets are used for each detected cosmic ray, to take into account the actual atmospheric conditions at the time of the measurement. Using the Global Data Assimilation System (GDAS), a global atmospheric model, we have implemented time-dependent, realistic atmospheric profiles in CORSIKA and CoREAS. We have produced realistic event-specific atmospheres for all air showers measured with LOFAR, an event set spanning several years and many different weather conditions. A complete re-analysis of our data set shows that for the majority of data, our previous correction factor performed rather well; we found only a small systematic shift of 2 g/cm$^2$ in the reconstructed $X_{\rm max}$. However, under extreme weather conditions, for example, very low air pressure, the shift can be up to 15 g/cm$^2$. We provide a correction formula to determine the shift in $X_{\rm max}$ resulting from a comparison of simulations done using the US-Std atmosphere and the GDAS-based atmosphere.

Cosmic rays are routinely measured at LOFAR, both with a dense array of antennas and an array of plastic scintillators. We compare the reconstruction of cosmic-ray energy using radio and particle techniques and show that they are in good agreement. We also present the radiation energy of air showers measured at LOFAR. This value scales quadratically with the electromagnetic energy in an air shower, which can in turn be related to the energy of the primary particle. Once the local magnetic field is accounted for, the radiation energy allows for a direct comparison between the LOFAR energy scale and that of the Pierre Auger Observatory, which are shown to be in agreement. This technique demonstrates the potential to cross-calibrate the energy scales of different experiments using radiation energy.

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

This is a collection of papers that have been contributed by the LOFAR Cosmic Ray Key Science Project (CRKSP) to the 36th International Cosmic Ray Conference held in Madison, Wisconsin, on July 24th to August 1st, 2019 (ICRC 2019). All papers contained here have been individually published in PoS(ICRC2019) with paper numbers 205, 362, 352, 416, and 363, in the order they appear in this collection. Minor modifications to the PoS versions have been applied where appropriate.

A Monte Carlo simulation program for the radio detection of Ultra High Energy (UHE) neutrino interactions in the Antarctic ice as viewed by the Antarctic Impulsive Transient Antenna (ANITA) is described in this article. The program, icemc, provides an input spectrum of UHE neutrinos, the parametrization of the Askaryan radiation generated by their interaction in the ice, and the propagation of the radiation through ice and air to a simulated model of the third and fourth ANITA flights. This paper provides an overview of the icemc simulation, descriptions of the physics models used and of the ANITA electronics processing chain, data/simulation comparisons to validate the predicted performance, and a summary of the impact of published results.

The LOFAR radio telescope is able to measure the radio emission from cosmic ray induced air showers with hundreds of individual antennas. This allows for precision testing of the emission mechanisms for the radio signal as well as determination of the depth of shower maximum $X_{\max}$, the shower observable most sensitive to the mass of the primary cosmic ray, to better than 20 g/cm$^2$. With a densely instrumented circular area of roughly 320 m$^2$, LOFAR is targeting for cosmic ray astrophysics in the energy range $10^{16}$ - $10^{18}$ eV. In this contribution we give an overview of the status, recent results, and future plans of cosmic ray detection with the LOFAR radio telescope.

Cosmic particles hitting Earth's moon produce radio emission via the Askaryan effect. If the resulting radio ns-pulse can be detected by radio telescopes, this technique potentially increases the available collective area for ZeV scale particles by several orders of magnitude compared to current experiments. The LOw Frequency ARray (LOFAR) is the largest radio telescope operating in the optimum frequency regime for this technique. In this contribution, we report on the status of the implementation of the lunar detection mode at LOFAR.

The LOw-Frequency ARray (LOFAR) is used to make precise measurements of radio emission from extensive air showers, yielding information about the primary cosmic ray. Interpreting the measured data requires an absolute and frequency-dependent calibration of the LOFAR system response. This is particularly important for spectral analyses, because the shape of the detected signal holds information about the shower development. We revisit the calibration of the LOFAR antennas in the range of 30 - 80 MHz. Using the Galactic emission and a detailed model of the LOFAR signal chain, we find an improved calibration that provides an absolute energy scale and allows for the study of frequency-dependent features in measured signals. With the new calibration, systematic uncertainties of 13% are reached, and comparisons of the spectral shape of calibrated data with simulations show promising agreement.

The ANtarctic Impulsive Transient Antenna (ANITA) NASA long-duration balloon payload completed its fourth flight in December 2016, after 28 days of flight time. ANITA is sensitive to impulsive broadband radio emission from interactions of ultra-high-energy neutrinos in polar ice (Askaryan emission). We present the results of two separate blind analyses searching for signals from Askaryan emission in the data from the fourth flight of ANITA. The more sensitive analysis, with a better expected limit, has a background estimate of $0.64^{+0.69}_{-0.45}$ and an analysis efficiency of $82\pm2\%$. The second analysis has a background estimate of $0.34^{+0.66}_{-0.16}$ and an analysis efficiency of $71\pm6\%$. Each analysis found one event in the signal region, consistent with the background estimate for each analysis. The resulting limit further tightens the constraints on the diffuse flux of ultra-high-energy neutrinos at energies above $10^{19.5}$ eV.

Recently, the ANITA collaboration reported on two upward-going extensive air shower events consistent with a primary particle that emerges from the surface of the ice. These events may be of $\nu_\tau$ origin, in which the neutrino interacts within the Earth to produce a $\tau$ lepton that emerges from the Earth, decays in the atmosphere, and initiates an extensive air shower. In this paper we estimate an upper bound on the ANITA acceptance to a diffuse $\nu_\tau$ flux detected via $\tau$-lepton-induced air showers within the bounds of Standard Model (SM) uncertainties. By comparing this estimate with the acceptance of Pierre Auger Observatory and IceCube and assuming SM interactions, we conclude that a $\nu_\tau$ origin of these events would imply a neutrino flux at least two orders of magnitude above current bounds.

These proceedings address a recent publication by the ANITA collaboration of four upward- pointing cosmic-ray-like events observed in the first flight of ANITA. Three of these events were consistent with stratospheric cosmic-ray air showers where the axis of propagation does not inter- sect the surface of the Earth. The fourth event was consistent with a primary particle that emerges from the surface of the ice suggesting a possible {\tau}-lepton decay as the origin of this event. These proceedings follow-up on the modeling and testing of the hypothesis that this event was of {\tau} neutrino origin.

The Antarctic Impulsive Transient Antenna (ANITA), a NASA long-duration balloon payload, searches for radio emission from interactions of ultra-high-energy neutrinos in polar ice. The third flight of ANITA (ANITA-III) was launched in December 2014 and completed a 22-day flight. We present the results of three analyses searching for Askaryan radio emission of neutrino origin. In the most sensitive of the analyses, we find one event in the signal region on an expected a priori background of $0.7^{+0.5}_{-0.3}$. Though consistent with the background estimate, the candidate event remains compatible with a neutrino hypothesis even after additional post-unblinding scrutiny.

We report on an upward traveling, radio-detected cosmic-ray-like impulsive event with characteristics closely matching an extensive air shower. This event, observed in the third flight of the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload, is consistent with a similar event reported in a previous flight. These events may be produced by the atmospheric decay of an upward-propagating $\tau$-lepton produced by a $\nu_{\tau}$ interaction, although their relatively steep arrival angles create tension with the standard model (SM) neutrino cross section. Each of the two events have $a~posteriori$ background estimates of $\lesssim 10^{-2}$ events. If these are generated by $\tau$-lepton decay, then either the charged-current $\nu_{\tau}$ cross section is suppressed at EeV energies, or the events arise at moments when the peak flux of a transient neutrino source was much larger than the typical expected cosmogenic background neutrinos.

The balloon-borne HiCal radio-frequency (RF) transmitter, in concert with the ANITA radio-frequency receiver array, is designed to measure the Antarctic surface reflectivity in the RF wavelength regime. The amplitude of surface-reflected transmissions from HiCal, registered as triggered events by ANITA, can be compared with the direct transmissions preceding them by O(10) microseconds, to infer the surface power reflection coefficient $\cal{R}$. The first HiCal mission (HiCal-1, Jan. 2015) yielded a sample of 100 such pairs, resulting in estimates of $\cal{R}$ at highly-glancing angles (i.e., zenith angles approaching $90^\circ$), with measured reflectivity for those events which exceeded extant calculations. The HiCal-2 experiment, flying from Dec., 2016-Jan., 2017, provided an improvement by nearly two orders of magnitude in our event statistics, allowing a considerably more precise mapping of the reflectivity over a wider range of incidence angles. We find general agreement between the HiCal-2 reflectivity results and those obtained with the earlier HiCal-1 mission, as well as estimates from Solar reflections in the radio-frequency regime. In parallel, our calculations of expected reflectivity have matured; herein, we use a plane-wave expansion to estimate the reflectivity R from both a flat, smooth surface (and, in so doing, recover the Fresnel reflectivity equations) and also a curved surface. Multiplying our flat-smooth reflectivity by improved Earth curvature and surface roughness corrections now provides significantly better agreement between theory and the HiCal 2a/2b measurements.

Contributions of the LOFAR Cosmic Ray Key Science Project to the 35th International Cosmic Ray Conference (ICRC 2017)

The Antarctic Impulsive Transient Antenna (ANITA) is a NASA long-duration balloon experiment with the primary goal of detecting ultra-high-energy ($>10^{18}\,\mbox{eV}$) neutrinos via the Askaryan Effect. The fourth ANITA mission, ANITA-IV, recently flew from Dec 2 to Dec 29, 2016. For the first time, the Tunable Universal Filter Frontend (TUFF) boards were deployed for mitigation of narrow-band, anthropogenic noise with tunable, switchable notch filters. The TUFF boards also performed second-stage amplification by approximately 45 dB to boost the $\sim\,\mu\mbox{V-level}$ radio frequency (RF) signals to $\sim$ mV-level for digitization, and supplied power via bias tees to the first-stage, antenna-mounted amplifiers. The other major change in signal processing in ANITA-IV is the resurrection of the $90^{\circ}$ hybrids deployed previously in ANITA-I, in the trigger system, although in this paper we focus on the TUFF boards. During the ANITA-IV mission, the TUFF boards were successfully operated throughout the flight. They contributed to a factor of 2.8 higher total instrument livetime on average in ANITA-IV compared to ANITA-III due to reduction of narrow-band, anthropogenic noise before a trigger decision is made.

The primary science goal of the NASA-sponsored ANITA project is measurement of ultra-high energy neutrinos and cosmic rays, observed via radio-frequency signals resulting from a neutrino- or cosmic ray- interaction with terrestrial matter (atmospheric or ice molecules, e.g.). Accurate inference of the energies of these cosmic rays requires understanding the transmission/reflection of radio wave signals across the ice-air boundary. Satellite-based measurements of Antarctic surface reflectivity, using a co-located transmitter and receiver, have been performed more-or-less continuously for the last few decades. Satellite-based reflectivity surveys, at frequencies ranging from 2--45 GHz and at near-normal incidence, yield generally consistent reflectivity maps across Antarctica. Using the Sun as an RF source, and the ANITA-3 balloon borne radio-frequency antenna array as the RF receiver, we have also measured the surface reflectivity over the interval 200-1000 MHz, at elevation angles of 12-30 degrees, finding agreement with the Fresnel equations within systematic errors. To probe low incidence angles, inaccessible to the Antarctic Solar technique and not probed by previous satellite surveys, a novel experimental approach ("HiCal-1") was devised. Unlike previous measurements, HiCal-ANITA constitute a bi-static transmitter-receiver pair separated by hundreds of kilometers. Data taken with HiCal, between 200--600 MHz shows a significant departure from the Fresnel equations, constant with frequency over that band, with the deficit increasing with obliquity of incidence, which we attribute to the combined effects of possible surface roughness, surface grain effects, radar clutter and/or shadowing of the reflection zone due to Earth curvature effects.

We observe a large fraction of circular polarization in radio emission from extensive air showers recorded during thunderstorms, much higher than in the emission from air showers measured during fair-weather circumstances. We show that the circular polarization of the air showers measured during thunderstorms can be explained by the change in the direction of the transverse current as a function of altitude induced by atmospheric electric fields. Thus by using the full set of Stokes parameters for these events, we obtain a good characterization of the electric fields in thunderclouds. We also measure a large horizontal component of the electric fields in the two events that we have analysed.

As of 2023, the low-frequency part of the Square Kilometre Array will go online in Australia. It will constitute the largest and most powerful low-frequency radio-astronomical observatory to date, and will facilitate a rich science programme in astronomy and astrophysics. With modest engineering changes, it will also be able to measure cosmic rays via the radio emission from extensive air showers. The extreme antenna density and the homogeneous coverage provided by more than 60,000 antennas within an area of one km$^2$ will push radio detection of cosmic rays in the energy range around 10$^{17}$ eV to ultimate precision, with superior capabilities in the reconstruction of arrival direction, energy, and an expected depth-of-shower-maximum resolution of 6~g/cm${^2}$.

For the interpretation of measurements of radio emission from extensive air showers, an important systematic uncertainty arises from natural variations of the atmospheric refractive index $n$. At a given altitude, the refractivity $N=10^6\, (n-1)$ can have relative variations on the order of $10 \%$ depending on temperature, humidity, and air pressure. Typical corrections to be applied to $N$ are about $4\%$. Using CoREAS simulations of radio emission from air showers, we have evaluated the effect of varying $N$ on measurements of the depth of shower maximum $X_{\rm max}$. For an observation band of 30 to 80 MHz, a difference of $4 \%$ in refractivity gives rise to a systematic error in the inferred $X_{\rm max}$ between 3.5 and 11 $\mathrm{g/cm^2}$, for proton showers with zenith angles ranging from 15 to 50 degrees. At higher frequencies, from 120 to 250 MHz, the offset ranges from 10 to 22 $\mathrm{g/cm^2}$. These offsets were found to be proportional to the geometric distance to $X_{\rm max}$. We have compared the results to a simple model based on the Cherenkov angle. For the 120 to 250 MHz band, the model is in qualitative agreement with the simulations. In typical circumstances, we find a slight decrease in $X_{\rm max}$ compared to the default refractivity treatment in CoREAS. While this is within commonly treated systematic uncertainties, accounting for it explicitly improves the accuracy of $X_{\rm max}$ measurements.

The low flux of the ultra-high energy cosmic rays (UHECR) at the highest energies provides a challenge to answer the long standing question about their origin and nature. Even lower fluxes of neutrinos with energies above $10^{22}$ eV are predicted in certain Grand-Unifying-Theories (GUTs) and e.g.\ models for super-heavy dark matter (SHDM). The significant increase in detector volume required to detect these particles can be achieved by searching for the nano-second radio pulses that are emitted when a particle interacts in Earth's moon with current and future radio telescopes. In this contribution we present the design of an online analysis and trigger pipeline for the detection of nano-second pulses with the LOFAR radio telescope. The most important steps of the processing pipeline are digital focusing of the antennas towards the Moon, correction of the signal for ionospheric dispersion, and synthesis of the time-domain signal from the polyphased-filtered signal in frequency domain. The implementation of the pipeline on a GPU/CPU cluster will be discussed together with the computing performance of the prototype.

We report here on a novel analysis of the complete set of four Stokes parameters that uniquely determine the linear and/or circular polarization of the radio signal for an extensive air shower. The observed dependency of the circular polarization on azimuth angle and distance to the shower axis is a clear signature of the interfering contributions from two different radiation mechanisms, a main contribution due to a geomagnetically-induced transverse current and a secondary component due to the build-up of excess charge at the shower front. The data, as measured at LOFAR, agree very well with a calculation from first principles. This opens the possibility to use circular polarization as an investigative tool in the analysis of air shower structure, such as for the determination of atmospheric electric fields.

The low flux of the ultra-high energy cosmic rays (UHECR) at the highest energies provides a challenge to answer the long standing question about their origin and nature. A significant increase in the number of detected UHECR is expected to be achieved by employing Earth's moon as detector, and search for short radio pulses that are emitted when a particle interacts in the lunar rock. Observation of these short pulses with current and future radio telescopes also allows to search for the even lower fluxes of neutrinos with energies above $10^{22}$ eV, that are predicted in certain Grand-Unifying-Theories (GUTs), and e.g. models for super-heavy dark matter (SHDM). In this contribution we present the initial design for such a search with the LOFAR radio telescope.

We report on four radio-detected cosmic-ray (CR) or CR-like events observed with the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload. Two of the four were previously identified as stratospheric CR air showers during the ANITA-I flight. A third stratospheric CR was detected during the ANITA-II flight. Here we report on characteristics these three unusual CR events, which develop nearly horizontally, 20-30~km above the surface of the Earth. In addition, we report on a fourth steeply upward-pointing ANITA-I CR-like radio event which has characteristics consistent with a primary that emerged from the surface of the ice. This suggests a possible $\tau$-lepton decay as the origin of this event, but such an interpretation would require significant suppression of the Standard Model ${\tau}$-neutrino cross section.

For fifty years, cosmic-ray air showers have been detected by their radio emission. We present the first laboratory measurements that validate electrodynamics simulations used in air shower modeling. An experiment at SLAC provides a beam test of radio-frequency (RF) radiation from charged particle cascades in the presence of a magnetic field, a model system of a cosmic-ray air shower. This experiment provides a suite of controlled laboratory measurements to compare to particle-level simulations of RF emission, which are relied upon in ultra-high-energy cosmic-ray air shower detection. We compare simulations to data for intensity, linearity with magnetic field, angular distribution, polarization, and spectral content. In particular, we confirm modern predictions that the magnetically induced emission in a dielectric forms a cone that peaks at the Cherenkov angle and show that the simulations reproduce the data within systematic uncertainties.

The first flight of the Antarctic Impulsive Transient Antenna (ANITA) experiment recorded 16 radio signals that were emitted by cosmic-ray induced air showers. For 14 of these events, this radiation was reflected from the ice. The dominant contribution to the radiation from the deflection of positrons and electrons in the geomagnetic field, which is beamed in the direction of motion of the air shower. This radiation is reflected from the ice and subsequently detected by the ANITA experiment at a flight altitude of 36km. In this paper, we estimate the energy of the 14 individual events and find that the mean energy of the cosmic-ray sample is 2.9 EeV. By simulating the ANITA flight, we calculate its exposure for ultra-high energy cosmic rays. We estimate for the first time the cosmic-ray flux derived only from radio observations. In addition, we find that the Monte Carlo simulation of the ANITA data set is in agreement with the total number of observed events and with the properties of those events.