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.

Large scale airshower simulations around the GZK cutoff are performed. An extensive analysis of the behaviour of the various subcomponents of the cascade is presented. We focus our investigation both on the study of total and partial multiplicities along the entire atmosphere and on the geometrical structure of the various cascades, in particular on the lateral distributions. The possibility of detecting new physics in Ultra High Energy Cosmic Rays (UHECR) at AUGER is also investigated. We try to disentangle effects due to standard statistical fluctuations in the first proton impact in the shower formation from the underlying interaction and comment on these points. We argue that theoretical models predicting large missing energy may have a chance to be identified, once the calibration errors in the energy measurements are resolved by the experimental collaborations, in measurements of inclusive multiplicities.

Stable Superstring Relics (SSR) provide some of the candidates for the possible origin of the Ultra High Energy Cosmic Rays (UHECR). After a brief overview of the motivations for introducing such relics, we address the question whether statistical fluctuations in the formation of the air showers generated by the primary spectrum of protons can be separated from a possible signal of new physics hidden in the first impact with the atmosphere. Our results are generated by using minimal modifications in the cross section of the primaries, and using available simulation codes used by the experimental collaborations. The results indicate that substantial increases in the cross section of the first impact, possibly due to new interactions, are unlikely to be detected in geometrical and/or variations of multiplicities in the cascade.