Our direct knowledge of solar eruptive events is limited to several decades and does not include extreme events, which can only be studied by the indirect proxy method over millennia, or by a large number of sun-like stars. There is a gap, spanning 1--2 orders of magnitude, in the strength of events between directly observed and reconstructed ones. Here, we study the proxy-method sensitivity to identify extreme solar particle events (SPEs). First, the strongest directly observed SPE (23-Feb-1956), used as a reference for proxy-based reconstructions, was revisited using the newly developed method. Next, the sensitivity of the cosmogenic-isotope method to detect a reference SPE was assessed against the precision and number of individual isotopic records, showing that it is too weak by a factor $\approx$30 to be reliably identified in a single record. Uncertainties of 10Be and 14C data are shown to be dominated by local/regional patterns and measurement errors, respectively. By combining several proxy records, an SPE 4--5 times stronger than the reference one can be potentially detected, increasing the present-day sensitivity by an order of magnitude. This will allow filling the observational gap in SPE strength distribution, thus enriching statistics of extreme events from 3--4 presently known ones to several tens. This will provide a solid basis for research in the field of extreme events, both for fundamental science, viz. solar and stellar physics, and practical applications, such as the risk assessments of severe space-based hazards for the modern technological society.

Using neutron monitor and space-borne data we performed an analysis of the second ground level enhancement of solar cycle 24, namely the event of 10 September 2017 (GLE 72) and derive the spectral and angular characteristics of GLE particles. We employ new neutron monitor yield function and a recently proposed model based on optimization procedure. The method consists of simulation of particle propagation in a model magnetosphere in order to derive the cut-off rigidity and neutron monitor asymptotic directions. Subsequently the rigidity spectrum and anisotropy of GLE particles are obtained in their dynamical evolution during the event on the basis of inverse problem solution. The derived angular distribution and spectra are briefly discussed.

The conventional definition of ground-level enhancement (GLE) events requires a detection of solar energetic particles (SEP) by at least two differently located neutron monitors. Some places are exceptionally well suitable for ground-based detection of SEP - high-elevation polar regions with negligible geomagnetic and reduced atmospheric energy/rigidity cutoffs. At present, there are two neutron-monitor stations in such locations on the Antarctic plateau: SOPO/SOPB (at Amundsen-Scott station, 2835 m elevation), and DOMC/DOMB (at Concordia station, 3233 m elevation). Since 2015, when the DOMC/DOMB station started continuous operation, a relatively weak SEP event that was not detected by sea-level neutron-monitor stations was registered by both SOPO/SOPB and DOMC/DOMB, and it was accordingly classified as a GLE. This would lead to a distortion of the homogeneity of the historic GLE list and the corresponding statistics. To address this issue, we propose to modify the GLE definition so that it maintains the homogeneity: A GLE event is registered when there are near-time coincident and statistically significant enhancements of the count rates of at least two differently located neutron monitors, including at least one neutron monitor near sea level and a corresponding enhancement in the proton flux measured by a space-borne instrument(s). Relatively weak SEP events registered only by high-altitude polar neutron monitors, but with no response from cosmic-ray stations at sea level, can be classified as sub-GLEs.

An important topic in the field of space weather is the precise assessment of the contribution of galactic cosmic rays and solar energetic particles on air crew exposure, specifically during eruptive events on the Sun. Here we present a model, a full chain analysis based on ground based measurements of cosmic rays with neutron monitors, subsequent derivation of particle spectral and angular characteristics and computation of dose rate. The model uses method for ground level enhancement analysis and newly numerically computed yield functions for conversion of secondary particle fluence to effective dose and/or the ambient dose equivalent. The precise an adequate information about the solar energetic particle spectra (SEPs) is the basis of the model. Since SEPs possess an essential isotropic part, specifically during the event onset, the angular characteristics should be also derived with good precision. This can be achieved using neutron monitor data during a special class of SEP events the ground level enhancements (GLEs). On the basis of the method representing a sequence of consecutive steps: computation of the NM asymptotic cones, NM rigidity cut-off and application of convenient optimization procedure, we derive the rigidity spectra and anisotropy characteristics of GLE particles. For the computation we use newly computed yield function of the standard sea-level 6NM64 neutron monitor for primary proton and alpha CR nuclei as well as 6NM64 yield function at altitudes ranging from the sea level up to 5000 m above the sea level. We derive the SEP spectra and pitch angle distributions in their dynamical development throughout the event. Subsequently on the basis of the derived spectra and angular characteristics and previously computed yield functions we calculate the effective dose and/or ambient dose equivalent during the GLE. Several examples are shown.

The possible effect of solar variability, accordingly cosmic rays variation on atmospheric physics and chemistry is highly debated over the last years. According to several recent models the induced by cosmic rays atmospheric ionization plays a key role in several different processes. At recent, an apparent effect on minor constituents and aerosols over polar regions during major solar proton events was observed. The ion production rate during ground level enhancements is a superposition of the contribution of cosmic rays with galactic and solar origin. The solar cycle 23 provided several strong ground level enhancements. The period of end October - beginning of November 2003 was characterized by a strong cosmic ray variability, namely a sequence of three GLEs was observed. In addition, there were several Forbush decreases, which led to a suppression of galactic cosmic ray flux. As a consequence the cosmic ray induced ion production in the atmosphere and the corresponding ionization effect were subject of dynamical changes. Here we compute the complex ion production due to cosmic rays during the Halloween sequence of ground level enhancements on October-November 2003 and we estimate the ionization effect. The spectral and angular characteristics of the solar protons are explicitly considered throughout the events as well their time evolution. The ionization effect during the period is computed at several altitudes above the sea level in a region with $R_{c}$ $\le$ 1 GV and $R_{c}$ $\le$ 2 GV.

The galactic cosmic rays are the main source of ionization in the Earth stratosphere and troposphere. The induced by primary cosmic ray particles ionization is important in various processes related to atmospheric physics and chemistry, specifically the minor constituents. The ion production in the atmosphere is enhanced compared to the average following major solar energetic particles events, specifically over the polar caps. During the solar cycle 23 we observed several strong ground level enhancements, one of the strongest among them been the Bastille day event on 14 July 2000. In the work presented here we apply a full Monte Carlo 3-D model in order to compute the cosmic ray induced ionization. The model is based on atmospheric shower simulation with the PLANETOCOSMICS code and the ion production rate is considered as a superposition of cosmic rays with galactic and solar origin. The ion production rate is computed as a function of the altitude above the sea level and the short, mid and long term ionization effect relative to the average due to galactic cosmic rays is computed.

A ground-based neutron monitor is a standard tool to measure cosmic ray variability near Earth, and it is crucially important to know its yield function for primary cosmic rays. Although there are several earlier theoretically calculated yield functions, none of them agrees with experimental data of latitude surveys of sea-level neutron monitors, thus suggesting for an inconsistency. A newly computed yield function of the standard sea-level 6NM64 neutron monitor is presented here separately for primary cosmic ray protons and $\alpha-$particles, the latter representing also heavier species of cosmic rays. The computations have been done using the GEANT-4 Planetocosmics Monte-Carlo tool and a realistic curved atmospheric model. For the first time, an effect of the geometrical correction of the neutron monitor effective area, related to the finite lateral expansion of the cosmic ray induced atmospheric cascade, is considered, that was neglected in the previous studies. This correction slightly enhances the relative impact of higher-energy cosmic rays (energy above 5--10 GeV/nucleon) in neutron monitor count rate. The new computation finally resolves the long-standing problem of disagreement between the theoretically calculated spatial variability of cosmic rays over the globe and experimental latitude surveys. The newly calculated yield function, corrected for this geometrical factor, appears fully consistent with the experimental latitude surveys of neutron monitors performed during three consecutive solar minima in 1976--77, 1986--87 and 1996--97. Thus, we provide a new yield function of the standard sea-level neutron monitor 6NM64 that is validated against experimental data.

New consistent and precise computations of the production of five cosmogenic radio-isotopes, 7Be, 10Be, 14C, 22Na and 36Cl, in the Earth's atmosphere by cosmic rays are presented in the form of tabulated yield functions. For the first time, a detailed set of the the altitude profiles of the production functions is provided which makes it possible to apply the results directly as input for atmospheric transport models. Good agreement with most of the earlier published works for columnar and global isotopic production rates is shown. Altitude profiles of the production are important, in particular for such tasks as studies of strong solar particle events in the past, precise reconstructions of solar activity on long-term scale, tracing air-mass dynamics using cosmogenic radio-isotopes, etc. As an example, computations of the $^{10}$Be deposition flux in the polar region are shown for the last decades and also for a period around 780 AD and confronted with the actual measurements in Greenland and Antarctic ice cores.

On the basis of neutron monitor data we estimate the energy spectrum, anisotropy axis direction and pitch-angle distribution of solar energetic particles during two major ground level enhancements (GLE 59 on 14 July 2000 and GLE 70 on 13 December 2006). For the analysis we use a newly computed neutron monitor yield function. The method consists of several consecutive steps: definition of the asymptotic viewing cones of neutron monitor stations considered for the data analysis by computations of cosmic ray particles propagation in a model magnetosphere with the MAGNETOCOSMICS code; computation of the neutron monitor model responses and derivation of the solar energetic particle characteristics on the basis of inverse problem solution. The pitch-angle distribution and rigidity spectrum of high-energy protons are obtained as function of time in the course of ground level enhancements. A comparison with previously reported results is performed and reasonable agreement is achieved. A discussion of the obtained results is included as well their possible application is included.

A new model CRAC:EPII (Cosmic Ray Atmospheric Cascade: Electron Precipitation Induced Ionization) is presented. The CRAC:EPII is based on Monte Carlo simulation of precipitating electrons propagation and interaction with matter in the Earth atmosphere. It explicitly considers energy deposit: ionization, pair production, Compton scattering, generation of Bremsstrahlung high energy photons, photo-ionization and annihilation of positrons, multiple scattering as physical processes accordingly. The propagation of precipitating electrons and their interactions with atmospheric molecules is carried out with the GEANT4 simulation tool PLANETOCOSMICS code using NRLMSISE 00 atmospheric model. The ionization yields is compared with an analytical parametrization for various energies of incident precipitating electron, using a flux of mono-energetic particles. A good agreement between the two models is achieved. Subsequently, on the basis of balloon-born measured spectra of precipitating electrons at 30.10.2002 and 07.01.2004, the ion production rate in the middle and upper atmosphere is estimated using the CRAC:EPII model

Recently an essential progress in development of physical models for cosmic ray induced ionization in the atmosphere is achieved. Basically, the models are full target, i.e. based on Monte Carlo simulation of an electromagnetic-muon-nucleon cascade in the atmosphere. In general, the contribution of proton nuclei in those models is highlighted, i.e. primary cosmic ray $\alpha$-particles and heavy nuclei are neglected or scaled to protons. The development of cosmic ray induced atmospheric cascade is sensitive to the energy and mass of the primary cosmic ray particle. The largest uncertainties in Monte Carlo simulations of a cascade in the Earth atmosphere are due to assumed hadron interaction models, the so-called hadron generators. In the work presented here we compare the ionization yield functions $Y$ for primary cosmic ray nuclei, such as protons, $\alpha$-particles, Oxygen and Iron nuclei, assuming different hadron interaction models. The computations are fulfilled with the CORSIKA 6.9 code using GHEISHA 2002, FLUKA 2011, UrQMD hadron generators for energy below 80 GeV/nucleon and QGSJET II for energy above 80 GeV/nucleon. The observed difference between hadron generators is widely discussed. The influence of different atmospheric parametrizations, namely US standard atmosphere, US standard atmosphere winter and summer profiles on ion production rate is studied. Assuming realistic primary cosmic ray mass composition, the ion production rate is obtained at several rigidity cut-offs - from 1 GV (high latitudes) to 15 GV (equatorial latitudes) using various hadron generators. The computations are compared with experimental data. A conclusion concerning the consistency of the hadron generators is stated.

We present the results of full new calculation of radiocarbon 14C production in the Earth atmosphere, using a numerical Monte-Carlo model. We provide, for the first time, a tabulated 14C yield function for the energy of primary cosmic ray particles ranging from 0.1 to 1000 GeV/nucleon. We have calculated the global production rate of 14C, which is 1.64 and 1.88 atoms/cm2/s for the modern time and for the pre-industrial epoch, respectively. This is close to the values obtained from the carbon cycle reservoir inventory. We argue that earlier models overestimated the global 14C production rate because of outdated spectra of cosmic ray heavier nuclei. The mean contribution of solar energetic particles to the global 14C is calculated as about 0.25% for the modern epoch. Our model provides a new tool to calculate the 14C production in the Earth's atmosphere, which can be applied, e.g., to reconstructions of solar activity in the past.

The general aim of this work is to obtain the lateral distribution of atmospheric Cherenkov light in extensive air showers produced by different primary particles precisely by. protons, Helium, Iron, Oxygen, Carbon, Nitrogen, Calcium, Silicon and gamma quanta in wide energy range at high mountain observation level of Chacaltaya cosmic ray station. The simulations are divided generally in two energy ranges 10GeV-10 TeV and 10 TeV-10 PeV. One large detector has been used for simulations, the aim being to reduce the statistical fluctuations of the obtained characteristics. The shape of the obtained lateral distributions of Cherenkov light in extensive air showers is discussed and the scientific potential for solution of different problems as well.

This work summarizes the results presented at 29th International Cosmic Ray Conference in Pune India. Generally the aim of this work is to obtain the lateral distribution of the atmospheric Cherenkov light in extensive air showers produced by different primary particles in wide energy range and at several observation levels and to fit the obtained lateral distributions. Using one large detector and partially modified CORSIKA code version are obtained the lateral distributions of Cherenkov light flux densities at several observation levels for different particle primaries precisely at 536 g/cm2 Chacaltaya, 700 g/cm2 Moussala and 875 g/cm2 Kartalska field observation levels for hadronic primaries and gamma quanta in the energy range 1011 eV-1016 eV. On the basis of the solution of over-determined inverse problem the approximation of these distributions is obtained. The same model function for all the primaries is used and for the different observation levels. The different model parameters for the different primaries and levels are obtained. The approximations are compared with polynomial approximation obtained with different method. Both approximations are used for detector efficiency estimation for the different experiments in preparation and estimation of the accuracy of the reconstruction techniques. At the same time inclined showers up to 30 degrees zenith angle are studied at Chacaltaya observation level. The obtained lateral distributions of vertical showers are compared with vertical showers model and the previously obtained approximation. This permits to adjust the reconstruction strategy and to study the model parameters behavior.

In this paper we present a new method for ground based gamma ray astronomy based only on atmospheric Cherenkov light flux analysis. The Cherenkov light flux densities in extensive air showers (EAS) initiated by different primaries are simulated, precisely primary protons in the energy range 100 GeV – 100 PeV, primary gamma quanta in the energy range 10 GeV – 10 PeV and primary iron nuclei in the energy range 10 TeV –100 PeV using the CORSIKA 6 code at high mountain observation level of 536 g/cm2. An approximation of lateral distribution of Cherenkov light in extensive air showers is obtained. The obtained approximation is a nonlinear fit such as Breit-Wigner with few parameters. A detailed study of the energy dependence of the proposed model function parameters is carried out and the fit of model parameters as a function of the primary energy is obtained as well. On the basis of the difference between the model parameters, precisely their behavior as a function of the energy, the strong no linearity of the model, we propose a method, which permits to make the distinction between gamma ray primaries from hadronic primaries. The possible backgrounds for ground based gamma ray astronomy are studied and the efficiency of the method is calculated. Different detector displacements are analyzed using the simulation of simplified mass spectrum of cosmic ray. The detector response is simulated taking into account the physical fluctuation of the processes, the statistical and possible systematic errors. The simulated and reconstructed events are compared and the accuracy in energy and primary mass reconstruction is obtained. Moreover the accuracy in shower axis determination is studied and criteria in shower axis position estimation are proposed.

A new method for energy and mass composition estimation of primary cosmic ray radiation based on atmospheric Cerenkov light flux in extensive air showers (EAS) analysis is proposed. The Cerenkov light flux in EAS initiated by primary protons and iron nuclei is simulated with CORSIKA 5.62 code for Chacaltaya observation level (536 g/cm2) in the energy range 10 TeV - 10 PeV. An adequate model, approximation of lateral distribution of Cerenkov light in showers is obtained. Using the proposed model and solution of overdetermined system of nonlinear equations based on Gauss Newton method with autoregularization, two different array detector arrangements are compared. The detector response for the detector sets is simulated. The accuracies in energy and shower axis determination are studied and the corresponding selection criteria are proposed. An approximation with nonlinear fit is obtained and the energy dependence of the proposed model function parameters is studied. The approximation of model parameters as function of the primary energy is carried out. This permits, taking into account the properties of the proposed method and model, to distinguish proton primaries from iron primaries. The detector response for the detector sets is simulated and the accuracies in energy determination are calculated. Moreover the accuracies in shower axis determination are studied and criteria in shower axis position estimation are proposed.