Monte-Carlo simulations of atmospheric muon production: Implication of the past martian environment

In order to evaluate the obliquity-driven atmospheric-density path length effect on nuclide production rate on Mars, we performed a Monte-Carlo simulation to produce the number of secondary particles such as muons, neutrons and protons in the martian atmosphere and to simulate that production of ¹⁰Be and ³⁶Cl in the martian regolith by muons and neutrons depends on how much atmosphere had been present for the past 10 million years. The vertical profile of the present martian atmosphere to generate secondary particles has been determined based on the data provided by the Viking missions. For other thickness profiles, we scaled Linsley's atmospheric model. Atmospheric shower has been generated with the SIBYLL 2.1 for high-energy hadronic interactions and EHSA for low energy photonuclear interactions. With increasing atmospheric thickness, more primary interactions occur in the atmosphere. Consequently the proton flux is reduced and the secondary cosmic ray flux, such as muons or energetic neutrons increases at surface. The result indicates that the muon production is more sensitive to obliquity-driven atmospheric variations than proton reduction. A thicker atmosphere would result in enhanced nuclide production at a place deeper than 5 m below the surface and the nuclides present in detectable concentrations. Application to the polar deposit is described.