SENECA is a hybrid air shower simulation written by H. Drescher that utilizes both Monte Carlo simulation and cascade equations. By using the cascade equations only in the high energy portion of the shower, where the shower is inherently one-dimensional, SENECA is able to utilize the advantages in speed from the cascade equations yet still produce complete, three dimensional particle distributions at ground level which capture the shower to shower variations coming from the early interactions. We present a comparison, on an event by event basis, of SENECA and CORSIKA, a well trusted MC simulation code. By using the same first interaction in both SENECA and CORSIKA, the effect of the cascade equations can be studied within a single shower, rather than averaged over many showers. Our study shows that for showers produced in this manner, SENECA agrees with CORSIKA to a very high accuracy with respect to densities, energies, and timing information for individual species of ground-level particles from both iron and proton primaries with energies between 1 EeV and 100 EeV. Used properly, SENECA produces ground particle distributions virtually indistinguishable from those of CORSIKA in a fraction of the time. For example, for a shower induced by a 10 EeV proton, SENECA is 10 times faster than CORSIKA, with comparable accuracy.

The color glass condensate approach describes successfully heavy ion collisions at RHIC. We investigate Iron-air collisions within this approach and compare results to event generators commonly used in air shower simulations. We estimate uncertainties in the extrapolation to GZK energies and discuss implications for air shower simulations.

The effects of classical gravitational radiation in models with large extra dimensions are investigated for ultra high energy cosmic rays (CRs). The cross sections are implemented into a simulation package (SENECA) for high energy hadron induced CR air showers. We predict that gravitational radiation from quasi-elastic scattering could be observed at incident CR energies above $10^9$ GeV for a setting with more than two extra dimensions. It is further shown that this gravitational energy loss can alter the energy reconstruction for CR energies $E_{\rm CR}\ge 5\cdot 10^9$ GeV.

The hadronic interaction model BBL implements the ideas of gluon saturation due to large densities. When approaching the black body limit at high energies, leading partons acquire large transverse momenta which breaks up their coherence. This leads to a suppression of forward scattering, and is therefore important for air showers. We discuss some general aspects of this new approach and their influence on air shower properties as seen by fluorescence and surface detectors: The position of the shower maximum is reduced due to stronger absorption in the atmosphere. The lateral distribution functions become flatter for the same reason. Muons are produced abundantly due to high multiplicities in the mid-rapidity region. The response of water Cherenkov detectors and comparisons to other interaction models are shown.

We give the first detailed study of the arrival time distribution of nucleons in UHECR air showers. We analyze in detail the influence of late arriving particles on the energy determination of the AGASA experiment, as well as how the arrival time distribution changes with distance from shower core. Our calculations are consistent with experimental observations of the AGASA group. Crucial to obtaining agreement, is the correct implementation of the energy loss for low-energy protons. We confirm AGASA's estimation of the error in their energy determination associated with late-arriving particles, assuming primary protons.

We discuss that hadron-induced atmospheric air showers from ultra-high energy cosmic rays are sensitive to QCD interactions at very small momentum fractions x where nonlinear effects should become important. The leading partons from the projectile acquire large random transverse momenta as they pass through the strong field of the target nucleus, which breaks up their coherence. This leads to a steeper x_F-distribution of leading hadrons as compared to low energy collisions, which in turn reduces the position of the shower maximum Xmax. We argue that high-energy hadronic interaction models should account for this effect, caused by the approach to the black-body limit, which may shift fits of the composition of the cosmic ray spectrum near the GZK cutoff towards lighter elements. We further show that present data on Xmax(E) exclude that the rapid ~ 1/x^0.3 growth of the saturation boundary (which is compatible with RHIC and HERA data) persists up to GZK cutoff energies. Measurements of pA collisions at LHC could further test the small-x regime and advance our understanding of high density QCD significantly.

We discuss particle production in the high energy limit of QCD. Due to a large gluon density, the interaction reaches the black disk limit and the projectile is resolved into its partonic structure at the saturation scale. This leads to suppression of forward particle production and hereby to a faster absorption of cosmic ray air showers. This property is most suitable for the distinction of evolution scenarios for the saturation scale, e.g. fixed and running coupling BFKL, the latter of which is favored by air shower measurements.

Air shower simulations are essential for interpreting data from cosmic ray experiments. At highest energies though, a microscopic treatment of a whole shower is not possible any more, since it would require a huge amount of CPU-time. We review hybrid approaches of air shower simulation which try to overcome this problem without giving rise to artificial fluctuations as generated by the thinning algorithm.

We discuss particle production in the high-energy, small-x limit of QCD where the gluon density of hadrons is expected to become nonperturbatively large. Strong modifications of the phase-space distribution of produced particles as compared to leading-twist models are predicted which reflect in the properties of cosmic ray induced air showers in the atmosphere. Assuming hadronic primaries, our results suggest a light composition near GZK cutoff energies. We also show that cosmic ray data discriminate among various QCD evolution scenarios for the rate of increase of the gluon density at small x, such as fixed-coupling and running-coupling BFKL evolution. There are clear indications for a slower growth of the gluon density as compared to RHIC and HERA, due e.g. to running-coupling effects.

The influence of high and low energy hadronic models on lateral distribution functions of cosmic ray air showers for AUGER energies is explored. A large variety of presently used high and low energy hadron interaction models are analysed and the resulting lateral distribution functions are compared. We show that the slope as well as the signal at 1000 m distance from the shower axis depend on both the high and low energy hadronic model used. The models are confronted with available hadron-nucleus data from accelerator experiments.

A new hybrid approach to air shower simulations is described. At highest energies, each particle is followed individually using the traditional Monte Carlo method; this initializes a system of cascade equations which are applicable for energies such that the shower is one-dimensional. The cascade equations are solved numerically down to energies at which lateral spreading becomes significant, then their output serves as a source function for a 3-dimensional Monte Carlo simulation of the final stage of the shower. This simulation procedure reproduces the natural fluctuations in the initial stages of the shower, gives accurate lateral distribution functions, and provides detailed information about all low energy particles on an event-by-event basis. It is quite efficient in computation time.

The SENECA model, a new hybrid approach to air shower simulations, is presented. It combines the use of efficient cascade equations in the energy range where a shower can be treated as one-dimensional, with a traditional Monte Carlo method which traces individual particles. This allows one to reproduce natural fluctuations of individual showers as well as the lateral spread of low energy particles. The model is quite efficient in computation time. As an application of the new approach, the influence of the low energy hadronic models on shower properties for AUGER energies is studied. We conclude that these models have a significant impact on the tails of lateral distribution functions, and deserve therefore more attention.

Practically all serious calculations of exclusive particle production in ultra-relativistic nuclear or hadronic interactions are performed in the framework of Gribov-Regge theory or the eikonalized parton model scheme. It is the purpose of this paper to point out serious inconsistencies in the above-mentioned approaches. We will demonstrate that requiring theoretical self-consistency reduces the freedom in modeling high energy nuclear scattering enormously. We will introduce a fully self-consistent formulation of the multiple-scattering scheme in the framework of a Gribov-Regge type effective theory. In addition, we develop new computational techniques which allow for the first time a satisfactory solution of the problem in the sense that calculation s of observable quantities can be done strictly within a self-consistent formalism.

The purpose of this paper is twofold: first we want to introduce a new type of hadronic interaction model (NEXUS), which has a much more solid theoretical basis as, for example, presently used models like QGSJET and VENUS, and ensures therefore a much more reliable extrapolation towards high energies. Secondly, we want to promote an extensive air shower (EAS) calculation scheme, based on cascade equations rather than explicit Monte Carlo simulations, which is very accurate in calculations of main EAS characteristics and extremely fast concerning computing time. We employ the NEXUS model to provide the necessary data on particle production in hadron-air collisions and present the average EAS characteristics for energies 10^14 - 10^17 eV. The experimental data of the casa-blanka group are analyzed in the framework of the new model.