Studies of the radiation environment on the Mars surface using the Geant4 toolkit

Studies of the radiation environment on the Mars surface using the Geant4 toolkit 1 corresp first-author chenjunliang_upc@163.com chenjunliang_upc@163.com chenjunliang_upc@163.com https://orcid.org/0000-0002-9065-2648 2 3 corresp tkdong@pmo.ac.cn tkdong@pmo.ac.cn tkdong@pmo.ac.cn 4 corresp zren@tongji.edu.cn zren@tongji.edu.cn zren@tongji.edu.cn 5 Key Laboratory of Modern Acoustics and Department of Physics, Nanjing University, Nanjing 210093, China Key Laboratory of Modern Acoustics and Department of Physics, Nanjing University, Nanjing 210093, China School of Electronic Engineering, Nanjing XiaoZhuang University, Hongjing Road, Nanjing 211171, China. School of Electronic Engineering, Nanjing XiaoZhuang University, Hongjing Road, Nanjing 211171, China. Key Laboratory of Dark Matter and Space Astronomy, Purple Mountain Observatory, CAS, 2 West Beijing Rood, Nanjing 210008, China Key Laboratory of Dark Matter and Space Astronomy, Purple Mountain Observatory, CAS, 2 West Beijing Rood, Nanjing 210008, China School of Physics Science and Engineering, Tongji University, Siping Rood 1239, Shanghai 200092, China School of Physics Science and Engineering, Tongji University, Siping Rood 1239, Shanghai 200092, China State Key Laboratory of Lunar and Planetary Sciences, Macau University of Science and Technology, Macao 999078, China State Key Laboratory of Lunar and Planetary Sciences, Macau University of Science and Technology, Macao 999078, China The radiation environment on the surface of Mars is a potential threat for future manned exploration missions to this planet. In this study, a simple geometrical model was built for simulating the radiation environment on the Mars surface caused by galactic cosmic rays (GCRs); the model was built and studied using the Geant4 toolkit. The simulation results were compared with the data reported by a radiation assessment detector (RAD). The simulated spectra of neutrons, photons, protons, α particles, and particle groups Z =3–5, Z =6–8, Z =9–13, and Z =14–24 were in a reasonable agreement with the RAD data. However, for deuterons, tritons, and 3 He, the simulations yielded much smaller values than for the corresponding RAD data. In addition, the particles’ spectra within the 90° zenith angle were also obtained. Based on these spectra, we calculated the radiation dose that would have been received by an average human body on Mars. The distribution of the dose throughout the human body was not uniform. The absorbed and equivalent doses for the brain were the highest among all of the organs, reaching 62.0±1.7 mGy/y and 234.1±8.0 mSv/y, respectively. The average absorbed and equivalent doses for the entire body were approximately 44 mGy/y and 153 mSv/y, respectively. Further analysis revealed that most of the radiation dose was owing to α particles, protons, and heavy ions. We then studied the shielding effect of the Mars soil with respect to the radiation. The body dose decreased significantly with increasing soil depth. At the depth of 1.5 m, the effective dose for the entire body was 17.9±2.4 mSv/y, lower than the dose limit for occupational exposure. At the depth of 3 m, the effective dose to the body was 2.7±1.0 mSv/y, still higher than the accepted dose limit. Galactic cosmic rays Radiation environment of the Mars Absorbed dose Equivalent dose Geant4 Introduction 1 Space exploration has developed rapidly, and significant progress has been made in the past decades. Mars has become a hotspot for deep-space exploration missions. Mars exploration programs have been carried out, to learn about the evolution of Mars and for determining whether life exists on Mars [ 1 1 ]. The manned mission to Mars has also been on schedule [ 2 2 ]. A key factor that must be considered for manned space exploration missions is the space radiation, which is caused by galactic cosmic rays (GCRs) and solar particle events (SPEs) [ 3 3 - 5 5 3-5 ]. In the space, energetic particles can penetrate spacecraft and harm the instruments and astronauts [ 6 6 ]. Radiation-related health risks constitute a serious challenge for long-term manned space exploration missions [ 7 7 , 8 8 ]. GCRs are mainly composed of charged nuclei, such as protons, α particles, and heavier nuclei with a broad range of energies, from sub-MeV to hundreds of TeV [ 9 9 , 10 10 ]. SPEs mainly consist of protons and α particles, with relatively lower energies compared with GCRs. Owing to the planet's thin atmospheric layer and the lack of a global magnetic field [ 11 11 ], the radiation intensity on Mars is much higher than that on Earth. Knowing the radiation environment on Mars and the radiation dose that is likely to be received by astronauts is very important for future manned missions to Mars [ 11 11 , 12 12 ]. This knowledge will allow decision-makers to accurately assess the corresponding health risks and is likely to provide a reference for designing radiation-protection gear [ 13 13 , 14 14 ]. The radiation environment on Mars is complex and is defined by several factors. The shielding effect of the Mars atmosphere is much weaker than that of the Earth atmosphere. As a result, energetic particles in GCRs and SEPs can penetrate the Mars atmosphere, reaching the planet’s surface. Owing to the irregularity and infrequency of SPEs, only GCRs were considered in this study. The interaction between GCRs and atmospheric molecules generates many secondary particles, which further interact with the planet’s atmosphere, resulting in the complex radiation spectrum on the Mars surface. Albedo particles that are produced on the planet’s surface also contribute to the radiation environment. The three main components of the radiation environment on the Mars surface were considered in this study. The Curiosity rover has carried out a long-term scientific exploration on Mars since it successfully landed on the planet in 2012. At the same time, an onboard radiation assessment detector (RAD) has been collecting radiation-related information on the Mars surface, including the particles' spectrum and the absorbed dose rate [ 15 15 ]. Some studies have shown that long-term exposure to the space radiation can harm humans [ 2 2 ]. Radiation can directly and indirectly affect DNA molecules in the human body, causing cell death or mutation and increasing the risk of cancer [ 16 16 - 18 18 16-18 ] and other diseases. The National Aeronautics and Space Administration (NASA) categorized the health risks owing to the space radiation into four groups: 1) degenerative tissue effects, 2) carcinogenesis, 3) acute and late central nervous system (CNS) effects, and 4) radiation-related syndromes [ 19 19 ]. Accurate calculation of the radiation dose received by astronauts is necessary for assessing the health risks associated with the space radiation. However, it is difficult to directly calculate the particles’ spectra on the Mars surface and the associated radiation dose. Monte Carlo simulations can be used for addressing this problem. Some transport codes have been used for studying the radiation environment in the space and on Mars, such as HZETRN [ 20 20 , 21 21 ], HETC–HEDS [ 22 22 ], FLUKA [ 23 23 ], PHITS [ 24 24 ], MCNP6 [ 25 25 , 26 26 ] and Geant4 [ 27 27 - 30 30 27-30 ]. In Ref. [ 20 20 - 26 26 20-26 ], Slaba et al. studied the radiation environment on the Mars surface using these tools. Most of these studies have mainly focused on simulating the radiation environment on the Mars surface. Our work addresses two aspects of this research area. First, we used the G4QMD model to deal with the inelastic scattering of GCRs in the Mars atmosphere and soil, and for obtaining the particles’ spectra on the Mars surface. Then, based on the simulated particles’ spectra, we calculated the radiation dose received by a typical human body. To our knowledge, this is the first study attempting to estimate the radiation dosage to individual organs of a typical human on the Mars surface. The results and conclusions of this study can be used when planning future manned Mars exploration missions. The radiation environment on Mars 1 Model setup 2 To reduce the simulation time, a simple geometrical model consisting of a cylinder (radius, 1 km) was built, to represent the Martian environment. The cylinder was divided into 12 layers to represent the Martian regolith and atmosphere. The first layer from the bottom corresponded to the Martian regolith (thickness, 3 m; density, 1.7 g/cm 3 ) [ 21 21 ]. The second layer was 4-m-thick, representing the atmospheric conditions on the Mar s surface. From the third layer to the twelfth layer, each layer was 1-m-thick, and the column densities of the different layers were 11.820, 5.915, 2.624, 1.069, 0.395, 0.131, 0.035, 0.009, 0.003, and 0.001 g/cm 2 , respectively. The overall column density of these layers was approximately 22 g/cm 2 [ 21 21 , 31 31 , 32 32 ], representing the Mars atmosphere. The detector was a cylinder (radius, 1 km; thickness, 3 mm) and was positioned1 m above the soil, for detecting and registering the Mars radiation environment. The planet’s soil and atmosphere components were derived from Ref. [ 21 21 ] and are listed in Table 1 Table 1 . Primary GCRs were incident from the top of the geometrical structure, in the isotropic fashion. Once a particle passed through the detector, the information about the event was recorded, including the particle's type, its kinetic energy, and its angle of incidence. 1 1   1 1 Martian atmosphere 1 1 Martian atmosphere 1 1 CO_2(95.7%), N (2.7%), Ar(1.6%) 1 1 2*Martian soil 1 1 SiO_2 (51.2%),Fe_2O_3 (9.3%),H_2O (7.4%) 1 1   1 1 Al_2CaK_2MgNa_2O_7 (32.1%) The radiation environment simulations were performed using the Geant4-10.5.1 package (version10.5.1). This code is a track structure code that was developed at CERN for simulating the passage of particles through matter [ 29 29 , 30 30 ]. It provides many models that can be chosen by the user according to the specific application, and it has been widely used in high-energy physics, medical physics, and radiobiology. The physical parameters that were used for simulating the interaction between GCRs and the Mars medium are listed in Table 2 Table 2 . The Emstandard_opt3 physics list was used for describing the electromagnetic interactions of photons, leptons, hadrons, and ions, with energies in the 1 keV–10 PeV range [ 33 33 ]. The G4HadronPhysicsQGSP_BIC_HP physics list [ 34 34 ] was used for describing the inelastic scattering of nucleons and mesons with nuclei, for energies in the 0 eV–100 TeV range. In this physics list, the quark gluon string model (QGS) was used for describing high-energy hadronic interactions (12 GeV–100 TeV range). The binary cascade model (BIC) was used for dealing with nucleon-nucleus interactions at intermediate energies, and high-precision neutron models (HPs) were used for dealing with neutron-nucleus interactions involving energies below 20 MeV. For the inelastic interactions of nuclei, the G4QMD model that was validated in our previous work was used [ 35 35 ]. In this model, all nucleons were considered as participant particles, and each nucleon was treated as a Gaussian wave packet. Equations of the particles’ motion were derived based on the time-dependent variation principle [ 36 36 ]. A detailed description of these models can be found in the Geant4 Physics Reference Manual[ 36 36 ]. 1 1 Models 1 1 Description 1 1 Emstandard_opt3 1 1 Electromagnetic model 1 1 G4HadronPhysicQGSP_BIC_HP 1 1 Hadronic inelastic model 1 1 G4IonQMDPhysics 1 1 Hadronic inelastic model for ions In this study, only GCRs were considered as contributing to the radiation environment of Mars. GCRs mainly consist of protons, α particles, and heavy ions. Although the proportion of heavy particles is very small, they are important for radiation protection, owing to their high relative biological effectiveness (RBE). In our simulations, particles with atomic numbers Z =1–26 were considered for GCRs. The energy spectra of these particles were calculated using the GCR model developed by Matthi ä et al. [ 10 10 ]. In that model, a single parameter W was used for describing the modulation effect of solar activity [ 10 10 ]. A part of the primary GCR spectra used in the simulation is shown in Fig. 1 Fig. 1 . The modulation parameter W was calculated according to the Oulu neutron monitor count rate (from http://cosmicrays.oulu.fi/) during the period from August 2012 to January 2013, and an average value ( W =56) was used in this work. Results 2 When GCRs pass through a planet's atmosphere, many neutrons and photons are induced owing to the interactions between the GCRs and atmospheric molecules. The simulated neutron and photon spectra on the Mars surface are plotted in Fig. 2 Fig. 2 , together with the experimental data [ 21 21 , 37 37 ]. As shown in this figure, the simulation results for neutrons were slightly higher than the corresponding experimental data, for energies below 100 MeV/n, while the simulation results for photons agreed well with the RAD data. For neutrons and photons, the average relative deviations of the simulation results from the corresponding experimental data were 60.8% and 56.0%, respectively. The simulation results for the charged particles’ spectra, for zenith angles smaller than 30°, are plotted in Figs. 3 Figs. 3 and 4 4 . Heavy particles can be divided into several groups, for proper comparison to the RAD data[ 15 15 ]. Evidently, the simulated spectra of protons, α particles, and particle groups Z =6-9, Z =9-13, and Z =14-24 were consistent with the corresponding RAD data, and the corresponding average relative deviations were 27.7%, 37.9%, 40.9%, 60.7%, and 40.1%. For the particle groups Z =3-5 and Z > 24, the simulation results were slightly lower than the measured data. However, the simulated spectra of deuterons, tritons, and 3 He were much lower than the corresponding RAD data, and the average relative deviation was approximately 90%. These particles are pure secondary particles generated in reactions between primary GCRs and atmospheric molecules, and the observed discrepancy might be owing to the lower generation cross-sections. Based on the previous work, the cross-sections for light fragments simulated by the G4QMD model were lower than the experimental data[ 35 35 ]. Heavy charged particles incident from the entire upper hemisphere were also recorded in the simulation, and the simulated spectra for zenith angles smaller than 90 are plotted in Fig. 5 Fig. 5 . Compared with the spectrum within the 30° zenith angle, for light particles, the flux averaged over the 90° zenith angle was close to the flux for zenith angles below 30°. However, for heavy particles, the flux for the zenith angles below 90° was slightly smaller than the flux for zenith angles below 30°. The spectra of electrons, muons, and pions were also calculated, owing to their significant contributions to the absorbed dose, and are plotted in Fig. 6 Fig. 6 . Radiation dose received by a typical human body 1 Accurate estimation of the dose received by astronauts is an essential precondition for successful manned deep-space exploration missions. In addition, with the wide use of nuclear technology in material physics and medical physics [ 38 38 - 46 46 38-46 ], it is also important to be able to calculate the radiation dose received by a typical human in the corresponding radiation environment. Based on the prior simulation results, we calculated the radiation dose that would have been received by astronauts on the Mars surface. According to the International Commission on Radiological Protection (ICRP) [ 47 47 ], the equivalent dose for a tissue, H T , is defined by <math display="block" id="M1"><mrow><msub><mi>H</mi><mi>T</mi></msub><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>R</mi></munder></mstyle><msub><mi>W</mi><mi>R</mi></msub><msub><mi>D</mi><mrow><mi>T</mi><mn>,</mn><mi>R</mi></mrow></msub></mrow></math> where D T,R is the absorbed dose for tissue T caused by radiation R , and W R is the radiation weighting factor for radiation R . The effective dose, E , is defined as <math display="block" id="M2"><mrow><mi>E</mi><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>T</mi></munder></mstyle><msub><mi>W</mi><mi>T</mi></msub><msub><mi>H</mi><mi>T</mi></msub></mrow></math> where W T is the tissue weighting factor for tissue T . Table 3 Table 3 and Table 4 Table 4 list the values of W T and W R , respectively. The ICRP recommends that the effective dose received during occupational exposure and in public should not exceed 20 mSv/y and 1 mSv/y [ 47 47 ], respectively. 1 1 Tissue 1 1 W T 1 1 Σ W T 1 1 Bone-marrow (red), colon, 1 1   1 1   1 1 Lung, stomach, breast, remainder tissues 1 1 0.12 1 1 0.72 1 1 Gonads 1 1 0.08 1 1 0.08 1 1 Bladder, oesophagus, liver, thyroid 1 1 0.04 1 1 0.16 1 1 Bone surface, brain, salivary glands, skin 1 1 0.01 1 1 0.04 1 1 Radiation type ( R ) 1 1 Radiation weighting factors ( W R ) 1 1 Photons 1 1 1 1 1 Electrons and muons 1 1 1 1 1 Protons and charged pions 1 1 2 1 1 α particles, fission fragments, heavy ions 1 1 20 1 1 Neutrons 1 1 2.5+18.2 exp {-[ln( E n ) 2 /6}, E n < 1 MeV 1 1   1 1 5.0+17.0 exp {- [ln (2 E n )] 2 /6}, 1 MeV ≤ E n ≤ 50 MeV 1 1   1 1 2.5+3.25 exp {-[ln(0.04 E n )] 2 /6}, E n > 50 MeV Model setup 2 The geometrical simulation model is shown in Fig. 7 Fig. 7 . The human phantom model [ 48 48 ] was borrowed from the Geant4 advanced example "human phantom" and placed on the Mars surface. A cubical regolith (10×10×6 m 3 ) was built, with the density of 1.7 g/cm 3 . Then, the simulated particles' spectra within the 90° zenith angle (in Figs. 5 Figs. 5 – 6 6 ) were used for simulating primary particles and projecting them onto the surface from the top of the geometrical model. For each step in the simulation, the energy deposited by different particles into different human body organs was recorded. After the simulation, the values of the total deposited energy, the absorbed dose, and the equivalent dose delivered by different particles to different organs were calculated. In the simulation, the total sampling number of each primary particle was set to 5 million. For each scenario, the simulation was repeated six times, and the uncertainties shown in this paper represent, in each case, one standard deviation from the corresponding mean. Results 2 Figure 8 Figure 8 shows the radiation dose received by different human body organs, for a typical human staying on the Mars surface for one year. The distribution of the deposited energies in the body organs is plotted in Fig. 8(a) Fig. 8(a) . Evidently, the energy deposited in the trunk was the highest, reaching (9.05±0.04)×10 12  MeV/y. The body part receiving the second highest amount of the radiation energy were legs, with the deposited energy per leg reaching (2.02±0.02)×10 12  MeV/y. This large deposition occurred because the mass and volume of the trunk and legs were larger than those of the other body parts. The amount of energy deposited in the head and skull was also high; however, the amounts of energy deposited in smaller organs, such as left/right adrenal and thyroid glands, wer among the lowest. Figure 8(b) Figure 8(b) shows the simulated absorbed dose distribution across the different body parts. The distribution of the absorbed dose differs significantly from the distribution of deposited energy. The absorbed dose was the highest for the brain, reaching 62.0±1.7 mGy/y, followed by the thyroid and skull. The absorbed doses for the thymus, upper spine, left/right lung, and left/right clavicle were also relatively high; however, the absorbed doses for the lower body organs were relatively low; for example, the absorbed doses for the left/right leg and leg bone were approximately 35 mGy/y. Evidently, the difference between the doses received by the different organs was significant, because primary particles were incident from above; consequently, some smaller-energy particles only affected the superficial layers of the human body. Thus, upper-body organs, such as the brain and the thyroid gland, received higher radiation doses than lower-body organs. The average absorbed dose for the entire body was 43.8±0.4 mGy/y. Taking the quality of different radiation sources into account, we calculated the equivalent doses for the different organs, and the results are shown in Fig. 8(c) Fig. 8(c) . The equivalent dose for the brain was the highest, reaching 234.1±8.0 mSv/y, agreeing well with the corresponding RAD data [ 31 31 ]. For the skull, adrenal, and upper spine, the equivalent doses were also relatively high, all above 200 mSv/y. For the leg bone, the equivalent dose was the lowest, just 116.4±1.7 mSv/y. The average equivalent dose for the entire body was approximately 153 mSv/y, much lower than the corresponding RAD data. We posit that this difference could be attributed to the fact that a typical human body is much larger than the RAD detector, and because the mass distribution in the body is not uniform. The lower part of the body and inner organs received l ess radiation, so that the average equivalent dose for the entire body was lower than the corresponding RAD data. To analyze the contributions of different particles to the radiation dose received by the human body, we extracted the deposited energy and the absorbed dose in the human body owing to electrons, photons, pions, muons, protons, α particles, and heavy ions ( Z ≥3), respectively. The results are plotted in Fig. 9 Fig. 9 , except for photons, because their contribution was far less than those of other particles. The energy deposited by protons was (1.08±0.01)×10 13 MeV/y, accounting for 56.1% of the total energy deposited in the human body. The second largest contribution was that of electrons; the corresponding deposited energy accounted for 22.6% of the total energy deposited in the body. The energies depositions by α particles, heavy ions, pions, and muons were smaller, accounting for 6.1%, 3.6%, 6.6%, and 3.9% of the total deposited energy, respectively. The contribution of neutrons to the total deposited energy was only 0.07%, much lower than that of protons. This was because the reaction cross-section of neutrons with matter was very small, and most of the deposited energy was owing to secondary particles. The contributions of the different particles to the absorbed dose were similar to those for the deposited energy. However, for the equivalent dose, the situation was different. α particles contributed mostly to the total equivalent dose, which reached 53.7±1.3 mSv/y, accounting for 35.1% of the total equivalent dose. The equivalent doses owing to protons and heavy ions were 49.1±0.6 mSv/y and 31.9±0.7 mSv/y, accounting for 32.1% and 20.9% of the total equivalent dose, respectively. The contributions of electrons, pions, muons, and neutrons were slightly smaller, accounting for only 6.5%, 3.8%, 1.4%, and 0.24% of the total equivalent dose, respectively. Although the deposited energies owing to α particles and heavy ions were not very high, they were the main contributors to the equivalent dose. This was because the radiation weighting factors of α particles and heavy ions were much larger than those of the other particles. Therefore, it is necessary to consider heavy ions in GCRs, when designing future manned deep-space exploration missions. It can be observed from Figs. 2 Figs. 2 , 3 3 and 4 4 that the simulated spectra of deuterons, tritons, and 3 He were much lower than those of the measured data. To assess the effect of the differences on the calculated doses to the human body, we multiplied the spectra by a scaling factor, for matching with the measured data, and recalculated the dose. The absorbed dose and the dose equivalents for the body organs, based on the scaled spectra, increased by approximately 22% on average, compared with the previous results. The average absorbed dose and the dose equivalent for the human body were 54.0 mGy/y and 191.6 mSv/y, respectively. The results show that these particles strongly affect the dose calculation, and the models in the Geant4 package should be improved for better describing nucleus-nucleus interactions. The shielding effect of the Martian soil 1 The simulation results suggest that the dose that would have been received by astronauts on the Mars surface is much higher than that received on the Earth surface. Therefore, appropriate radiation protection measures are necessary. When astronauts will carry out long-term exploration missions on the Mars surface, it is likely that they will use Martian soil to build shelters, for protecting themselves from radiation. Thus, the shielding effect of Martian soil is worth investigating. Model setup 2 A hollow cube (5×5×3 m 3 ) was built and placed in the Mars soil, and the human phantom was placed in the cube. The thickness of the soil above the human phantom, denoted by d, was adjustable. In the simulation, we considered 9 levels of the soil thickness, namely d =0.1 m, 0.3 m, 0.5 m, 0.7 m, 1.0 m, 1.5 m, 2.0 m, 2.5 m, and 3 m; for each one of these scenarios we computed the dose received by the human body. The primary particles were projected from the top of the soil according to the simulated particles' spectra on the Mars surface. The tissue energy deposition owing to different particles was recorded. Results 2 The radiation doses on the human body at different depths under the Mars surface are plotted in Fig. 10 Fig. 10 . The deposited energy and dose rate on the Mars surface were also included for comparison. Evidently, with increasing the soil depth, the deposited energy, absorbed dose, and equivalent dose in the organs all decreased rapidly. At the depths of 0.1 m, 0.5 m, and 1.0 m, the absorbed dose was lower by 7.6%, 41.3%, and 71.1%, respectively, compared with the dose on the Mars surface. The reduction rate of the absorbed dose decreased as the soil thickness continued to increase. As shown in Fig. 10(c) Fig. 10(c) , the equivalent dose in the organs dropped faster than the absorbed dose at small depths. For example, at the depth of 0.1 m, the equivalent doses on the brain and whole body decreased by 30.6% and 28.1%, respectively. We attribute this to the fact that energetic heavy ions collided with the Mars soil and broke into lighter fragments, with smaller radiation weighting factors. The variation of the effective dose with the soil depth was similar to that of the dose equivalent and at the soil depths of 0.1 m, 0.5 m, and 1.0 m, the effective dose on the entire body was 121.0±9.7 mSv/y, 67.6 ±7.9 mSv/y, and 35.0±2.4 mSv/y, respectively. At the depth of 1.5 m, the average equivalent dose and effective dose on the entire body were 17.4±0.5 mSv/y and 17.9±2.4 mSv/y, respectively, lower than the dose limit of 20 mSv/y for occupational exposure. At the depth of 3 m, the equivalent doses received by different organs were between 2.0 mSv/y and 3.8 mSv/y, and the effective dose on the entire body was 2.7±1.0 mSv/y, still higher than the dose limit for public exposure. To confine the effective dose below 1 mSv/y, the thickness of the soil needs to be larger than 3 m. Conclusions 1 In this study, we first simulated the radiation environment on the Mars surface using the Geant4 toolkit. In the simulation, a simple geometrical model and G4QMD physics model were used. The simulation results for neutrons, photons, protons, α particles, and particle groups Z =6–8, Z =9–13, and Z =14–24 were in a reasonable agreement with the corresponding RAD data, with average deviations of 60.8%, 56.0%, 27.7%, 37.9%, 40.9%, 60.7%, and 40.1%, respectively. However, for deuterons, tritons, and 3 He, which are pure secondary particles, the simulation results were much smaller than the corresponding RAD data, which might have been owing to the underestimation of the production cross-sections in the G4QMD model. In addition, the particles’ spectrum from the entire upper hemisphere was also obtained. Based on the simulated radiation environment, we calculated the radiation dose received by the human body on the Mars surface. The results showed that the doses absorbed by different organs differed significantly from each other. The absorbed doses for the brain, thyroid, and skull were relatively high, but the doses for the legs and male genitalia were much lower. The average absorbed dose and equivalent dose of the entire body were 43.8±0.4 mGy/y and 152.9±2.0 mSv/y, respectively. Further analysis showed that most of the dose was induced by α particles, protons, and heavy ions. The shielding effect of the Mars soil on radiation was also studied. With increasing soil depth, the deposited energy, absorbed dose, and equivalent dose for the organs decreased rapidly. At the depth of 1.5 m, the effective dose for the entire body was 17.9±2.4 mSv/y, lower than the dose limit recommended by the ICRP for occupational exposure. However, for confining the effective dose to the safe range for public exposure, the soil thickness should exceed 3 m. References G.B. Yu , E.H. Liu , G.L. Liu , et al. , Moderate resolution imaging camera (MoRIC) of China’s first mars mission Tianwen-1 . Earth Planet. Phys. 4 ( 4 ), 364 - 370 ( 2020 ). doi: 10.26464/epp2020056 http://doi.org/10.26464/epp2020056 Moderate resolution imaging camera (MoRIC) of China’s first mars mission Tianwen-1 M. Durante , Space radiation protection: Destination mars . Life Sci. Space Res. (Amst). 1 2 - 9 ( 2014 ). doi: 10.1016/j.lssr.2014.01.002 http://doi.org/10.1016/j.lssr.2014.01.002 Space radiation protection: Destination mars F.A. Cucinotta , M. Durante , Cancer risk from exposure to galactic cosmic rays: implications for space exploration by human beings . The Lancet Oncology 7 , 431 - 435 ( 2006 ). doi: 10.1016/S1470-2045(06)70695-7 http://doi.org/10.1016/S1470-2045(06)70695-7 Cancer risk from exposure to galactic cosmic rays: implications for space exploration by human beings N.A. Schwadron , L. Townsend , K. Kozarev et al. , Earth-moon-mars radiation environment module framework . Space Weather 8 , S00E02 ( 2010 ). doi: 10.1029/2009SW000523 http://doi.org/10.1029/2009SW000523 Earth-moon-mars radiation environment module framework Á. Vicente-Retortillo , F. Valero , L. Vázquez et al. , A model to calculate solar radiation fluxes on the Martian surface . J. Space Weather Space Clim. 5 , A33 ( 2015 ). doi: 10.1051/swsc/2015035 http://doi.org/10.1051/swsc/2015035 A model to calculate solar radiation fluxes on the Martian surface F.F. Badavi , D.O. Adams , J.W. Wilson , On the validity of the aluminum equivalent approximation in space radiation shielding applications . Adv. Space Research 46 , 719 - 727 ( 2010 ). doi: 10.1016/j.asr.2010.04.006 http://doi.org/10.1016/j.asr.2010.04.006 On the validity of the aluminum equivalent approximation in space radiation shielding applications M. Durante , F.A. Cucinotta , Physical basis of radiation protection in space travel . Rev. Mod. Phys. 83 , 1245 - 1281 ( 2011 ). doi: 10.1103/RevModPhys.83.1245 http://doi.org/10.1103/RevModPhys.83.1245 Physical basis of radiation protection in space travel M. Vuolo , G. Baiocco , S. Barbieri et al. , Exploring innovative radiation shielding approaches in space: A material and design study for a wearable radiation protection spacesuit . Life Sci. Space Res. (Amst) 15 , 69 - 78 ( 2017 ). doi: 10.1016/j.lssr.2017.08.003 http://doi.org/10.1016/j.lssr.2017.08.003 Exploring innovative radiation shielding approaches in space: A material and design study for a wearable radiation protection spacesuit A.I. Mrigakshi , D. Matthiä , T. Berger , et al. , How Galactic Cosmic Ray models affect the estimation of radiation exposure in space . Adv. Space Res. 51 , 825 - 834 ( 2013 ). doi: 10.1016/j.asr.2012.10.017 http://doi.org/10.1016/j.asr.2012.10.017 How Galactic Cosmic Ray models affect the estimation of radiation exposure in space D. Matthiä , T. Berger , A.I. Mrigakshi , et al. , A ready-to-use galactic cosmic ray model . Adv. Space Research 51 , 329 - 338 ( 2013 ). doi: 10.1016/j.asr.2012.09.022 http://doi.org/10.1016/j.asr.2012.09.022 A ready-to-use galactic cosmic ray model B. Ehresmann , C.J. Zeitlin , D.M. Hassler , et al. , The charged particle radiation environment on Mars measured by MSL/RAD from November 15, 2015 to January 15, 2016 . Life Sci. Space Res. (Amst) 14 , 3 - 11 ( 2017 ). doi: 10.1016/j.lssr.2017.07.004 http://doi.org/10.1016/j.lssr.2017.07.004 The charged particle radiation environment on Mars measured by MSL/RAD from November 15, 2015 to January 15, 2016 T. Sato , K. Niita , V.A. Shurshakov , et al. , Evaluation of dose rate reduction in a spacecraft compartment due to additional water shield . Cosmic Research 49 , 319 - 324 ( 2011 ). doi: 10.1134/S0010952511040083 http://doi.org/10.1134/S0010952511040083 Evaluation of dose rate reduction in a spacecraft compartment due to additional water shield F.A. Cucinotta , S. Hu , N.A. Schwadron et al. , Space radiation risk limits and Earth-Moon-Mars environmental models . Space Weather 8 , S00E09 ( 2010 ). doi: 10.1029/2010SW000572 http://doi.org/10.1029/2010SW000572 Space radiation risk limits and Earth-Moon-Mars environmental models S. Trovati , F. Ballarini , G. Battistoni et al. , Human exposure to space radiation: role of primary and secondary particles . Radiat. Prot. Dosimetry 122 , 362 - 366 ( 2006 ). doi: 10.1093/rpd/ncl438 http://doi.org/10.1093/rpd/ncl438 Human exposure to space radiation: role of primary and secondary particles B. Ehresmann , C. Zeitlin , D.M. Hassler et al. , Charged particle spectra obtained with the Mars Science Laboratory Radiation Assessment Detector (MSL/RAD) on the surface of Mars . J. Geophys. Res. Planets 119 , 468 - 479 ( 2014 ). doi: 10.1002/2013JE004547 http://doi.org/10.1002/2013JE004547 Charged particle spectra obtained with the Mars Science Laboratory Radiation Assessment Detector (MSL/RAD) on the surface of Mars L. de la Fuente Rosales , S. Incerti , Z. Francis et al. , Accounting for radiation-induced indirect damage on DNA with the Geant 4-DNA code . Phys. Med. 51 , 108 - 116 ( 2018 ). doi: 10.1016/j.ejmp.2017.09.048 http://doi.org/10.1016/j.ejmp.2017.09.048 Accounting for radiation-induced indirect damage on DNA with the Geant 4-DNA code N. Lampe , M. Karamitros , V. Breton et al. , Mechanistic DNA damage simulations in Geant4-DNA part 1: A parameter study in a simplified geometry . Phys. Med. 48 , 135 - 145 ( 2018 ). doi: 10.1016/j.ejmp.2018.02.011 http://doi.org/10.1016/j.ejmp.2018.02.011 Mechanistic DNA damage simulations in Geant4-DNA part 1: A parameter study in a simplified geometry J. Chen , S. Yun , T. Dong et al. , Investigate the radiation-induced damage on an atomistic DNA model by using Geant4-DNA toolkit . Nucl. Instrum. Meth. Sect. B 494-495 , 59 - 67 ( 2021 ). doi: 10.1016/j.nimb.2021.03.010 http://doi.org/10.1016/j.nimb.2021.03.010 Investigate the radiation-induced damage on an atomistic DNA model by using Geant4-DNA toolkit J.C. Chancellor , G.B. Scott , J.P. Sutton , Space radiation: the number one risk to astronaut health beyond low earth orbit . Life (Basel) , 4 , 491 - 510 ( 2014 ). doi: 10.3390/life4030491 http://doi.org/10.3390/life4030491 Space radiation: the number one risk to astronaut health beyond low earth orbit T.C. Slaba , N.N. Stoffle , Evaluation of HZETRN on the Martian surface: Sensitivity tests and model results , Life Sci. Space Res. (Amst) , 14 , 29 - 35 ( 2017 ). doi: 10.1016/j.lssr.2017.03.001 http://doi.org/10.1016/j.lssr.2017.03.001 Evaluation of HZETRN on the Martian surface: Sensitivity tests and model results D. Matthiä , B. Ehresmann , H. Lohf , et al. , The Martian surface radiation environment—a comparison of models and MSL/RAD measurements . J. Space Weather Space Clim. 6 , A13 ( 2016 ). doi: 10.1051/swsc/2016008 http://doi.org/10.1051/swsc/2016008 The Martian surface radiation environment—a comparison of models and MSL/RAD measurements W.C. de Wet , L.W. Townsend , A calculation of the radiation environment on the Martian surface . Life Sci. Space Res. (Amst) 14 , 51 - 56 ( 2017 ). doi: 10.1016/j.lssr.2017.07.008 http://doi.org/10.1016/j.lssr.2017.07.008 A calculation of the radiation environment on the Martian surface J.W. Norbury , T.C. Slaba , N. Sobolevsky et al. , Comparing HZETRN, SHIELD, FLUKA and GEANT transport codes . Life Sci. Space Res. (Amst) 14 , 64 - 73 ( 2017 ). doi: 10.1016/j.lssr.2017.04.001 http://doi.org/10.1016/j.lssr.2017.04.001 Comparing HZETRN, SHIELD, FLUKA and GEANT transport codes J. Flores-McLaughlin , Radiation transport simulation of the Martian GCR surface flux and dose estimation using spherical geometry in PHITS compared to MSL-RAD measurements . Life Sci. Space Res. (Amst) 14 , 36 - 42 ( 2017 ). doi: 10.1016/j.lssr.2017.07.007 http://doi.org/10.1016/j.lssr.2017.07.007 Radiation transport simulation of the Martian GCR surface flux and dose estimation using spherical geometry in PHITS compared to MSL-RAD measurements D. Matthiä , D.M. Hassler , W. de Wet , et al. , The radiation environment on the surface of Mars—Summary of model calculations and comparison to RAD data . Life Sci. Space Res. (Amst) 14 , 18 - 28 ( 2017 ). doi: 10.1016/j.lssr.2017.06.003 http://doi.org/10.1016/j.lssr.2017.06.003 The radiation environment on the surface of Mars—Summary of model calculations and comparison to RAD data H.N. Ratliff , M.B.R. Smith , L. Heilbronn , Simulation of the GCR spectrum in the Mars curiosity rover's RAD detector using MCNP6 . Life Sci. Space Res (Amst) , 14 , 43 - 50 ( 2017 ). doi: 10.1016/j.lssr.2017.07.003 http://doi.org/10.1016/j.lssr.2017.07.003 Simulation of the GCR spectrum in the Mars curiosity rover's RAD detector using MCNP6 D. Matthiä , T. Berger , The radiation environment on the surface of Mars—Numerical calculations of the galactic component with GEANT4/PLANETOCOSMICS . Life Sci. Space Res. 14 , 57 - 63 ( 2017 ). doi: 10.1016/j.lssr.2017.03.005 http://doi.org/10.1016/j.lssr.2017.03.005 The radiation environment on the surface of Mars—Numerical calculations of the galactic component with GEANT4/PLANETOCOSMICS S. Agostinelli , J. Allison , K. Amako et al. , Geant4—a simulation toolkit . Nucl. Instrum. Meth. Sect. A 506 , 250 - 303 ( 2003 ). doi: 10.1016/S0168-9002(03)01368-8 http://doi.org/10.1016/S0168-9002(03)01368-8 Geant4—a simulation toolkit J. Allison , K. Amako , J. Apostolakis et al. , Geant4 developments and applications . IEEE T. Nucl. Sci. 53 , 270 - 278 ( 2006 ). doi: 10.1109/TNS.2006.869826 http://doi.org/10.1109/TNS.2006.869826 Geant4 developments and applications J. Allison , K. Amako , J. Apostolakis et al. , Recent developments in Geant4 . Nucl. Instrum. Meth. Sect. A 835 , 186 - 225 ( 2016 ). doi: 10.1016/j.anucene.2014.08.021 http://doi.org/10.1016/j.anucene.2014.08.021 Recent developments in Geant4 D.M. Hassler , C. Zeitlin , R.F. Wimmer-Schweingruber , et al. , Mars’ surface radiation environment measured with the Mars science laboratory's curiosity rover . Science , 343 , 1244797 ( 2014 ). doi: 10.1126/science.1244797 http://doi.org/10.1126/science.1244797 Mars’ surface radiation environment measured with the Mars science laboratory's curiosity rover C.G. Justus , D.L. Johnson , Mars Global Reference Atmospheric Model 2001 Version (Mars-GRAM 2001) Users Guide , Marshall Space Flight Center , NASA/TM-2001-210961 ( 2001 ). Mars Global Reference Atmospheric Model 2001 Version (Mars-GRAM 2001) Users Guide H. Burkhardt , V. M. Grichine , P. Gumplinger et al. , Geant4 standard electromagnetic package for HEP applications . IEEE Symposium Conference Record Nuclear Science 2004 , 3 , 1907 - 1910 ( 2004 ). doi: 10.1109/NSSMIC.2004.1462617 http://doi.org/10.1109/NSSMIC.2004.1462617 Geant4 standard electromagnetic package for HEP applications A. Schumann , J. Petzoldt , P. Dendooven et al. , Simulation and experimental verification of prompt gamma-ray emissions during proton irradiation . Phys. Med. Biol. 60 4197 - 4207 ( 2015 ). doi: 10.1088/0031-9155/60/10/4197 http://doi.org/10.1088/0031-9155/60/10/4197 Simulation and experimental verification of prompt gamma-ray emissions during proton irradiation J. Chen , S. Yun , T. Dong , et al. , Validation of Geant4 physics models for nuclear beams in extended media . Nucl. Instrum. Meth. Sect. B 434 , 113 - 119 ( 2018 ). doi: 10.1016/j.nimb.2018.08.022 http://doi.org/10.1016/j.nimb.2018.08.022 Validation of Geant4 physics models for nuclear beams in extended media Physics Reference Manual . https://geant4-userdoc.web.cern.ch/UsersGuides/PhysicsReferenceManual/html/index.html https://geant4-userdoc.web.cern.ch/UsersGuides/PhysicsReferenceManual/html/index.html Physics Reference Manual J. Köhler , C. Zeitlin , B. Ehresmann , et al. , Measurements of the neutron spectrum on the Martian surface with MSL/RAD . J. Geophys. Research Planets 119 , 594 - 603 ( 2014 ). doi: 10.1002/2013JE004539 http://doi.org/10.1002/2013JE004539 Measurements of the neutron spectrum on the Martian surface with MSL/RAD W.M. Sun , N. Du , WD. Tian et al. , Secondary and activated X( γ ) radiation of SPHIC particle therapy facility , Nucl. Sci. Tech. 32 ( 4 ), 1 - 8 ( 2021 ). doi: 10.1007/s41365-021-00870-6 http://doi.org/10.1007/s41365-021-00870-6 Secondary and activated X(γ) radiation of SPHIC particle therapy facility CY. Li , XB. Xia , J. Cai et al. , Radiation dose distribution of liquid fueled thorium molten salt reactor . Nucl. Sci. Tech. 32 : 22 ( 2021 ). doi: 10.1007/s41365-021-00857-3 http://doi.org/10.1007/s41365-021-00857-3 Radiation dose distribution of liquid fueled thorium molten salt reactor C.M. Niu , M.H. Li , J.R. Dai , Cage-like radiotherapy system for noncoplanar radiotherapy . Nucl. Sci. Tech. 32 : 12 ( 2021 ). doi: 10.1007/s41365-021-00848-4 http://doi.org/10.1007/s41365-021-00848-4 Cage-like radiotherapy system for noncoplanar radiotherapy P. Huang , S.Q. Li , R.D. Yang et al. , Photon energy response optimization using few-channel spectroscopy dose method for Si-PIN photodetector applied in personal dose equivalent measurements . Nucl. Sci. Tech. 31 , 36 ( 2020 ). doi: 10.1007/s41365-020-0748-2 http://doi.org/10.1007/s41365-020-0748-2 Photon energy response optimization using few-channel spectroscopy dose method for Si-PIN photodetector applied in personal dose equivalent measurements X.J. Ding , M. Yu , X. Zheng et al. , Stability study of Disperse Blue 79 under ionizing radiation . Nucl. Sci. Tech. 31 : 21 ( 2020 ). doi: 10.1007/s41365-020-0724-x http://doi.org/10.1007/s41365-020-0724-x Stability study of Disperse Blue 79 under ionizing radiation Z. Ahmadi Ganjeh , M. Eslami-Kalantari , A.A. Mowlavi , Dosimetry calculations of involved and noninvolved organs in proton therapy of liver cancer: a simulation study . Nucl. Sci. Tech. 30 , 173 ( 2019 ). doi: 10.1007/s41365-019-0698-8 http://doi.org/10.1007/s41365-019-0698-8 Dosimetry calculations of involved and noninvolved organs in proton therapy of liver cancer: a simulation study C.Q. Chen , X.B. Xia , Z.H. Zhang et al. , Radiological environmental impact analysis of a 2-MW thorium molten salt reactor during an accident . Nucl. Sci. Tech. 30 , 78 ( 2019 ). doi: 10.1007/s41365-019-0605-3 http://doi.org/10.1007/s41365-019-0605-3 Radiological environmental impact analysis of a 2-MW thorium molten salt reactor during an accident S. Abdullahi , A.F. Ismail , S. Mohd Fadzil , et al . Assessment of the long-term possible radiological risk from the use of ceramic tiles in Malaysia . Nucl. Sci. Tech. 30 , 19 ( 2019 ). doi: 10.1007/s41365-019-0558-6 http://doi.org/10.1007/s41365-019-0558-6 Assessment of the long-term possible radiological risk from the use of ceramic tiles in Malaysia Y.H. Wang , Q. Li , L. Chen et al. , Simulation study of the dose and energy responses of FNTD personal neutron dosimetry . Nucl. Sci. Tech. 30 , 32 ( 2019 ). doi: 10.1007/s41365-019-0546-x http://doi.org/10.1007/s41365-019-0546-x Simulation study of the dose and energy responses of FNTD personal neutron dosimetry A.D. Wrixon , New ICRP recommendations . J. Radiol. Prot. 28 161 - 168 ( 2008 ). doi: 10.1088/0952-4746/28/2/r02 http://doi.org/10.1088/0952-4746/28/2/r02 New ICRP recommendations W.S. Snyder , Jr H.L., Fisher , M.R. Ford et al. , Estimates of absorbed fractions for monoenergetic photon sources uniformly distributed in various organs of a heterogeneous phantom . J. Nucl. Med. 10 Suppl 3 : 7 - 52 ( 1969 ). Estimates of absorbed fractions for monoenergetic photon sources uniformly distributed in various organs of a heterogeneous phantom All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Jun-Liang Chen, Su-Jun Yun, Tie-Kuang Dong, Zhong-Zhou Ren and Xiao-Ping Zhang. The first draft of the manuscript was written by Jun-Liang Chen and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript. 10.1007/s41365-022-00987-2 This work was supported by the National Natural Science Foundation of China (Nos.12035011, 11535004, 11905103, 11947211, 11975167, 11761161001, 11565010, 11961141003, 11805103, 11673075, 11303107, 11120101005, and 11235001), the National Key R&D Program of China (Nos. 2018YFA0404403 and 2016YFE0129300), the Science and Technology Development Fund of Macau (No. 008/2017/AFJ), the Fundamental Research Funds for the Central Universities (Nos. 22120210138 and 22120200101), and by the China Post- doctoral Science Foundation (Nos. 2019M660095 and 2020T130478). 2021-07-20 2021-11-19 2021-11-23 2022-01-25 2022-02-16 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2022 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-022-00987-2 http://dx.doi.org/10.1007/s41365-022-00987-2 . 2022 Nuclear Science and Techniques 01 33