We present the first hybrid measurement of the average muon number in air showers at ultra-high energies, initiated by cosmic rays with zenith angles between $62^\circ$ and $80^\circ$. 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^\circ$, which arrives at the Surface Detector array at an altitude of 1450 m above sea level, contains on average $(2.68 \pm 0.04 \pm 0.48\,\mathrm{(sys.)}) \times 10^{7}$ muons with energies larger than 0.3 GeV. The logarithmic gain $\mathrm{d}\ln{N_\mu} / \mathrm{d}\ln{E}$ of muons with increasing energy between $4\times 10^{18}$ eV and $5 \times 10^{19}$ eV is measured to be $(1.029\, \pm\, 0.024\, \pm 0.030\,\mathrm{(sys.)})$.

We analyze the distribution of arrival directions of ultra-high energy cosmic rays recorded at the Pierre Auger Observatory in 10 years of operation. The data set, about three times larger than that used in earlier studies, includes arrival directions with zenith angles up to $80^\circ$, thus covering from $-90^\circ$ to $+45^\circ$ in declination. After updating the fraction of events correlating with the active galactic nuclei (AGNs) in the V\'eron-Cetty and V\'eron catalog, we subject the arrival directions of the data with energies in excess of 40 EeV to different tests for anisotropy. We search for localized excess fluxes and for self-clustering of event directions at angular scales up to $30^\circ$ and for different threshold energies between 40~EeV and 80~EeV. We then look for correlations of cosmic rays with celestial structures both in the Galaxy (the Galactic Center and Galactic Plane) and in the local Universe (the Super-Galactic Plane). We also examine their correlation with different populations of nearby extragalactic objects: galaxies in the 2MRS catalog, AGNs detected by Swift-BAT, radio galaxies with jets and the Centaurus~A galaxy. None of the tests shows a statistically significant evidence of anisotropy. The strongest departures from isotropy (post-trial probability ${\sim}1.4$\%) are obtained for cosmic rays with $E>58$~EeV in rather large windows around Swift AGNs closer than 130~Mpc and brighter than $10^{44}$~erg/s (18$^\circ$ radius) and around the direction of Centaurus~A (15$^\circ$ radius).

The surface detector array of the Pierre Auger Observatory provides information about the longitudinal development of the muonic component of extensive air showers. Using the timing information from the flash analog-to-digital converter traces of surface detectors far from the shower core, it is possible to reconstruct a muon production depth distribution. We characterize the goodness of this reconstruction for zenith angles around 60 deg. and different energies of the primary particle. From these distributions we define X(mu)max as the depth along the shower axis where the production of muons reaches maximum. We explore the potentiality of X(mu)max as a useful observable to infer the mass composition of ultrahigh-energy cosmic rays. Likewise, we assess its ability to constrain hadronic interaction models.

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.

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.

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85{\deg} to +20{\deg}, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

The Pierre Auger Observatory is making significant contributions towards understanding the nature and origin of ultra-high energy cosmic rays. One of its main challenges is the monitoring of the atmosphere, both in terms of its state variables and its optical properties. The aim of this work is to analyze aerosol optical depth $\tau_{\rm a}(z)$ values measured from 2004 to 2012 at the observatory, which is located in a remote and relatively unstudied area of the Pampa Amarilla, Argentina. The aerosol optical depth is in average quite low - annual mean $\tau_{\rm a}(3.5~{\rm km})\sim 0.04$ - and shows a seasonal trend with a winter minimum - $\tau_{\rm a}(3.5~{\rm km})\sim 0.03$ -, and a summer maximum - $\tau_{\rm a}(3.5~{\rm km})\sim 0.06$ -, and an unexpected increase from August to September - $\tau_{\rm a}(3.5~{\rm km})\sim 0.055$). We computed backward trajectories for the years 2005 to 2012 to interpret the air mass origin. Winter nights with low aerosol concentrations show air masses originating from the Pacific Ocean. Average concentrations are affected by continental sources (wind-blown dust and urban pollution), while the peak observed in September and October could be linked to biomass burning in the northern part of Argentina or air pollution coming from surrounding urban areas.

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 which 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.

We describe a new method of identifying night-time clouds over the Pierre Auger Observatory using infrared data from the Imager instruments on the GOES-12 and GOES-13 satellites. We compare cloud identifications resulting from our method to those obtained by the Central Laser Facility of the Auger Observatory. Using our new method we can now develop cloud probability maps for the 3000 km^2 of the Pierre Auger Observatory twice per hour with a spatial resolution of ~2.4 km by ~5.5 km. Our method could also be applied to monitor cloud cover for other ground-based observatories and for space-based observatories.

The observation of ultrahigh energy (UHE) neutrinos has become a priority in experimental astroparticle physics. UHE neutrinos can be detected with a variety of techniques. In particular, neutrinos can interact in the atmosphere (downward-going neutrinos) or in the Earth crust (Earth-skimming neutrinos), producing air showers that can be observed with arrays of detectors at the ground. With the Surface Detector Array of the Pierre Auger Observatory we can detect these types of cascades. The distinguishing signature for neutrino events is the presence of very inclined showers produced close to the ground (i.e. after having traversed a large amount of atmosphere). In this work we review the procedure and criteria established to search for UHE neutrinos in the data collected with the ground array of the Pierre Auger Observatory. This includes Earth-skimming as well as downward-going neutrinos. No neutrino candidates have been found, which allows us to place competitive limits to the diffuse flux of UHE neutrinos in the EeV range and above.

The Pierre Auger Observatory in Malarg\"ue, Argentina, is designed to study the properties of ultra-high energy cosmic rays with energies above 1018 eV. It is a hybrid facility that employs a Fluorescence Detector to perform nearly calorimetric measurements of Extensive Air Shower energies. To obtain reliable calorimetric information from the FD, the atmospheric conditions at the observatory need to be continuously monitored during data acquisition. In particular, light attenuation due to aerosols is an important atmospheric correction. The aerosol concentration is highly variable, so that the aerosol attenuation needs to be evaluated hourly. We use light from the Central Laser Facility, located near the center of the observatory site, having an optical signature comparable to that of the highest energy showers detected by the FD. This paper presents two procedures developed to retrieve the aerosol attenuation of fluorescence light from CLF laser shots. Cross checks between the two methods demonstrate that results from both analyses are compatible, and that the uncertainties are well understood. The measurements of the aerosol attenuation provided by the two procedures are currently used at the Pierre Auger Observatory to reconstruct air shower data.

To interpret the mean depth of cosmic ray air shower maximum and its dispersion, we parametrize those two observables as functions of the first two moments of the $\ln A$ distribution. We examine the goodness of this simple method through simulations of test mass distributions. The application of the parameterization to Pierre Auger Observatory data allows one to study the energy dependence of the mean $\ln A$ and of its variance under the assumption of selected hadronic interaction models. We discuss possible implications of these dependences in term of interaction models and astrophysical cosmic ray sources.

A thorough search for large scale anisotropies in the distribution of arrival directions of cosmic rays detected above $10^{18}$ eV at the Pierre Auger Observatory is reported. For the first time, these large scale anisotropy searches are performed as a function of both the right ascension and the declination and expressed in terms of dipole and quadrupole moments. Within the systematic uncertainties, no significant deviation from isotropy is revealed. Upper limits on dipole and quadrupole amplitudes are derived under the hypothesis that any cosmic ray anisotropy is dominated by such moments in this energy range. These upper limits provide constraints on the production of cosmic rays above $10^{18}$ eV, since they allow us to challenge an origin from stationary galactic sources densely distributed in the galactic disk and emitting predominantly light particles in all directions.

A thorough search for large scale anisotropies in the distribution of arrival directions of cosmic rays detected above $10^{18}$ eV at the Pierre Auger Observatory is presented. This search is performed as a function of both declination and right ascension in several energy ranges above $10^{18}$ eV, and reported in terms of dipolar and quadrupolar coefficients. Within the systematic uncertainties, no significant deviation from isotropy is revealed. Assuming that any cosmic ray anisotropy is dominated by dipole and quadrupole moments in this energy range, upper limits on their amplitudes are derived. These upper limits allow us to challenge an origin of cosmic rays above $10^{18}$ eV from stationary galactic sources densely distributed in the galactic disk and emitting predominantly light particles in all directions.

Observations of cosmic ray arrival directions made with the Pierre Auger Observatory have previously provided evidence of anisotropy at the 99% CL using the correlation of ultra high energy cosmic rays (UHECRs) with objects drawn from the Veron-Cetty Veron catalog. In this paper we report on the use of three catalog independent methods to search for anisotropy. The 2pt-L, 2pt+ and 3pt methods, each giving a different measure of self-clustering in arrival directions, were tested on mock cosmic ray data sets to study the impacts of sample size and magnetic smearing on their results, accounting for both angular and energy resolutions. If the sources of UHECRs follow the same large scale structure as ordinary galaxies in the local Universe and if UHECRs are deflected no more than a few degrees, a study of mock maps suggests that these three methods can efficiently respond to the resulting anisotropy with a P-value = 1.0% or smaller with data sets as few as 100 events. Using data taken from January 1, 2004 to July 31, 2010 we examined the 20, 30,..., 110 highest energy events with a corresponding minimum energy threshold of about 49.3 EeV. The minimum P-values found were 13.5% using the 2pt-L method, 1.0% using the 2pt+ method and 1.1% using the 3pt method for the highest 100 energy events. In view of the multiple (correlated) scans performed on the data set, these catalog-independent methods do not yield strong evidence of anisotropy in the highest energy cosmic rays.

The Surface Detector array of the Pierre Auger Observatory can detect neutrinos with energy between 10^17 eV and 10^20 eV from point-like sources across the sky south of +55 deg and north of -65 deg declinations. A search has been performed for highly inclined extensive air showers produced by the interaction of neutrinos of all flavours in the atmosphere (downward-going neutrinos), and by the decay of tau leptons originating from tau neutrinos interactions in the Earth's crust (Earth-skimming neutrinos). No candidate neutrinos have been found in data up to 2010 May 31. This corresponds to an equivalent exposure of ~3.5 years of a full surface detector array for the Earth-skimming channel and ~2 years for the downward-going channel. An improved upper limit on the diffuse flux of tau neutrinos has been derived. Upper limits on the neutrino flux from point-like sources have been derived as a function of the source declination. Assuming a differential neutrino flux k_PS E^-2 from a point-like source, 90% C.L. upper limits for k_PS at the level of ~5 x 10^-7 and 2.5 x 10^-6 GeV cm^-2 s^-1 have been obtained over a broad range of declinations from the searches of Earth-skimming and downward-going neutrinos, respectively.

We report a measurement of the proton-air cross-section for particle production at the center-of-mass energy per nucleon of 57 TeV. This is derived from the distribution of the depths of shower maxima observed with the Pierre Auger Observatory: systematic uncertainties are studied in detail. Analysing the tail of the distribution of the shower maxima, a proton-air cross-section of $[505\pm22(stat)^{+28}_{-36}(sys)]$ mb is found.

The Pierre Auger Observatory is planned to be upgraded so that the energy spectrum of cosmic rays can be studied down to 0.1 EeV and the muon component of showers can be determined. The former will lead to a spectrum measured by one technique from 0.1 EeV to beyond 100 EeV while the latter will aid identification of the primary particles. These enhancements consist of three high elevation telescopes (HEAT) and an infilled area having both surface detectors and underground muon counters (AMIGA). The surface array of the Auger Observatory will be enhanced over a 23.5 km2 area by 85 detector pairs laid out as a graded array of water-Cherenkov detectors and 30 m2 buried muon scintillator counters. The spacings in the array will be 433 and 750 m. The muon detectors will comprise highly segmented scintillators with optical fibres ending on multi-anode phototubes. The AMIGA complex will be centred 6.0 km away from the fluorescence detector installation at Coihueco and will be overlooked by the HEAT telescopes. We describe the design features of the AMIGA enhancement.

An upper limit of 16% (at 95% c.l.) is derived for the photon fraction in cosmic rays with energies greater than 10^19 eV, based on observations of the depth of shower maximum performed with the hybrid detector of the Pierre Auger Observatory. This is the first such limit on photons obtained by observing the fluorescence light profile of air showers. This upper limit confirms and improves on previous results from the Haverah Park and AGASA surface arrays. Additional data recorded with the Auger surface detectors for a subset of the event sample, support the conclusion that a photon origin of the observed events is not favored.

The stability of the fluorescence telescopes of the Pierre Auger Observatory is monitored with the optical relative calibration setup. Optical fibers distribute light pulses to three different diffuser groups within the optical system. The total charge per pulse is measured for each pixel and compared with reference calibration measurements. This allows monitoring the short and long term stability with respect of the relative timing between pixels and the relative gain for each pixel. The designs of the LED calibration unit (LCU) and of the Xenon flash lamp used for relative calibration, are described and their capabilities to monitor the stability of the telescope performances are studied. We report the analysis of relative calibration data recorded during 2004. Fluctuations in the relative calibration constants provide a measure of the stability of the FD.

We discuss the statistical and systematic uncertainties in the event reconstruction (core location, and determination of S(1000), i.e., the signal at a distance of 1000 m from the shower core) by the Pierre Auger surface detector for showers with zenith angle less than 60 degrees. The method is based on a maximum likelihood method where the reference lateral distribution function is obtained through the experimental data. We also discuss S(1000) as primary energy estimator.