Abstract
Future space travel to the earth’s moon or the planet Mars will likely lead to the selection of experienced International Space Station (ISS) or lunar crew persons for subsequent lunar or mars missions. The major risk for space travel is the galactic cosmic rays (GCR) risks of cancer, circulatory diseases and detriments in cognition. However large uncertainties in risk prediction occur due to the quantitative and qualitative differences in heavy ion microscopic energy deposition leading to differences in biological effects compared to low LET radiation. In addition, there are sparse radiobiology data and absence of epidemiology data for heavy ions and other high LET radiation. Non-targeted effects (NTEs) are found in radiobiology studies to increase the biological effectiveness of high let radiation at low dose for cancer related endpoints. In this paper the most recent version of the NASA Space Cancer Risk model (NSCR-2022) is used to predict mission risks while considering NTEs in solid cancer risk predictions. I discuss predictions of space radiation risks of cancer and circulatory disease mortality for US Whites and US Asian-Pacific Islander (API) populations for 6-month ISS, 80-day lunar missions, and combined ISS-lunar mission. Results predict NTE increase cancer risks by about ∼2.3 fold over a model that ignores NTEs. US API are predicted to have a lower cancer risks of about 30% compared to US Whites. Cancer risks are slightly less than additive for multiple missions, which is due to the decease of risk with age of exposure and the increased competition with background risks as radiation risks increase. The inclusion of circulatory risks increases mortality estimates about 25% and 37% for females and males, respectively in the model ignoring NTEs, and 20% and 30% when NTEs are assumed to modify solid cancer risk. The predictions made here for combined ISS and lunar missions suggest risks are within risk limit recommendations by the NCRP for such missions.
1. Introduction
The International Space Station (ISS) has been continuously occupied since the year 2000 by crew sizes of 3 to 7 astronauts whom typically enjoy 6-month missions. In the past astronaut flight careers have transcended more than one space program, including astronauts from the Mercury and Gemini programs on Apollo and Skylab missions, Space Transportation System (STS) (space shuttle) program astronauts participating on NASA-MIR program and ISS missions, and an astronaut participating in Gemini, Apollo, and STS missions. The large number of ISS veterans are likely candidates for lunar and possibly Mars or other deep space missions in the future. The high occupational radiation exposures for long duration missions presents the challenge of estimating the cumulative risk, including uncertainties in the predictions over an astronaut’s career of several long duration missions.
Space radiation risk assessments carry large uncertainties due to lack of knowledge on the radiobiology of high LET radiation such as heavy ions and secondary neutrons [1–3]. The major concern for space travel is galactic cosmic rays (GCR) risks of cancer, circulatory diseases and detriments in cognition. Solar particle events are also a concern however are readily mitigated by passive shielding [4], while shielding of GCR is difficult due to their high energies. High LET radiation produces both quantitative and qualitative differences in microscopic energy deposition leading to differences in biological effects compared to low LET radiation. In addition, there are sparse radiobiology data and the absence of epidemiology data for heavy ions and other high LET radiation such as secondary neutrons and low energy protons and helium ions. Non-targeted effects (NTEs) are found in low dose radiobiology studies to increase biological effectiveness for low doses of high LET radiation [5–10] for surrogate endpoints of cancer risk. Non-targeted effects include bystander effects and changes to the tissue micro-environment, including altered biochemical signaling and the induction of chronic inflammation.
In conventional radiation protection on Earth radiation weighting factors or LET dependent quality factors (QF) are used to estimate an effective dose [11,12]. However, this approach is severely limited for higher levels of exposure and risks for long duration (>3 months) space missions, which are well above a recommended effective dose limit of 50 mSv per year [13]. The large uncertainties in radiation quality effects for heavy ions and high LET radiation suggest uncertainties in risk estimates should be evaluated and reported in radiation protection programs [1,11,14–17]. The NSCR model [15,18–20] uses track structure concepts to form a radiation quality function with parameters fit to available proton, fission neutron, helium and heavy ion data for tumor induction in animals and surrogate endpoints in cell culture models related to cancer risk. NSCR evaluates the overall uncertainty in risk models using Monte-Carlo propagation of uncertainties to combined uncertainties in epidemiology data, organ doses and fluence spectra, and dose-rate and radiation quality effects. The NSCR model was extended to include the impact of NTE’s on radiation quality descriptors in recent work [21] and in this paper, we further investigate NTE role on risk predictions.
In this paper I apply the most recent version of NSCR, denoted NSCR-2022 [21–24] to report on risk predictions of cancer and circulatory diseases for average populations of different ages that would potentially participate in multiple long duration space missions, essentially ISS and lunar missions. Risk predictions vary with the background risks for the model population and here we make predictions for representative highest (Whites) and lowest (Asian-Pacific Islanders (API)) risk groups amongst US racial and ethnic groups based on a recent assessment [23]. Because NTE’s are shown to have a dose response that deviates from linearity, we investigate the possibility of deviation from linearity on a combined exposure using 35-year old persons at first mission and 40-year old at second mission in our study.
2. Cancer Risk Evaluation
The instantaneous cancer incidence or mortality rates, λI and λM, respectively, are modeled as functions of the tissue averaged absorbed dose DT, or dose-rate DTr, sex, age at exposure aE, and attained age a or latency L, which is the time after exposure L=a-aE. The λI (or λM) is a sum over rates for each tissue that contributes to cancer risk, λIT (or λMT). These dependencies vary for each cancer type that could be increased by radiation exposure. The total risk of exposure induced cancer (REIC) is calculated by folding the instantaneous radiation cancer incidence-rate with the probability of surviving to time t, which is given by the survival function S0(t) for the background population times the probability for radiation cancer death at previous time, summing over one or more space mission exposures, and then integrating over the remainder of a lifetime [17,20]: where z is the dummy integration variable. In Eq. (1), Nm is the number of missions (exposures), and for each exposure, j, there is a minimum latency of 5-years for solid cancers, 2-years for leukemia, and 5-years for circulatory diseases assumed. The upper limit of Eq. (1) is set at 100 year in the present calculations. Tissue specific REIC estimates are similar to Eq. (1) using the single term from λI of interest. The equation for risk of exposure induced death (REID) estimates is similar to Eq. (1) with the incidence rate replaced by the mortality rate (defined below) in the outer integral of Eq. (1).
The tissue-specific cancer incidence rate for an organ absorbed dose, DT, in multiplicative transfer models, denoted after adjustment for radiation quality and low dose and dose-rates through introduction of a scaling factor RQF [17,20] where λ0IT is the tissue-specific cancer incidence rate in the reference population, and where ERRT and is the tissue specific excess relative risk per Sievert, respectively. In previous reports we also considered a mixture model using a weighted average of the multiplicative and additive risk models. The choice of the multiplicative model in the present paper is discussed below. The tissue specific rates for cancer mortality λMT are modeled following the BEIR VII report [25] whereby the incidence rate of equation (2) is scaled by the age, sex, and tissue specific ratio of rates for mortality to incidence in the population under study:
Risks predictions for circulatory disease mortality were made in the same manner as a previous report [26], however using the recent meta-analysis of ERR for mortality after low LET radiation reported by Little et al. [27,28]. This approach uses ERR from mortality data after low LET radiation exposure. Circulatory disease incidence data following radiation exposure are sparse and preclude considerations herein. There is an uncertainty related to the dependence on possible treatments reducing circulatory disease mortality with time period, host country and other factors. Such factors are likely distinct for astronauts in recent or future time periods compared to the various historical epidemiology data considered by Little et al. [28]. A dose-rate reduction effectiveness factor is not applied because the meta-analysis results used are based on low dose (<0.5 Sv) or chronic exposures. For circulatory disease risks because a relative biological effectiveness (RBE) factor is distinct from the quality factor (QF), organ dose equivalents are estimated by the NCRP for deterministic effects (tissue reactions) [11,29].
3. Space Radiation Quality Descriptor
The RQF is estimated from RBE’s determined from low dose and dose-rate particle data relative to acute γ-ray exposures for doses of about 0.5–3 Gy, which we denote as RBEγAcute to distinguish from RBEmax, which is based on less accurate initial slope estimates for γ-rays. The basic approach is illustrated in Figure 1, which shows a stochastic track for an 56Fe ion with energy of 1 GeV/u [30]. The track is described as consisting of a so-called track core of ultra-high ionizations within 100 nm of the ions path, and a penumbra of ionizations by single energetic δ-rays extending many microns from the ions path. The quality function, RQF, considers these two contributions using a parametric approach motivated by the Katz track structure model [31,32]. However, allowing the penumbra term to follow the dose response estimated from γ-ray epidemiology studies, which is generally assumed to be linear over the relevant dose range for space missions. In addition, since we are scaling directly to the acute γ-ray responses, the dose-rate modifier is assumed to influenced by results from experimental models used in RBE determinations. The dose-rate modifier (assumed to be similar to estimate of a dose and dose-rate reduction effectiveness factor (DDREF)) adjusts the penumbra like but not the core like term of ultra-high ionizations. I note here that for chronic exposures in space it is unlikely that large regions of a genome or even individual cells would be traversed by two ion cores for time-periods of less than 1 month [33].
In this approach the scaling factor includes terms representing the penumbra-like and core-like of low and high ionization densities, respectively based on the following functions: where with the parametric function [20] where E is the particles kinetic energy per nucleon, L is the LET in keV/μm, Z is the particles charge number, Z* the effective charge number, and β the particles speed relative the speed of light. The LET for protons and helium in tissue are calculated using the National Institute of Standards [34] data. Effective charge is used to scale LET of heavy ions to protons [35].
The functional form of Eq. (5) reflects predictions from both radial dose and stochastic Monte-Carlo track structure descriptions of biological action cross sections dependence on ion charge and kinetic energy [36,37], which supports its usage in the parametric NSCR model to represent available experimental data. An ancillary condition is used to correlate the values of the parameter κ as a function of m as [20] where κ0 is value for the most likely value m=3. This constraint fixes the peak effectiveness with kinetic energy for each heavy ion charge group in the model to be consistent with results from experiments [15]. In Monte-Carlo uncertainty analysis conditional sampling is used where m is selected from a cumulative distribution function (CDF) followed by selection of κ(m), which is then distributed with a normal distribution with standard deviation (SD) shown in Table 1.
The three model parameters (Σ0/αγ, κ and m) in Eq. (4) are fit to radiobiology data for tumors in mice or surrogate cancer endpoints as described previously [20,38]. Distinct parameters are used for estimating solid cancer and leukemia risks based on estimates of smaller RBEs for acute myeloid leukemia and thymic lymphoma in mice compared to those found for solid cancers [20]. In addition, it is noted that RBE’s for male hepatocellular carcinoma in mice are much larger than for female mice and for most other solid cancers for both male and female mice [20].
The space radiation quality model corresponds to a pseudo-biological action cross section in units of µm2 given by,
Figure 2 compares RQF from Eq. (4) for solid cancers to the ICRP LET dependent quality factor divided by DDREF=2 as a function of kinetic energy for the major GCR components of protons, 16O, 28Si and 56Fe. The most important differences are at the lowest energies where protons are shown to be more effective and heavy ions less effective in NSCR-2022 compared to ICRP QF [12], and at the highest energies for heavy ions where the large penumbra of δ-rays contribution to microscopy energy deposition is predicted to reduce effectiveness, which is not described in the ICRP function. The ICRP function does not represent the branching of RBE with charge number for ions of identical LET and instead uses LET alone. In addition, the ICRP provides no uncertainty assessment for the QF, while uncertainty assessment is a main focus of NSCR model.
4. Non-Targeted Effects Quality Factor Model
In the NTE model we assume the targeted effects (TE) contribution is valid with a linear response to the lowest dose or fluence considered, while an additional NTE contribution occurs [21]. where Dcr is found as the dose (or fluence) where the NTE and TE are equivalent that is given by
This corresponds to a pseudo-biological action cross section is given by, where F is the particle fluence (in units of 1/µm2) and the η function represents the NTE contribution, which is parameterized as a function of x=Z*2/β2as:
In Eq. (11) the area, Abys, reflects the region of bystander cells surrounding a cell traversed directly by a particle that receives an oncogenic signal. The parameters η0/αγ and η1 are estimated from low dose radiobiology experiments for mouse Harderian gland tumor induction [8,39] and chromosomal aberrations [9]. The second factor on the right-hand side of Eq. (11) describes the “turning on” of NTE at very low doses.
The parameter Abys is difficult to estimate from reported heavy ion experiments because of lack of low dose responses (<0.01 Gy) and it is value is likely correlated with estimates of η0. However, several experiments with α-particles provide data to make estimates [40,41]. Figure 3 (panels A, C, and E) shows results of fitting the model equations to experiments with α-particle induced chromosomal and chromatid aberration in wild-type CHO cells [41], and neoplastic transformation of C3H10T1/2 [40]. In Figure 3 (panels B, D, and F) these same data are shown with the y-axis divided by the mean number of ions that traversed the cell versus the number of ions that traversed. The mean number of ions traversed is estimated as the ion fluence times the cross-sectional area of the cell nucleus (i.e. H=FxAcell), which is 62.7 and 250 μm2 for CHO and C3H10T1/2 cells, respectively. Using this scaling a linear response would appear as a horizontal line. There is a clear increase above linearity as the number of ion traversals decreases below unity that is reflective of a bystander effect. The magnitude of the increase is impacted by the control values found in the experiment. Here we note a null result is found for the controls and the lowest dose tested (0.0017 Gy) for chromosomal aberrations [41], while a non-zero result is found for chromatid breaks [41] and the transformation experiment [40]. The parameter Abys in Eq. (11) is estimated to vary between 900 for CHO cells and 1600 μm2 for the larger C3H10T1/2 cells. We note that larger distances for bystander effects are suggested in experiments using a micro-beam to target specific areas within or outside of cell nuclei [42].
In Figure 4 we show a similar scaled response for Harderian gland tumors where the cell area is reported as 24 μm2 indicating a nuclear radius of 2.76 μm. A similar pattern is found for tumor induction as a function of heavy ion traversals per cell as that for the α-particle endpoints in Figure 3 whereby a rapid increase occurs as the number of ion traversals decreases below unity. In Figure 4 we also show the experimental data with background subtracted, (P-P0)/Hits which suggest the increasing ratio (P/H) with decreasing number of traversals is not due to the approach towards a non-zero background. For the tumor prevalence the impact of cell survival is considered, however over the low dose range (<0.2 Gy) increases in target cell survival would only increase P/H by about 20% based on cell survival experiments with heavy ions [43]. Other possibilities include radiation induced senescence or modulation of DNA damage responses reducing cancer susceptibility at higher doses to a larger extent than at lower doses. These observations which suggest bystander or other NTEs are not surprising as very few cells will survive a direct traversal by a heavy ion with LET > 100 keV/ μm of a cell nucleus. This is illustrated in Figure 6, which shows amorphous track model [44,45] and Figure 7 for stochastic radiation track structures of heavy ions of high LET [30]. In Figure 7 the 56Fe nuclei at LET=450 keV/μm is similar to the track of the Nb beam used in the Alpen et al. experiments [39]. Clearly a large portion of a cell nucleus receives an extremely large number of ionizations for direct cell nucleus traversals of the very high LET ions (>200 keV/μm) used in the Harderian Gland experiments of Alpen et al. [39] leaving very few surviving cells. In addition, the dose-response for Harderian gland tumor induction by γ-rays does not support a high tumor prevalence from cells traversed from δ-rays produced by passing heavy ions when the ion does not traverse the cell nucleus.
5. Space Radiation Exposures Organ Exposures
GCR exposures include primary and secondary H, He and HZE particles, and secondary neutrons, mesons, electrons, and γ-rays produced over a wide energy range. We used the HZE particle transport computer code (HZETRN) [46] with quantum fragmentation model nuclear interaction cross sections [47] and GCR environmental model [15] to estimate particle energy spectra for particle type j, φj(Z,E) as described previously [15,20]. For the radiation quality model described above, a mixed-field pseudo-action cross section is formed by weighting the particle flux spectra, φj(E) for particle species, j, contributing to GCR exposure evaluated with the HZETRN code with the pseudo-biological action cross section for mono-energetic particles and summing over all particles and kinetic energies:
Equations for the mixed-field pseudo-action cross section in the NTE model as folded with particle specific energy spectra as [21]:
6. Parameter Values and Uncertainty Analysis
For the various parameters that enter into the model equations for REIC and REID, probability distribution functions (PDFs) are formulated based on experimental data to represent plausible ranges of values. The physics models for dose and dose equivalent agree with spaceflight data in low Earth orbit [13,15], in transit to Mars [48] and on the Mars surface [49] to within 15% in most comparisons as illustrated in Figure 5. Similar good agreement was shown for comparisons organ doses to a phantom torso experiment flown on an STS mission [13,15]. Albedo particles on the lunar or Mars surface are estimated to increase the organ dose equivalents by 10 to 20% largely due to neutrons and protons, and we use a mean increase of 15% [49]. The physics uncertainties are estimated using a normal distribution with standard deviations of 0.35 for Z<4 ions and 0.25 for Z>4 ions. Larger deviations may occur for specific energy values of fluence spectra, however are constrained by overall uncertainty in integral quantities such as dose and dose equivalent.
Models for tissue specific ERR for cancer incidence are taking directly from the published reports from the Radiation Effects Research Foundation (RERF) for site-specific cancer incidence [50–59]. Statistical and dosimetry errors in the model parameters from these reports are used directly in our Monte-Carlo uncertainty analysis. Because additive risk models were not reported for several of the tissues in these reports [50–59] and to be consistent with the multiplicative model use for circulatory risk estimates [27] we report on the multiplicative model results in this paper.
A Bayesian approach is used to combine the DDREF estimate from the atomic-bomb survivor data with estimates from mouse tumor induction [18,20]. We studied the BEIR VII estimates and one by Hoel [60] and combined these with mouse data for solid cancers in mice and found small differences on their impact to GCR risk predictions due to the dominance of the core like terms [20]. For NSCR-2022 we use the Bayesian derived PDF with a mean DDREF=2.
Solid tumor data in mice are used directly to model the value and PDF for the parameter Ʃ0/αγ, where the PDF is represented by a three-parameter logistic function [20]. We use values derived excluding the Harderian Gland (HG) tumors as this organ does not appear in humans, and hepatocellular carcinoma in male mice because it was found to be an outlier from values for other cancers. HG tumor response data are used in estimation of other RQF parameters because it is the most detailed data set for radiation quality dependence for tumor induction [8,39]. The parameter values and uncertainties for the parameters m and κ are estimated from experiments on detailed radiation quality dependences with 5 ion’s or more in cell culture for chromosomal aberrations, neoplastic transformation, and HPRT mutations and the HG mouse tumor induction [20,21,38]. Using the constraint of Eq. (6) shows that the value of m has small influence on results for the range of 2 to 4. The value of the κ parameter is set at a higher value for light ions (κ=1000) based on experiments with fission neutrons and low energy light ions (<10 MeV/u). This difference in the κ parameter value between light and heavy ions is attributed to the influence of the varying effectiveness of δ-rays of different energies with these energies increasing with ion velocity. Electrons with energies below ∼20 keV are observed to have increased biological effectiveness [61,62], which plays a larger role for lower energy ions (<10 MeV/u) compared to relativistic energy ions. This energy region makes a large contribution for tracks of low energy light ions, which are abundant from slowing down of GCR primaries and secondaries from neutrons and other ions. In contrast, heavy ions have low abundance and short ranges below kinetic energies of 10 MeV/u and therefore are not a major contributor to space radiation risks.
7. Predictions for Space Missions
We considered astronauts of age 35-years at first missions and age 40-years at second mission. In a recent publication [23] it was shown that US Whites have the highest risks due to their higher background cancer rates and life-table, while US Blacks having a modestly lower radiation risk compared to US Whites. US Asian-Pacific Islanders (API) were found to be at the lowest radiation risks in the multiplicative model, with US Hispanics slightly higher. We therefore studied US Whites and US API groups as representative of the most and least sensitive model populations. GCR exposures are slowly varying for shielding below 100 g/cm2 of aluminum and we used 30 g/cm2 for an ISS exposure and 20 g/cm2 for lunar segments. We assumed 180-day ISS missions and considered GCR and trapped protons as described previously [13]. For possible lunar missions we used a slow Earth to moon transit time of 10 days and a 60-day surface time based on reported scenarios [63]. Average solar cycle conditions are used for ISS segment while an average solar minimum is used for the lunar segment.
Table 2 shows predictions of REIC and Table 3 for REID with the inclusion of circulatory disease risk predictions. Results predict NTE increase REIC by about ∼2.3 fold over the model that ignored NTEs. The predictions with NTEs show the lack of REID <3% with a 95% confidence interval for multiple ISS and lunar missions. Comparison of US Whites to US API show a decrease in cancer risks of about 30% compared to US whites. Cancer risks for multiple missions are slightly less than additive, which is due to declining risk with age at exposure and increased competition with background risks as radiation risks increase. Females have higher risks compared to males largely due to breast and ovarian cancers, and higher risk for lung cancer as predicted by low LET epidemiology data. For males the use of a higher RBE for male liver cancers based on mouse studies would narrow this difference to some extent as liver cancers makeup about ∼10% of the overall cancer risk. This difference in high LET biological effectiveness for liver cancer between males and females requires further study.
The inclusion of circulatory risks increases mortality estimates about 25% and 37% for females and males, respectively in the model ignoring NTEs, and 20% and 30% when NTEs are assumed to modify RBE’s for solid cancers. Circulatory disease risks for US Whites and US API are predicted to be nearly the same. At this time RBE’s for circulatory diseases or their possible surrogate endpoints are extremely sparse, which leads to an uncertainty that is difficult to quantify. Using the NCRP recommended RBE model for non-cancer risks [11,29] suggests the proportion of circulatory disease risk in the overall risk is lower for space radiation compared to low LET exposures on Earth. A healthy worker effect for circulator risk found for astronauts would also suggest a lower proportion [15]. In addition, we note that differences in treatment of circulatory diseases with time period and exposed group leads to uncertainty when applying historical epidemiology data to astronauts in current or future missions, while the impact of such differences for cancer risk has been estimated using the BEIR VII [25] recommended approach shown in Eq. (3).
8. Summary
Cancer risk is the dominant risk established for long-term space missions. The predictions made here for combined ISS and lunar missions suggest risks are within recommendations of the NCRP for such missions [11,17] whose value for radiation protection of astronauts was described in a recent report [64]. Participation in a 3rd long duration mission could possibly exceed the NCRP recommendations for ISS and lunar operations, however would be dependent on solar cycle, time between missions, age, sex, and possibly other factors. An important issue for cancer risk assessment is a possible deviation from a linear response for high LET radiation at low doses due to NTEs, which impacts RBE estimates. NTEs include bystander effects where cells traversed by heavy ions transmit oncogenic signals to nearby cells, genomic instability in the progeny of irradiated cells, and tissue microenvironment changes related to cancer development [5–10].
Mechanistic studies that have used low doses (<0.1 Gy), where less than one particle traverses a cell nucleus, suggests NTE dominate low dose risk from high LET radiation. A concern is the potential limitation of experimental models to quantify NTEs, however we note these are similar to the models used in assigning radiation quality factors in the past. The effects of NTEs are shown to significantly increase risk predictions for ISS and lunar missions in the present paper and for a Mars mission in a previous report [21]. A major concern is that most radiobiology studies continue to employ heavy ion doses corresponding to more than one particle per cell nucleus (>∼0.1 Gy) that will not be sensitive to NTEs, while GCR exposures occur where NTE’s are predicted to dominate risks. Organ absorbed doses from heavy ions are less than 0.1 Gy for a Mars mission, and less than 0.05 Gy for multiple ISS and lunar mission scenarios. Clearly the development of sensitive experimental models and data sets at low doses of heavy ions (<0.1 Gy) is needed to improve risk assessments. In addition, it is recommended the use of relativistic super heavy ions (A>58) such as high-energy Nb, La or Pb beams at low doses should be investigated as an optimal tool to investigate NTEs for space radiobiology applications.
Data Availability
All data produced in the present work are contained in the manuscrip