We discuss three of the known puzzles of the cosmic ray anisotropy in the PeV and sub-PeV energy region. They are: (1) the so called inverse anisotropy, (2) the irregularity in the energy dependence of the amplitude and phase of the first harmonic and (3) the contribution of the single source.

We give a short review of processes of stochastic acceleration in the Galaxy. We discuss: how to estimate correctly the number of accelerated particles, and at which condition the stochastic mechanism is able to generate power-law nonthermal spectra. We present an analysis of stochastic acceleration in the Galactic halo and discuss whether this mechanism can be responsible for production of high energy electrons there, which emit gamma-ray and microwave emission from the giant Fermi bubbles. Lastly, we discuss whether the effects of stochastic acceleration can explain the CR distribution in the Galactic disk (CR gradient).

Evidence for a local 'Single Source' of cosmic rays is amassing by way of the recent precise measurements of various cosmic ray energy spectra from the AMS-02 instrument. To observations of individual cosmic ray nuclei, electrons, positrons and antiprotons must now be added the determination of the boron-to-carbon ratio and the energy spectrum of lithium to 2000 GV with high precision. Our analysis leads us to claim that, with certain assumptions about propagation in the Galaxy, the results confirm our arguments regarding the presence of a local single source, perhaps, a supernova remnant (SNR). An attempt is made to determine some of the properties of this SNR and its progenitor star.

An analysis of the energy spectra of cosmic rays and particularly the precise data from the AMS-02 experiment support the view about the important role of the Local Bubble in the nearby interstellar medium. It is suggested that the bulk of cosmic rays below about 200 GV of rigidity (momentum/charge ratio) comes from the modest number of supernova remnants in the Local Bubble which appear to have occurred some $10^6$ years ago and contributed to its formation. At higher rigidities the contribution from a 'Local Source', a single supernova remnant generated some $10^5$ years ago seems to dominate up to, at least 1000 GV.

Evidence for a local 'Single Source' of cosmic rays is amassing by way of the recent precise measurements of various cosmic ray energy spectra from the AMS-02 instrument. To observations of individual cosmic ray nuclei, electrons, positrons and antiprotons must now be added the determination of the boron-to-carbon ratio and the energy spectrum of lithium to 2000 GV with high precision. Our analysis leads us to claim that, with certain assumptions about the propagation in the Galaxy, the results confirm our arguments regarding the presence of a local single source, perhaps, a supernova remnant (SNR). An attempt is made to determine some of the properties of this SNR and its progenitor star.

The dependence of the cosmic ray intensity on Galactocentric distance is known to be much less rapid than that to be thought-to-be sources: supernova remnants. This is an old problem ('the radial gradient problem') which has led to a number of possible 'scenarios'. Here, we use recent data on the supernova's radial distribution and correlate it with the measured HII electron temperature ({\em T}). We examined two models of cosmic ray injection and acceleration and in both of them the injection efficiency increases with increasing ambient temperature {\em T}. The increase is expected to vary as a high power of {\em T} in view of the strong temperature dependence of the tail of the Maxwell-Boltzmann distribution of particle energies. Writing the efficiency as proportional to $T^n$ we find $n\approx 8.4$. There is thus, yet another possible explanation of the radial gradient problem.

In view of the fact that the AMS-02 instrument has recently been used to make preliminary observations of the ratio of the antiprotons (\-{P}) to protons (P) in the primary cosmic radiation we have returned to our idea of signatures of a local recent supernova. We find that at the present level of accuracy there is no inconsistency between our predictions for the \-{P}/P ratio to some hundreds of GeV using the preliminary observations.

The methodical properties of the original difference method for the search of the anisotropy at the knee region of the primary cosmic radiation energy spectrum are analyzed. The main feature of the suggested method is a study of the difference in the EAS characteristics in different directions but not their intensity. It is shown that the method is stable to the random experimental errors and allows to separate the anomalies related to the laboratory coordinate system from the anomalies in the celestial coordinates. The method uses multiple scattering of the charge particles in the Galaxy magnetic fields to study the whole celestial sphere including the regions outside of the line of sight of the installation.

We present preliminary results obtained by a novel difference method for the study of the nature of the knee in the energy spectrum of the primary cosmic radiation. We have applied this method to data from the GAMMA experiment in Armenia. The analysis provides evidence for the possible existence of a nearby source of primary cosmic rays in the Southern hemisphere.

The existence of the cosmic ray Halo in our Galaxy has been discussed for more than half a century. If it is real it could help to explain some puzzling features of the cosmic ray flux: its small radial gradient, nearly perfect isotropy and the low level of the fine structure in the energy spectra of the various particles. All these features could be understood if: (a) the Halo has a big size (b) cosmic rays in the Halo have a unform spatial distribution and (c) the cosmic ray density in the Halo is comparable or even higher than that in the Galactic Disk. The main topic of the paper concerns the present status of the anisotropy and a model for its formation. In our model the extremely small amplitude of the dipole anisotropy is due to the dilution of the anisotropy in the Disk by the dominating isotropic cosmic rays from the Halo. Some minor deviations from complete isotropy in the sub-PeV and PeV energy regions point out to the possible contribution of the Single Source with the phase of its first harmonic opposite to the phase produced by the Disk.

Ultra-high energy cosmic rays interacting with the radiation fields in the universe cause electromagnetic cascades resulting in a flux of extragalactic gamma rays, detectable to some 100 GeV. Recent precise measurements of the extragalactic gamma ray flux by Fermi-LAT, coupled with estimates of the background from active galactic nuclei of various types, allows limits to be set on the cascade component. By comparison with prediction and, making various assumptions, ie taking a particular model, limits can be set on the maximum energy to which ultra-high energy particle can be accelerated. If our model is correct, it is unlikely that the maximum energy is above 100 EeV, in turn, the apparent 'GZK' cut-off in the measured ultra-high energy spectrum could instead be due to a fall-off in the intrinsic emergent particle spectrum. However, it is not plausible to be dogmatic at the present time because of uncertainty in many of the parameters involved. We have used recent estimates of the range of parameters and have found that although our model has parameters in the allowable ranges the uncertainties are so large that our result is not unique, although the method is satisfactory. The result must thus, so far, be taken as an indication only.

Recent work by Aplin and Lockwood [1] was interpreted by them as showing that there is a multiplying ratio of order 10$^{12}$ for the infra-red energy absorbed in the ionization produced by cosmic rays in the atmosphere to the energy content of the cosmic rays themselves. We argue here that the interpretation of the result in terms of infra-red absorption by ionization is incorrect and that the result is therefore most likely due to a technical artefact

All the components of cosmic rays have 'structure' in their energy spectra at some level, ie deviations from a simple power law, and their examination is relevant to the origin of the particles. Emphasis, here, is placed on the large-scale structures in the spectra of nuclei (the 'knee' at about 3 PeV), that of electrons and positrons (a shallow 'upturn' at about 100 GeV) and the positron to electron plus positron ratio (an upturn starting at about 5 GeV). Fine structure is defined as deviations from the smooth spectra which already allow for the large-scale structure. Search for the fine structure has been performed in the precise data on positron to electron plus positron ratio measured by the AMS-02 experiment. Although no fine structure is indicated, it could in fact be present at the few percent level.

We argue that the cosmic ray positron excess observed in ATIC-2, Fermi LAT, PAMELA, HESS and recently in the precision AMS-02 experiment can be attributed to the production in a local, middle-aged supernova remnant (SNR). Using the prediction of our model of cosmic ray acceleration in SNR we estimate that the SNR responsible for the observed positron excess is located between 250 and 320pc from the Sun and is 170-380 kyear old. The most probable candidate for such a source is the SNR which gave birth to the well-known Geminga pulsar, but is no longer visible. Other contenders are also discussed.

Several large cosmic ray (CR) detectors have recently provided data on the arrival directions of CR, which taken together with previous data recorded over many decades allow the amplitgude and phase of the first harmonic to be derived with reasonable precision and up to higher energies. We find a high degree of consistency amongst the various measurements. The new data indicate that at an energy above ~0.1 PeV a change of the CR anisotropy sets in. The amplitude of the first harmonic, which rises to 3 TeV, then diminishes and begins to rise again. The direction of the phase also changes to the opposite one. A measure of understanding follows from the use of two-dimensional maps of CR excesses over the mean background. When the energy of CR approaches the PeV region, the excess of CR moves from the Galactic Anti-Centre to the opposite direction of the Galactic Centre. The possible role of such potential CR sources as the supernovae Monogem Ring and Vela, which could help to explain some of the observed results, is discussed.

Recent measurements of cosmic gamma ray intensities up to TeV energies have been used to estimate the spectral shape of the parent cosmic ray particles present in the interstellar medium. The case is made for the particle spectrum in the Inner Galaxy being flatter than locally and in the Outer Galaxy. Of various possible explanations we make the case for the propagation of the particles being different in the more turbulent interstellar medium of the Inner Galaxy. The characteristic parameter alpha for the so called anomalous diffusion is expected to be less in the Inner Galaxy than that locally and the corresponding power law spectral exponent of the cosmic ray particles gamma will differ from that locally. Arguments are presented favouring a value of alpha lower than the local one by about Delta alpha = 0.2, the consequence is that Delta gamma = 0.1 for the parent particles.

Following on previous work by others, which gave evidence for few-day changes in the European Diurnal Temperature Range (DTR) apparently correlated with Cosmic Ray Forbush Decreases, we have made an independent study. We find no positive evidence. An analysis has also been made of the Fourier components of the time series of the DTR value (taken as deviations from a +/-10 day running mean). Evidence for a number of interesting periods is found, including one at about 27 days, albeit with a variability with time. The same period of solar irradiance (particularly in the UV) is favoured as the explanation.

The energy spectrum of the primary cosmic radiation in the energy range 1 - 100 PeV and the extensive air shower (EAS) characteristics obtained on the basis of the expanded data bank of the GAMMA experiment (Mt. Aragats, Armenia) are presented. With increased statistics we confirm our previous results on the energy spectrum. The spectral index above the knee is about -3.1, but at energies beyond 20 PeV a flattening of the spectrum is observed. The existence of the 'bump' at about 70 PeV is confirmed with a significance of more than 4{\sigma}. In the energy range of 10 - 100 PeV the shower age becomes energy independent and we observe a direct proportionality of the EAS size to the primary energy. This suggests an approximately constant depth of the EAS maximum in this energy range. This is evidence in favour of an increasing average mass of primary particles at energies above 20 PeV. The additional source scenario, which is a possible explanation of the 'bump' in the spectrum, also leads to the conclusion of increasing mass of the primary cosmic rays. A comparison with the data of other experiments is presented.

An analysis is made of the fine structure in the cosmic ray energy spectrum: new facets of present observations and their interpretation and the next step. It is argued that less than about 10% of the intensity of the helium `peak' at the knee at $\simeq 5PeV$ is due to just a few sources (SNR) other than the single source. The apparent concavity in the rigidity spectra of protons and helium nuclei which have maximum curvature at about 200 GV is confirmed by a joint analysis of the PAMELA, CREAM and ATIC experiments. The spectra of heavier nuclei also show remarkable structure in the form of `ankles' at several hundred GeV/nucleon. Possible mechanisms are discussed. The search for `pulsar peaks' has not yet proved successful.

The present status of the GAMMA facility consisting of an enlarged surface EAS array (116 of 1 m^2 scintillation detectors) and underground muon carpet (150 m^2 detectors) is described. The recent results on mass composition and energy spectrum at the energy region above the knee obtained on the basis of the GAMMA experimental data are presented. It is shown that the power law after the knee is not invariable like -3.1. The slope of the energy spectrum becomes more flat at E0>20 PeV. The strong irregularities of the energy spectrum at about 70-80 PeV are discussed in comparison with other experiments. The bump can be described by a two-component model of primary cosmic ray origin, where additional (pulsar) Fe components are included with a very flat power law energy spectrum.

It has recently been suggested that neutron stars inside the shells of young supernova remnants (SNR) are the sources of PeV cosmic rays and that the interaction of the particles with the radiation field in the SNR causes electron pair production, which has relevance to recent observations of 'high' positron fluxes. Furthermore, the character of the interaction is such that the well-known knee in the cosmic ray energy spectrum can be explained. Our examination of the mechanism leads us to believe that the required parameters of SN and pulses are so uncommon that the knee and positron fraction can only be explained if a single, local and recent SN - and associated pulsar - are concerned.

The inter-relation of clouds, solar irradiance and surface temperature is complex and subject to different interpretations. Here, we continue our recent work, which related mainly to the period from 1960 to the present, back to 1900 with further, but less detailed, analysis of the last 1000 years. The last 20 years is examined especially. Attention is given to the mean surface temperature, solar irradiance correlation, which appears to be present (with decadal smoothing) with a 22-year period; it is stronger than the 11-year cycle correlation with one year resolution. UV in the solar radiation is a likely cause. Cloud data are taken from synoptic observations back to 1952 and, again, there appears to be a correlation - with opposite phase for high and low clouds - at the 20-30y level. Particular attention is devoted to answering the question, 'what fraction of the observed increase in mean Global temperature (~0.7^oC) can be attributed to solar, as distinct from man-made, effects?' We conclude that a best estimate is 'essentially' all from 1900 to 1956 and <14% from 1956 to the present.

An analysis is made of the manner in which the cosmic ray intensity at Earth has varied over its existence and its possible relevance to both the origin and the evolution of life. Much of the analysis relates to the 'high energy' cosmic rays ($E>10^{14}eV;=0.1PeV$) and their variability due to the changing proximity of the solar system to supernova remnants which are generally believed to be responsible for most cosmic rays up to PeV energies. It is pointed out that, on a statistical basis, there will have been considerable variations in the likely 100 My between the Earth's biosphere reaching reasonable stability and the onset of very elementary life. Interestingly, there is the increasingly strong possibility that PeV cosmic rays are responsible for the initiation of terrestrial lightning strokes and the possibility arises of considerable increases in the frequency of lightnings and thereby the formation of some of the complex molecules which are the 'building blocks of life'. Attention is also given to the well known generation of the oxides of nitrogen by lightning strokes which are poisonous to animal life but helpful to plant growth; here, too, the violent swings of cosmic ray intensities may have had relevance to evolutionary changes. A particular variant of the cosmic ray acceleration model, put forward by us, predicts an increase in lightning rate in the past and this has been sought in Korean historical records. Finally, the time dependence of the overall cosmic ray intensity, which manifests itself mainly at sub-10 GeV energies, has been examined. The relevance of cosmic rays to the 'global electrical circuit' points to the importance of this concept.

Our model involving cosmic ray acceleration in supernova remnants has been used to predict cosmic ray intensities over long periods of time on a statistical basis. If, as is highly probable, extensive air showers caused by PeV cosmic rays are needed to initiate terrestrial lightning then past dramatic changes in PeV intensities may have had important biological effects. The model has been used to estimate the manner in which the PeV cosmic ray intensity at Earth has varied over the past tens of thousand years.

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.

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.

An update of the status of the knee in the cosmic ray energy spectrum at 3-4 PeV is presented. We argue that the evidence in favour of the presence of a 'single source' is even stronger than before.

We present an upper limit of galactic diffuse gamma-ray flux, measured with the GAMMA experiment at energy about 175 TeV. The results were obtained using selection of muon poor extensive air showers at mountain level (700 g/cm^2, Mt. Aragats, Armenia) for 5 GeV muon energy threshold. An upper limit for the differential flux at energy E=175+/-25 TeV is equal to (5.8-7.0)x10^-12 [erg m^2 s sr]^-1 for 95% confidence level.

The variation with time from 1956-2002 of the globally averaged rate of ionization produced by cosmic rays in the atmosphere is deduced and shown to have a cyclic component of period roughly twice the 11 year solar cycle period. Long term variations in the global average surface temperature as a function of time since 1956 are found to have a similar cyclic component. The cyclic variations are also observed in the solar irradiance and in the mean daily sun spot number. The cyclic variation in the cosmic ray rate is observed to be delayed by 2-4 years relative to the temperature, the solar irradiance and daily sun spot variations suggesting that the origin of the correlation is more likely to be direct solar activity than cosmic rays. Assuming that the correlation is caused by such solar activity, we deduce that the maximum recent increase in the mean surface temperature of the Earth which can be ascribed to this activity is $\lesssim14%$ of the observed global warming.

We present all-particle primary cosmic-ray energy spectrum in the 3-200 PeV energy range obtained by a multi-parametric event-by-event evaluation of the primary energy. The results are obtained on the basis of an expanded EAS data set detected at mountain level (700 g/cm^2) by the GAMMA experiment. The energy evaluation method has been developed using the EAS simulation with the SIBYLL interaction model taking into account the response of GAMMA detectors and reconstruction uncertainties of EAS parameters. Nearly unbiased (<5%) energy estimations regardless of a primary nuclear mass with an accuracy of about 15-10% in the 3-200 PeV energy range respectively are attained. An irregularity ('bump') in the spectrum is observed at primary energies of ~74 PeV. This bump exceeds a smooth power-law fit to the data by about 4 standard deviations. Not rejecting stochastic nature of the bump completely, we examined the systematic uncertainties of our methods and conclude that they cannot be responsible for the observed feature.

Recent results from the HESS gamma ray telescope have shown the presence of both a diffuse, extended, flux of gamma rays above ~0.4 TeV and discrete sources in and near the Galactic Centre. Here, we put forward a possible explanation in terms of the diffusion of cosmic ray protons from a succession of supernova remnants (SNR) in the SgrA* region of the Galaxy plus a contribution from SNR in the rest of the Galactic Centre Region, to be called the Galactic Centre Ridge (GCR). Protons are favoured over electrons because the mG magnetic fields in the Region will attenuate energetic electrons severely. Prominent features are the need for 'anomalous diffusion' of the protons in the whole region and the adoption of low efficiency for SNR acceleration in the high density regions. The latter is related by us to the well-known low 'cosmic ray gradient' in the Galaxy. A corroborating feature is the close correlation of inferred cosmic ray intensity with the smoothed intensity of 5 GHZ radio radiation. We attribute this to the presence of the SNR in the GCR.

The EGRET results for gamma ray intensities in and near the Galactic Plane have been analysed in some detail. Attention has been concentrated on energies above 1 GeV and the individual intensities in a $4^{\circ}$ longitude bin have been determined and compared with the large scale mean found from a nine-degree polynomial fit. Comparison has been made of the observed standard deviation for the ratio of these intensities with that expected from variants of our model. The basic model adopts cosmic ray origin from supernova remnants, the particles then diffusing through the Galaxy with our usual `anomalous diffusion'. The variants involve the clustering of SN, a frequency distribution for supernova explosion energies, and 'normal', rather than 'anomalous' diffusion. It is found that for supernovae of unique energy, and our usual anomalous diffusion, clustering is necessary, particularly in the Inner Galaxy. An alternative, and preferred, situation is to adopt the model with a frequency distribution of supernova energies. The results for the Outer Galaxy are such that no clustering is required. If their explosion energies are distributed then supernovae can be the origin of UHECR.

The EGRET results for gamma-ray intensities in and near the Galactic Plane have been analysed in some detail. Attention has been concentrated on energies above 1 GeV and the individual intensities in a $4^\circ$ longitude bin have been determined and compared with the large scale mean found from a nine-degree polynomial fit. Comparison has been made of the observed standard deviation for the ratio of these intensities with that expected from variants of our model. The basic model adopts cosmic ray origin from supernova remnants, the particles then diffusing through the Galaxy with our usual 'anomalous diffusion'. The variants involve the clustering of SN, a frequency distribution for supernova explosion energies, and 'normal', rather than 'anomalous' diffusion. It is found that for supernovae of unique energy, and our usual anomalous diffusion, clustering is necessary, particularly in the Inner Galaxy. An alternative, and preferred, situation is to adopt the model with a frequency distribution of supernova energies. The results for the Outer Galaxy are such that no clustering is required.

We study the effect of the stochastic character of supernova explosions on the anisotropy of galactic cosmic rays below the knee. We conclude that if the bulk of cosmic rays are produced in supernova explosions the observed small and nearly energy independent amplitude of the anisotropy and its phase are to the large extent determined by the history of these explosions in the vicinity of the solar system, namely by the location and the age of the supernova remnants, within a few kpc, which give the highest contribution to the total intensity at the present epoch. Among the most important factors which result in the small magnitude and the energy independence of the anisotropy amplitude are the mixed primary mass composition, the effect of the Single Source and the Galactic Halo. Special attention is given to the phase of the anisotropy. It is shown that the excessive flux from the Outer Galaxy can be due to the location of the solar system at the inner edge of the Orion Arm which has the enhanced density and rate of supernova explosions.

In view of recent developments attention is directed again at two aspects of the well known 'knee' in the cosmic ray energy spectrum at 3 PeV: the mass of the predominant particles at this energy and their source. It is inevitable in a subject such as this that ideas - and conclusions - evolve. Earlier, we had used a particular acceleration model and the nature of the local ISM to infer that the particles are mainly oxygen nuclei. Initially no specific source was identified. More recently, however, we have specified the Monogem Ring supernova remnant (SNR) as the likely source; this is at just the right distance and age and the energies are reasonable. Concerning the mass composition at the knee, a quantity more difficult to determine, recent direct measurements, which extend to higher energies than hitherto, show a likely flattening in the spectrum above 10^4 GeV/nucleon for He-nuclei, a flattening which, if extrapolated to higher energies, would meet the measured spectrum in the knee region. The other nuclei do not show this feature. He-nuclei in the knee region would also be marginally consistent with KASCADE EAS data, although there are serious problems with EAS mass estimates in that experiment. Concerning the acceleration, recent models applied to the Monogem Ring SNR allow a satisfactory explanation in terms of either oxygen or helium, but with the latter being a distinct possibility and perhaps more likely.

It is generally regarded that the bulk of cosmic rays originate in the Galaxy and that those below the 'knee' (the rapid steepening in the energy spectrum) at a few PeV come from Galactic supernovae, the particles being accelerated by the shocks in the supernova remnants. At higher energies, there are problems in that conventional SNR - which surely constitute the bulk of the sources - have a natural limit at a few tens of PeV (for iron nuclei). The question of the origin of particles above this limit is thus an open one. Here we examine a number of possibilities: a variety of supernovae and hypernovae, pulsars, a Giant Galactic Halo and an Extragalactic origin. A relevant property of any model is the extent to which it can provide the lack of significant irregularity of the energy spectrum. Although it is appreciated that spectral measurements are subject to systematic as well as random errors we consider that contemporary data are good enough to allow at least some progress in this field. In the search for origin above PeV energies we conclude that shocks in the Galactic Halo, whatever their source (Galactic wind, relativistic plasmoids - 'cannonballs', multiple shocks from supernovae etc.) are most likely, pulsrs such as B0656+14 and hypernovae come a close second although such a suggestion is not without its difficulties. What is most important is that trapping of particles in the Halo is needed to reduce irregularities of the energy spectra both below and above the 'knee' caused by the stochastic nature of supernova explosions and other potential (discrete) Galactic sources. We argue that precise experimental studies of spectral 'irregularities' will provide considerable help in the search for cosmic ray origin.

The determination of the primary cosmic ray mass composition from the longitudinal development of atmospheric cascades is still a debated issue. In this work we discuss several data analysis methods and show that if the entire information contained in the longitudinal profile is exploited, reliable results may be obtained. Among the proposed methods FCC ('Fit of the Cascade Curve'), MTA ('Multiparametric Topological Analysis') and NNA ('Neural Net Analysis') with conjugate gradient optimization algorithm give the best accuracy.

The 'Sling effect' appears when a fragment of a projectile nucleus emitted after its peripheral collision with a target nucleus is caused to rotate with a high spin. The spinning fragment has a deformed shape and looks like an oblate ellipsoid. Due to the virtual non-compressibility of nuclear matter, and the polarization of the spin in the plane transverse to the input momentum of the projectile nucleus such an ellipsoid has a reduced mean interaction cross-section compared with a non-spinning fragment which has a spherical shape. Purely geometrical arguments dictate that such an ellipsoidal nucleus should have additional fluctuations of cross-section even at a fixed impact parameter dependent on the orientation angle between the axis of the ellipsoid and the vector connecting the centers of projectile and target nucleus. The number of 'wounded nucleons' in the projectile nucleus participating in the interaction correlates strongly with the interaction cross-section. All these effects lead to a non-exponential attenuation of fragments and increased probability for a fragment to penetrate down to a larger depth in the absorber, than normal. If the sling effect appears in the interaction of a primary cosmic ray nucleus with nuclei in the atmosphere the induced atmospheric cascade will have a slower attenuation, and thereby can help to reduce some important inconsistencies in the interpretation of the existing experimental data on extensive air showers, observed in the lower half of the atmosphere. The paperr gives numerical estimates of the sling effect.

It is evident that models of the knee should match the observational phenomenology. In this talk I discuss a few aspects of phenomenology, which are important not only for the understanding of the knee origin, but also for the general problem of the origin of cosmic rays. Among them are the shape of the energy spectrum, its irregularity, the sharpness of the knee and its fine structure. The classification of models is given and some examples of the most recent models are discussed. The most probable conclusion deduced from this examination is that the knee has an astrophysical origin and the so called 'source' models of the knee are most likely among them.

The detemination of the primary cosmic ray mass composition from the longitudinal development of atmospheric cascades is still an open problem. In this work we propose a new method of the multiparametric topological analysis and show that if both X$_{max}$ - the depth of shower maximum and N$_{max}$ - the number of charged particles in the shower maximum are used, reliable results can be obtained.

The determination of the primary cosmic ray mass composition from the characteristics of extensive air showers (EAS), obtained at an observation level in the lower half of the atmosphere, is still an open problem. In this work we propose a new method of the Multiparametric Topological Analysis and show its applicability for the determination of the mass composition of the primary cosmic rays at the PeV energy region.

Recent findings indicate that the Monogem Ring supernova remnant and the associated pulsar PSR B0656+14 may be the 'Single Source' responsible for the formation of the sharp knee in the cosmic ray energy spectrum at $\sim$3PeV. We estimate the contribution of the pulsar B0656+14 to the cosmic rays in the PeV region and conclude that the pulsar cannot contribute more than 15% to the cosmic ray intensity at the knee. Therefore it cannot be the dominant source there and an SNR is still needed. We also examine the possibility of the pulsar giving the peak of the extensive air shower intensity observed from the region inside the Monogem Ring. The estimates of the gamma-ray flux produced by cosmic ray particles from this pulsar indicate that it can be the source of the observed peak, if the particles were confined within the SNR during a considerable fraction of its total age. We also estimate the contribution of Geminga and Vela pulsars to cosmic rays at the knee.

Recent findings indicate that the Monogem Ring supernova remnant (SNR) and the associated pulsar B0656+14 may be the 'Single Source' responsible for the knee in the cosmic ray (CR) energy spectrum at ~3 PeV. We estimate the contribution of this pulsar to CR in the PeV region. We conclude that although the pulsar can contribute to the formation of the knee, it cannot be the domimant source and a SNR is still needed. We also examine the possibility of the pulsar giving the peak of the extensive air shower (EAS) intensity observed from the region inside the Monogem Ring. If the experimental EAS results concerning a narrow source are confirmed, they can be important, since they give evidence: (i) for the acceleration of protons and heavier nuclei by the pulsar; (ii) for the existence of the confinement mechanism in SNR; (iii) that CR produced by the Monogem Ring SNR and associated pulsar B0656+14 were released recently giving rise to the formation of the knee and the observed narrow peak in the EAS intensity; (iv) for the Monogem Ring and the associated pulsar B0656+14 being identified as the Single Source proposed in our Single Source Model of the knee. A number of predictions of the examined scenario are made.

Some six years ago, we (Erlykin and Wolfendale, 1997) proposed the `single source model' in which a local, recent supernova remnant (SNR) was responsible for the `knee' in the cosmic ray (CR) energy spectrum at ~3PeV. Stimulated by the paper by Bhadra (2002), which drew attention to a possible gamma ray signature of this local remnant, we now study the situation for the local source and we conclude that, in contrast to Bhadra's conclusion, the non-observation of this remnant is understandable - at least using our SNR model. It is due to the fact that this SNR, being local, develops in the local hot interstellar medium (HISM) with its low density of gas and also being nearby it will be an extended source occupying up to 40$^\circ$ of the sky and thus indistinguishable from the background.

The idea that the shocks from supernova remnants (SNR) expanding into the interstellar medium (ISM) accelerate CR is still popular one but a number of authors have drawn attention to the fact that the experimental evidence for the presence of gamma rays from the expected interaction of CR with gas in the remnant is poor. Indeed, it is claimed that many SNR are not `seen' in GeV or TeV gamma rays, whereas `they should have been'. We have looked at this problem and we conclude that the idea of CR production in SNR cannot be faulted in this way, if the evacuation of ambient gas by the stellar wind of the progenitor star and, frequently, by associated earlier close-by SN, is taken into account; such phenomena are expected for the important Type II SN which result from very massive stars. Other SNR have, apparently, been seen and the interaction of SNR-accelerated particles with adjacent molecular clouds has been deemed responsible. However, we worry about this interpretation because of the slow progress of the SNR shock through such clouds, although electron-effects may, indeed, contribute. This paper is devoted mainly to the visibility of nearby SNR in gamma rays although many of the arguments also relate to remote SNR. For the nearby SNR another problem enters the scene: the large angular spread of the remnant. It is especially important for the old SNR, where cosmic rays have already diffused to a large distance from the SNR center. We have also examined the effect of the 'anomalous diffusion' of CR propagation in the non-homogeneous interstellar medium on the visibility of SNR for gamma rays of different energies.

An analysis has been made of the present situation with respect to the high energy hadron-nucleus and nucleus-nucleus interaction models as applied to cosmic rays. As is already known, there are inconsistencies in the interpretation of experimental data on the primary mass composition, which appear when different EAS components are used for the analyses, even for the same experiment. In the absence of obvious experimental defects, there is a clear need for an improvement to the existing models; we argue that the most promising way is to enlist two effects which should be present in nucleus- nucleus collisions but have not been allowed for before. These are: a few percent energy transfer into the EAS electromagnetic component due to electron- positron pair production or electromagnetic radiation of the quark-gluon plasma and a small slow-down of the cascading process in its initial stages associated with the extended lifetime of excited nuclear fragments. The latter process displaces the shower maximum deeper into the atmosphere.

An analysis is made of the masses and spectral features for cosmic rays in the PeV region, insofar as they have a bearing on the problem of the interaction of cosmic ray particles. In our Single Source Model we identified two 'peaks' seen in a summary of the world's data on primary spectra, and claimed that they are probably due to oxygen and iron nuclei from a local, recent supernova. In the present work we examine other possible mass assignments. We conclude that of the other possibilities only Helium and Oxygen (instead of O and Fe) has much chance of success; the original suggestion is still preferred, however. Concerning our location with respect to the SNR shell, the analysis suggests that we are close to it - probably just inside.

The problem of the origin of cosmic rays in the shocks produced by supernova explosions at energies below the so called 'knee' (at ~3*10$^6$ GeV) in the energy spectrum is addressed, with special attention to the propagation of the particles through the inhomogenious interstellar medium and the need to explain recent anisotropy results, [1]. It is shown that the fractal character of the matter density and magnetic field distribution leads to the likelihood of a substantial increase of spatial fluctuations in the cosmic ray energy spectra. While the spatial distribution of cosmic rays in the vicinity of their sources (eg. inside the Galactic disk) does not depend much on the character of propagation and is largely determined by the distribution of their sources, the distribution at large distances from the Galactic disk depends strongly on the character of the propagation. In particular, the fractal character of the ISM leads to what is known as 'anomalous diffusion' and such diffusion helps us to understand the formation of Cosmic Ray Halo. Anomalous diffusion allows an explanation of the recent important result from the Chacaltaya extensive air shower experiment [1], viz. a Galactic Plane Enhancement of cosmic ray intensity in the Outer Galaxy, which is otherwise absent for the case of the so-called 'normal' diffusion. All these effects are for just one reason: anomalous diffusion emphasizes the role of local phenomena in the formation of cosmic ray characteristics in our Galaxy and elsewhere.

An analysis has been made of the present situation with the high energy hadron-nucleus and nucleus-nucleus interaction models. As is already known there are inconsistencies in the interpretation of experimental data on the primary mass composition, which appear when different EAS components are used for the analyses, even for the same experiment. In the absence of obvious experimental defects, there is a clear need for an improvement to the existing models; we argue that the most promising way is to introduce two effects which should be present in nucleus-nucleus collisions and have not been allowed for before. These are: a few percent energy transfer into the EAS electromagnetic component due to electron-positron pair production or electromagnetic radiation of quark-gluon plasma and a small slow-down of the cascading process in its initial stages associated with the extended lifetime of excited nuclear fragments. The latter process displaces the shower maximum deeper into the atmosphere.

A sudden steepening of the cosmic-ray energy spectrum (the knee) is observed at an energy of about 3 PeV (1 PeV = 10$^{15}$eV). The recent results on extensive air showers allow us to conclude that: a) the knee has an astrophysical origin; b) the 'sharpness' and the fine structure of the knee rule out 'Galactic Modulation' as the origin of the knee; c) most likely the knee is the result of the explosion of a single, recent, nearby supernova.

A sudden steepening of the cosmic-ray energy spectrum (the knee) is observed at an energy of about 3 PeV (1 PeV = $10^{15}$ eV). The experimental study of the PeV cosmic rays has intensified greatly during the last 3 years. The recent results on extensive air showers allow us to conclude that: a) the knee has an astrophysical origin; b) the `sharpness' and the fine structure of the knee rule out `Galactic Modulation' as the origin of the knee; c) most likely the knee is the result of the explosion of a single, recent, nearby supernova.

The recently published precise spectrum of cosmic ray protons from the Alpha Magnetic Spectrometer has been examined in some detail from the standpoint of a search for deviations from a smooth, simple, power law. We find a significant excess ($\simeq$ 10% over $\sim$ 0.3 interval in logE) centered on 50 GeV in the published data. It is possible that the 'unfolding technique' adopted by the experimenters causes an overestimate of the excess but it is difficult to reduce it much below about (5 $\pm$ 3)%. We have examined other recent data, too. There is also evidence here, for an excess in the same energy region, although of only (1.6 $\pm$ 0.9)%. A value of (3 $\pm$ 1.5)% would be consistent with all the data. There are hints of similar excesses for heavier nuclei and electrons. Possible explanations are put forward for an excess, should it prove to be genuine.