We present a study of the multiplicities, of the lateral distributions and of the ratio of the electromagnetic to the hadronic components in the air showers, generated by the collision in the atmosphere of an incoming high energy cosmic ray and mediated by the formation of a mini black hole, predicted in TeV scale gravity models with large extra dimensions. The analysis is performed via a large scale simulation of the resulting cascades over the entire range ($10^{15}-10^{19}$eV) of ultra high initial energies, for several values of the number of large extra dimensions, for a variety of altitudes of the initial interaction and with the energy losses in the bulk taken into account. The results are compared with a representative of the standard events, namely the shower due to the collision of a primary proton with a nucleon in the atmosphere. Both the multiplicities and the lateral distribution of the showers show important differences between the two cases and, consequently, may be useful for the observational characterization of the events. The electromagnetic/hadronic ratio is strongly fluctuating and, thus, less decisive for the altitudes considered.

We present a study of the decay of an intermediate mini black hole at the first impact of a cosmic ray particle with the atmosphere, in the context of D-brane world scenarios with TeV scale gravity and large extra dimensions. We model the decay of the black hole using the semiclassical approximation and include the corrections coming from energy loss into the bulk. Extensive simulations show that mini black hole events are characterized by essentially different multiplicities and wider lateral distributions of the air showers as a function of the energy of the incoming primary, as compared to standard events. Implications for their detection and some open issues on their possible discovery are also briefly addressed.

The supersymmetric evolution of the fragmentation functions (or timelike evolution) within N=1 $QCD$ is discussed and predictions for the fragmentation functions of the theory (into final protons) are given. We use a backward running of the supersymmetric DGLAP equations, using a method developed in previous works. We start from the usual QCD parameterizations at low energy and run the DGLAP back, up to an intermediate scale -assumed to be supersymmetric- where we switch-on supersymmetry. From there on we assume the applicability of an N=1 supersymmetric evolution (ESAP). We elaborate on possible application of these results to High Energy Cosmic Rays near the GZK cutoff.