Based on extensive air shower simulations it is shown that the electron distributions with respect to the two angles, determining electron direction at a given shower age, for a fixed electron energy and lateral distance, are universal. It means that the distributions do not depend on the primary particle energy or mass (thus, neither on the interaction model), shower zenith angle or shower to shower fluctuations, if they are taken at the same shower age. Together with previous work showing the universality of the distributions of the electron energy, of the lateral distance (integrated over angles) and of the angle (integrated over lateral distance) for fixed electron energy this paper completes a full universal description of the electron states at various shower ages. Analytical parametrizations of the full electron states are given. It is also shown that some distributions can be described by a smaller than five numbers of variables, the new ones being products of the old ones raised to some powers. The accuracy of the present parametrization is sufficiently good for applying to showers with the primary energy uncertainty of 14% (as it is at the Pierre Auger Observatory). The shower fluctuations in the chosen bins of the multidimensional variable space are about 6%, determining the minimum uncertainty needed for parametrization of the universal distributions. An analytical way of estimation of the effect of the geomagnetic field is given. Thanks to the universality of the electron distributions in any shower a new method of shower reconstruction can be worked out from data of the observatories using the fluorescence technique. The light fluxes (both fluorescence and Cherenkov) for any shower age can be exactly predicted for a shower with any primary energy and shower maximum depth, so that the two quantities can be obtained by best fitting the predictions to the measurements

It is shown that the shape of any electron distribution in a high energy air shower is the same in all such showers, if taken at the same age, independently of the primary energy, mass and thus, of the interaction model.

The Pierre Auger Observatory is the world's largest cosmic ray observatory. Our current exposure reaches nearly 40,000 km$^2$ str and provides us with an unprecedented quality data set. The performance and stability of the detectors and their enhancements are described. Data analyses have led to a number of major breakthroughs. Among these we discuss the energy spectrum and the searches for large-scale anisotropies. We present analyses of our X$_{max}$ data and show how it can be interpreted in terms of mass composition. We also describe some new analyses that extract mass sensitive parameters from the 100% duty cycle SD data. A coherent interpretation of all these recent results opens new directions. The consequences regarding the cosmic ray composition and the properties of UHECR sources are briefly discussed.

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

One of the methods for studying the highest energy cosmic rays is to measure the fluorescence light emitted by the extensive air showers induced by them. To reconstruct a shower cascade curve from measurements of the number of photons arriving from the subsequent shower track elements it is necessary to take into account the multiple scatterings that photons undergo on their way from the shower to the detector. In contrast to the earlier Monte-Carlo work, we present here an analytical method to treat the Rayleigh and Mie scatterings in the atmosphere. The method consists in considering separately the consecutive 'generations' of the scattered light. Starting with a point light source in a uniform medium, we then examine a source in a real atmosphere and finally - a moving source (shower) in it. We calculate the angular distributions of the scattered light superimposed on the not scattered light registered from a shower at a given time. The analytical solutions (although approximate) show how the exact numerical results should be parametrised what we do for the first two generations (the contribution of the higher ones being small). Not allowing for the considered effect may lead to an overestimation of shower primary energy by ~15% and to an underestimation of the primary particle mass.

Ultra High Energy (UHE) particles coming from discrete extragalactic sources are potential candidates for EAS events above a few tens of EeV. In particular, galaxies with huge infrared luminosity triggered by collision and merging processes are possible sites of UHECR acceleration. Using the PSCz catalogue of IR galaxies we calculate a large scale anisotropy of UHE protons originating in the population of the luminous infrared galaxies (LIRGs). Small angle particle scattering in weak irregular extragalactic magnetic fields as well as deflection by regular Galactic field are taken into account. We give analytical formulae for deflection angles with included energy losses on cosmic microwave background (CMB). The hypotheses of the anisotropic and isotropic distributions of the experimental data above 40 EeV from AGASA are checked, using various statistical tests. We show that on the basis of the small scale clustering analysis there is a much better correlation of the UHECRs data below GZK cut-off with the predictions of the LIRG origin than with those of isotropy. We derive analytical formulae for a probability of a given number of doublets, triplets and quadruplets for any density distribution of independent events on the sky. The famous AGASA UHE triple event is found to be very well correlated on the sky with the brightest extragalactic infrared source within 70 Mpc - merger galaxies Arp 299 (NGC 3690 + IC 694).