Compilation of papers presented by the JEM-EUSO Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin.

Context. Recent observations of unexpected structures in the Galactic Cosmic Ray (GCR) spectrum and composition, as well as growing evidence for episodes of intense dynamical activity in the inner regions of the Galaxy, call for an evaluation of the high-energy particle acceleration associated with such activity and its potential impact on the global GCR phenomenology. Aims. We investigate whether particles accelerated during high-power episodes around the Galactic center can account for a significant fraction of the observed GCRs, or conversely what constraints can be derived regarding their Galactic transport if their contributions are negligible. Methods. We address these questions by studying the contribution of a continuous source of energetic particles at the Galactic center to the local GCRs. Particle transport in the Galaxy is described with a two-zone analytical model. We solve for the contribution of a Galactic Center Cosmic-Ray (GCCR) source using Green functions and Bessel expansion, and discuss the required injection power for these GCCRs to influence the global GCR phenomenology at Earth. Results. We find that, with standard parameters for particle propagation in the galactic disk and halo, the GCCRs can make a significant or even dominant contribution to the total CR flux observed at Earth. Depending on the parameters, such a source can account for both the observed proton flux and B-to-C ratio (in the case of a Kraichnan-like scaling of the diffusion coefficient), or potentially produce spectral and composition features. Conclusions. Our results show that the contribution of GCCRs cannot be neglected a priori, and that they can influence the global GCR phenomenology significantly, thereby calling for a reassessement of the standard inferences from a scenario where GCRs are entirely dominated by a single type of sources distributed throughout the Galactic disk.

EUSO-Balloon is a pathfinder mission for the Extreme Universe Space Observatory onboard the Japanese Experiment Module (JEM-EUSO). It was launched on the moonless night of the 25$^{th}$ of August 2014 from Timmins, Canada. The flight ended successfully after maintaining the target altitude of 38 km for five hours. One part of the mission was a 2.5 hour underflight using a helicopter equipped with three UV light sources (LED, xenon flasher and laser) to perform an inflight calibration and examine the detectors capability to measure tracks moving at the speed of light. We describe the helicopter laser system and details of the underflight as well as how the laser tracks were recorded and found in the data. These are the first recorded laser tracks measured from a fluorescence detector looking down on the atmosphere. Finally, we present a first reconstruction of the direction of the laser tracks relative to the detector.

The Telescope Array (TA) shows a 20$^{\circ}$ hotspot as well as an excess of UHECRs above 50~EeV when compared with the Auger spectrum. We consider the possibility that both the TA excess and hotspot are due to a dominant source in the Northern sky. We carry out detailed simulations of UHECR propagation in both the intergalactic medium and the Galaxy, using different values for the intergalactic magnetic field. We consider two general classes of sources: transients and steady, adopting a mixed UHECR composition that is consistent with the one found by Auger. The spatial location of the sources is draw randomly. We generate Auger-like and TA-like data sets from which we determine the spectrum, the sky maps and the level of anisotropy. We find that, while steady sources are favored over transients, it is unlikely to account for all the currently available observational data. While we reproduce fairly well the Auger spectrum for the vast majority of the simulated data sets, most of the simulated data sets with a spectrum compatible with that of TA (at most a few percent depending on density model tested) show a much stronger anisotropy than the one observed. We find that the rare cases in which both the spectrum and the anisotropy are consistent require a steady source within $\sim 10$ Mpc, to account for the flux excess, and a strong extragalactic magnetic field $\sim 10$ nG, to reduce the excessive anisotropy.

GeV-TeV gamma-ray and PeV-EeV neutrino backgrounds provide a unique window on the nature of the ultra-high-energy cosmic-rays (UHECRs). We discuss the implications of the recent Fermi-LAT data regarding the extragalactic gamma-ray background (EGB) and related estimates of the contribution of point sources as well as IceCube neutrino data on the origin of the UHECRs. We calculate the diffuse flux of cosmogenic $\gamma$-rays and neutrinos produced during the UHECRs propagation and derive constraints on the possible cosmological evolution of UHECR sources. In particular, we show that the mixed-composition scenario which is in agreement with both (i) Auger measurements of the energy spectrum and composition up to the highest energies and (ii) the ankle-like feature in the light component detected by KASCADE-Grande, is compatible with both the Fermi-LAT measurements and with current IceCube limits.

GeV-TeV gamma-rays and PeV-EeV neutrino backgrounds provide a unique window on the nature of the ultra-high-energy cosmic-rays (UHECRs). We discuss the implications of the recent Fermi-LAT data regarding the extragalactic gamma-ray background (EGB) and related estimates of the contribution of point sources as well as IceCube neutrino data on the origin of the UHECRs. We calculate the diffuse flux of cosmogenic $\gamma$-rays and neutrinos produced by the UHECRs and derive constraints on the possible cosmological evolution of UHECR sources. In particular, we show that the mixed-composition scenario considered in \citet{Globus2015b}, which is in agreement with both (i) Auger measurements of the energy spectrum and composition up to the highest energies and (ii) the ankle-like feature in the light component detected by KASCADE-Grande, is compatible with both the Fermi-LAT measurements and with current IceCube limits. We also discuss the possibility for future experiments to detect associated cosmogenic neutrinos and further constrain the UHECR models, including possible subdominant UHECR proton sources.

The Telescope Array (TA) shows a 20$^{\circ}$ hotspot as well as an excess of UHECRs above 50~EeV when compared with the Auger spectrum. We consider the possibility that both the TA excess and hotspot are due to a dominant source in the Northern sky. We carry out detailed simulations of UHECR propagation in both the intergalactic medium and the Galaxy, using different values for the intergalactic magnetic field. We consider two general classes of sources: transients and steady, adopting a mixed UHECR composition that is consistent with the one found by Auger. The spatial location of the sources is draw randomly. We generate Auger-like and TA-like data sets from which we determine the spectrum, the sky maps and the level of anisotropy. We find that, while steady sources are favored over transients, it is unlikely to account for all the currently available observational data. While we reproduce fairly well the Auger spectrum for the vast majority of the simulated data sets, most of the simulated data sets with a spectrum compatible with that of TA (at most a few percent depending on density model tested) show a much stronger anisotropy than the one observed. We find that the rare cases in which both the spectrum and the anisotropy are consistent require a steady source within $\sim 10$ Mpc, to account for the flux excess, and a strong extragalactic magnetic field $\sim 10$ nG, to reduce the excessive anisotropy.

We present a complete phenomenological model accounting for the evolution of the cosmic-ray spectrum and composition with energy, based on the available data over the entire spectrum. We show that there is no need to postulate any additional component, other than one single Galactic component depending on rigidity alone, and one extragalactic component, whose characteristics are similar to those derived from a study of particle acceleration at mildly relativistic shocks in a GRB environment (Globus et al., 2015). In particular, we show that the resulting cosmic ray spectrum and composition satisfy the various constraints derived from the current data in the Galactic/extragalactic transition region, notably from the measurements of KASCADE Grande and Auger. Finally, we derive some generic features that a working phenomenological scenario may exhibit to give a global account of the cosmic ray data with a minimum number of free parameters.

We study the acceleration of CR protons and nuclei at GRB internal shocks. Physical quantities and their time evolution are estimated using the internal shock modeling implemented by Daigne & Mochkovitch 1998. We consider different hypotheses about the way the energy dissipated at internal shocks is shared between accelerated CR, e- and B field. We model CR acceleration at mildly relativistic shocks, including all the significant energy loss processes. We calculate CR and neutrino release from single GRBs, assuming that nuclei heavier than protons are present in the relativistic wind. Protons can only reach maximum energies of ~ 10^19.5 eV, while intermediate and heavy nuclei are able to reach values of ~ 10^20 eV and above. The spectra of nuclei escaping from the acceleration site are found to be very hard while the combined spectrum of protons and neutrons is much softer. We calculate the diffuse UHECR flux expected on Earth using the GRB luminosity function from Wanderman & Piran 2010. Only the models assuming that the prompt emission represent a very small fraction of the energy dissipated at internal shocks, and that most of this dissipated energy is communicated to accelerated CR, are able to reproduce the magnitude of the UHECR flux observed. For these models, the observed shape of the UHECR spectrum can be well reproduced and the evolution of the composition is compatible with the trend suggested by Auger. We discuss implications of the softer proton component for the GCR to EGCR transition in the light of the recent composition analyses (KASCADE-Grande experiment). The associated secondary particle diffuse fluxes do not upset any current observational limit. Diffuse neutrino flux from GRB sources should however be detected with the lifetime of neutrino observatories.

UHECRs have attracted a lot of attention due to their challengingly high energies and their potential value to constrain physical processes and astrophysical parameters in the most energetic sources of the universe. Current detectors have failed to detect significant anisotropies which had been expected to allow source identification. Some indications about the UHECR composition, which may become heavier at the highest energies, has even put into question the possibility that such a goal could be achieved soon. We investigate the potential value of a new-generation detector, with 10 times larger exposure, to overcome the current situation and make significant progress in the detection of anisotropies and thus in the study of UHECRs. We take as an example the expected performances of the JEM-EUSO, assuming a uniform full-sky coverage with a total exposure of 300,000 km2 sr yr. We simulate realistic UHECR sky maps for a wide range of possible astrophysical scenarios allowed by the current constraints, taking into account the energy losses and photo-dissociation of the UHECRs, as well as their deflections by magnetic fields. These sky maps, built for the expected statistics of JEM-EUSO as well as for the current Auger statistics, as a reference, are analyzed from the point of view of their intrinsic anisotropies, using the two-point correlation function. A statistical study of the resulting anisotropies is performed for each astrophysical scenario, varying the UHECR source composition and spectrum as well as the source density. We find that significant anisotropies are expected to be detected by a next-generation UHECR detector, for essentially all the astrophysical scenarios studied, and give precise, quantitative meaning to this statement. Our results show that a gain of one order of magnitude in exposure would make a significant difference compared to the existing detectors.

We evaluate the exposure during nadir observations with JEM-EUSO, the Extreme Universe Space Observatory, on-board the Japanese Experiment Module of the International Space Station. Designed as a mission to explore the extreme energy Universe from space, JEM-EUSO will monitor the Earth's nighttime atmosphere to record the ultraviolet light from tracks generated by extensive air showers initiated by ultra-high energy cosmic rays. In the present work, we discuss the particularities of space-based observation and we compute the annual exposure in nadir observation. The results are based on studies of the expected trigger aperture and observational duty cycle, as well as, on the investigations of the effects of clouds and different types of background light. We show that the annual exposure is about one order of magnitude higher than those of the presently operating ground-based observatories.

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

Contributions of the JEM-EUSO Collaboration to the 32nd International Cosmic Ray Conference, Beijing, August, 2011.

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.

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malarg\"ue and averaged monthly models, the utility of the GDAS data is shown.

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.

We study the production of cosmogenic neutrinos and photons during the extragalactic propagation of ultra-high-energy cosmic rays (UHECRs). For a wide range of models in cosmological evolution of source luminosity, composition and maximum energy we calculate the expected flux of cosmogenic secondaries by normalizing our cosmic ray output to experimental spectra and comparing the diffuse flux of GeV-TeV gamma-rays to the experimental one measured by the Fermi satellite. Most of these models yield significant neutrino fluxes for current experiments like IceCube or Pierre Auger. Furthermore, we discuss the possibilities of signing the presence of UHE proton sources either within or outside the cosmic ray horizon using neutrinos or photons observations even if the cosmic ray composition becomes heavier at the highest energies. We discuss the possible constraints that could be brought on the UHECR origin from the different messengers and energy ranges.

In this paper we review the extragalactic propagation of ultrahigh energy cosmic-rays (UHECR). We present the different energy loss processes of protons and nuclei, and their expected influence on energy evolution of the UHECR spectrum and composition. We discuss the possible implications of the recent composition analyses provided by the Pierre Auger Observatory. The influence of extragalactic magnetic fields and possible departures from the rectilinear case are also mentioned as well as the production of secondary cosmogenic neutrinos and photons and the constraints their observation would imply for the UHECRs origin. Finally, we conclude by briefly discussing the relevance of a multi messenger approach for solving the mystery of UHECRs.

We present the results of an analysis of data recorded at the Pierre Auger Observatory in which we search for groups of directionally-aligned events (or `multiplets') which exhibit a correlation between arrival direction and the inverse of the energy. These signatures are expected from sets of events coming from the same source after having been deflected by intervening coherent magnetic fields. The observation of several events from the same source would open the possibility to accurately reconstruct the position of the source and also measure the integral of the component of the magnetic field orthogonal to the trajectory of the cosmic rays. We describe the largest multiplets found and compute the probability that they appeared by chance from an isotropic distribution. We find no statistically significant evidence for the presence of multiplets arising from magnetic deflections in the present data.

Studies of the cosmic ray energy spectrum at the highest energies with the Pierre Auger Observatory

Technical reports on operations and monitoring of the Pierre Auger Observatory

Studies of the composition of the highest energy cosmic rays with the Pierre Auger Observatory, including examination of hadronic physics effects on the structure of extensive air showers.

Ongoing and planned enhancements of the Pierre Auger Observatory