Issue 
J. Space Weather Space Clim.
Volume 9, 2019
Planetary Space Weather



Article Number  A2  
Number of page(s)  19  
DOI  https://doi.org/10.1051/swsc/2018051  
Published online  14 January 2019 
Research Article
Implementation and validation of the GEANT4/AtRIS code to model the radiation environment at Mars
Institute of Experimental and Applied Physics (IEAP), University of Kiel, Leibnitzstr. 11, 24118 Kiel, Germany
^{*} Corresponding authors: astro.jingnan@gmail.com; banjac@physik.unikiel.de
Received:
9
October
2018
Accepted:
22
December
2018
A new GEANT4 particle transport model – the Atmospheric Radiation Interaction Simulator (AtRIS, Banjac et al. 2018. J Geophys Res Space Phys 123. https://doi.org/10.1029/2018JA026042) – has been recently developed in order to model the interaction of radiation with planets. The upcoming instrumentational advancements in the exoplanetary science, in particular transit spectroscopy capabilities of missions like JWST and EELT, have motivated the development of a particle transport code with a focus on providing the necessary flexibility in planet specification (atmosphere and soil geometry and composition, tidal locking, oceans, clouds, etc.) for the modeling of radiation environment for exoplanets. Since there are no factors limiting the applicability of AtRIS to Mars and Venus, AtRIS’ unique flexibility opens possibilities for new studies.
Following the successful validation against Earth measurements (Banjac et al. 2018. J Geophys Res Space Phys 123. https://doi.org/10.1029/2018JA026042), this work applies AtRIS with a specific implementation of the Martian atmospheric and regolith structure to model the radiation environment at Mars. We benchmark these first modeling results based on different GEANT4 physics lists with the energetic particle spectra recently measured by the Radiation Assessment Detector (RAD) on the surface of Mars. The good agreement between AtRIS and the actual measurement provides one of the first and sound validations of AtRIS and the preferred physics list which could be recommended for predicting the radiation field of other conceivable (exo)planets with an atmospheric environment similar to Mars.
Key words: Particle radiation in space / Particle transport model / Martian exploration / Planetary space weather
© J. Guo et al., Published by EDP Sciences 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
There are mainly two types of energetic particles in the heliosphere that may impose risks for deep space and planetary missions: Galactic Cosmic Rays (GCRs) and Solar Energetic Particles (SEPs). GCRs are energetic charged particles, comprised of 2% electrons and 98% atomic nuclei with the later contributed by 87% protons, 12% helium, and about 1% heavier nuclei (Z ≥ 3) (Simpson, 1983). They have energies from less than 1 MeV nuc^{−1} up to hundreds of TeVs with a powerlaw distribution at energies above ~GeV nuc^{−1} (Allkofer, 1975). SEPs on the other hand are energetic particles (mainly protons and electrons) emitted from the Sun and accelerated by solar flares and/or Coronal Mass Ejection (CME) associated shocks. SEPs generally have energies from a few keV up to several hundreds MeV (occasionally even reaching 1 or 2 GeV or further above) and can reach significantly higher fluxes at these energies compared to background GCRs.
It is likely that the first humanvisited planet will be our neighbour planet Mars. It has a thin atmosphere with its surface pressure less than 1% of that at Earth’s surface making it much easier for high energy particles to reach the Martian surface. Therefore, the assessment of the Martian radiation environment is necessary and fundamental for (a) mitigating radiation risks for nearfuture robotic and crewed missions and (b) better understanding the impact of energetic particles on the preservation of organic biosignatures on Mars. Primary GCRs and SEPs passing through the Martian atmosphere may undergo inelastic interactions with the ambient atomic nuclei losing their energies and also creating secondary particles via spallation and fragmentation processes. These secondary particles may further interact with the atmosphere as they propagate downwards and even with the Martian regolith, finally resulting in very complex spectra including both primaries and secondaries at the surface of Mars (e.g., Saganti et al., 2002).
There are various particle transport codes such as High charge (Z) and Energy TRaNsport (HZETRN, Slaba et al., 2016; Wilson et al., 2016), Particle and Heavy Ion Transport code System (PHITS, Sato et al., 2013) and GEometry And Tracking (GEANT4, Agostinelli et al., 2003; Allison et al., 2016) which can be employed for studying the particle spectra and radiation dose at Mars. Various studies have also combined these particle transport codes with different GCR and/or SEP spectra for estimating the radiation exposure on the surface of Mars (e.g., Simonsen et al., 1990; Simonsen & Nealy, 1993; Saganti et al., 2004; Keating et al., 2005; De Angelis et al., 2006; Ehresmann et al., 2011; McKennaLawlor et al., 2012; Guo et al., 2017a). Articles collected in a special issue at Life Sciences in Space Research (Hassler et al., 2017) have included most recent studies to model the radiation environment at Mars which is also compared with in situ measurement at the surface of Mars (Ehresmann et al., 2017; Guo et al., 2017b) by the Radiation Assessment Detector (RAD, Hassler et al., 2012) on board the Curiosity rover belonging to the Mars Science Laboratory (MSL, Grotzinger et al., 2012).
In particular, PLANETOCOSMICS (http://cosray.unibe.ch/~laurent/planetocosmics/) is a toolkit based on GEANT4 with a specific application purpose to simulate particle transport in planetary magnetic fields as well as interactions when passing through planetary environments and creating secondary particles (Desorgher, 2005). Modeling the radiation environment on the surface of Mars using PLANETOCOSMICS has been carried out by previous researchers (e.g., Dartnell et al., 2007; Ehresmann et al., 2011; Gronoff et al., 2015; Matthiä et al., 2016; Guo et al., 2018c) and has been validated by Matthiä et al. (2016) when compared to energetic charged and neutral particle spectra on the surface of Mars measured by MSL/RAD.
Since the landing of MSL at the Gale Crater of Mars on August 6, 2012, RAD has been providing the first in situ detection of the radiation environment at the surface of Mars (Hassler et al., 2014). These novel measurements provide evaluations of the Martian radiation level at an atmospheric depth around 22 g cm^{−2}. It has been found that the measured GCR radiation dose rate is anticorrelated with the surface pressure (proportional to the atmospheric column depth) which changes both daily and seasonally (Rafkin et al., 2014; Guo et al., 2015, 2017a) up to ±25%. Meanwhile the varying heliospheric conditions modulate the GCR fluxes and thus RAD measured radiation is anticorrelated with solar activity both in the long term evolution (Guo et al., 2015) and in the short term due to e.g., interplanetary CMEs and their associated shocks passing Mars (Witasse et al., 2017; von Forstner et al., 2018; Winslow et al., 2018; Guo et al., 2018b). RAD also measures the flux spectra of energetic charged particles, such as protons, deuterons, tritons, ^{3}He and ^{4}He, carbon, nitrogen, oxygen and iron ions which penetrate downwards and stop within the detector set with energies up to about 100 MeV nuc^{−1} (Ehresmann et al., 2014, 2017). With an inversion method exploiting the response matrix of the neutral particle detection efficiency in the scintillators, RAD can also measure spectra of neutral particles on Mars from ~10 to 900 MeV (Kohler et al., 2014; Guo et al., 2017b).
The last solar cycle (no. 24) has been rather quiet and RAD has only detected a few SEP events at the Martian surface and most of them are rather insignificant apart from the September 10, 2017 event (Ehresmann et al., 2018; Zeitlin et al., 2018; Guo et al., 2018a). This was the first Ground Level Enhancement (GLE) resulting from a very intense SEP detected at the surface of two different planets: Earth and Mars. In general, SEP events are sporadic, often rather impulsive and could be extremely hazardous especially when the observer has a direct magnetic connection to the particle acceleration region and injection site at the Sun. In the case of Mars, there is no global magnetosphere that could shield the atmosphere from energetic particles. In addition, its thin atmosphere can only stop charged particles below ~150 MeV nuc^{−1} (Guo et al., 2018c). Therefore it is important to reliably model SEP events and their induced radiation in order to give immediate and precise alerts for future human missions to Mars. To do so, it is essential to well understand the interactions of atmospheric molecules with incoming particles.
Recently, a new GEANT4 particle transport model – the Atmospheric Radiation Interaction Simulator (AtRIS, Sect. 2.2) – has been developed by Banjac et al. (2018) in order to model the interaction of radiation with planets. The upcoming instrumentational advancements in the exoplanetary science, in particular transit spectroscopy capabilities of missions like JWST and EELT, have motivated the development of this particle transport code with a focus on providing the necessary flexibility in planet specification (atmosphere and soil geometry and composition, tidal locking, oceans, clouds, etc.) for the modeling of exoplanets. The application of AtRIS to model the radiation environment on Mars has been realized and validated for the first time in this study.
2 Model description and implementation
2.1 The Mars Climate Database (MCD)
The Mars Climate Database (MCD, http://wwwmars.lmd.jussieu.fr) offers the possibility to access Martian atmospheric properties, such as temperature, density and composition, for different altitudes, seasons and even the time of the day on Mars. MCD has been developed using different Martian atmospheric circulation models which are further compared and modified by the observation results from past and current Mars missions (Lewis et al., 1999). Therefore, it provides a Martian atmospheric environment which can be implemented into the planetary particle transport toolkit for the purpose of simulating high energetic particles interacting with the Martian atmospheric atoms.
In our atmospheric setup for the AtRIS particle transport model, we use the composition, density and temperature profiles from MCD at Gale Crater on Mars where MSL’s rover Curiosity landed on August 6, 2012 (coordinate: 4.5° S, 137.4° E). The elemental composition of the Martian atmosphere consists of C, O, N, Ar, and H with more than 95% of the molecules being CO_{2}. The atmospheric condition is set to be “clim aveEUV” for climatology scenario with average solar EUV radiation. The Martian solar longitude is set to be 200.5° when the surface pressure is close to the annual average pressure at Gale Crater measured by MSL ~ 840 Pa (e.g., Guo et al., 2015). Figure 1 shows the atmospheric structure, i.e., density and pressure versus the atmospheric altitude used for the AtRIS simulations.
Fig. 1 Martian atmospheric structure at Gale Crater implemented in the AtRIS simulations as imported from the MCD model. 
2.2 Model description: AtRIS
GEANT4 is a Monte Carlo approach widely used for simulating the interactions of particles as they traverse matter (Agostinelli et al., 2003).
The Atmospheric Radiation Interaction Simulator, AtRIS, is a GEANT4 based particle transport code developed to simulate the propagation of energetic particles through planetary atmosphere and regolith. AtRIS allows rather flexible geometry and composition definitions of the planet. This flexibility and ease of use has been achieved through a custom interface called the planet specification format, that has been optimized for the specification of planets and their atmospheres. For Earth, an interface based on the socalled NRLMSISE00 (Picone et al., 2002) model is provided within AtRIS. Similarly, for Mars, the MCD interface has been implemented (see Sect. 2.1).
AtRIS can calculate ion and electron pair production rates, secondary particle distributions (as a function of energy, directionality, planet altitude and so on), as well as absorbed and equivalent dose rates for a 30 cm diameter ICRU sphere phantom composed of water. The tracking of charged particles through magnetic fields is not implemented in AtRIS. As Mars lacks a global magnetic field, it serves as a good validation object for particle propagation modeling of AtRIS. Control over hadronic and electromagnetic processes is provided via the standard GEANT4 messenger and compounded physics list naming scheme (Allison et al., 2016; Geant4_Collaboration, 2017). Different physics lists tested and compared in this study will be discussed in detail in Section 3.1.
The main features of AtRIS are (i) the planet specification format as explained above, (ii) the Atmospheric Response Matrices (ARMs), quantifying the relation between primary energy and altitudedependent ionization as well as equivalent/absorbed dose rate, and finally, (iii) the spectrum folding procedure used to calculate net quantities like the electron–ion pair production rate, by implementing a convolution of a measured spectrum and the ARM. A more detailed description is given in Banjac et al. (2018, b), where the results of AtRIS were compared with and validated against different kinds of Earth measurements (i.e., ion pair production and secondary particle fluxes, absorbed dose rate and dose equivalent).
In the current AtRIS setup for the Martian environment as shown in Table 1, we use a sphere with a radius of 3390 km (approximately the average radius of Mars) representing Mars. The soil composition is approximated as 50%Si, 40%O, and 10%Fe (mass fractions) and the crust (soil) sheet is 100 m thick. The maximum height of the atmosphere from the MCD model is 100 km which is divided into 500 m thick layers. The accumulated column depth at the surface is about 22 g cm^{−2} corresponding to a surface pressure of 830 Pa (Fig. 1). Table 1 provides the main features of the atmospheric and regolith structures used for the AtRIS simulations.
Martian atmospheric and regolith properties used in AtRIS and four different physics lists of GEANT4 (more details in the text) tested in this study.
2.3 The matrix realization of AtRIS
Recently, Guo et al. (2018c) have developed a generalized approach based on the GEANT4/PLANETOCOSMICS transport code and the MCD Martian environment setups to quickly model the Martian surface radiation level of any given incoming proton/helium ion spectra. Such an approach called “Planetomatrix” can be mathematically described by ARMs and visually illustrated as 2D histograms (Figs. 1 and 2 of Guo et al., 2018c).
Similar to this Planetomatrix approach, we construct the “response function” of the Martian atmosphere based on simulated particle spectra from AtRIS. We can describe the statistical transformation of the atomic and nuclear interaction process for a particle spectrum (of particle type i as a function of energy E_{0}) above the Martian atmosphere resulting in a particle spectrum (of particle type j as a function of energy E) on the Martian surface in a matrix _{ij}(E_{0}, E). As particles could also interact with the Martian atmosphere and regolith and produce albedo particles contributing to the upward fluxes, we have also considered the generation of such upward particles in our simulations. Since the energy spectra of upward and downwardtraveling particles are dissimilar, we have separately constructed the upward and downward directed matrices for each primary–secondary case. Once such a matrix is constructed, it can be folded with any incoming SEP or GCR spectra within the energy range of the simulated particles for calculating the surface spectra without rerunning the particle transport simulations. The upward and downward secondary particle spectra generated by different types of primary particles can be combined. The following steps are used for constructing the matrix and folding it with a given incoming spectra to obtain the surface spectra:
Simulate primary particles (type i) with energy in the range of 1 and 10^{5} MeV through the Mars AtRIS model described in Section 2.2. A flat spectrum over 50 energy bins distributed uniformly in logarithmic scale is used. There are N_{i}(E_{0}) particles simulated in each energy bin.
Based on the simulation results create a matrix M_{ij}(E_{0}, E_{j}) for each secondary particle type j with certain directions, e.g., downwarddirected protons on the surface of Mars generated by primary protons. Each column of the matrix is a histogram of secondary particles H_{i}(E_{j}) representing the number of particles at energy E_{j} created by N_{i}(E_{0}) primary particles with energy E_{0}.
Divide each histogram H_{i}(E_{j}) in M_{ij}(E_{0}, E_{j}) by the number of simulated particles located in each incoming energy bin N_{i}(E_{0}) to generate the normalized histogram due to a single primary particle with energy E_{0}. The normalized histograms constitute the normalized matrix which is .
For a given GCR/SEP spectrum f_{i} which is often an isotropic flux in deep space (e.g., with units of particles/sec/m^{2}/sr/MeV), first interpolate this spectrum using the incoming energies (E_{0}) of the matrix M_{ij}(E_{0}, E_{j}) so that f_{i} = f_{i}(E_{0}). Then multiply it with W(E_{0}) which is the incoming bin width (MeV) of the matrix such that F_{i}(E_{0}) = W(E_{0})f_{i}(E_{0}) with the unit of particles/sec/m^{2}/sr.
Calculate the total primary particle count rate (particles/sec) taking into account of the integrated geometric factor of particles arriving at the planet by multiplying the above F_{i}(E_{0}) with Ω_{in}A_{in}. Here A_{in} is the area of which primary particles were fed into the simulation. When using a sphere covering the top of the planetary atmosphere as the source sphere, A_{in} is 4πR^{2} with R (meter) being the radius of this sphere (i.e., R_{top} in Table 1). Ω_{in} is the integration of the cosine of zenith angle of the source particles and is π when they were injected isotropically inward from the source sphere.
Fold such scaled incoming particle count rate (energy dependent) with the normalized matrix to obtain the histogram of secondaries (particles/sec) corresponding to each column of the matrix (at certain incoming primary energy).
Scale the secondary particle histogram at each column into a differential spectra per geometric factor by dividing it with Ω_{out}A_{out}W(E_{j}). Here W(E_{j}) is the energy bin width (MeV) of the secondary particle histogram. A_{out} is the surface area of the planet where the output secondaries are counted and it is 4πR^{2}_{surf} with R_{surf} (see Table 1) in centimeter. And the Ω_{out} is the integration of the cosine of zenith angle of secondaries crossing the surface plane.
In this study, we consider three cases of secondary particle directions from AtRIS simulations: (a) all secondaries propagating downwards with the fluxweighted solid angle Ω_{out} of π, (b) all secondaries propagating upwards with Ω_{out} = π and (c) secondary downward directed particles with a zenith angle within 36°, i.e., Q_{out} = 1.1 sr. This corresponds to the inner cone angle used for counting downward propagating charged particles stopping inside the RAD detector (Hassler et al., 2012).
Finally the surface secondary spectra F_{ij}(E_{j}) induced by primary particle spectrum F_{i}(E_{0}) obtained through step 1–7 has unit of particles/sec/cm^{2}/sr/MeV. The above procedure can be also mathematically summarized as:(1)where is the atmospheric response matrix which takes into account the scaling of the energy bin widths of the output histograms and the geometric factors used in the simulation. The surface spectrum of particle type j resulting from primary particle type i is the multiplication of the matrix with the primary particle spectrum and this operation is essentially the sum product over the second axis (E_{0}) of and F_{i}, i.e.,(2)
The simulation is often set up with the areas A_{in} and A_{out} equal to each other where the planet is approximated as a rectangular box which has the same size on top of the atmosphere and on the regolith surface. In the current AtRIS/MCD setup, A_{in} is only slightly bigger than A_{out} as R_{top} is the radius of the planet including the atmospheric layer and it is slightly larger than R_{surf}. Given three cases of solid angle of secondaries detected on the surface, i.e., downward, upward and within 36° zenith angle (RAD inner cone), a matrix for each case can be constructed for certain secondary particles induced by given primary particle type.
Although the construction of each matrix is timeconsuming, the multiplication of different input spectra F_{i}(E_{0}) with such a matrix to generate different surface spectra F_{j}(E_{j}) is very much simplified. The spectra of the surface secondary particle are the combined spectra resulting from different types of primary particles which arrive at the top the Martian atmosphere:(3)(4)(5)here , and represent the energy spectra of secondary particle type j, resulting from various primary particle types, averaged in the downward, upward and within 36° zenith angle respectively. As a first step of the model verification, we have simulated protons and ^{4}He ions as primary particles because they constitute the majority of GCR particles. We considered the secondary charged particles (type j) whose spectra are also measured by RAD on the surface of Mars, such as protons, deuterons, tritons, ^{4}He and ^{3}He ions as output particles.
3 Validation of AtRIS
GEANT4 (version 10.4.p02 used here) offers a wide variety of models for handling physical processes within different energy ranges. It is not yet entirely clear what is the most accurate and efficient physics list describing high energy (from hundreds of MeV to tens of GeV) particles integrating with the Martian atmosphere (Matthiä et al., 2016). One of the goals of MSL/RAD is to help validate the appropriate physics list which could precisely model the high energetic cosmic ray interaction with the Mars atmosphere (Hassler et al., 2012).
In this study, we employ four different physics lists (Table 1): QGSP_BIC_HP (model A), QGSP_BERT_HP (model B), FTFP_BERT_HP (model C) and also FTFP_INCLXX_HP (model D) when applying AtRIS to model the Martian radiation environment. These physics lists are named based on the combination of various models in different energy ranges (Geant4_Collaboration, 2017):
QGS stands for the Quark Gluon String model for high energy particles (>~20 GeV).
FTF represents the Fritiof model for particles >~5 GeV.
P stands for the Precompound model used for deexcitation process for nucleoninduced reactions below 1–2 MeV.
BIC uses the Binary Cascade Model for particles <~10 GeV.
BERT uses the Bertini Cascade Model in the middle energy range of <~10 GeV.
INCL represents the Liege Intranuclear Cascade model in the middle energy range.
HP option switches on the high precision neutron elastic and inelastic scattering model for neutrons below 20 MeV.
At energies where different physics models are overlapping, a weighted combination of various models is considered.
Specifically speaking, model A) QGSP_BIC_HP and model B) QGSP_BERT_HP use the same Quark Gluon String (QGS) model for the high energy range while different cascade models for the lower energy range (<~10 GeV). FTFP_BERT_HP in model C is recommended by Geant4 collaboration “for cosmic ray applications where good treatment of very high energy particles is required” (Geant4_Collaboration, 2017). It contains all standard electromagnetic (EM) processes by default. It uses Bertinistyle cascade for hadrons <5 GeV and the FTF (Andersson et al., 1987; NilssonAlmqvist & Stenlund, 1987) model for simulating the interaction of mesons, nucleons and hyperons in the 3 GeV–100 TeV energy range. Finally and in model D, the Liege Intranuclear Cascade model INCL++ has been recently extended in GEANT4 to handle reactions between 3 and 15 GeV incident energies (Mancusi et al., 2014) and extensive benchmarks have shown the INCL++ model has a very good predictive power for the particles related to neutron production in spallation reactions (Leray et al., 2011). A map of physics models in different energy ranges used for the INCLXX physics list can be found here http://irfu.cea.fr/dphn/Spallation/physlist.html.
We have validated the accuracy of the AtRIS approach when applied with different physics lists (named model A, B, C and D) to the Martian environment, especially for the generation of charged particle spectra by primary protons and ^{4}He ions. We first compare the ARMs of each primary–secondary pair when using different physics lists as shown in Section 3.1. Then we compare such modeled surface secondary spectra based on different physics lists with the MSL/RAD measurements as explained in Section 3.2.
3.1 Matrices constructed via AtRIS
Using the matrix approach described in Section 2.3, we compare the ARMs derived from above four different models: A) QGSP_BIC_HP, B) QGSP_BERT_HP, C) FTFP_BERT_HP and D) FTFP_INCLXX_HP. To ease the comparison between different simulation results, normalized matrices (step 3 in Sect. 2.3) have been plotted and discussed.
The interaction of primary protons with the Martian atmosphere has been modeled and ARMs have been constructed for five different secondary types including Hydrogen isotopes (proton, deuteron and triton) and Helium isotopes (^{3}He and ^{4}He). Each secondary type is differentiated as downward and upwarddirected with the respective matrix. A total of ten matrices have been generated for primary proton interactions. Similarly ten other matrices represent the primary ^{4}He ion interaction with the Martian atmosphere. Although there are other secondaries such as electrons, muons, positrons which also contribute to the radiation environment on Mars, we focus on the five hydrogen and helium isotopes as their spectra have been well measured on the surface by MSL/RAD (Ehresmann et al., 2014) and a direct modelobservation comparison is possible.
3.1.1 Primary proton matrices
Figure 2 shows the matrices which describe the surfaced downward secondary protons which result from the primary protons under an atmospheric configuration as described in Section 2.1. The general statistics of the matrices constructed from four different physics lists are rather similar as shown in the top row of the figure. We have also accounted for the Poisson uncertainty in the Monte Carlo simulations in each bin of the matrix and scaled it with the statistics in the corresponding bin to obtain the normalized uncertainties as plotted in the second row of the figure. As shown, the normalized uncertainties are relatively small (below 20%) in most bins apart from the low energy output bins where the bin widths are smaller (energies are binned into logarithmic scale) with lower statistics.
Fig. 2 Matrices which describe the surface downward protons created by the primary protons. First row: the normalized Martian atmospheric matrices of primary protons generating surface downward directed protons from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 
To better quantify the differences between models, we have calculated the (relative) image difference between matrices of different models as shown in the last two rows of Figure 2, which we define as Image Differencing Matrices (IDM). It is visible that there are some systematic differences between different models and boundaries of patterns are formed at energies where physics processes are switching in the model. For instance, A–B IDM reflects the difference between BIC and BERT models for incoming particles below ~10 GeV and the consistency between two physics lists at higher energies where QGSP is used for both models. Alternatively, B–C IDM shows the difference between QGSP and FTFP at high energies and the agreement of BERT at lower energies. In comparison, C–D shows a good agreement between FTFP_BERT_HP and FTFP JNCLXX_HP apart from low output energy ranges where BERT seems to have a higher efficiency. This is because when simulating p–p interactions at energies below 20 MeV, FTFPJNCLXX_HP uses the models from the HP library. INCL is turned on only when one of the involved particles has the atomic number larger than one.
As shown in the bottom row, the mean relative differences are about 20.9% between A and B, 23.2% between B and C and 17.7% between C and D which are not much larger in comparison to the mean relative uncertainties of each model (second row). Therefore the relative IDM indicates a general good agreement between different physics lists for protons generating secondary protons on the surface of Mars. Besides, a 2D Gaussian filter has also been applied to the IDM of models C and D, shown in the last column of the last two rows, in order to smooth out features resulting from fine bins with low statistics.
Figures 3 and 4 show the matrices of primary protons generating secondary deuterium and tritium hydrogen isotopes detected as downward particles on the surface of Mars. Compared to the generation of secondary protons, the p–d and p–t reaction probability is much lower as can be seen in the matrices of the top panels which have the same color bar scale as the proton–proton matrices. The deuterium and tritium matrices lack the high energy diagonal component, while most secondaries are due to spallation reactions of high energy protons (above a few GeV) with the atmosphere nucleus. Model D is clearly more efficient in generating deuterium and tritium particles via such reactions as is readily seen in the model D matrices. Model A in the range with the binary cascade model is least efficient in such productions. Indicated by the IDM between different models, in the energy range below ~10 GeV, deuteron and triton generation efficiency of different models follow: BIC < BERT < INCL; at higher energies, FTFP has a slightly higher efficiency than QGSP.
Fig. 3 Matrices which describe the surface downward deuterons created by the primary protons. First row: the normalized Martian atmospheric matrices of primary protons generating surface downward directed ^{2}H ions from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 
Fig. 4 Matrices which describe the surface downward tritons created by the primary protons. First row: the normalized Martian atmospheric matrices of primary protons generating surface downward directed ^{3}H ions from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 
Figure 5 shows the matrices of primary protons resulting into surface downward ^{4}He ions. These reactions are even more rare than the above proton–deuteron and proton–triton cases. Despite of the low statistics of such reactions, the IDM between A and B indicates a slightly higher efficiency of BERT and the IDM between B and C reflects a higher production rate in FTFP compared to QGSP. And it is also noticeable from the IDM between C and D that INCL predicts a higher probability of such spallation interactions in the relevant energy range.
Fig. 5 Matrices which describe the surface downward ^{4}He particles created by the primary protons. First row: the normalized Martian atmospheric matrices of primary protons generating surface downward directed ^{4}He ions from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 
Finally, Figure 6 shows surface downward ^{3}He ions induced by primary protons. The features of matrices in this case are very similar to those of deuterium, tritium and ^{4}He ions in Figures 3–5 indicating similar cascading processes in generating these different types of secondaries. However ^{3}He ions have a slightly lower statistics than ^{4}He ions which may suggest some of the ^{3}He ions are from fragmented ^{4}He ions first generated in the atmosphere.
Fig. 6 Matrices which describe the surface downward ^{3}He particles created by the primary protons. First row: the normalized Martian atmospheric matrices of primary protons generating surface downward directed ^{3}He ions from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and a is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 
Primary protons and their secondaries generated in the atmosphere may reach the Martian surface and generate upward particles in the regolith. As an example, we only show the matrices constructed for the case of upward protons as shown in Figure 7. Similar to the downward proton case in Figure 2 the IDM of A–B and B–C show higher efficiency of BERT in the middle energy range and FTFP in the high energy part. The IDM of C–D shows features that are more similar to those of downward protons rather than those of deuterium, tritium and ^{3}He particles, indicating that these are mostly downward propagating protons reaching the soil and get redirected upwards rather than protons produced via spallation in the regolith. This is also suggested by the similar atmospheric cutoff energies of primary protons when comparing the downward and upward secondary proton matrices.
Fig. 7 Matrices which describe the surface upward protons created by the primary protons. First row: the normalized Martian atmospheric matrices of primary protons generating surface upward directed protons from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 
The atmospheric cutoff energy is about 160 and 200 MeV in Figures 2 and 7 respectively which approximately corresponds to the minimum energy a proton needs to traverse through ~22 g cm^{−2} of atmosphere.
3.1.2 Primary ^{4}He ion matrices
Similar to the ARMs of primary protons, the primary ^{4}He ion interaction with the Martian atmosphere and the generation of secondaries have been mapped to matrices with some examples shown in Figures 8–11.
Fig. 8 Matrices which describe the surface downward protons created by the primary ^{4}He particles. First row: the normalized Martian atmospheric matrices of primary ^{4}He generating surface downward directed protons from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins 
Figure 8 displays the matrices for ^{4}He generated surface protons. The general shape from four physics lists looks rather similar and also very comparable to the proton–proton matrices in Figure 2, indicating that a substantial amount of ^{4}He ions fragment into protons in the atmosphere and follow the same energy loss processes as in the proton–proton case. It is noticeable that model D predicts a slightly higher efficiency in the primary energy range from a few GeV to ~20 GeV corresponding to the INCL model in the physics list of model D.
Figure 9 shows the ^{4}He–^{4}He matrices from four physics lists which have similar structures. It is visible that there are three components in these matrices. The first is represented by the topright diagonal component showing some ^{4}He particles with energies ≥~2–3 GeV arriving at the Martian surface as ^{4}He ions with similar energies, i.e., these relativistic particles propagate through the 22 g cm^{−2} of atmosphere without any interactions. The second component is shown as a vertical branch connecting to the lowerleft part of the diagonal structure. These are primary ions with energies between ~650 MeV and ~3–4 GeV which traverse through the atmosphere with substantial ionization energy loss. The third component is the blob structure at the lowerright part of the matrices (similar to the pd, pt matrices in Figs. 3 and 4) which represent high energy primary particles interacting with the atmosphere and generating secondary ^{4}He ions via spallation process. There is a gap beneath the diagonal component while there is an enhanced production of ^{4}He produced protons at this range (Fig. 8). This suggests that ^{4}He particles easily fragment into protons as they traverse through the atmosphere.
Fig. 9 Matrices which describe the surface downward ^{4}He created by the primary ^{4}He particles. First row: the normalized Martian atmospheric matrices of primary ^{4}He generating surface downward directed ^{4}He from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 
As shown by the IDM in the lower panel, all four models result in very small differences especially for the first component. The difference between different models is mostly due to statistics as the IDM shows rather similar features compared to the uncertainty matrices and the mean values of IDM are comparable to the mean matrix uncertainties. Model D seems to predict a larger effect of the third component, especially in the energy range where the INCL model is applied.
Figure 10 represents the ^{4}He–deuteron case and also shows the threecomponent structure in four models. Here the first and second components correspond to ^{4}He ions fragmenting into ^{2}H particles as they propagate through the atmosphere. It is also shown by the IDM of A–B that Bertini cascade is more efficient and Binary cascade. The third component predicted by model D is much more enhanced compared to model A and B and the transition between the second and the third components in model D is much smoother. The average difference between model D and C is significant and as large as 100%. Similar features of the matrices are also shown in Figure 11 for the ^{4}He–^{3}He case with model D predicating an enhancement of ^{3}He production especially at the range between the second and third components. This might be due to ^{4}He being more efficient in fragmenting into ^{3}He particles.
Fig. 10 Matrices which describe the surface downward deuterons created by the primary ^{4}He particles. First row: the normalized Martian atmospheric matrices of primary ^{4}He generating surface downward directed ^{2}H from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 
Fig. 11 Matrices which describe the surface downward ^{3}He created by the primary ^{4}He particles. First row: the normalized Martian atmospheric matrices of primary ^{4}He generating surface downward directed ^{3}He from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 
3.2 Comparison of modeled spectra with RAD proton measurements
In order to generate surface secondary particle spectra, we fold the above matrices with the primary proton and helium GCR spectra calculated from a standard GCR model – the BadhwarO’Neill 2010 (BON10, O’Neill, 2010). The BON10 spectra have been computed using the spherically symmetric Fokker–Planck equation, with the Local Interstellar Spectrum (LIS) at the boundary of the heliosphere (~100 AU) as a boundary condition. The model uses an input parameter of the solar modulation Φ, in units of volt, as derived from the International Sunspot Number. The modulation parameter Φ is a depiction of solar activity and its modulation magnitude on the GCR particle fluxes. The modeled spectra at Earth have been fitted and adjusted to measurements obtained by a multitude of instruments, including all available observations of GCR particle energy spectra from 1955 to 2010. Note that we use the nearEarth GCR spectra to approximate those at Mars with a distance of ~1.5 AU to the Sun. The radial gradient of GCRs in the inner heliosphere has been estimated to be about 3% per AU for protons with energies around 1 GeV (e.g., Gieseler & Heber, 2016). This translates into about 1.5% of increase of GCR flux from Earth and Mars. However, this correction is much smaller compared to other systematic and statistical uncertainties in the models employed in the current work and therefore we omit the radial gradient of the heliospheric modulation of GCRs.
The solar modulation parameter used here as input for the BON10 GCR model is equal to the average value, 550 mV, during the periods of the selected RAD measurements (Ehresmann et al., 2017). Primary protons and ^{4}He ions from the GCR model have been folded with the matrices to generate surface secondary particles. The resulting secondary spectra of protons, ^{4}He and ^{3}He ions, deuterons, tritons in the energy range of 10–100 MeV nuc^{−1} in the RAD view cone have been compared with MSL/RAD surface measurements as shown in Figure 12. Although RAD is not always positioned perpendicular to the surface as the rover body is inclined occasionally, WimmerSchweingruber et al. (2015) have found from the measurement that the variability of downward particle fluxes with the zenith angle is small. Therefore a separate set of matrices with surface particles constrained within the RAD view cone (<36° of downward zenith angle) has been constructed. Besides, we have also propagated the Poisson uncertainties from the simulations through the matrix multiplication for obtaining the uncertainties of the modeled surface spectra (Eq. (3)).
Fig. 12 Martian surface spectra of H (panel a) and He (panel b) isotopes induced by primary GCR proton and ^{4}He particles modeled via four different physics lists (Table 1). The spectra are averaged within 36° of downward zenith angle corresponding to the RAD view cone. (a) Input: GCR H + ^{4}He, output: downward H isotopes, (b) input: GCR H + ^{4}He, output: downward He isotopes. 
The proton spectra modeled by four different physics lists all agree very well with each other and also the RAD measurement. The modeled fluxes are slightly lower than the observed spectra which is reasonable as GCR particles heavier than ^{4}He have not been considered in the model. Despite the large uncertainties due to low statistics, it is clearly shown that the fluxes for both deuterium and tritium particles from model A, B and C are considerably lower than the measured spectra while those predicted by model D match well with the data. As explained in Section 3.1 and shown in Figures 3, 4, and 10, model D has a higher efficiency in generating secondary particles via spallation process (or the third component) which are located at the energy range RAD measures.
For Helium isotopes, all four models show a good agreement of ^{4}He ions, although with a slightly higher flux, compared to the RAD data. For ^{3}He case, model D again has a better prediction in comparison to the RAD measurement in the energy range of 10–100 MeV/nuc. The higher production rate of ^{3}He in model D is clearly shown in Figures 6 and 11 as explained in Section 3.1.
We have also compared the modeled spectra in a larger range beyond the energy range of the RAD data as shown in Figure 13(a) and (d) for H and He isotopes respectively. For protons and ^{4}He ions, all models agree with each other at different energy ranges. However for deuterons, tritons and ^{3}He ions, model D shows significant enhancement of the flux compared to other models for particles up to ~500 MeV/nuc. At higher energies the discrepancy between different models is much smaller as also shown in the deuterium, tritium and ^{3}He matrices. The RAD measurement locates within the energy range where predictions from different physics lists differ the most and the current comparison strongly supports the FTFP JNCLXX_HP model for GCR particles transported through the Martian environment.
Fig. 13 Martian surface spectra of H (a, b, c) and He (d, e, f) isotopes induced by primary GCR proton and ^{4}He particles modeled via four different physics lists (Table 1). The spectra are averaged within different solid angles: (a) and (d) for the RAD view cone, (b) and (e) for the downward direction and (c) and (f) for the upward direction. 
The angular dependence of the secondary particles on the surface of Mars is also very important for future human exploration of Mars and understanding radiation effects on Mars. The downward spectra within the RAD view cone shown in Figure 13(a) and those averaged over the whole downwarddirected solid angle shown in panel (b) have very similar levels of flux. This is consistent with the results found by WimmerSchweingruber et al. (2015) from the RAD measurement that within 15° of the Rover tilt angles, the radiation field is mostly isotropic.
However the upward spectra are rather dissimilar from the downward spectra and all upward spectra lack the high energy component (>~GeV) compared to downward spectra. Particle fluxes of protons and ^{4}He ions are particularly lower than in the downward direction where more primary GCR particles contribute to the surface downward radiation. As the upward flux, mostly produced in the Martian regolith, is much lower than the downward flux, this suggests the effectiveness of shielding using the Martian regolith for potential human habitat on Mars.
4 Summary and conclusion
We have applied the novel tool AtRIS (Banjac et al., 2018) with a specific implementation of the Martian atmospheric and regolith structure to model the radiation environment at Mars.
We have validated the accuracy of the AtRIS model when applied with four different physics lists to the Martian environment, especially for the generation of charged particle spectra by primary protons and ^{4}He ions. We first visualized the atmospheric response matrix of each primary–secondary chain via a 2D histogram where the energydependent efficiency of primary particles generating secondaries is nicely shown. We compared these matrices obtained via four different physics lists of GEANT4 by calculating and visualizing their IDMs. The IDM results suggest that in general the Bertini cascade model is more efficient than the Binary cascade model while INCL model is most productive especially for deuterium, tritium and ^{3}He particles. For two primary particle types of protons and ^{4}He ions, four different physics lists agree with each other reasonably well. For other hydrogen and helium isotopes, model D FTFP_INCLXX_HP generates significantly more secondaries due to spallation processes of primary particles in the energy range of ~GeV and 20 GeV where the Liege Intranuclear Cascade model is applied, as shown in Figures 3, 4, 6, 10 and 11.
To benchmark the AtRIS model based on different GEANT4 physics lists when applied to the Martian environment, we compared the AtRIS predicted spectra with the energetic particle spectra recently measured by MSL/RAD on the surface of Mars (Ehresmann et al., 2017). We folded the above ARMs with BON GCR spectra of protons and ^{4}He particles to obtain the secondary particle spectra of Hydrogen and Helium isotopes: protium, deuterium, tritium, ^{3}He and ^{4}He particles. This process has summed up the same secondary type generated by both proton and ^{4}He primaries. But we have ignored the contribution by other heavier GCR particles. This may have caused the slightly smaller flux of proton in all four models in comparison with the MSL/RAD measurement. The comparison shows that all models have an agreed prediction of the proton and helium flux in comparison to the RAD measurement in the energy range of 10–100 MeV/nuc. However, model D has a higher prediction of deuterium, tritium and ^{3}He flux and a better match with the RAD data. This is due to the same reason we have already observed in the ARMs.
In general, the good agreement between AtRIS and the actual measurement extends the validation of AtRIS to nonEarth systems. For the specific case of a very thin atmosphere, we have analyzed how different physics lists may impact the generation of secondary particles through the Martian atmosphere. While we have seen that FTFPJNCLXX_HP provides a better agreement with the MSL/RAD data for the specific cases that were examined, a further study is needed to validate this for heavier primary ions and other secondary particles as well as for planets with thicker atmospheres and different atmospheric compositions (Earthlike and Venuslike). For instance, secondary muons are a major surface radiation component at the largest moon of Saturn–Titan. Heavier ions and neutrons have more enhanced biological effectiveness and should also be carefully examined and included when using AtRIS to predict the equivalent radiation dose.
In order to make AtRIS easily accessible by the community, we are preparing an online documentation including examples for AtRIS in the form of a wiki page. In 2019, AtRIS will be published under a GNU GPL licence.
Acknowledgments
Part of this work was supported by the German Deutsche Forschungs Gemeinschaft, DFG project 282759267. The work is also supported by DLR and DLR’s Space Administration grant numbers 50QM0501, 50QM1201 and 50QM1701 to the Christian Albrechts University, Kiel. We thank Bent Ehresmann for the helpful discussions and support. J.G. and R.F.W.S. acknowledge the International Space Science Institute, which enabled part of the collaborations in this paper through the ISSI International Team 353 Radiation Interactions at Planetary Bodies. The authors and the editor thank Rami Vainio and an anonymous referee for their assistance in evaluating this paper.
References
 Agostinelli S, Allison J, Amako KA, Apostolakis J, Araujo H, et al. 2003. GEANT4 – A simulation toolkit. Nucl Instrum Methods Phys Res Sect A: Accel Spectrom Detect Assoc Equip 506(3): 250–303. [Google Scholar]
 Allison J, Amako K, Apostolakis J, Arce P, Asai M, et al. 2016. Recent developments in Geant4. Nucl Instrum Methods Phys Res Sect A: Accel Spectrom Detect Assoc Equip 835: 186–225. [Google Scholar]
 Allkofer OC. 1975. Introduction to cosmic radiation. NASA STI/Recon Tech Rep A 75. [Google Scholar]
 Andersson B, Gustafson G, NilssonAlmqvist B. 1987. A model for lowpT hadronic reactions with generalizations to hadronnucleus and nucleusnucleus collisions. Nucl Phys B 281(1–2): 289–309. [Google Scholar]
 Banjac S, Herbst K, Heber B. 2018a. The Atmospheric Radiation Interaction Simulator (AtRIS). J Geophys Res Space Phys 123. DOI: 10.1029/2018JA026042. [Google Scholar]
 Dartnell L, Desorgher L, Ward J, Coates A. 2007. Modelling the surface and subsurface Martian radiation environment: Implications for astrobiology. Geophys Res Lett 34: L02207. [CrossRef] [Google Scholar]
 De Angelis G, Wilson J, Clowdsley M, Qualls G, Singleterry R. 2006. Modeling of the Martian environment for radiation analysis. Rad Meas 41(9): 1097–1102. [CrossRef] [Google Scholar]
 Desorgher L. 2005. PLANETOCOSMICS software user manual. Accessible from the GEANT4/PLANETOCOSMICS web page. [Google Scholar]
 Ehresmann B, Burmeister S, WimmerSchweingruber R, Reitz G. 2011. Influence of higher atmospheric pressure on the Martian radiation environment: Implications for possible habitability in the Noachian epoch. J Geophys Res (Space Phys) 116(A15): 10106. [Google Scholar]
 Ehresmann B, Hassler D, Zeitlin C, Guo J, WimmerSchweingruber R, et al. 2018. Energetic particle radiation environment observed by RAD on the surface of Mars during the September 2017 event. Geophys Res Lett 45(11): 5305–5311. [CrossRef] [Google Scholar]
 Ehresmann B, Zeitlin C, Hassler DM, WimmerSchweingruber RF, Bohm E, et al. 2014. Charged particle spectra obtained with the Mars Science Laboratory Radiation Assessment Detector (MSL/RAD) on the surface of Mars. J Geophys Res Planet 119(3): 468–479. [Google Scholar]
 Ehresmann B, Zeitlin CJ, Hassler DM, Matthiä D, Guo J, et al. 2017. The charged particle radiation environment on Mars measured by MSL/RAD from November 15, 2015 to January 15, 2016. Life Sci Space Res 14: 3–11. [CrossRef] [Google Scholar]
 Geant4_Collaboration. 2017. Geant4 physics reference manual 10.4. Accessible from the GEANT4 web page. URL http://geant4userdoc.web.cern.ch/geant4userdoc/UsersGuides/PhysicsListGuide/html/index.html. [Google Scholar]
 Gieseler J, Heber B. 2016. Spatial gradients of GCR protons in the inner heliosphere derived from Ulysses COSPIN/KET and PAMELA measurements. A&A 589: A32. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gronoff G, Norman RB, Mertens CJ. 2015. Computation of cosmic ray ionization and dose at Mars. I: A comparison of HZETRN and planetocosmics for proton and alpha particles. Adv Space Res 55(7): 1799–1805. [CrossRef] [Google Scholar]
 Grotzinger JP, Crisp J, Vasavada AR, Anderson RC, Baker CJ, et al. 2012. Mars Science Laboratory mission and science investigation. Space Sci Rev 170(1–4): 5–56. [NASA ADS] [CrossRef] [Google Scholar]
 Guo J, Dumbovic M, WimmerSchweingruber RF, Temmer M, Lohf H, et al. 2018a. Modeling the evolution and propagation of the 2017 September 9th and 10th CMEs and SEPs arriving at Mars constrained by remotesensing and insitu measurement. Space Weather 16: 1156–1169. [NASA ADS] [CrossRef] [Google Scholar]
 Guo J, Lillis R, WimmerSchweingruber RF, Zeitlin C, Simonson P, et al. 2018b. Measurements of Forbush decreases at Mars: Both by MSL on ground and by MAVEN in orbit. A&A 611: A79. [CrossRef] [EDP Sciences] [Google Scholar]
 Guo J, Slaba TC, Zeitlin C, WimmerSchweingruber RF, Badavi FF, et al. 2017a. Dependence of the Martian radiation environment on atmospheric depth: Modelling and measurement. J Geophys Res Planet Sci 122: 329–341, DOI: 10.1002/2016JE005206. [Google Scholar]
 Guo J, Zeitlin C, WimmerSchweingruber R, Hassler DM, Kohler J, Ehresmann B, Bottcher S, Bohm E, Brinza DE. 2017b. Measurements of the neutral particle spectra on Mars by MSL/RAD from 20151115 to 20160115. Life Sci Space Res 14: 12–17. [CrossRef] [Google Scholar]
 Guo J, Zeitlin C, WimmerSchweingruber RF, McDole T, Kiihl P, Appel JC, Matthiä D, Krauss J, Kohler J. 2018c. A generalized approach to model the spectra and radiation dose rate of solar particle events on the surface of Mars. Astron J 155(1): 49. [NASA ADS] [CrossRef] [Google Scholar]
 Guo J, Zeitlin C, WimmerSchweingruber RF, Rafkin S, Hassler DM, et al. 2015. Modeling the variations of dose rate measured by RAD during the First MSL Martian Year: 2012–2014. Astrophys J 810(1): 24. [NASA ADS] [CrossRef] [Google Scholar]
 Hassler D, Zeitlin C, WimmerSchweingruber R, Bottcher S, Martin C, et al. 2012. The radiation assessment detector (RAD) investigation. Space Sci Rev 170(1–4): 503–558. [NASA ADS] [CrossRef] [Google Scholar]
 Hassler DM, Norbury JW, Reitz G. 2017. Mars Science Laboratory radiation assessment detector (MSL/RAD) modeling workshop proceedings. Life Sci Space Res 14: 1–2. [CrossRef] [Google Scholar]
 Hassler DM, Zeitlin C, WimmerSchweingruber RF, Ehresmann B, Rafkin S, et al. 2014. Mars surface radiation environment measured with the Mars Science Laboratory’s curiosity rover. Science 343(6169): 1244797. [CrossRef] [Google Scholar]
 Keating A, Mohammadzadeh A, Nieminen P, Maia D, Coutinho S, Evans H, Pimenta M, Huot JP, Daly E. 2005. A model for Mars radiation environment characterization. IEEE Trans Nucl Sci 52(6): 2287–2293. [Google Scholar]
 Kohler J, Zeitlin C, Ehresmann B, WimmerSchweingruber R, Hassler D, et al. 2014. Measurements of the neutron spectrum on the Martian surface with MSL/RAD. J Geophys Res Planet 119(3): 594–603. [Google Scholar]
 Leray S, David J, Khandaker M, Mank G, Mengoni A, Otsuka N, Filges D, Gallmeier F, Konobeyev A, Michel R. 2011. Results from the IAEA benchmark of spallation models. J Korean Phys Soc 59(2): 791–796. [CrossRef] [Google Scholar]
 Lewis SR, Collins M, Read PL, Forget F, Hourdin F, Fournier R, Hourdin C, Talagrand O, Huot JP. 1999. A climate database for Mars. J Geophys Res Planet (1991–2012) 104(E10): 24177–24194. [Google Scholar]
 Mancusi D, Boudard A, Cugnon J, David JC, Kaitaniemi P, Leray S. 2014. Extension of the Liege intranuclearcascade model to reactions induced by light nuclei. Phys Rev C 90(5): 054602. [Google Scholar]
 Matthiä D, Ehresmann B, Lohf H, Kohler J, Zeitlin C, et al. 2016. The Martian surface radiation environment – A comparison of models and MSL/RAD measurements. J Space Weather Space Clim 6(13): 1–17. DOI: 10.1051/swsc/2016008. [CrossRef] [Google Scholar]
 McKennaLawlor S, Gongalves P, Keating A, Morgado B, Heynderickx D, et al. 2012. Characterization of the particle radiation environment at three potential landing sites on Mars using ESAs MEREM models. Icarus 218(1): 723–734. [CrossRef] [Google Scholar]
 NilssonAlmqvist B, Stenlund E. 1987. Interactions between hadrons and nuclei: The Lund Monte Carlo – FRITIOF version 1.6. Comput Phys Commun 43(3): 387–397. [Google Scholar]
 O’Neill PM. 2010. BadhwarO’Neill 2010 galactic cosmic ray flux model revised. IEEE Trans Nucl Sci 6(57): 3148–3153. [Google Scholar]
 Picone J, Hedin A, Drob DP, Aikin A. 2002. NRLMSISE00 empirical model of the atmosphere: Statistical comparisons and scientific issues. J Geophys Res Space Phys 107(A12): 1468. [Google Scholar]
 Rafkin SC, Zeitlin C, Ehresmann B, Hassler D, Guo J, et al. 2014. Diurnal variations of energetic particle radiation at the surface of Mars as observed by the Mars Science Laboratory Radiation Assessment Detector. J Geophys Res (Planet) 119: 1345–1358. [Google Scholar]
 Saganti PB, Cucinotta FA, Wilson JW, Schimmerling W. 2002. Visualization of particle flux in the human body on the surface of Mars. J Rad Res 43(Suppl): S119–S124. [CrossRef] [Google Scholar]
 Saganti PB, Cucinotta FA, Wilson JW, Simonsen LC, Zeitlin C. 2004. Radiation climate map for analyzing risks to astronauts on the Mars surface from galactic cosmic rays. Space Sci Rev 110(1–2): 143–156. [CrossRef] [Google Scholar]
 Sato T, Niita K, Matsuda N, Hashimoto S, Iwamoto Y, et al. 2013. Particle and heavy ion transport code system, PHITS, version 2.52. J Nucl Sci Technol 50(9): 913–923. [Google Scholar]
 Simonsen L, Nealy J, Townsend L, Wilson J. 1990. Radiation exposure for manned mars surface missions, NASA technical paper series. NASATP2979, L16708. MARCH 1990. PP. 24. Document, Vol. 1. [Google Scholar]
 Simonsen LC, Nealy JE. 1993. Mars surface radiation exposure for solar maximum conditions and 1989 solar proton events. NASA technical paper series, 3300. [Google Scholar]
 Simpson J. 1983. Elemental and isotopic composition of the galactic cosmic rays. Ann Rev Nucl Part Sci 33(1): 323–382. [Google Scholar]
 Slaba TC, Wilson JW, Badavi FF, Reddell BD, Bahadori AA. 2016. Solar proton exposure of an ICRU sphere within a complex structure part II: Raytrace geometry. Life Sci Space Res 9: 77–83. [CrossRef] [Google Scholar]
 von Forstner JLF, Guo J, WimmerSchweingruber RF, Hassler DM, Temmer M, et al. 2018. Using forbush decreases to derive the transit time of ICMEs propagating from 1 AU to Mars. J Geophys Res Space Phys 123(1): 39–56. [CrossRef] [Google Scholar]
 Wilson JW, Slaba TC, Badavi FF, Reddell BD, Bahadori AA. 2016. Solar proton exposure of an ICRU sphere within a complex structure: Combinatorial geometry. Life Sci Space Res 9: 69–76. [CrossRef] [Google Scholar]
 WimmerSchweingruber RF, Kohler J, Hassler DM, Guo J, Appel JK, et al. 2015. On determining the zenith angle dependence of the Martian radiation environment at Gale Crater altitudes. Geophys Res Lett 42(24): 10557–10564 [CrossRef] [Google Scholar]
 Winslow RM, Schwadron NA, Lugaz N, Guo J, Joyce CJ, et al. 2018. Opening a window on ICMEdriven GCR modulation in the inner solar system. Astrophys J 856(2): 139. [CrossRef] [Google Scholar]
 Witasse O, SánchezCano B, Mays M, Kajdic P, Opgenoorth H, et al. 2017. Interplanetary coronal mass ejection observed at STEREOA, Mars, comet 67P/ChuryumovGerasimenko, Saturn, and New Horizons enroute to Pluto. Comparison of its Forbush decreases at 1.4, 3.1 and 9.9 AU. J Geophys Res Space Phys 122: 7865–7890. [Google Scholar]
 Zeitlin C, Hassler D, Guo J, Ehresmann B, WimmerSchweingruber R, et al. 2018. Analysis of the radiation hazard observed by RAD on the surface of mars during the September 2017 solar particle event. Geophys Res Lett 45(12): 5845–5851. [Google Scholar]
Cite this article as: Guo j, Banjac S, Röostel L, Terasa J.C, Herbst K, et al. 2019. Implementation and validation of the GEANT4/AtRIS code to model the radiation environment at Mars. J. Space Weather Space Clim. 9, A2.
All Tables
Martian atmospheric and regolith properties used in AtRIS and four different physics lists of GEANT4 (more details in the text) tested in this study.
All Figures
Fig. 1 Martian atmospheric structure at Gale Crater implemented in the AtRIS simulations as imported from the MCD model. 

In the text 
Fig. 2 Matrices which describe the surface downward protons created by the primary protons. First row: the normalized Martian atmospheric matrices of primary protons generating surface downward directed protons from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 

In the text 
Fig. 3 Matrices which describe the surface downward deuterons created by the primary protons. First row: the normalized Martian atmospheric matrices of primary protons generating surface downward directed ^{2}H ions from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 

In the text 
Fig. 4 Matrices which describe the surface downward tritons created by the primary protons. First row: the normalized Martian atmospheric matrices of primary protons generating surface downward directed ^{3}H ions from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 

In the text 
Fig. 5 Matrices which describe the surface downward ^{4}He particles created by the primary protons. First row: the normalized Martian atmospheric matrices of primary protons generating surface downward directed ^{4}He ions from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 

In the text 
Fig. 6 Matrices which describe the surface downward ^{3}He particles created by the primary protons. First row: the normalized Martian atmospheric matrices of primary protons generating surface downward directed ^{3}He ions from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and a is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 

In the text 
Fig. 7 Matrices which describe the surface upward protons created by the primary protons. First row: the normalized Martian atmospheric matrices of primary protons generating surface upward directed protons from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 

In the text 
Fig. 8 Matrices which describe the surface downward protons created by the primary ^{4}He particles. First row: the normalized Martian atmospheric matrices of primary ^{4}He generating surface downward directed protons from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins 

In the text 
Fig. 9 Matrices which describe the surface downward ^{4}He created by the primary ^{4}He particles. First row: the normalized Martian atmospheric matrices of primary ^{4}He generating surface downward directed ^{4}He from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 

In the text 
Fig. 10 Matrices which describe the surface downward deuterons created by the primary ^{4}He particles. First row: the normalized Martian atmospheric matrices of primary ^{4}He generating surface downward directed ^{2}H from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 

In the text 
Fig. 11 Matrices which describe the surface downward ^{3}He created by the primary ^{4}He particles. First row: the normalized Martian atmospheric matrices of primary ^{4}He generating surface downward directed ^{3}He from physics lists A, B, C and D (Table 1). Second row: Normalized statistical uncertainty of the matrices. Third row: Image Differencing Matrices (IDM) of model A and B (1st column), B and C (2nd column), C and D (3rd column) which is also visualized with a 2D Gaussian filter applied (4th column). Fourth row: the normalized IDM. represents the average value in each matrix (while bins with zero statistics have been excluded) and σ is the standard deviation. All matrices are shown in the energy range from 1 MeV to 100 GeV (in logarithmic scale) for both input bins and output bins. 

In the text 
Fig. 12 Martian surface spectra of H (panel a) and He (panel b) isotopes induced by primary GCR proton and ^{4}He particles modeled via four different physics lists (Table 1). The spectra are averaged within 36° of downward zenith angle corresponding to the RAD view cone. (a) Input: GCR H + ^{4}He, output: downward H isotopes, (b) input: GCR H + ^{4}He, output: downward He isotopes. 

In the text 
Fig. 13 Martian surface spectra of H (a, b, c) and He (d, e, f) isotopes induced by primary GCR proton and ^{4}He particles modeled via four different physics lists (Table 1). The spectra are averaged within different solid angles: (a) and (d) for the RAD view cone, (b) and (e) for the downward direction and (c) and (f) for the upward direction. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.