To exploit the full potential of radio measurements of cosmic-ray air showers at MHz frequencies, a detector timing synchronization within 1 ns is needed. Large distributed radio detector arrays such as the Auger Engineering Radio Array (AERA) rely on timing via the Global Positioning System (GPS) for the synchronization of individual detector station clocks. Unfortunately, GPS timing is expected to have an accuracy no better than about 5 ns. In practice, in particular in AERA, the GPS clocks exhibit drifts on the order of tens of ns. We developed a technique to correct for the GPS drifts, and an independent method is used to cross-check that indeed we reach a nanosecond-scale timing accuracy by this correction. First, we operate a ``beacon transmitter'' which emits defined sine waves detected by AERA antennas recorded within the physics data. The relative phasing of these sine waves can be used to correct for GPS clock drifts. In addition to this, we observe radio pulses emitted by commercial airplanes, the position of which we determine in real time from Automatic Dependent Surveillance Broadcasts intercepted with a software-defined radio. From the known source location and the measured arrival times of the pulses we determine relative timing offsets between radio detector stations. We demonstrate with a combined analysis that the two methods give a consistent timing calibration with an accuracy of 2 ns or better. Consequently, the beacon method alone can be used in the future to continuously determine and correct for GPS clock drifts in each individual event measured by AERA.

Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to extend its range of detection and to directly measure the muon content of the particle showers. It consists of an infill of surface water-Cherenkov detectors accompanied by buried scintillator detectors used for muon counting. The main objectives of the AMIGA engineering array, referred to as the Unitary Cell, are to identify and resolve all engineering issues as well as to understand the muon-number counting uncertainties related to the design of the detector. The mechanical design, fabrication and deployment processes of the muon counters of the Unitary Cell are described in this document. These muon counters modules comprise sealed PVC casings containing plastic scintillation bars, wavelength-shifter optical fibers, 64 pixel photomultiplier tubes, and acquisition electronics. The modules are buried approximately 2.25 m below ground level in order to minimize contamination from electromagnetic shower particles. The mechanical setup, which allows access to the electronics for maintenance, is also described in addition to tests of the modules' response and integrity. The completed Unitary Cell has measured a number of air showers of which a first analysis of a sample event is included here.

Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E-CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +/- 0.16 (1.61 +/- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec \theta)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec \theta)_\mathrm{max}$.

The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.

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.

We present the first hybrid measurement of the average muon number in air showers at ultrahigh energies, initiated by cosmic rays with zenith angles between 62 degrees and 80 degrees. The measurement is based on 174 hybrid events recorded simultaneously with the surface detector array and the fluorescence detector of the Pierre Auger Observatory. The muon number for each shower is derived by scaling a simulated reference profile of the lateral muon density distribution at the ground until it fits the data. A 10(19) eV shower with a zenith angle of 67 degrees, which arrives at the surface detector array at an altitude of 1450 m above sea level, contains on average (2.68 +/- 0.04 +/- 0.48(sys)) x 10(7) muons with energies larger than 0.3 GeV. The logarithmic gain d lnN mu/d ln E of muons with increasing energy between 4 x 10(18) eV and 5 x 10(19) eV is measured to be (1.029 +/- 0.024) 0.030(sys)).

The energy spectrum above 4×10 eV is presented, obtained from 5936 events with zenith angles exceeding 62◦ collected by the Surface Detector of the Pierre Auger Observatory from 1 January 2004 to 31 December 2010. Showers with such large zenith angles are muon-dominated at ground level and the radial symmetry around the shower axis is broken due to geomagnetic deflections. They are analysed separately from showers with smaller zenith angles using two-dimensional models of the muon density at ground, allowing one to reconstruct a global muon number for every event. The conversion of the muon number to energy is obtained using the sub-sample of events detected simultaneously with both the Surface and the Fluorescence Detector. The spectrum obtained displays suppression near 4 × 10 eV compatible with the analysis that uses less inclined events.

The optical calibration of the fluorescence telescopes is a significant contribution to the overall uncertainty of energy measurements made by the Pierre Auger Observatory and Telescope Array Project. Some sources of uncertainty, such as the fluorescence yield in air, affect both experiments similarly. However, the optical calibration of the fluorescence telescopes is a source of independent uncertainty. The Pierre Auger and Telescope Array collaborations have taken initial steps to establish a relative end-to-end optical calibration of the fluorescence telescopes. An Octocopter carrying a portable light source has been flown in front of fluorescence telescopes at both Pierre Auger and Telescope Array sites. Laboratory calibration measurements of the light source before and after the flights provide a common baseline for the relative end-to-end calibration. We expect this system will lead to a common photonic calibration for both experiments. After giving a brief description of the UV light source and the Octocopter used for the measurements, the parameters and the calibration procedures for the light source will be discussed. First results on telescope images of the light source will be presented.

We present the results of an analysis of the large angular scale distribution of the arrival directions of cosmic rays with energy above 4 EeV detected at the Pierre Auger Observatory including for the first time events with zenith angle between 60° and 80°. We perform two Rayleigh analyses, one in the right ascension and one in the azimuth angle distributions, that are sensitive to modulations in right ascension and declination, respectively. The largest departure from isotropy appears in the E > 8 ?> EeV energy bin, with an amplitude for the first harmonic in right ascension r 1 &agr; = ( 4.4 ± 1.0 ) × 10 − 2 ?> , that has a chance probability P ( ≥ r 1 &agr; ) = 6.4 × 10 − 5 ?> , reinforcing the hint previously reported with vertical events alone.

This paper briefly reviews the current understanding we have of the highest observed energy cosmic rays, being particles with energy in excess of 5×1019eV, which hit us from unknown locations in the cosmos, at very low fluxes. I then describe the Pierre Auger project, whose goal is to find the nature and source of these macroscopic-energy particles.

Joint contributions of the Pierre Auger Collaboration and the Telescope Array Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than 60 detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

A thorough search of the sky exposed at the Pierre Auger Cosmic Ray Observatory reveals no statistically significant excess of events in any small solid angle that would be indicative of a flux of neutral particles from a discrete source. The search covers from −90° to +15° in declination using four different energy ranges above 1 EeV (1018 eV). The method used in this search is more sensitive to neutrons than to photons. The upper limit on a neutron flux is derived for a dense grid of directions for each of the four energy ranges. These results constrain scenarios for the production of ultrahigh energy cosmic rays in the Galaxy.

Using the data taken at the Pierre Auger Observatory between December 2004 and December 2012, we have examined the implications of the distributions of depths of atmospheric shower maximum (Xmax), using a hybrid technique, for composition and hadronic interaction models. We do this by fitting the distributions with predictions from a variety of hadronic interaction models for variations in the composition of the primary cosmic rays and examining the quality of the fit. Regardless of what interaction model is assumed, we find that our data are not well described by a mix of protons and iron nuclei over most of the energy range. Acceptable fits can be obtained when intermediate masses are included, and when this is done consistent results for the proton and iron-nuclei contributions can be found using the available models. We observe a strong energy dependence of the resulting proton fractions, and find no support from any of the models for a significant contribution from iron nuclei. However, we also observe a significant disagreement between the models with respect to the relative contributions of the intermediate components.

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

We report a study of the distributions of the depth of maximum, Xmax, of extensive air-shower profiles with energies above 10^17.8 eV as observed with the fluorescence telescopes of the Pierre Auger Observatory. The analysis method for selecting a data sample with minimal sampling bias is described in detail as well as the experimental cross-checks and systematic uncertainties. Furthermore, we discuss the detector acceptance and the resolution of the Xmax measurement and provide parameterizations thereof as a function of energy. The energy dependence of the mean and standard deviation of the Xmax-distributions are compared to air-shower simulations for different nuclear primaries and interpreted in terms of the mean and variance of the logarithmic mass distribution at the top of the atmosphere.

The emission of radio waves from air showers has been attributed to the so-called geomagnetic emission process. At frequencies around 50 MHz this process leads to coherent radiation which can be observed with rather simple setups. The direction of the electric field induced by this emission process depends only on the local magnetic field vector and on the incoming direction of the air shower. We report on measurements of the electric field vector where, in addition to this geomagnetic component, another component has been observed that cannot be described by the geomagnetic emission process. The data provide strong evidence that the other electric field component is polarized radially with respect to the shower axis, in agreement with predictions made by Askaryan who described radio emission from particle showers due to a negative charge excess in the front of the shower. Our results are compared to calculations which include the radiation mechanism induced by this charge-excess process.