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.

EUSO-TA is a ground-based florescence detector built to validate the design of an ultra-high energy cosmic ray fluorescence detector to be operated in space. EUSO-TA detected the first air shower events with the technology developed within the JEM-EUSO program. It operates at the Telescope Array (TA) site in Utah, USA. With the external trigger provided by the Black Rock Mesa fluorescence detectors of Telescope Array (TA-FDs), EUSO-TA observed nine ultra-high energy cosmic ray events and several laser events from the Central Laser Facility of Telescope Array and portable lasers like the JEM-EUSO Global Light System prototype. The reconstruction parameters of the cosmic ray events which crossed the EUSO-TA field of view (both detected and not detected by EUSO-TA), were provided by the Telescope Array Collaboration. As the TA-FDs have a wider field of view than EUSO-TA ($\sim$30 times larger), they allow the cosmic ray energy reconstruction based on the observation of most of the extensive air-shower profiles, including the shower maximum, while EUSO-TA only observes a portion of the showers, usually far from the maximum. For this reason, the energy of the cosmic rays corresponding to the EUSO-TA signals appear lower than the actual ones. In this contribution, the analysis of the cosmic-ray events detected with EUSO-TA is discussed.

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 development of an analytically iterative method for solving steady-state as well as unsteady-state problems of cosmic-ray (CR) modulation is proposed. Iterations for obtaining the solutions are constructed for the spherically symmetric form of the CR propagation equation. The main solution of the considered problem consists of the zero-order solution that is obtained during the initial iteration and amendments that may be obtained by subsequent iterations. The finding of the zero-order solution is based on the CR isotropy during propagation in the space, whereas the anisotropy is taken into account when finding the next amendments. To begin with, the method is applied to solve the problem of CR modulation where the diffusion coefficient $\kappa$ and the solar wind speed $u$ are constants with an Local Interstellar Spectra (LIS) spectrum. The solution obtained with two iterations was compared with an analytical solution and with numerical solutions. Finally, solutions that have only one iteration for two problems of CR modulation with $u = constant$ and the same form of LIS spectrum were obtained and tested against numerical solutions. For the first problem, $\kappa$ is proportional to the momentum of the particle $p$, so it has the form $\kappa=k_0\eta$, where $\eta=\frac{p}{m_0c}$. For the second problem, the diffusion coefficient is given in the form $\kappa=k_0\beta\eta$, where $\beta=\frac{v}{c}$ is the particle speed relative to the speed of light. There was a good matching of the obtained solutions with the numerical solutions as well as with the analytical solution for the problem where $\kappa = constant$.

The cosmic rays modulation inside the heliosphere is well described by a transport equation introduced by Parker in 1965. To solve this equation several approaches were followed in the past. Recently the Monte Carlo approach becomes widely used in force of his advantages with respect to other numerical methods. In the Monte Carlo approach, the transport equation is associated to a fully equivalent set of Stochastic Differential Equations. This set is used to describe the stochastic path of a quasi-particle from a source, e.g., the interstellar medium, to a specific target, e.g., a detector at Earth. In this work, we present both the Forward-in-Time and Backward-in-Time Monte Carlo solutions. We present an implementation of both algorithms in the framework of HelMod Code showing that the difference between the two approach is below 5\% that can be quoted as the systematic uncertain of the Method itself.

We implemented a website to deal with main effects on Cosmic Ray access to the Earth, i.e. the Solar Modulation and the Geomagnetic Field effect. In helmod.org the end user can easily access a web interface to results catalog of the HelMod Monte Carlo Code. This Model uses a Monte Carlo Approach to solves the Parker Transport Equation, obtaining a modulated proton flux for a period (monthly average) between January 1990 and december 2007. geomagsphere.org is instead based on GeoMag Backtracing Code, that solves the Lorentz equation with a Runge-Kutta method of 6th order, and, reversing charge sign and velocity, reconstruct particle trajectories in the Earth Magnetosphere back in time. We use last models of internal (IGRF-11) and external (Tsyganenko 1996 -T96- and 2005 -T05-) field components valid up to 2015. Particles are backtraced to the outer (magnetopause) or inner boundary to separate Primary (allowed trajectory) from Secondary (forbidden) Cosmic Rays. This code has been used both for reproducing known effects as East-West effect and rigidity cutoff calculations. In geomagsphere.org the user can choose the external field model from Tsyganenko (T96 or T05) and obtain for a fixed position and date from 1st Jan. 1968 (T96) and 1st Jan. 1995 (T05) respectively till 31$^{st}$ Dec 2012, the vertical rigidity cutoff estimation obtained with the backtracing technique with a rigidity step of 0.1 GV. For a more precise calculation (0.01 GV), requiring more CPU time, results are sent to the user by email (mail model)