Recently the CERN ALICE experiment, in its dedicated cosmic ray run, observed muon bundles of very high multiplicities, thereby confirming similar findings from the LEP era at CERN (in the CosmoLEP project). Originally it was argued that they apparently stem from the primary cosmic rays with a heavy masses. We propose an alternative possibility arguing that muonic bundles of highest multiplicity are produced by strangelets, hypothetical stable lumps of strange quark matter infiltrating our Universe. We also address the possibility of additionally deducing their directionality which could be of astrophysical interest. Significant evidence for anisotropy of arrival directions of the observed high multiplicity muonic bundles is found. Estimated directionality suggests their possible extragalactic provenance.

It is known from some time that to explain some apparently {\it exotic events} observed in cosmic ray experiments, a deeply penetrating component in cosmic rays is needed, which, in turn, leads frequently to the nonexponential distributions of some observables. We argue that such unexpected behavior of some cosmic rays data can be most naturally explained by attributing them to the possible fluctuations of cross sections.

We present arguments aiming to reconcile the apparently contradictory results concerning the chemical composition of cosmic rays of highest energy, coming recently from Auger and HiRes collaborations. In particular, we argue that the energy dependence of the mean value and root mean square fluctuation of shower maxima distributions observed by the Auger experiment are not necessarily caused by the change of nuclear composition of primary cosmic rays.

We present arguments aiming at reconciling apparently contradictory results concerning the chemical composition of cosmic rays of highest energy, coming recently from the Auger and HiRes collaborations. In particular, we argue that the energy dependence of the mean value and root mean square fluctuation of shower maxima distributions observed by the Auger experiment are not necessarily caused by the change of nuclear composition of primary cosmic rays. They could also be caused by the change of distribution of the first interaction point in the cascade. A new observable, in which this influence is strongly suppressed, is proposed and tested.

Energy spectrum of cosmic rays (CR) exhibits power-like behavior with very characteristic "knee" structure. We consider a generalized statistical model for the production process of cosmic rays which accounts for such behavior in a natural way either by assuming existence of temperature fluctuations in the source of CR, or by assuming specific temperature distribution of CR sources. Both possibilities yield the so called Tsallis statistics and lead to the power-like distribution. We argue that the "knee" structure arises as result of abrupt change of fluctuations in the source of CR. Its possible origin is briefly discussed.

We derive the constitutive equations of causal relativistic dissipative hydrodynamics ($d$-hydrodynamics) from perfect nonextensive hydrodynamics ($q$-hydrodynamics) using the nonextensive/dissipative correspondence (NexDC) proposed by us recently. The $q$-hydrodynamics can be thus regarded as a possible model for the $d$-hydrodynamics facilitating its application to high energy multiparticle production processes. As an example we have shown that applying the NexDC to the perfect 1+1 $q$-hydrodynamics, one obtains a proper time evolution of the bulk pressure and the Reynolds number.

The showers of cosmic rays entering the Earth's atmosphere are main sources of information on cosmic rays and are also believed to provide information on elementary interactions at energies not accessible to accelerators. In this context we would like first to remind the role of inelasticity K and elementary cross section $\sigma$ and then argue that similar in importance are fluctuations of different observables. The later will be illustrated by multiplicity fluctuations in hadronic and nuclear collisions.

Recently new data from the Cosmo-LEP project appeared, this time from DELPHI detector. They essentially confirm the findings reported some time ago by ALEPH, namely the appearence of bundles of muons with unexpectedly high multiplicities, which so far cannot be accounted by present day models. We argue, using arguments presented by us some time ago, that this phenomenon could be regarded as one more candidate for the presence in the flux of cosmic rays entering the Earth's atmosphere from outer space nuggets od Strange Quark Matter (SQM) in form of so called strangelets.

Assuming that cosmic rays entering the Earth's atmosphere contain a small admixture of nuggets of strange quark matter in form of strangelets one can explain a number of apparently "strange" effects observed in different cosmic rays experiments. We shall demonstrate here that the mass spectrum of such strangelets filles the "nuclear desert" gap existing between the heaviest elements observed in Universe and the next "nuclear-like objects" represented by neutron and strange stars.

The apparently observed power-like abundance of elements in Universe is discussed in more detail with special emphasis put on the strangelets.

We provide description of rapidity spectra of particles produced in p-anti-p collisions using anomalous diffusion approach to account for their non-equilibrium character. In particular, we exhibit connection between multiproduction processes and anomalous diffusion described through the nonlinear Focker-Planck equation with nonlinearity given by the nonextensivity parameter q describing the underlying Tsallis q-statistics and demonstrate how it leads to the Feynman scaling violation in these collisions. The q parameter obtained this way turns out to be closely connected to parameter 1/k converting the original poissonian multiplicity distribution to its observed Negative Binomial form. The inelasticity of reaction has been also calculated and found to slightly decrease with the increasing energy of reaction. Keywords: inelasticity, rapidity spectra, nonextensivity

We propose simple Monte Carlo method for calculating parton distribution in nuclei. The EMC effect is automatically included by means of two parameters, which characterize the change of nuclear pion field. Good agreement with experimental data in the broad range of variable x is obtained.

We discuss the possible imprints of Strange Quark Matter (SQM) in cosmic ray data. In particular, we investigate the propagation of SQM through the atmosphere and discuss: i) direct candidates for strangelets, ii) exotic events interpreted as signals of SQM and iii) muon bundles and delayed neutrons in Extensive Air Showers. The physics and astrophysics of SQM is shortly reviewed. We point out the possibility that extreme energy cosmic rays are the results of the decay of unstable primordial objects. Finally, the abundance of possible candidates for strangelets and their mass spectrum are estimated and compared with the astrophysical limits and prospects of the possible observation of SQM in accelerator experiments are outlined.

The hypothesis is discussed that muon bundles of extremaly high multiplicity observed recently by ALEPH detector (in the dedicated cosmic-ray run) can originate from the strangelets colliding with the atmosphere.

We argue that atmospheric cascades can be regarded as example of the self-organized criticality and studied by using Levy flights and nonextensive approach. It allows us to understand the scale-invariant energy fluctuations inside cascades in a natural way.

The requirements for observing Centauro-like phenomena and their role as possible probes of deeply penetrating component in Cosmic Rays are discussed.