A MC approach to simulate up- and down-going neutrino showers including local topographic conditions

The extragalactic flux of protons is predicted to be suppressed above the famous Greisen–Zatsepin–Kuzmin cut-off at about EGZK ≃ 50 EeV due to the resonant photo-pion production with the cosmic microwave background. Current cosmic ray data do not give a conclusive confirmation of the GZK cut-off and the quest about the origin and the chemical composition of the highest energy cosmic rays is still open. Amongst other particles neutrinos are expected to add to the composition of the cosmic radiation at highest energies. We present an approach to simulate neutrino induced air showers by a full Monte Carlo simulation chain. Starting with neutrinos at the top of the atmosphere, the performed simulations take into account the details of the neutrino propagation inside the Earth and atmosphere as well as inside the surrounding mountains. The products of the interactions are input for air shower simulations. The mountains are modelled by means of a digital elevation map. To exemplify the potential and features of the developed tools we study the possibility to detect neutrino induced extensive air showers with the fluorescence detector set-up of the Pierre Auger Observatory. Both, down-going neutrinos and up-going neutrinos are simulated and their rates are determined. To evaluate the sensitivity, as a function of the incoming direction, the aperture, the acceptance and the total observable event rates are calculated for the Waxman–Bahcall (WB) bound.