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.

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.

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$\sigma$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $\alpha_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

We present an in-depth and systematic variability study of a sample of 20 powerful blazars, including 12 BL Lacs and 8 flat spectrum radio quasars, applying various analysis tools such as flux distribution, symmetry analysis, and time series analysis on the decade-long Fermi/LAT observations. The results show that blazars with steeper \gama-ray spectral indexes are found to be more variable; and the \gama-ray flux distribution closely resembles lognormal probability distribution function. The statistical variability properties of the sources as studied by power spectral density analysis are consistent with \emph{flicker noise} ($P(\nu)\propto1/\nu$) -- an indication of long-memory processes at work. Statistical analysis of the distribution of flux rise and decay rates in the light curves of the sources, aimed at distinguishing between particle acceleration and energy dissipation timescales, counter-intuitively suggests that both kinds of rates follow a similar distribution and the derived mean variability timescales are in the order of a few weeks. The corresponding emission region size is used to constrain location of \gama-ray production sites in the sources to be a few parsecs. Additionally, using Lomb-Scargle periodogram and weighted wavelet z-transform methods and extensive Monte Carlo simulations, we detected year timescale quasi-periodic oscillations in the sources S5 0716+714, Mrk 421, ON +325, PKS 1424-418 and PKS 2155-304. The detection significance was computed taking proper account of the red-noise and other artifacts inherent in the observations. We explain the results in the light of current blazar models with relativistic shocks propagating down the jet viewed close to the line of sight.

With the Surface Detector array (SD) of the Pierre Auger Observatory we can detect neutrinos with energy between $10^{17}\,$eV and $10^{20}\,$eV from point-like sources across the sky, from close to the Southern Celestial Pole up to $60^\circ$ in declination, with peak sensitivities at declinations around $\sim -53^\circ$ and $\sim+55^\circ$, and an unmatched sensitivity for arrival directions in the Northern hemisphere. A search has been performed for highly-inclined air showers induced by neutrinos of all flavours with no candidate events found in data taken between 1 Jan 2004 and 31 Aug 2018. Upper limits on the neutrino flux from point-like steady sources have been derived as a function of source declination. An unrivaled sensitivity is achieved in searches for transient sources with emission lasting over an hour or less, if they occur within the field of view corresponding to the zenith angle range between $60^\circ$ and $~95^\circ$ where the SD of the Pierre Auger Observatory is most sensitive to neutrinos.

The determination of the primary energy of extensive air showers using the fluorescence detection technique requires an estimation of the energy carried away by particles that do not deposit all their energy in the atmosphere. This estimation is typically made using Monte Carlo simulations and thus depends on the assumed primary particle mass and on model predictions for neutrino and muon production. In this work we present a new method to obtain the invisible energy from events detected by the Pierre Auger Observatory. The method uses measurements of the muon number at ground level, and it allows us to reduce significantly the systematic uncertainties related to the mass composition and the high energy hadronic interaction models, and consequently to improve the estimation of the energy scale of the Observatory.

Neutrinos with energies above $10^{17}$ eV are detectable with the Surface Detector Array of the Pierre Auger Observatory. The identification is efficiently performed for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for Earth-skimming $\tau$ neutrinos with nearly tangential trajectories relative to the earth. No neutrino candidates were found in $\sim\,14.7$ years of data taken up to 31 August 2018. This leads to restrictive upper bounds on their flux. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high-energy neutrinos with an $E_\nu^{-2}$ spectrum in the energy range $1.0 \times 10^{17}~{\rm eV} - 2.5 \times 10^{19}~{\rm eV}$ is $E^2 {\rm d}N_\nu/{\rm d}E_\nu < 4.4 \times 10^{-9}~{\rm GeV~cm^{-2}~s^{-1}~sr^{-1}}$, placing strong constraints on several models of neutrino production at EeV energies and on the properties of the sources of ultra-high-energy cosmic rays.

Contributions of the Pierre Auger Collaboration to the 36th International Cosmic Ray Conference (ICRC 2019), 24 July - 1 August 2019, Madison, Wisconsin, USA.

The search for Cosmic-Ray Ensembles (CRE), groups of correlated cosmic rays that might be distributed over very large areas, even of the size of the planet, requires a globally spread and dense network of detectors, as proposed by the Cosmic-Ray Extremely Distributed Observatory (CREDO) Collaboration. This proposal motivates an effort towards exploring the potential of using even very much diversified detection technologies within one system, with detection units located even in hard-to-reach places, where, nevertheless, the sensors could work independently - without human intervention. For these reasons we have developed a dedicated communication solution enabling the connection of many different types of detectors, in a range of environments. The proposed data transmission system uses radio waves as an information carrier on the 169MHz frequency band in a contrast to the typical commercially used frequencies in a IoT systems (868MHz). The connectivity within the system is based on the star topology, which ensures the least energy consumption. The solution is now being prepared to being implemented using the prototype detection system based on the CosmicWatch open hardware design: a portable, pocket size, and economy particle detector using the scintillation technique. Our prototype detector is equipped with a dedicated software that integrates it with the already operational CREDO server system.

Ultra-high energy photons constitute one of the most important pieces of the astroparticle physics problems. Their observation may provide new insight on several phenomena such as supermassive particle annihilation or the GZK effect. Because of the absence of any significant photon identification by a leading experiments such as the Pierre Auger Observatory, we consider a screening phenomenon called preshower effect which could efficiently affect ultra-high energy photon propagation. This effect is a consequence of photon interactions with the geomagnetic field and results in large electromagnetic cascade of particles several thousands kilometers above the atmosphere. This collection of particles, called cosmic-ray ensembles (CRE), may reach the atmosphere and produce the well-known air showers. In this paper we propose to use gamma-ray telescopes to look for air showers induced by CRE. Possible sources of ultra-high energy photons include the GZK effect and Super Heavy Dark Matter particles. Simulations involving the preshower effect and detectors response are performed and properties of these peculiar air showers are investigated. The use of boosted decision trees to obtain the best cosmic-ray ensemble/hadron separation, the aperture and event rate predictions for a few models of photon production are also presented.

The Cosmic-Ray Extremely Distributed Observatory (CREDO) is a project dedicated to global studies of extremely extended cosmic-ray phenomena, the cosmic-ray ensembles (CRE), beyond the capabilities of existing detectors and observatories. Up to date cosmic-ray research has been focused on detecting single air showers, while the search for ensembles of cosmic-rays, which may spread over a significant fraction of the Earth, is a scientific terra incognita. The key idea of CREDO is to combine existing cosmic-ray detectors (large professional arrays, educational instruments, individual detectors, such as smartphones, etc.) into a worldwide network, thus enabling a global analysis. The second goal of CREDO involves a large number of participants (citizen science!), assuring the geographical spread of the detectors and managing manpower necessary to deal with vast amount of data to search for evidence for cosmic-ray ensembles. In this paper the status and perspectives of the project are presented.

Propagation of ultra-high energy photons in the galactic and intergalactic space gives rise to cascades comprising thousands of photons. Using Monte Carlo simulations, we investigate the development of such cascades in the solar magnetosphere, and find that the photons in the cascades are distributed over hundreds of kilometers as they arrive at the top of the Earth's atmosphere. We also perform similar study for cascades starting as far as 10 Mpc away from us using relevant magnetic field models. A few photons correlated in time are expected to arrive at the Earth from the latter type of cascade. We present our simulation results and discuss the prospects for detection of these cascades with the Cosmic-Ray Extremely Distributed Observatory.

Propagation of ultra-high energy photons in the solar magnetosphere gives rise to cascades comprising thousands of photons. We study the cascade development using Monte Carlo simulations and find that the photons in the cascades are spatially extended over hundreds of kilometers as they arrive at the top of the Earth's atmosphere. We compare results from simulations which use two models of the solar magnetic field, and show that although signatures of such cascades are different for the models used, for practical detection purpose in the ground-based detectors, they are similar.

The profile of the longitudinal development of showers produced by ultra-high energy cosmic rays carries information related to the interaction properties of the primary particles with atmospheric nuclei. In this work, we present the first measurement of the average shower profile in traversed atmospheric depth at the Pierre Auger Observatory. The shapes of profiles are well reproduced by the Gaisser-Hillas parametrization within the range studied, for E > 10^{17.8} eV. A detailed analysis of the systematic uncertainties is performed using 10 years of data and a full detector simulation. The average shape is quantified using two variables related to the width and asymmetry of the profile, and the results are compared with predictions of hadronic interaction models for different primary particles.

An overview of the multi-messenger capabilities of the Pierre Auger Observatory is presented. The techniques and performance of searching for Ultra-High Energy neutrinos, photons and neutrons are described. Some of the most relevant results are reviewed, such as stringent upper bounds that were placed to a flux of diffuse cosmogenic neutrinos and photons, bounds placed on neutrinos emitted from compact binary mergers that were detected by LIGO and Virgo during their first and second observing runs, as well as searches for high energy photons and neutrons from the Galactic center that constrain the properties of the putative Galactic PeVatron, observed by the H.E.S.S.\ collaboration. The observation of directional correlations between ultra-high energy cosmic rays and either high energy astrophysical neutrinos or specific source populations, weighted by their electromagnetic radiation, are also discussed. They constitute additional multi-messenger approaches aimed at identifying the sources of high energy cosmic rays.

We present a detailed study of the large-scale anisotropies of cosmic rays with energies above 4 EeV measured using the Pierre Auger Observatory. For the energy bins [4,8] EeV and $E\geq 8$ EeV, the most significant signal is a dipolar modulation in right ascension at energies above 8 EeV, as previously reported. In this paper we further scrutinize the highest-energy bin by splitting it into three energy ranges. We find that the amplitude of the dipole increases with energy above 4 EeV. The growth can be fitted with a power law with index $\beta=0.79\pm 0.19$. The directions of the dipoles are consistent with an extragalactic origin of these anisotropies at all the energies considered. Additionally we have estimated the quadrupolar components of the anisotropy: they are not statistically significant. We discuss the results in the context of the predictions from different models for the distribution of ultrahigh-energy sources and cosmic magnetic fields.

One of the main objectives of cosmic-ray studies are precise measurements of energy and chemical composition of particles with extreme energies. Large and sophisticated detectors are used to find events seen as showers starting in the Earth's atmosphere with recorded energies larger than 100 EeV. However, a Cosmic-Ray Ensemble (CRE) developing before reaching the Earth as a bunch of correlated particles may spread over larger areas and requires an extended set of detectors to be discovered. The Cosmic-Ray Extremely Distributed Observatory (CREDO) is a solution to find such phenomena. Even simple detectors measuring the particle arrival time only are useful in this approach, as they are sufficient both to provide candidate CRE events and to determine the direction from which they are arriving.

The Cosmic-Ray Extremely Distributed Observatory (CREDO) is a project dedicated to global studies of extremely extended cosmic-ray phenomena, the cosmic-ray ensembles (CRE), beyond the capabilities of existing detectors and observatories. Up to date cosmic-ray research has been focused on detecting single air showers, while the search for ensembles of cosmic-rays, which may overspread a significant fraction of the Earth, is a scientific terra incognita. Instead of developing and commissioning a completely new global detector infrastructure, CREDO proposes approaching the global cosmic-ray analysis objectives with all types of available detectors, from professional to pocket size, merged into a worldwide network. With such a network it is possible to search for evidences of correlated cosmic-ray ensembles. One of the observables that can be investigated in CREDO is a number of spatially isolated events collected in a small time window which could shed light on fundamental physics issues. The CREDO mission and strategy requires active engagement of a large number of participants, also non-experts, who will contribute to the project by using common electronic devices (e.g. smartphones). In this note the status and perspectives of the project is presented.

In the past few years, cosmic-rays beyond the GZK cut-off ($E > 5 \times 10^{19}$ eV) have been detected by leading collaborations such as Pierre Auger Observatory. Such observations raise many questions as to how such energies can be reached and what source can possibly produce them. Although at lower energies, mechanisms such as Fermi acceleration in supernovae front shocks seem to be favored, top-down scenarios have been proposed to explain the existence of ultra-high energy cosmic-rays: the decay of super-massive long-lived particles produced in the early Universe may yield to a flux of ultra-high energy photons. Such photons might be presently generating so called super-preshowers, an extended cosmic-ray shower with a spatial distribution that can be as wide as the Earth diameter. The Cosmic Ray Extremely Distributed Observatory (CREDO) mission is to find such events by means of a network of detectors spread around the globe. CREDO's strategy is to connect existing detectors and create a worldwide network of cosmic-ray observatories. Moreover, citizen-science constitutes an important pillar of our approach. By helping our algorithms to recognize detection patterns and by using smartphones as individual cosmic-ray detectors, non-scientists can participate in scientific discoveries and help unravel some of the deepest mysteries in physics.

With the Auger Engineering Radio Array (AERA) of the Pierre Auger Observatory, we have observed the radio emission from 561 extensive air showers with zenith angles between 60$^\circ$ and 84$^\circ$. In contrast to air showers with more vertical incidence, these inclined air showers illuminate large ground areas of several km$^2$ with radio signals detectable in the 30 to 80\,MHz band. A comparison of the measured radio-signal amplitudes with Monte Carlo simulations of a subset of 50 events for which we reconstruct the energy using the Auger surface detector shows agreement within the uncertainties of the current analysis. As expected for forward-beamed radio emission undergoing no significant absorption or scattering in the atmosphere, the area illuminated by radio signals grows with the zenith angle of the air shower. Inclined air showers with EeV energies are thus measurable with sparse radio-antenna arrays with grid sizes of a km or more. This is particularly attractive as radio detection provides direct access to the energy in the electromagnetic cascade of an air shower, which in case of inclined air showers is not accessible by arrays of particle detectors on the ground.

Although the photon structure is most efficiently studied with the accelerator instruments, there is also a scientifically complementary potential in investigations on photons produced in the outer space. This potential is already being explored with gamma ray telescopes, ultra-high energy cosmic ray observatories and, since very recently, by the Cosmic-Ray Extremely Distributed Observatory (CREDO). Unlike the former instruments focused on detection of single photons, CREDO aims at the detection of cascades (ensembles) of photons originating even at astrophysical distances. If at least a part of such a cascade reaches Earth, it might produce a unique pattern composed of a number of air showers observable by an appropriately dense array of standard detectors. If the energies of air showers constituting the pattern are relatively low and if the typical distances between the neighbors are large, the ensemble character of the whole phenomenon might remain uncovered, unless the CREDO strategy is implemented.

Any considerations on propagation of particles through the Universe must involve particle interactions: processes leading to production of particle cascades. While one expects existence of such cascades, the state of the art cosmic-ray research is oriented purely on a detection of single particles, gamma rays or associated extensive air showers. The natural extension of the cosmic-ray research with the studies on ensembles of particles and air showers is being proposed by the CREDO Collaboration. Within the CREDO strategy the focus is put on generalized super-preshowers (SPS): spatially and/or temporally extended cascades of particles originated above the Earth atmosphere, possibly even at astrophysical distances. With CREDO we want to find out whether SPS can be at least partially observed by a network of terrestrial and/or satellite detectors receiving primary or secondary cosmic-ray signal. This paper addresses electromagnetic SPS, e.g. initiated by VHE photons interacting with the cosmic microwave background, and the SPS signatures that can be seen by gamma-ray telescopes, exploring the exampleof Cherenkov Telescope Array. The energy spectrum of secondary electrons and photons in an electromagnetic super-preshower might be extended over awide range of energy, down to TeV or even lower, as it is evident from the simulation results. This means that electromagnetic showers induced by such particles in the Earth atmosphere could be observed by imaging atmospheric Cherenkov telescopes. We present preliminary results from the study of response of the Cherenkov Telescope Array to SPS events, including the analysis of the simulated shower images on the camera focal plane and implementedgeneric reconstruction chains based on the Hillas parameters.

The main objective of the Cosmic-Ray Extremely Distributed Observatory (CREDO) is the detection and analysis of extended cosmic ray phenomena, so-called super-preshowers (SPS), using existing as well as new infrastructure (cosmic-ray observatories, educational detectors, single detectors etc.). The search for ensembles of cosmic ray events initiated by SPS is yet an untouched ground, in contrast to the current state-of-the-art analysis, which is focused on the detection of single cosmic ray events. Theoretical explanation of SPS could be given either within classical (e.g., photon-photon interaction) or exotic (e.g., Super Heavy Dark Matter decay or annihilation) scenarios, thus detection of SPS would provide a better understanding of particle physics, high energy astrophysics and cosmology. The ensembles of cosmic rays can be classified based on the spatial and temporal extent of particles constituting the ensemble. Some classes of SPS are predicted to have huge spatial distribution, a unique signature detectable only with a facility of the global size. Since development and commissioning of a completely new facility with such requirements is economically unwarranted and time-consuming, the global analysis goals are achievable when all types of existing detectors are merged into a worldwide network. The idea to use the instruments in operation is based on a novel trigger algorithm: in parallel to looking for neighbour surface detectors receiving the signal simultaneously, one should also look for spatially isolated stations clustered in a small time window.

The Cosmic-Ray Extremely Distributed Observatory (CREDO) is an infrastructure for global analysis of extremely extended cosmic-ray phenomena, so-called super-preshowers, beyond the capabilities of existing, discrete, detectors and observatories. To date cosmic-ray research has been focused on detecting single air showers, while the search for ensembles of cosmic-ray events induced by super-preshowers is a scientific terra incognita - CREDO explores this uncharted realm. Positive detection of super-preshowers would have an impact on ultra-high energy astrophysics, cosmology and the physics of fundamental particle interactions as they can theoretically be formed within both classical (photon-photon interactions) and exotic (Super Heavy Dark Matter particle decay and interaction) scenarios. Some super-preshowers are predicted to have a significant spatial extent - a unique signature only detectable with the existing cosmic-ray infrastructure taken as a global network. An obvious, although yet unprobed, super-preshower 'detection limit' would be located somewhere between an air shower, induced by a super-preshower composed of tightly collimated particles, and a super-preshower composed of particles spread so widely that only few of them can reach the Earth. CREDO will probe this detection limit, leading to either an observation of an as yet unseen physical phenomenon, or the setting upper limits to the existence of large extraterrestrial cascades which would constrain fundamental physics models. While CREDO's focus is on testing physics at energies close to the Grand Unified Theories range, the broader phenomena are expected to be composed of particles with energies ranging from GeV to ZeV. This motivates our advertising of this concept across the astroparticle physics community.

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.

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 \pm 0.7 (stat) \pm 6.7 (sys) 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-principle 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.

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

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.