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 have used recently published data from the Pierre Auger Observatory and the Telescope Array Project to conclude some inferences concerning the origin and composition of ultra high-energy cosmic rays (UHECR). We advocate for model the UHECR flux as a combination of the Galactic and Extragalactic components exchanging the dominant role at the famous feature of the energy spectrum called 'the ankle'. We put a special emphasis on the individual experiment energy resolution and possible biases it introduces. We have eventually combined the data from both big experiments to the one 'world average' UHECR spectrum and we used it, supported by the analysis of the distribution of the Extensive Air Shower (EAS) development maximum depth to find the UHECR mass composition of both: Galactic and Extragalactic components. We studied the dependence of the on mass composition on the cosmic ray source distribution in space. We present results for some models existing in the literature. We found the evidence of the deficit of Extragalactic sources in the vicinity of the Galaxy of the the UHECR, what is in concordance with the small scale anisotropy analysis results.

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.

We have updated our analysis of the 9-year WMAP data using the collection of polarization maps looking for the presence of additional evidence for a finite 'cosmic ray foreground' for the CMB. We have given special attention to high Galactic latitudes, where the recent BICEP2 findings were reported. The method of examining the correlation with the observed gamma ray flux proposed in our earlier papers and applied to the polarization data shows that the foreground related to cosmic rays is still observed even at high Galactic altitudes and conclusions about gravitational waves are not yet secure. Theory has it that there is important information about inflationary gravitational waves in the fine structure of the CMB polarization properties (polarization vector and angle) and it is necessary to examine further the conclusions that can be gained from studies of the CMB maps, in view of the disturbing foreground effects.

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.

The Galactic cosmic ray (GCR) intensity has been postulated by others to vary cyclically with a peak to valley ratio of ~3:1, as the Solar System moves from the Spiral Arm to the Inter-Arm regions of the Galaxy. These intensities have been correlated with global temperatures and used to support the hypothesis of GCR induced climate change. In this paper we show that the model used to deduce such a large ratio of Arm to Interarm GCR intensity requires unlikely values of some of the GCR parameters, particularly the diffusion length in the interstellar medium, if as seems likely to be the case, the diffusion is homogeneous. Comparison is made with the existing gamma ray astronomy data and this also indicates that the ratio is not large. The variation in the intensity is probably of order 10-20% and should be no more than 30% as the Solar System moves between these two regions, unless the conventional parameters of the GCR are incorrect. In addition we show that the variation of the GCR intensity, as the trajectory of the Solar System oscillates about the Galactic Plane, is too small to account for the extinctions of species as has been postulated unless, again, conventional assumptions about the GCR parameters are not correct.

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.

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 origin of the prominent 'knee' in the cosmic ray energy spectrum at an energy of several PeV is still uncertain. A recent mechanism has shown promise, however; this involves particles from a very young pulsar interacting with the radiation field from a very young supernova remnant. The ensuing nuclear reaction of the particles with the photons by way of e+e- production then causes the characteristic knee. In an earlier paper we argued that the mechanism would imply only one source of a very rare type - if it were to explain the spectral shape. Here we examine the mechanism in more detail and conclude that for even a single source to work its characteristics would need to be so unusual that the mechanism would not be possible for any known type of pulsar-supernova combination.

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.

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.

A decrease in the globally averaged low level cloud cover, deduced from the ISCCP infra red data, as the cosmic ray intensity decreased during the solar cycle 22 was observed by two groups. The groups went on to hypothesise that the decrease in ionization due to cosmic rays causes the decrease in cloud cover, thereby explaining a large part of the presently observed global warming. We have examined this hypothesis to look for evidence to corroborate it. None has been found and so our conclusions are to doubt it. From the absence of corroborative evidence, we estimate that less than 23%, at the 95% confidence level, of the 11-year cycle change in the globally averaged cloud cover observed in solar cycle 22 is due to the change in the rate of ionization from the solar modulation of cosmic rays.

A very recent observation by the Auger Observatory group claims strong evidence for cosmic rays above 56 EeV being protons from Active Galactic Nuclei. If, as would be expected, the particles above the ankle at about 2 EeV are almost all of extragalactic origin then it follows that the characteristics of the nuclear interactions of such particles would need to be very different from conventional expectation -- a result that follows from the measured positions of 'shower maximum' in the Auger' work. Our own analysis gives a different result, viz that the detected particles are still 'massive' specifically with a mean value of <ln A> = 2.2 +- 0.8. The need for a dramatic change in the nuclear physics disappears.

It has been claimed by others that observed temporal correlations of terrestrial cloud cover with `the cosmic ray intensity' are causal. The possibility arises, therefore, of a connection between cosmic rays and Global Warming. If true, the implications would be very great. We have examined this claim to look for evidence to corroborate it. So far we have not found any and so our tentative conclusions are to doubt it. Such correlations as appear are more likely to be due to the small variations in solar irradiance, which, of course, correlate with cosmic rays. We estimate that less than 15% of the 11-year cycle warming variations are due to cosmic rays and less than 2% of the warming over the last 35 years is due to this cause.

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.

A detailed search has been made for evidence of foreground contributions to the Cosmic Microwave Background (CMB) WMAP, such foregrounds being related to our Galaxy. We find remarkable results. On the largest angular scales we find significant differences between the power in the fluctuations for positive and negative Galactic latitudes and between the four Galactic Quadrants. There are also differences between the power spectrum at latitudes within 10 deg of the Plane and at higher latitudes. The `slope' of the power spectrum has similar variations. In an effort to find the origin of these Galactic-style effects we have examined the evidence from Galactic gamma rays, specifically from the EGRET instrument. A correlation is found between gamma ray intensities and the CMB and other cosmic ray indicators. Most of the large scale Galactic asymmetries (eg north, south difference and Quadrant variations) have analogues in cosmic ray asymmetries and also in some other Galactic properties, such as the column density of gas. Thus, it is possible to hypothesise on direct cosmic ray-induced contributions, although it may be that cosmic rays are simply the indicators of Galactic `conditions' which are influencing the residual CMB fluctuations. Irrespective of the actual cause of the correlations we have endeavoured to extrapolate to the situation where the residual foreground is minimised. The effect on the usually-derived cosmological properties - universal matter content, ionization, the tensor contribution, tilt and curvature - is briefly examined.

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.

We have argued [J. Szabelski et al. (2002)] that the well-known `ankle' in the cosmic ray energy spectrum, at logE (eV) ~ 18.7-19.0, marks the transition from mainly Galactic sources at lower energies to mainly extragalactic above. Recently, however, there have been claims for lower transitional energies, specifically from logE (eV) ~ 17.0 [G. Thompson et al. (2004)] via 17.2-17.8 [V.S. Berezinsky et al. (2004)] to 18.0 [A.M. Hillas (2004)]. In our model the ankle arises naturally from the sum of simple power law-spectra with slopes differing by ~ 1.8; from differential slope -3.8 for Galactic particles (near logE = 19) to ~ -2.0 for extragalactic sources. In the other models, on the other hand, the ankle is intrinsic to the extragalactic component alone, and arises from the shape of the rate of energy loss versus energy for the (assumed) protons interacting with the cosmic microwave background (CMB). Our detailed analysis of the world's data on the ultra-high energy spectrum shows that taken together, or separately, the resulting mean sharpness of the ankle (second difference of the log(intensity*E^3) with respect to logE) is consistent with our `mixed' model. For explanation in terms of extragalactic particles alone, however, the ankle will be at the wrong energy - for reasonable production models and of insufficient magnitude if, as seems likely, there is still a significant fraction of heavy nuclei at the ankle energy.

It is often claimed that there should be a 'GZK cut-off" in the flux of extragalactic cosmic rays, arising from interactions between the cosmic rays and the cosmic micro-wave background photons (e.g. 1 and 2). Some experiments (3 and 4) show particles of even higher energy than this value and this has led to claims for exotic processes (e.g. 5 and 6). We contend that such claims are unnecessary - there is no predicted cut-off, rather a continuation of the injection spectrum at reduced intensity. We have combined the world's data and shown that the prediction for a rather flat injection spectrum (exponent: 1.9 - 2.2) in the case of universal particle injection provides a reasonable fit to the data. Conventional forms for the particle attenuation in the intergalactic medium (e.g. 7 and 8) have been assumed. Either protons or iron nuclei (or a mixture) will suffice. Attention is drawn to another aspect, too, that of the losses on the infra-red radiation which may be intense near to strong sources and for sources in galaxy clusters. The attendant magnetic fields near the sources leads to significantly long diffusion times through the strong infra-red fields. Two 'case histories' are considered.

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.

We have cross-correlated the WMAP data with several surveys of extragalactic sources and find evidence for temperature decrements associated with galaxy clusters and groups detected in the APM Galaxy Survey survey and the ACO catalogue. We interpret this as evidence for the thermal Sunyaev-Zeldovich (SZ) effect from the clusters. Most interestingly, the signal may extend to ~1 deg around both groups and clusters and we suggest that this may be due to hot `supercluster' gas. We have further cross-correlated the WMAP data with clusters identified in the 2MASS galaxy catalogue and also find evidence for temperature decrements there. From the APM group data we estimate the mean Compton parameter as y(z<0.2)=7x10^(-7). We have further estimated the gas mass associated with the galaxy group and cluster haloes. Assuming temperatures of 5 keV for ACO clusters and 1 keV for APM groups and clusters, we derive average gas masses. Using the space density of APM groups we then estimate Omega_gas. For an SZ extent of theta_max=20', kT=1 keV and h=0.7, this value of Omega_gas=0.04 is consistent with the standard value of Omega_baryon=0.044 but if the indications we have found for a more extended SZ effect out to theta=60' are confirmed, then higher values of Omega_gas will be implied. Finally, the contribution to the WMAP temperature power spectrum from the extended SZ effect around the z<0.2 APM+ACO groups and clusters is 1-2 orders of magnitude lower than the l=220 first acoustic peak. But if a similar SZ effect arises from more distant clusters then this contribution could increase by a factor >10 and then could seriously affect the WMAP cosmological fits.

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.