We present a method to analyse arrival directions of ultra-high-energy cosmic rays (UHECRs) using a classifier defined by a deep convolutional neural network trained on a HEALPix grid. To illustrate a high effectiveness of the method, we employ it to estimate prospects of detecting a large-scale anisotropy of UHECRs induced by a nearby source with an (orbital) detector having a uniform exposure of the celestial sphere and compare the results with our earlier calculations based on the angular power spectrum. A minimal model for extragalactic cosmic rays and neutrinos by Kachelrie{\ss}, Kalashev, Ostapchenko and Semikoz (2017) is assumed for definiteness and nearby active galactic nuclei Centaurus A, M82, NGC253, M87 and Fornax A are considered as possible sources of UHECRs. We demonstrate that the proposed method drastically improves sensitivity of an experiment by decreasing the minimal required amount of detected UHECRs or the minimal detectable fraction of from-source events several times compared to the approach based on the angular power spectrum. We also test robustness of the neural networks against different models of the large-scale Galactic magnetic fields and variations of the mass composition of UHECRs, and consider situations when there are two nearby sources or the dominating source is not known a priori. In all cases, the neural networks demonstrate good performance unless the test models strongly deviate from those used for training. The method can be readily applied to the analysis of data of the Telescope Array, the Pierre Auger Observatory and other cosmic ray experiments.

TUS (Tracking Ultraviolet Set-up) is the world's first orbital detector of ultra-high-energy cosmic rays (UHECRs). It was launched into orbit on 28th April 2016 as a part of the scientific payload of the Lomonosov satellite. The main aim of the mission was to test the technique of measuring the ultraviolet fluorescence and Cherenkov radiation of extensive air showers generated by primary cosmic rays with energies above ~100 EeV in the Earth atmosphere from space. During its operation for 1.5 years, TUS registered almost 80,000 events with a few of them satisfying conditions anticipated for extensive air showers (EASs) initiated by UHECRs. Here we discuss an event registered on 3rd October 2016. The event was measured in perfect observation conditions as an ultraviolet track in the nocturnal atmosphere of the Earth, with the kinematics and the light curve similar to those expected from an EAS. A reconstruction of parameters of a primary particle gave the zenith angle around 44$^\circ$ but an extreme energy not compatible with the cosmic ray energy spectrum obtained with ground-based experiments. We discuss in details all conditions of registering the event, explain the reconstruction procedure and its limitations and comment on possible sources of the signal, both of anthropogenic and astrophysical origin. We believe this detection represents a significant milestone in the space-based observation of UHECRs because it proves the capability of an orbital telescope to detect light signals with the apparent motion and light shape similar to what are expected from EASs. This is important for the on-going development of the future missions KLYPVE-EUSO and POEMMA, aimed for studying UHECRs from space.

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.

JEM-EUSO is an international program for the development of space based Ultra High Cosmic Ray observatories. The program consists of a series of missions which are either under development or in the data analysis phase. The instruments consist of wide field of view telescopes, which operate in the UV range, designed to detect the fluorescence light emitted by extended air showers in the atmosphere. We describe the simulation software ESAF and show the physical assumptions done in it. We present here the implementation of the POEMMA, K-EUSO, TUS, Mini-EUSO, EUSO-SPB and EUSO-TA detectors in the simulation software ESAF.

KLYPVE-EUSO (K-EUSO) is a planned orbital detector of ultra-high-energy cosmic rays (UHECRs), which is to be deployed on board the International Space Station. K-EUSO is expected to have a uniform exposure over the celestial sphere and register from 120 to 500 UHECRs at energies above 57 EeV in a 2-year mission. We employed the TransportCR and CRPropa 3 packages to estimate prospects of detecting a large-scale anisotropy of ultra-high-energy cosmic rays from a nearby source with K-EUSO. Nearby active galactic nuclei Centaurus A, M82, NGC 253, M87 and Fornax A were considered as possible sources of UHECRs. A minimal model for extragalactic cosmic rays and neutrinos by Kachelriess, Kalashev, Ostapchenko and Semikoz (2017) was chosen for definiteness. We demonstrate that an observation of $\gtrsim300$ events will allow detecting a large-scale anisotropy with a high confidence level providing the fraction of from-source events is $\simeq$10-15%, depending on a particular source. The threshold fraction decreases with an increasing sample size. We also discuss if an overdensity originating from a nearby source can be observed at around the ankle in case a similar anisotropy is found beyond 57 EeV. The results are generic and hold for other future experiments with a uniform exposure of the celestial sphere.

TUS (Tracking Ultraviolet Set-up), the first orbital telescope of ultra-high energy cosmic rays (UHECRs), has demonstrated that instruments of this kind have much broader capabilities and can also detect meteors, transient luminous events, anthropogenic glow and other processes taking place in Earth atmosphere in the UV frequency range. In this short letter, we address the question if an orbital UHECR detector can also register gamma-ray bursts (GRBs) via the fluorescent glow of irradiated nocturnal atmosphere. An analysis of the Fermi GBM catalog of GRBs and properties of several active and perspective instruments reveals that a detector with parameters of the optical system similar to those of the KLYPVE-EUSO (K-EUSO) or POEMMA telescopes and an appropriate \slow-mode" trigger will be able to observe only two GRBs in average every year of operation. POEMMA will register up to 3-4 GRBs in a year of operation in the stereo mode.

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 origin and nature of extreme energy cosmic rays (EECRs), which have energies above the 50 EeV, the Greisen-Zatsepin-Kuzmin (GZK) energy limit, is one of the most interesting and complicated problems in modern cosmic-ray physics. Existing ground-based detectors have helped to obtain remarkable results in studying cosmic rays before and after the GZK limit, but have also produced some contradictions in our understanding of cosmic ray mass composition. Moreover, each of these detectors covers only a part of the celestial sphere, which poses problems for studying the arrival directions of EECRs and identifying their sources. As a new generation of EECR space detectors, TUS (Tracking Ultraviolet Set-up), KLYPVE and JEM-EUSO, are intended to study the most energetic cosmic-ray particles, providing larger, uniform exposures of the entire celestial sphere. The TUS detector, launched on board the Lomonosov satellite on April 28, 2016, from Vostochny Cosmodrome in Russia, is the first of these. It employs a single-mirror optical system and a photomultiplier tube matrix as a photo-detector and will test the fluorescent method of measuring EECRs from space. Utilizing the Earth's atmosphere as a huge calorimeter, it is expected to detect EECRs with energies above 100 EeV. It will also be able to register slower atmospheric transient events: atmospheric fluorescence in electrical discharges of various types including precipitating electrons escaping the magnetosphere and from the radiation of meteors passing through the atmosphere. We describe the design of the TUS detector and present results of different ground-based tests and simulations.

TUS (Tracking Ultraviolet Set-up), the first orbital detector of extreme energy cosmic rays (EECRs), those with energies above 50 EeV, was launched into orbit on April 28, 2016, as a part of the Lomonosov satellite scientific payload. The main aim of the mission is to test a technique of registering fluorescent and Cherenkov radiation of extensive air showers generated by EECRs in the atmosphere with a space telescope. We present preliminary results of its operation in a mode dedicated to registering extensive air showers in the period from August 16, 2016, to November 4, 2016. No EECRs have been conclusively identified in the data yet, but the diversity of ultraviolet emission in the atmosphere was found to be unexpectedly rich. We discuss typical examples of data obtained with TUS and their possible origin. The data is important for obtaining more accurate estimates of the nocturnal ultraviolet glow of the atmosphere, necessary for successful development of more advanced orbital EECR detectors including those of the KLYPVE (K-EUSO) and JEM-EUSO missions.

TUS (Track Ultraviolet Setup) is the first space experiment aimed to check the possibility of registering extreme energy cosmic rays (EECRs) at E>50 EeV by measuring the fluorescence signal of extensive air showers in the atmosphere. The detector operates as a part of the scientific payload of the Lomonosov satellite for more than a year. We describe an algorithm of searching for EECR events in the TUS data and briefly discuss a number of candidates selected by formal criteria.

TUS is the world's first orbital detector of extreme energy cosmic rays (EECRs), which operates as a part of the scientific payload of the Lomonosov satellite since May 19, 2016. TUS employs the nocturnal atmosphere of the Earth to register ultraviolet (UV) fluorescence and Cherenkov radiation from extensive air showers generated by EECRs as well as UV radiation from lightning strikes and transient luminous events, micro-meteors and space debris. The first months of its operation in orbit have demonstrated an unexpectedly rich variety of UV radiation in the atmosphere. We briefly review the design of TUS and present a few examples of events recorded in a mode dedicated to registering EECRs.

On April 28, 2016, the Lomonosov satellite, equipped with a number of scientific instruments, was launched into orbit. Here we present briefly some of the results obtained with the first orbital telescope of extreme energy cosmic rays TUS and by a group of detectors aimed at multi-messenger observations of gamma-ray bursts.

In the present-day cosmic ray data, the strongest indication of anisotropy of the ultrahigh energy cosmic rays is the 20-degree hotspot observed by the Telescope Array with the statistical significance of 3.4\sigma. In this work, we study the possibility of detecting such a spot by space-based all-sky observatories. We show that if the detected luminosity of the hotspot is attributed to a physical effect and not a statistical fluctuation, the KLYPVE and JEM-EUSO experiments would need to collect ~300 events with E>57 EeV in order to detect the hotspot at the 5\sigma\ confidence level with the 68% probability. We also study the dependence of the detection prospects on the hotspot luminosity.

Two types of orbital detectors of extreme energy cosmic rays are being developed nowadays: (i) TUS and KLYPVE with reflecting optical systems (mirrors) and (ii) JEM-EUSO with high-transmittance Fresnel lenses. They will cover much larger areas than existing ground-based arrays and almost uniformly monitor the celestial sphere. The TUS detector is the pioneering mission developed in SINP MSU in cooperation with several Russian and foreign institutions. It has relatively small field of view (+/-4.5 deg), which corresponds to a ground area of 6.4x10^3 sq.km. The telescope consists of a Fresnel-type mirror-concentrator (~2 sq.m) and a photo receiver (a matrix of 16x16 photomultiplier tubes). It is to be deployed on the Lomonosov satellite, and is currently at the final stage of preflight tests. Recently, SINP MSU began the KLYPVE project to be installed on board of the Russian segment of the ISS. The optical system of this detector contains a larger primary mirror (10 sq.m), which allows decreasing the energy threshold. The total effective field of view will be at least +/-14 degrees to exceed the annual exposure of the existing ground-based experiments. Several configurations of the detector are being currently considered. Finally, JEM-EUSO is a wide field of view (+/-30 deg) detector. The optics is composed of two curved double-sided Fresnel lenses with 2.65 m external diameter, a precision diffractive middle lens and a pupil. The ultraviolet photons are focused onto the focal surface, which consists of nearly 5000 multi-anode photomultipliers. It is developed by a large international collaboration. All three orbital detectors have multi-purpose character due to continuous monitoring of various atmospheric phenomena. The present status of development of the TUS and KLYPVE missions is reported, and a brief comparison of the projects with JEM-EUSO is given.

Results of a blind search for localized regions of an excessive flux of cosmic rays in the energy range from 50 TeV to 20 PeV with the data of the FIAN KLARA-Chronotron experiment, the EAS MSU array and the Prototype of the EAS-1000 array are presented. A number of regions with a significant excess of the registered flux over an expected isotropic background are found. Some of the regions are present in at least two of the data sets considered.

We study anisotropy of cosmic rays in the energy range 0.2-1.4 EeV at intermediate angular scales using the public data set of the Pierre Auger Observatory. At certain scales, the analysis reveals a number of deviations from the isotropic distribution with the statistical significance above three standard deviations. It also demonstrates that the anisotropy evolves with energy. If confirmed with the complete Auger or Telescope Array data sets, the result can shed new light on the structure of galactic magnetic fields and the problem of transition from galactic to extragalactic cosmic rays.

A search for anisotropy at different scales in the arrival directions of cosmic rays with energies above 1 PeV is performed with data of the Tunka-133 array. A number of regions with angular sizes up to approximately 30 degrees that deviate from the uniform background with a statistical significance greater than three standard deviations are found using the shuffling technique both in the complete data set and its subsets in different energy ranges. Some regions with an excess of cosmic rays correlate with the distribution of possible cosmic ray sources. A comparison of the results with a similar analysis of the legacy data of the EAS MSU array is presented.

The origin of ultra-high-energy cosmic rays (UHECRs) is one of the most intriguing problems of modern cosmic ray physics. We briefly review the main astrophysical models of their origin and the forthcoming orbital experiments TUS and JEM-EUSO, and discuss how the new data can help one solve the long-standing puzzle.

Recent results of Milagro, Tibet, ARGO-YBJ and IceCube experiments on the small-scale anisotropy of Galactic cosmic rays (CRs) with energies from units up to a few hundred TeV arise a question on a possible nature of the observed phenomenon, as well as on the anisotropy of CRs at higher energies. An analysis of a small-scale anisotropy of CRs with energies at around PeV registered with the EAS MSU array presented in the article, reveals a number of regions with an excessive flux. A typical size of the regions varies from 3 up to 12 degrees. We study correlation of these regions with positions of potential astrophysical sources of CRs and discuss a possible origin of the observed anisotropy.

Three regions of excessive flux of cosmic rays with energies of the order of PeV are found in the experimental data of the EAS MSU array at a confidence level greater than 4\sigma. For two of them, there are similar regions in the experimental data of the EAS-1000 Prototype array. One of the interesting features of the regions is the absence of supernova remnants in their vicinities, traditionally considered as the main sources of Galactic cosmic rays, but the presence of isolated pulsars, some of which are able to accelerate heavy nuclei up to energies close to PeV. In our opinion, this favors the assumption that isolated pulsars are able to contribute to the flux of Galactic cosmic rays more than is usually assumed.

An analysis of arrival directions of extensive air showers (EAS) registered with the EAS MSU and EAS-1000 prototype arrays has revealed a region of excessive flux of PeV cosmic rays in the direction toward pulsars PSR J1840+5640 and LAT PSR J1836+5925 at significance level up to 4.5sigma. The first of the pulsars was discovered almost 30 years ago and is a well-studied old radio pulsar located at the distance of 1.7pc from the Solar system. The second pulsar belongs to a new type of pulsars, discovered by the space gamma-ray observatory Fermi, pulsations of which are not observed in optical and radio wavelengths but only in the gamma-ray range of energies (gamma-ray-only pulsars). In our opinion, the existence of the region of excessive flux of cosmic rays registered with two different arrays provides a strong evidence that isolated pulsars can give a noticeable contribution to the flux of Galactic cosmic rays in the PeV energy range.

Angular distributions of extensive air showers with different number of charged particles in the range 2.5x10^5--4x10^7 are derived using the experimental data obtained with the EAS MSU array. Possible approximations of the obtained distributions with different empiric functions available in literature, are analysed. It is shown that the exponential function provides the best approximation of the angular distributions in the sense of the chi-squared criterion.

An analysis of arrival directions of extensive air showers generated by cosmic rays in the PeV energy range and registered with the EAS MSU and EAS-1000 Prototype arrays reveals a considerable number of regions of excessive flux of cosmic rays. We present results of comparative analysis of regions found in the two data sets, estimate probabilities of their appearance, and discuss correlation of their locations with coordinates of possible astrophysical sources of PeV cosmic rays.

New results of an analysis of arrival directions of extensive air showers registered with the EAS--1000 Prototype array from August 1997 till February 1999 are presented. The method of Alexandreas et al., which has been used for analysis of data registered with CYGNUS, Milagrito, HEGRA AIROBICC, KASCADE and a number of other experiments, is employed. The existence of zones of excessive flux of cosmic rays with energies in the region of the knee is confirmed, as well as closeness of the zones to coordinates of possible astrophysical cosmic ray sources.

We study arrival directions of 1.4x10^6 extensive air showers (EAS) registered with the EAS--1000 Prototype Array in the energy range 0.1--10 PeV. By applying an iterative algorithm that provides uniform distribution of the data with respect to sidereal time and azimuthal angles, we find a number of zones with excessive flux of cosmic rays (CRs) at >=3 sigma level. We compare locations of the zones with positions of galactic supernova remnants (SNRs), pulsars, open star clusters, and regions of ionized hydrogen and find remarkable coincidences, which may witness in favour of the hypothesis that certain objects of these types, including the SNRs Cassiopeia A, the Crab Nebula, the Monogem Ring and some other, provide a noticeable contribution to the flux of CRs in the PeV range of energies. In addition, we find certain signs of a contribution from the M33 galaxy and a number of comparatively nearby groups of active galactic nuclei and interacting galaxies, in particular those in the Virgo cluster of galaxies. The results also provide some hints for a search of possible sources of ultra-high energy (UHE) cosmic rays and support an earlier idea that a part of both UHE and PeV CRs may originate from the same astrophysical objects.

Bursts of the count rate of extensive air showers (EAS) lead to the appearance of clusters in time series that represent EAS arrival times. We apply methods of nonlinear time series analysis to twenty EAS cluster events found in the data set obtained with the EAS-1000 prototype array. In particular, we use the Grassberger-Procaccia algorithm to compute the correlation dimension of the time series in the vicinity of the clusters. We find that four cluster events produce signs of chaos in the corresponding time series. By applying a number of supplementary methods we assess that the nature of the observed behaviour of the correlation dimension is likely to be deterministic. We suggest a simple qualitative model that might explain an origin of clusters in general and "possibly chaotic" clusters in particular. Finally, we compare our conclusions with the results of similar investigations performed by the EAS-TOP and LAAS groups.

We perform cluster analysis of arrival times of extensive air showers (EAS) registered with the EAS-1000 Prototype array in the period from August, 1997 till February, 1999. We present twenty cluster events each consisting of one or several EAS clusters, study the dynamics of their development, and analyze the angular distribution of EAS in clusters. We find that there may be certain correlation between EAS clusters with mean electron number $\sim10^5$ on the one hand and gamma-ray bursts and ultrahigh energy cosmic rays on the other.

We apply certain methods of nonlinear time series analysis to the extensive air shower clusters found earlier in the data set obtained with the EAS-1000 Prototype array. In particular, we use the Grassberger-Procaccia algorithm to compute the correlation dimension of samples in the vicinity of the clusters. The validity of the results is checked by surrogate data tests and some additional quantities. We compare our conclusions with the results of similar investigations performed by the EAS-TOP and LAAS groups.

We have analysed arrival times of extensive air showers (EAS) registered with the EAS-1000 prototype array during the period from August, 1997 till February, 1999. Our analysis has revealed that though the vast majority of samples of consecutive time intervals between EAS arrival times obey the exponential distribution, there are sequences of showers that have another distribution and thus violate the homogeneity hypothesis. The search for correlation between such events and clusters of showers and events with big delays between arrival times was also carried out.