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.

Non-recurrent short term variations of the galactic cosmic-ray (GCR) flux above 70 MeV n$^{-1}$ were observed between 2016 February 18 and 2017 July 3 aboard the European Space Agency LISA Pathfinder (LPF) mission orbiting around the Lagrange point L1 at 1.5$\times$10$^6$ km from Earth. The energy dependence of three Forbush decreases (FDs) is studied and reported here. A comparison of these observations with others carried out in space down to the energy of a few tens of MeV n$^{-1}$ shows that the same GCR flux parameterization applies to events of different intensity during the main phase. FD observations in L1 with LPF and geomagnetic storm occurrence is also presented. Finally, the characteristics of GCR flux non-recurrent variations (peaks and depressions) of duration $<$ 2 days and their association with interplanetary structures are investigated. It is found that, most likely, plasma compression regions between subsequent corotating high-speed streams cause peaks, while heliospheric current sheet crossing cause the majority of the depressions.

Appreciable hard X-ray (HXR) and gamma-ray emissions in the 0.04 to 150 MeV energy range associated with the 2003 October 29 solar flare (X10/3B) were observed at 20:38 to 20:58 UT by the SONG instrument onboard the CORONAS-F mission. To restore flare gamma-ray spectra we fitted the SONG energy loss spectra with a three-component model of the incident spectrum: (1) a power law in energy, assumed to be due to electron bremsstrahlung; (2) a broad continuum produced by prompt nuclear de-excitation gamma-lines; and (3) a broad gamma-line generated from pion decay. We also restored spectra from the RHESSI data, compared them with the SONG spectra and found a reasonable agreement between these spectra in the 0.1 to 10 MeV energy range. The pion-decay emission was observed from 20:44:20 UT and had its maximum at 20:48 to 20:51 UT. The power-law spectral index of accelerated protons estimated from the ratio between intensities of different components of gamma rays changed with time. The hardest spectrum with a power-law index S = - 3.5 to - 3.6 was observed at 20:48 to 20:51 UT. Time histories of the pion-decay emission and proton spectrum were compared with changes of the locations of flare energy release as shown by RHESSI hard X-ray images and remote H alpha brightening. An apparent temporal correlation between processes of particle acceleration and restructuring of flare magnetic field was found. In particular, the protons were accelerated to sub-relativistic energies after radical change of the character of footpoint motion from a converging motion to a separation motion.

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.

We present a short overview of the event associated with the recent strong solar flare on 10 September 2017 (X 8.2) based on the available data both from satellite GOES-13 and from selected neutron monitors. The onset time of SPE/GLE at 1 AU was found between 16:06 and 16:08 UT. The GLE effect was anisotropic with a maximum increase of 6 %. The maximum energy of accelerated protons was about 6 GeV. We estimated the release time of sub-relativistic protons into open field lines as 15:53 to 15:55 UT.

It is shown that the hypothesis of the axion mechanism of Sun luminosity suggesting that the solar axion particles are born in the core of the Sun and may be efficiently converted back into $\gamma$-quanta in the magnetic field of the solar overshoot tachocline is physically relevant. As a result, it is also shown that the intensity variations of the $\gamma$-quanta of axion origin, induced by the magnetic field variations in the tachocline, directly cause the Sun luminosity %and total solar irradiance (TSI) variations and eventually characterize the active and quiet states of the Sun. Within the framework of this mechanism estimations of the strength of the axion coupling to a photon ($g_{a \gamma} = 3.6 \cdot 10^{-11} GeV^{-1}$) and the hadronic axion particle mass ($m_a \sim 2.3 \cdot 10^{-2} eV$) have been obtained. It is also shown that the claimed axion parameters do not contradict any known experimental and theoretical model-independent limitations.

Energetic (suprathermal) solar particles, accelerated in the interplanetary medium, contribute to the solar wind pressure, in particular during high solar activity periods. We estimated the effect of the increase of solar wind pressure due to suprathermal particles on magnetospheric transmissivity of galactic cosmic rays in the case of one recent solar event.

The Cosmic Rays propagation was studied in details using the HelMod-2D Monte Carlo code, that includes a general description of the diffusion tensor, and polar magnetic-field. The Numerical Approach used in this work is based on a set of Stochastic Differential Equations fully equivalent to the well know Parker Equation for the transport of Cosmic Rays. In our approach the Diffusion tensor in the frame of the magnetic field turbolence does not depends explicitly by Solar Latitude but varies with time using a diffusion parameter obtained by Neutron Monitors. The parameters of the Model were tuned using data during the solar Cycle 23 and Ulysses latitudinal Fast Scan in 1995. The actual parametrization is able to well reproduce the observed latitudinal gradient of protons and the southward shift of the minimum of latitudinal intensity. The description of the model is also available online at website www.helmod.org. The model was then applied on Pamela/Ulysses proton intensity from 2006 up to 2009. The model during this 4-year continous period agree well with both PAMELA (at 1 AU) and Ulysses data (at various solar distance and solar latitude). The agreement improves when considering the ratio between this data. Studies done also with particles with different charge (e.g. electrons) allow us to explain the presence (or not) of protons and electrons latitudinal gradients observed by Ulysses during the Latitudinal Fast Scan in 1995 and 2007.

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)

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.

Low energy galactic cosmic rays as well as particles accelerated to high energies either at the solar surface or in the interplanetary medium have access to the atmosphere above a given position on the Earth depending upon the state of the magnetosphere. The interpretation of the cosmic ray anisotropy, deduced from the neutron monitor (NM) network, must assume the variability of the magnetospheric configuration. Along with a short review of changes of the geomagnetic cutoffs in the disturbed magnetosphere reported in the earlier papers, we present the results of computations of transmissivity function and asymptotic directions for selected points on the ground and for a low altitude polar orbiting satellite as well. The computations, based on different available models of geomagnetic field of external sources are performed for quiet time periods and for strong geomagnetic disturbances occurred in 2003 and 2004.

The planned JEM-EUSO (Extreme Universe Space Observatory onboard the ISS Japanese Experimental Module) will measure the energy spectra of cosmic rays up to the range of 1000 EeV and will search for direction to their sources. It will observe the extensive air showers generated in the atmosphere by high energy cosmic ray primary particle from the space. The instantaneous aperture of the telescope will exceed by one order the aperture of the largest ground based detectors. JEM-EUSO apparatus is a large telescope with a diameter of 2.5 m with fast UV camera. Slovakia is responsible for the determination of the UV background, which influences the operational efficiency of the experiment and for the analysis of fake trigger events.

The cosmic rays differential intensity inside the heliosphere, for energy below 30 GeV/nuc, depends on solar activity and interplanetary magnetic field polarity. This variation, termed solar modulation, is described using a 2-D (radius and colatitude) Monte Carlo approach for solving the Parker transport equation that includes diffusion, convection, magnetic drift and adiabatic energy loss. Since the whole transport is strongly related to the interplanetary magnetic field (IMF) structure, a better understanding of his description is needed in order to reproduce the cosmic rays intensity at the Earth, as well as outside the ecliptic plane. In this work an interplanetary magnetic field model including the standard description on ecliptic region and a polar correction is presented. This treatment of the IMF, implemented in the HelMod Monte Carlo code (version 2.0), was used to determine the effects on the differential intensity of Proton at 1\,AU and allowed one to investigate how latitudinal gradients of proton intensities, observed in the inner heliosphere with the Ulysses spacecraft during 1995, can be affected by the modification of the IMF in the polar regions.

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

The differential intensities of Cosmic Rays at Earth were calculated using a 2D stochastic Montecarlo diffusion code and compared with observation data. We evaluated the effect of stretched and compressed heliospheres on the Cosmic Ray intensities at the Earth. This was studied introducing a dependence of the diffusion parameter on the heliospherical size. Then, we found that the optimum value of the heliospherical radius better accounting for experimental data. We also found that the obtained values depends on solar activity. Our results are compatible with Voyager observations and with models of heliospherical size modulation.

A propagation model of galactic cosmic protons through the Heliosphere was implemented using a 2-D Monte Carlo approach to determine the differential intensities of protons during the solar cycle 23. The model includes the effects due to the variation of solar activity during the propagation of cosmic rays from the boundary of the heliopause down to Earth's position. Drift effects are also accounted for. The simulated spectra were found in agreement with those obtained with experimental observations carried out by BESS, AMS and PAMELA collaborations. In addition, the modulated spectrum determined with the present code for the year 1995 exhibits the latitudinal gradient and equatorial southward offset minimum found by Ulysses fast scan in 1995.

Spectra of Galactic Cosmic Rays (GCRs) measured at the Earth are the combination of several processes: sources production and acceleration, propagation in the interstellar medium and propagation in the heliosphere. Inside the solar cavity the flux of GCRs is reduced due to the solar modulation, the interaction which they have with the interplanetary medium. We realized a 2D stochastic simulation of solar modulation to reproduce CR spectra at the Earth, and evaluated the importance in our results of the Local Interstellar Spectrum (LIS) model and its agreement with data at high energy. We show a good agreement between our model and the data taken by AMS-01 and BESS experiments during periods with different solar activity conditions. Furthermore we made a prediction for the flux which will be measured by AMS-02 experiment.

We implemented a quasi time-dependent 2D stochastic model of solar modulation describing the transport of cosmic rays (CR) in the heliosphere. Our code can modulate the Local Interstellar Spectrum (LIS) of a generic charged particle (light cosmic ions and electrons), calculating the spectrum at 1AU. Several measurements of CR antiparticles have been performed. Here we focused our attention on the CR antiproton component and the antiproton over proton ratio. We show that our model, using the same heliospheric parameters for both particles, fit the observed anti-p/p ratio. We show a good agreement with BESS-97 and PAMELA data and make a prediction for the AMS-02 experiment.

Galactic Cosmic Rays (GCRs) are mainly protons confined in the galactic magnetic field to form an isotropic flux inside the galaxy. Before reaching the Earth orbit they enter the Heliosphere and undergo diffusion, convection, magnetic drift and adiabatic energy loss. The result is a reduction of particles flux at low energy (below 10 GeV), called solar modulation. We realized a quasi time-dependent 2D Stochastic Simulation of Solar Modulation that is able to reproduce CR spectra once known the Local Interstellar Spectrum (LIS). We were able to estimate the different behaviors associated to the polarity dependence of the Heliospheric modulation for particles as well as for antiparticles. We show a good agreement with the antiproton/proton ratio measured by AMS-01, Pamela, BESS, Heat and Caprice and we performed a prediction for the AMS-02 Experiment.

The solar modulation, a combination of diffusion, convection, magnetic drift and energy loss inside the heliosphere is usually seen as a depletion in the Galactic cosmic ray (CR) flux at low energy (less than 10 GeV/nuc). Antiparticles such as antiprotons or positrons undergo the same processes of respective particles but with a different magnitude depending on the Solar magnetic field polarity. For electrons and positrons, due to the small mass, energy loss mechanisms as inverse compton, synchrotron, bremsstrahlung and ionization have to be taken into account, together with the typical adiabatic losses considered in the heliosphere. We developed a Monte Carlo stochastic simulation with the aim to compare the solar modulation of particles and antiparticles in the same observation period. We are able to estimate the different behaviours associated to the charge sign dependent processes of the heliospheric modulation. We compared the simulated positron fraction with measurements performed by AMS-01 and PAMELA. We also present the prediction for the AMS-02 experiment.

We show existence of strong negative correlation between the temporal variations of magnetic field toroidal component of the solar tachocline (the bottom of convective zone) and the Earth magnetic field (Y-component). The possibility that hypothetical solar axions, which can transform into photons in external electric or magnetic fields (the inverse Primakoff effect), can be the instrument by which the magnetic field of convective zone of the Sun modulates the magnetic field of the Earth is considered. We propose the axion mechanism of Sun luminosity and "solar dynamo - geodynamo" connection, where an energy of solar axions emitted in M1 transition in 57Fe nuclei is modulated at first by the magnetic field of the solar tachocline zone (due to the inverse coherent Primakoff effect) and after that is resonance absorbed in the core of the Earth, thereby playing the role of an energy source and a modulator of the Earth magnetic field. Within the framework of this mechanism estimations of the strength of an axion coupling to a photon (ga_gamma about 1.64*10^-9 GeV^-1), the axion-nucleon coupling (ga_gamma_eff about 10^-5) and the axion mass (ma about 30 eV) have been obtained.

The 28 October 2003 flare gave us the unique opportunity to compare the acceleration time of high-energy protons with the escaping time of those particles which have been measured onboard spacecraft and by neutron monitors network as GLE event. High-energy emission time scale and shock wave height and velocity time dependencies were also studied.

In this report we present two complex eruptive solar events and the associated Cosmic Ray effects (Forbush decrease). We use combined recordings from a number of Earthbound Receivers, Space Experiments and data archives (such as the ARTEMIS-IV Radio spectrograph, the Athens NEUTRON MONITOR, the LASCO CME Lists, the SONG of the {CORONAS-F} satellite, etc.). The influence of solar transients on the interplanetary medium conditions and the cosmic ray flux is analysed and discussed. The observed time sequence of events of this time period indicates that the initiation of CMEs is closely related to the appearance of type II and IV radio bursts and strong solar flares. Their effects extend from the lower corona to the near Earth vicinity affecting Cosmic Ray measurements and space weather. As regards the Forbush decrease our data indicate significant amplification at the presence of a MHD shock.

We show existence of strong negative correlation between the temporal variations of magnetic field toroidal component of the solar tachocline (the bottom of convective zone) and the Earth magnetic field (Y-component). The possibility that hypothetical solar axions, which can transform into photons in external electric or magnetic fields (the inverse Primakoff effect), can be the instrument by which the magnetic field of convective zone of the Sun modulates the magnetic field of the Earth is considered. We propose the axion mechanism of "solar dynamo-geodynamo" connection, where an energy of axions, which form in the Sun core, is modulated at first by the magnetic field of the solar tachocline zone (due to the inverse coherent Primakoff effect) and after that is absorbed in the liquid core of the Earth under influence of the terrestrial magnetic field, thereby playing the role of an energy source and a modulator of the Earth magnetic field. Within the framework of this mechanism new estimations of the strength of an axion coupling to a photon (ga_gamma about 5*10^-9 GeV^-1) and the axion mass (ma ~ 30 eV) have been obtained.

Various aspects of the connection between cloud cover (CC) and cosmic rays (CR) are analysed. Many features of this connection indicate that there is no direct causal connection between low cloud cover (LCC) and CR in spite of the evident long-term correlation between them. However, most of these features are indirect. If only some part of the LCC is connected and varies with CR, then its value, obtained from the joint analysis of their 11-year variations, and averaged over the globe, should be most likely less than 20%. The most significant argument against a causal connection of CR and LCC is the anticorrelation between LCC and the medium cloud cover (MCC). The scenario of the parallel influence of the solar activity on the global temperature and CC on one side and CR on the other, which can lead to the observed correlations, is discussed and advocated.

Various aspects of the connection between cloud cover (CC) and cosmic rays (CR) are analysed. We argue that the anticorrelation between the temporal behaviour of low (LCC) and middle (MCC) clouds evidences against causal connection between them and CR. Nevertheless, if a part of low clouds (LCC) is connected and varies with CR, then its most likely value averaged over the Globe should not exceed 20% at the two standard deviation level.

Recent research has shown that distinct physical regions in the Venusian induced magnetosphere are recognizable from the variations of strength and of wave/fluctuation activity of the magnetic field. In this paper the statistical properties of magnetic fluctuations are investigated in the Venusian magnetosheath, terminator, and wake regions. The latter two regions were not visited by previous missions. We found 1/f fluctuations in the magnetosheath, large-scale structures near the terminator and more developed turbulence further downstream in the wake. Location independent short-tailed non-Gaussian statistics was observed.

During the period when the Sun was intensely active on October-November 2003, two remarkable solar neutron events were observed by the ground-based neutron monitors. On October 28, 2003, in association with an X17.2 large flare, solar neutrons were detected with high statistical significance (6.4 sigma) by the neutron monitor at Tsumeb, Namibia. On November 4, 2003, in association with an X28 class flare, relativistic solar neutrons were observed by the neutron monitors at Haleakala in Hawaii and Mexico City, and by the solar neutron telescope at Mauna Kea in Hawaii simultaneously. Clear excesses were observed at the same time by these detectors, with the significance calculated as 7.5 sigma for Haleakala, and 5.2 sigma for Mexico City. The detector onboard the INTEGRAL satellite observed a high flux of hard X-rays and gamma-rays at the same time in these events. By using the time profiles of the gamma-ray lines, we can explain the time profile of the neutron monitor. It appears that neutrons were produced at the same time as the gamma-ray emission.