{'fraction_non_alphanumeric': 0.036617451999391155, 'fraction_numerical': 0.0277675096218661, 'mean_word_length': 4.644329896907217, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 10, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Revisiting the cosmic-ray induced Venusian radiation dose in the context of habitability

The Atmospheric Radiation Interaction Simulator (AtRIS) was used to model the altitude-dependent Venusian absorbed dose and the Venusian dose equivalent. For the first time, we modeled the dose rates for different shape-, size-, and composition-mimicking detectors (phantoms): a CO$_2$-based phantom, a water-based microbial cell, and a phantom mimicking human tissue. Based on a new model approach, we give a reliable estimate of the altitude-dependent Venusian radiation dose in water-based microorganisms here for the first time. These microorganisms are representative of known terrestrial life. We also present a detailed analysis of the influence of the strongest ground-level enhancements measured at the Earth's surface, and of the impact of two historic extreme solar events on the Venusian radiation dose. Our study shows that because a phantom based on Venusian air was used, and because furthermore, the quality factors of different radiation types were not taken into account, previous model efforts have underestimated the radiation hazard for any putative Venusian cloud-based life by up to a factor of five. However, because we furthermore show that even the strongest events would not have had a hazardous effect on putative microorganisms within the potentially habitable zone (51 km - 62 km), these differences may play only a minor role.


Introduction
As discussed in the literature, planetary atmospheres within our Solar System are exposed to a harsh radiation environment consisting of, for example, the solar wind, ultraviolet-and X-ray radiation, and cosmic rays (CRs) of Galactic or solar origin. The atmospheres may also include magnetospheric particle precipitation.
According to Reames (1999), solar energetic particles (SEPs) are accelerated in solar flares, coronal mass ejections, or at interplanetary shocks. Galactic cosmic rays (GCRs), however, are most likely produced by diffusive shock acceleration at the shocks of supernova remnants (see, e.g., Hillas 2005;Büsching et al. 2005). Thus, because of their different origins and acceleration mechanisms, the energy distribution of these two components differs significantly. While SEPs exceed the GCR flux by four orders of magnitude at lower energies (several keV), a steep decrease at energies above several MeV occurs (see, e.g., Herbst et al. 2015). Only in case of extreme SEP events that are detectable on the Earth's surface, socalled ground-level enhancements (GLEs), the spectrum exceeds energies of hundreds of MeV to a few GeV. In addition, differences in shape and composition of these GLE events have been reported (Schmelz et al. 2012). Although GCRs with energies of up to several TeV enter the heliosphere omnidirectionally, the heliospheric magnetic field (HMF), which is carried outward with the solar wind, modulates the GCR flux over a solar cycle.
When they are present, CRs encounter yet another magnetic filter in the form of the planetary magnetic field before they impinge on planetary atmospheres. There, CRs can be deflected by the Lorentz force, which leads further modulates the CR flux. However, CR interactions with atmospheric neutrals are the primary driver of the planetary for nearly two billion years (see, e.g., Donahue et al. 1982;Donahue & Hodges 1992). This could have been long enough for life (as we know it from Earth) to have been evolved. As suggested by Grinspoon & Bullock (2013), with increasing atmospheric temperature, the surface water may have evaporated, which may have forced any putative microorganisms to migrate into the still habitable Venusian cloud layer. Consequently, these microorganisms may have been able to adapt to the limited-water and low-pH environments of the Venusian clouds (see, e.g., Schulze-Makuch et al. 2004;Grinspoon & Bullock 2013). According to Limaye et al. (2018), the atmospheric movement of life-essential nutrients, water, and organics is, among others, most likely regulated by diurnal cycles and atmospheric storms. The VeGa 1 and VeGa 2 balloon measurements, for example, showed that the typical vertical motion within the Venusian atmosphere is in the order of 1 to 2 m/s (Blamont et al. 1986). Thus, Venusian cloud particles are assumed to have atmospheric residence times comparable with the two-to three-month terrestrial Hadley circulation. Based on recent measurements of stationary gravity waves at the top of the Venusian cloud layer, which most likely are the result of the surface topography (see Fukuhara et al. 2017), it becomes evident that vertical motions within the cloud layers are also possible. Because of the stable Venusian cloud lapse rates, therefore, airborne (nutrient-rich) particles could be sustained within the cloud layer for extended time intervals (Limaye et al. 2018). Dartnell et al. (2015) investigated whether GCRs and SEPs, including extreme SEP events from the historical record, could represent a constraint on the possibility of cloud-based life on Venus. This study, however, did not include modeling on water-based shape-, size-, and composition-mimicking detectors (phantoms), which more realistically describes the radiation effects in living organisms.
In this study, we revisit the altitude-dependent Venusian radiation environment by studying the cosmic-ray induced atmospheric absorbed and effective dose rates. To study the impact of CRs from the astrobiological point of view, for the first time we study the radiation dose that different phantoms, in particular, a microbial cell and a phantom mimicking human tissue (based on the International Commission on Radiation Protection (ICRP) phantom (see, e.g., ICRP 2007) would suffer within the Venusian atmosphere. The cosmic-ray induced radiation dose is modeled with the newly developed Geant4-based Atmospheric Radiation Interaction Simulator (AtRIS, Banjac et al. 2019b), a 3D simulation code developed to model the radiation dose of exoplanets, which varies from hot Jupiters to Earth-like planets. In addition to the GCR-induced atmospheric radiation due to primary particles with Z = 1 − 28, we investigate the influence of 67 recent GLE events that occurred on Earth between 1942 and 2017 (see, e.g., Raukunen et al. 2018). From historical records of so-called cosmogenic radionuclides, it is moreover known that much stronger SEP events occurred in the past (see, e.g., Mekhaldi et al. 2015). Two events therefore are of particular interest: the Carrington event, and the AD775 event.
Based on assumptions about the GLE event spectrum (see Dartnell et al. 2015), we also investigate the influence of such intense events on the Venusian atmospheric radiation and their influence on the habitability within the Venusian HZ in the cloud layer.

Atmospheric Radiation Interaction Simulator (AtRIS)
AtRIS, a Geant4-based (Monte Carlo) particle transport code, was recently developed by Banjac et al. (2019b) in order to study, among others, the propagation and interaction of energetic particles within diverse (exo)planetary atmospheres and surfaces. As particular input, the user has to define an interface for the planetary atmospheric model. By now, interfaces for the terrestrial NRLMSISE-00 (Picone et al. 2002) model, the Mars Climate Database (MCD), and the dayside model of the Venus International Reference Atmosphere (VIRA, see, e.g., Kliore et al. 1992) have been implemented. Thereby, atmospheric altitude-dependent ion and electron pair production rates, secondary particle distributions, and absorbed dose rates and dose equivalent rates can be calculated. Thus, it is rather easy to implement GCR and SEP event spectra. Furthermore, the user can decide on any hadronic and electromagnetic physics lists provided through the standard Geant4 message and compounded physics list naming scheme (Allison et al. 2016;Geant4_Collaboration 2018).
The main features of AtRIS are i) the planet specification format (PSF), ii) the atmospheric response matrices (ARMs), which quantify the relationship between primary energy and altitude-dependent ionization, absorbed dose rate and dose equivalent, and iii) the spectrum-folding procedure used to calculate net quantities such as the electron-ion pair production rate by implementing a convolution of a spectrum and the ARM. A detailed description of the main features is given in Banjac et al. (2019b). At this point, AtRIS has successfully been validated for the terrestrial, the Martian, and the Venusian atmosphere (see Banjac et al. 2019b,a;Guo et al. 2019;Herbst et al. 2019a, respectively). Herbst et al. (2019a) calculated the atmospheric ionization within the Venusian atmosphere. To do this, the mean energy required to produce an ion-electron pair had to be taken into account. However, different values exist in the literature: while Borucki et al. (1982) gave a value of 33.5 eV, Simon Wedlund et al. (2011) suggested that a value of 28.7 ± 4.3 eV might be more reasonable in CO 2 -dominated atmospheres, which leads to differences in the ionization rates in the order of 15% (see Herbst et al. 2019a). These differences are also expected when the atmospheric radiation dose is modeled. However, for comparative reasons, all results discussed in this study are based on the "traditional" value proposed by Borucki et al. (1982). Nevertheless, as previously discussed in Nordheim et al. (2015) and Herbst et al. (2019a), dividing the presented ionization and radiation profiles by a factor of 1.17 results in values corresponding to the use of the ionization potential by Simon Wedlund et al. (2011).
One of the main AtRIS outputs is the so-called ARMs, which reflect the altitude-dependent ionization, absorbed dose, and dose rate equivalent as a function of the primary particle energy. According to this, the absorbed dose of a specific phantom is derived using the following approach: 1. Definition of the particle type j and its kinetic energy E i . 2. Using the AtRIS-specific precalculated relative ionization efficiency I R, j , which is given by where E d represents the average ionization energy a particle causes within a well-defined phantom. 3. Calculation of the average absorbed dose a phantom experiences when particles of a given energy pass through it.
We note that the absorbed dose corresponds to the ionization energy per mass of the phantom m(r) = ρ · 4 3 π · r 3 . Thus, . (2) 4. If required, multiplication with the geometric factor g = πA(r) = 4π 2 · r 2 (?), which is proportional to the particle flux going through the phantom, represents the normalized average dose D j,n , We like to point out that the shape and size of different phantoms can vary. Different shapes, however, account for different geometric factors: g = π · 2a 2 for a slab detector, g = π · 2r 2 for a disk detector, and a spherical detector, for example, a planetary atmospheric layer, has a geometrical factor of π · 4πr 2 (see, e.g., Banjac et al. 2019a). By placing phantoms of different sizes and shapes within the same radiation field, the number of particles that can penetrate the phantom and thereby deposit energy within, therefore decreases as the square of the phantom radius. 5. Convolution with the primary particle spectrum and summation over all energy bins and particles.
In this way, the dosimetric quantities of a reference phantom can be calculated without the need of further time-consuming computations. This is helpful for any parametric study like the one presented in Sec. 4.7.1, where the altitude-dependent absorbed dose rate profiles for 67 GLE events are computed.
The Venusian absorbed dose rates have currently only been computed for a Venusian air phantom (see, e.g. Dartnell et al. 2015). However, life as we know it is based on water rather than CO 2 . According to climate modeling efforts of Way et al. (2016), Venus may early in its planetary evolution have had a liquid-water ocean and surface temperatures that may have allowed for habitable conditions. This and the fact that with increasing temperatures the Venusian water content may have evaporated and, therewith, microorganisms may have been transported in the cloud layer (Grinspoon & Bullock 2013) lead us to study the absorbed dose in H 2 O-based phantoms to better describe the impact of CRs in the Venusian atmosphere on life as we know it (see Sec. 4.3).
We also further calculate the dose equivalent H, which is more relevant for estimating the radiation hazard to H 2 Obased (human-like) phantoms because it takes into account the varying biological effectiveness of different radiation types. With this, the dose equivalent H is calculated by where w R, j is the radiation weighting factor defined by ICRP 60 (1991) and most recently modified in ICRP 103 (2007). Thereby, w R, j account for the different biological effectiveness of different types of radiation and for all particles w R, j ≥ 1, as listed in , for E > 50 MeV , where E is the kinetic energy of the neutron in MeV. For all other particles, w R is equal to 1. Because w R is dimensionless, the dose equivalent H has the same dimension as the absorbed dose. A more advanced approximation of w R, j exists in the literature. However, this approximation relies on the determination of the linear energy transfer (LET), which is difficult and somewhat complex. Although the approximation can correct for overestimating the influence of relativistic heavy minimally ionizing particles (MIPS), for example, 100 GeV α particles, it is not suitable for cases in which the LET spectrum cannot be measured. Moreovoer, at column depths below 200 g/cm 2 (corresponding to an altitude of ∼ 62.5 km within the Venusian atmosphere), the flux of relativistic heavy particles is relatively low, so that Eq. (5) is a valid approximation. Furthermore, we note that these weighting factors are tailored explicitly for biological studies with a focus on human radiation health. Thus, they may not be directly applicable to astrobiology studies with a focus on the survival of microorganisms, for example.

Simulation setup
As discussed in Herbst et al. (2019a), the Venusian environment was modeled for spherical geometry using a core with a radius of ∼ 6052 km, a soil composition of 50% Si, 40% O, and 10% Fe, and a crust (soil) sheet with a thickness of 100 m. The upper layer of the Venusian atmosphere was set at an altitude of 150 km, and the atmosphere was divided into 500 m thick layers. Furthermore, the same atmospheric composition, source altitude, and particle input spectra were used.
The atmospheric setup was thus based on the VIRA model (see Kliore et al. 1992), which uses the day-side atmospheric parameter set by Seiff et al. (1985) at low-to mid-atmospheric layers (surface to 100 km) for low latitudes (φ < 30 • ) and those proposed by Keating et al. (1985) for altitudes up to 150 km at low latitudes (φ < 16 • ). The atmospheric composition was chosen to consist of 96.5% CO 2 and 3.5% N 2 .
To simulate the hadronic and electromagnetic interactions, and thus the evolution of secondary particle cascades as well as the atmospheric radiation dose, Dartnell et al. (2015) based their studies on the QGSP_BIC_HP model using the quark-gluon-string (QGS) model for high-energy particle interactions and the binary cascade model (BIC) for particle energies below 10 GeV. As discussed in Herbst et al. (2019a), the more accurate model to use is the FTFP_BERT_HP model, which is explicitly recommended for radiation protection and shielding applications. Therefore, our computations use the Bertini-style cascade for hadrons with energies below 5 GeV and the Fritiof (FTF) model (see, e.g., Nilsson-Almqvist & Stenlund 1987) for particle interactions of mesons, nucleons, and hyperons with energies between 3 GeV and 100 TeV. In addition, the standard em constructor (see, e.g., Allison et al. 2016) was used to model electromagnetic interactions. We note that differences due to the use of these two different hadronic interaction models only occur at altitudes above 120 km and below 5 km (see Fig. 4 in Herbst et al. 2019a, for more details). Therefore, in the context of the habitability of the Venusian cloud layer, the computed physical radiation doses of our study are comparable to those of Dartnell et al. (2015).
As discussed in Herbst et al. (2019a), we used the CREME2009 model 1 (see, e.g. Tylka et al. 1997) as input for the time-dependent GCR flux at the top of the Venusian atmosphere. This model provides the differential intensities for all particles from protons (Z = 1) to nickel (Z = 28) for energies between 1 MeV/nuc up to 100 GeV/nuc derived for a solar distance of 1 AU. In this study, we investigate the atmospheric radiation dose during quiet solar minimum and maximum conditions. As discussed in Herbst et al. (2019a), the high-energy tail of the particle spectrum (in particular up to 10 TeV/nuc) may have a much higher impact on the atmospheric radiation dose as previously thought. Therefore, an element-dependent power-law extension of the spectra was applied. Although Venus is about 30% closer to the Sun than Earth, a rescaling of the modeled GCR spectra is not necessary because the radial gradients of GCR particles in the inner Solar System are negligible (see, e.g., Heber et al. 1996;Morales-Olivares & Caballero-Lopez 2010;Gieseler & Heber 2016).
Furthermore, the spectra of particularly strong SEP events are provided by the CREME2009 model. Because of their acceleration at shock fronts of solar flares (SFs) and coronal mass ejections (CMEs), SEPs have energies of only up to some hundred MeV and reach the Solar System planets anisotropically. In contrast to GCRs, the flux of SEP events strongly depends on the orbital distance. Thus, a 1/R 2 scaling was applied to the spectra. Following the investigations by Herbst et al. (2019a), we studied the influence of 67 GLE events measured within the neutron monitor era (Raukunen et al. 2018). The strongest GLE event ever measured in situ occurred on 23 February 1956 (hereafter GLE05). However, based on observations of the so-called Carrington event, and from the cosmogenic radionuclide records of 10 Be, 14 C, and 36 Cl covering the past ∼ 10,000 years, today we now know that even much stronger SEP events have occurred throughout the Sun's history, in particular around BC 660, AD774-775, and AD993-994 (see O'Hare et al. 2019;Mekhaldi et al. 2015;Miyake et al. 2012, respectively). Based on model efforts, Kovaltsov & Usoskin (2014) and Herbst et al. (2015) found that the event of AD775 was most likely well over one order of magnitude stronger than GLE05.
We here investigate the impact of different GCRs nuclei, their upper energy limit, and the impact of solar modulation on the Venusian atmospheric radiation dose, and we highlight the differences between absorbed dose and dose equivalent. We also study the influence of 67 terrestrial GLE events on both quantities and discuss the impact of potential super-events like the Carrington-and the AD774/775 event. However, because the spectral shape and the upper event energy of the latter is not known, the influence of the two historical events was modeled based on measured event spectra of other historical events.

Contribution of different GCR nuclei on the atmospheric ionization
As shown in Herbst et al. (2019a), the atmosphere primarily is ionized through GCR protons (Z=1) and helium nuclei (Z=2). The left panel of Fig. 2 shows a more detailed study of the particle-dependent GCR-induced Venusian ionization of Z = 1 (solid red line), Z = 2 (solid blue line), and Z = 3-28 (solid green line) particles. Information on the relative contribution of the different GCR species is also given as dashed lines on the left-hand side. The corresponding values are given in the upper x-axis. Primary protons are the main source of the CR-induced ionization at all altitudes, with increasing significance toward lower altitudes (up to 80%). An anticorrelated behavior can be found for Z= 3-28 primary Table 2. Modeled absorbed dose-rate values at the lower and upper limit of the HZ (51 km and 62 km, respectively) proposed by Dartnell et al. (2015). Here, we list the values corresponding to the three scenarios.

Contribution of different GCR nuclei to the absorbed dose and the dose equivalent
By applying the approach discussed in Sec. 2, we modeled the altitude-dependent absorbed dose rates and the equivalent dose rates. We show them for example in Fig. 2. Although the atmospheric ionization (left panel) and the absorbed dose rates (right panel) are closely related, considerable differences in the two altitude-dependent profiles are clear. The total absorbed dose rate (solid black line) and the relative contributions of protons (in red), alpha particles (in blue), and Z>2 nuclei (in green) are constant in between 140 km and 85 km. Here, protons account for 60% of the total absorbed dose, while alpha particles and Z>2 nuclei account for about 20% each. Below 85 km, this picture changes. Here, the influence of Z>2 nuclei rapidly decreases, while the influence of protons increases up to 80%. However, down to 10 km, the contribution of alpha particles is always within 20%, but it increases to 40 % at altitudes below. Within the potential Venusian habitable zone (highlighted in green, 51 km -62 km), for example, protons contribute more than 70% to the total absorbed dose rate, while Z> 2 particles contribute less than 10% and can therefore be neglected.
In order to give a more realistic estimate of the radiation hazard to potential microorganisms in the Venusian clouds, it is more convenient to investigate the dose equivalent that accounts for the biological effectiveness of different radiation types (see Section 4.6 for more details), as shown in the lower panel of Fig. 2. Here the impact of primary particles with Z>2 between 140 km and 85 km becomes evident because with a contribution of more than 35% to the total equivalent dose, they are the dominant source. However, although they have a much higher biological effectiveness than the primary protons and alpha particles, they do not deliver the majority of the equivalent dose at all altitudes. Below 85 km, their contribution clearly rapidly decreases below 10 % within the potential HZ. In order to model the absorbed dose rates and the dose equivalent rates from the upper cloud layer down, it is therefore sufficient to take only the interactions of GCR protons and alpha particles into account.

Investigating absorbed dose rates in different phantoms
A commonly used approach to determine the (Venusian) atmospheric radiation dose is to divide the calculated ionelectron production rates by the density of the ambient planetary atmosphere (see, e.g., Dartnell et al. 2015). In this way, Dartnell et al. (2015) analyzed the radiation constraints on a potential Venusian aerial biosphere. They only considered the radiation dose on the gas, however, and therefore did not explicitly study the radiation impact on biological systems such as cells, microorganisms, or organic tissue. To study the exposure of life as we know it, phantoms that better reflect biological systems can therefore be used. The ICRU proposed the so-called ICRU sphere, a tissue-equivalent sphere with a radius of 15 cm and a density of 1 g/cm 3 . With a composition of 76.2% oxygen, 11.1% carbon, 10.1% hydrogen, and 2.6% nitrogen, the ICRU phantom is designed to represent and reflect life as we know it from Earth (McNair 1981).
However, because it is not clear if (possible) life on Venus has formed and evolved in comparable prebiotic conditions, using this phantom to model the radiation hazard of Venus may lead to false conclusions. Because life is in general expected to be water based (see, e.g., Ball 2017), a more realistic approach needs to be taken in order to give a more advanced estimate of the altitude-dependent radiation hazard to potential Venusian life. In this section, we therefore investigate the impact of CRs on different phantoms: phantom 1: based on Venusian air (as discussed in Dartnell et al. 2015) phantom 2: an ICRU-sphere equivalent consisting of 100% liquid water with a radius of 15 cm phantom 3: a water-based phantom corresponding to the size of a microbial cell of spherical shape and a radius equal to 50 µm The GCR-induced absorbed dose rate profiles for phantom 1 (in blue), phantom 2 (in black), and phantom 3 (in red) are shown in Fig. 3. Different approaches lead to different estimates of the altitude-dependent radiation dose. Consequently, the modeled absorbed dose rates depend on the size and the composition of the assumed phantom. While the absorbed dose rates above 100 km are much higher for phantom 1, the dose rates below 70 km are up to 60% lower than the absorbed dose phantom 2 would suffer within the Venusian atmosphere (see the dashed lines for the relative differences with respect to phantom 2). A comparison of the absorbed dose rate values at the upper and lower limit of the HZ is given in Table 2. Within the HZ, the absorbed dose rates of phantom 2 (phantom 3) are up to four (three) times higher than   For the dose in the air plot (blue line), the relative difference includes dotted lines that correspond to maximum and minimum differences, as described in detail in the text.
those of the Venusian air (phantom 1). We note that, due to its composition, phantom 2 creates an appreciable amount of self-irradiation, for example, due to the production of secondary particles within the 15 cm water-sphere, which further leads to an additional energy deposition. However, compared to phantom 3, the 50µm water-based microbial cell, this effect only accounts for differences in the order of a factor of 1.7 and does not explain the differences, in particular, those between phantom 1 and phantom 3. Consequently, the previously introduced approach by Dartnell et al. (2015) systematically underestimates the absorbed dose of potential water-based microorganisms below 70 km, and in particular, within the potential Venusian HZ, by up to 45%.
We note that the absorbed dose rates of phantom 2 have a systematical error at altitudes above ∼ 100 km that results from the use of the VIRA model with the two atmospheric models by Seiff et al. (1985) and Keating et al. (1985), which are discontinuous at their interface. As discussed in Section 2, the absorbed dose rates of phantom 2 were computed by dividing with the atmospheric density ρ of the ambient air, which then results in a jump of the absorbed dose rate values. This jump is not visible in Dartnell et al. (2015), therefore we assume that a smooth transition between the two atmospheric models was implemented.
Liquid water may have been present on the Venusian surface. Model studies by Grinspoon & Bullock (2013) suggest that this water may have evaporated and in addition, transported microorganisms into the Venusian cloud layer. Our further investigations therefore use the water-based 50µm microbial cell phantom (if not stated otherwise). Herbst et al. (2019a) showed that it is crucial to consider primary protons up to energies of 10 TeV in order to accurately reflect the Venusian atmospheric ion-pair production rate at altitudes below 40 km. In Fig. 4 we follow up on this investigation by studying the altitude profile of the induced absorbed dose rates as a function of the maximum primary proton energy. Although substantial differences between the upper limit of 100 GeV (solid blue line), 1 TeV (solid red line), and habitable zone 0% 10 TeV (solid black line) can be observed at altitudes below 40 km, only differences of up to 20 % can be found between an upper limit of 100 GeV and 10 TeV within the potential HZ. While studies of the radiation dose profiles of the entire altitude should include primaries with energies up to 10 TeV, this is therefore not mandatory when the habitability of a possible aerial biosphere is investigated. An upper limit of 1 TeV, as assumed by Dartnell et al. (2015), is sufficient here.

Investigating the influence of solar activity
We used the CREME model by Tylka et al. (1997) to calculate the primary particle spectra of Z= 1-28 nuclei as a function of the solar modulation. In order to study the influence of the modulation of GCR particles upon the Venusian atmospheric radiation dose, the absorbed dose rates during solar minimum and maximum conditions were computed. As shown in Fig. 5, the altitude-dependent values during low solar activity are given in black, and those during solar maximum conditions are displayed in red. The relative differences on the left show with respect to the induced radiation doses during solar minimum conditions (dashed lines) that the modulation of GCRs due to the solar activity only affects altitudes above 50 km and also the potential HZ. While the dose rates induced during high solar activity are about 50 % lower than the absorbed dose during solar minimum at altitudes above 80 km, the differences within the Venusian HZ still account for up to 20% and therefore are about 1.4 times higher during solar minimum conditions.

Investigating the dose equivalent
According to the approach we discussed in Section 2, the deposited ionization energy of each Z = 1-28 nuclei was scaled with the corresponding radiation weighting factor w R, j given in Table 1 and Eq. (5). Fig. 6 shows the corresponding modeled altitude-dependent absorbed dose rate profile (solid black line) and the dose equivalent profile (solid red line) that is the result of primary protons and helium nuclei alone during solar minimum conditions. As discussed in Sect. 4.1, this is a valid approximation because Z>2 nuclei only contribute with less than 20% to the Venusian radiation dose. The altitude-dependent ratio between the two measures therefore gives us the mean radiation weighting factor, which is comparable to the average quality factor Q (ICRP 60 1991). At the Martian surface, for example, the average quality factor is around three (see Hassler et al. 2014).
According to our results, the average quality factor at the upper Venusian cloud layer (∼72 km) is almost ten and slightly decreases toward a second peak within the HZ at around 57 km. At the lower boundary of the HZ and the lower cloud layer, in particular between 57 km and 47 km, Q rapidly drops to negligible values. This drop-off indicates that neither neutrons nor protons or nuclear fragments with the highest w R, j values contribute to the absorbed dose rates below the Venusian cloud layers. We therefore conclude that the atmospheric ionization and the absorbed dose rates at lower altitudes are solely induced by secondary muons.

Investigating the solar cosmic-ray-induced atmospheric radiation
From observations, we know that the spectral shape of GLE events and GCRs differs strongly. Although the energy spectrum of such strong SEPs is orders of magnitudes higher in the low-energy range, their energy spectra often do not exceed energies of a few GeV/nuc. In addition, different particle acceleration and transport mechanisms exist for the two populations. This means that the energy spectra of GLE events and their temporal evolution differ significantly from one another.
Although remarkable efforts revealing the temporal evolution of SEP events can be found in the literature (see, e.g., Belov et al. 2005;Bombardieri et al. 2007;Plainaki et al. 2007;Matthiä et al. 2009;Miroshnichenko 2018, and references therein), the model efforts to study the GLE-induced atmospheric effects are often based on a mean energy spectrum. As previously discussed in Nordheim et al. (2015), Plainaki et al. (2016), and Herbst et al. (2019a), the Venusian atmospheric ionization is significantly affected by these strong SEP events. We investigated their influence on the atmospheric absorbed dose rates.
In order to study the mean influence of strong SEP events on the atmospheric Venusian radiation dose, in a first step, we investigated one of the strongest GLE events of the modern era that occurred in October 1989. Based on the studies October 1989, and two specific cases were investigated: scenario 1, in which the 5-minute average event peak flux based on the GOES measurements was taken into account, and scenario 2, in which the 180-hour average of the particle flux was used.
The results are displayed in the upper panel of Fig. 7. To convey an idea of the impact strength of the event on the Venusian atmosphere, we also plot the GCR-induced radiation dose taking into account primary protons and alpha particles during solar minimum conditions. Even though most GLE events are limited to energies below 10 GeV, strong SEP events like the October 1989 event can significantly increase the Venusian radiation dose rate at altitudes above ∼ 50 km, and thus also well within the potential HZ.
The hazards of GLE events can better be estimated by studying the dose equivalent altitude profile, which is shown for the example of scenario 2 in the lower panel. The ratio between the two measures, again reflecting the average quality factor Q, is plotted on the left. While GCR particles induce a radiation field with Q ≈ 10 (see discussion in Section 4.6 and Fig. 6) at the top of the cloud layer, the 180-hour average of the GLE in October 1989 shows an average quality factor of ≈ 20, indicating that mostly neutrons are produced within the atmospheric layer. This shows that extreme SEP events can cause significant short-term increases in the radiation hazard in the cloud region above 50 km.
However, other differences with respect to the GCR-induced radiation dose become apparent: in the case of the GLEinduced radiation field,Q is drastically reduced below 50 km. Furthermore, below 20 km, no more secondary particles are produced, while Q is below 4 at altitudes above 80 km, which is much smaller than the GCR-inducedQ, as we show in the lower panel of Fig. 7.

Total absorbed dose rate induced by modern strong SEP events
In a second step, we studied the influence of 67 out of the 72 GLE events observed throughout the ground-based instrumental era on the Venusian atmosphere. Because Venus is closer to the Sun, the GLE spectra were scaled to the Venusian orbit in order to reflect a reliable radiation hazard they might have caused within the Venusian atmosphere. We assumed that all of the investigated GLE events measured at Earth also affected the Venusian atmosphere. Because SEPs are accelerated along the heliospheric magnetic field, a good connection between the Sun and the planet is mandatory in order for the event to arrive at Venus. This might not have been the case, however. This investigation should therefore be considered as a case study for the effect of strong SEP events on the Venusian radiation dose. The corresponding total absorbed dose rates, the sum of both the GCR-(solid black line) and GLE-induced absorbed dose rates, are displayed as colored lines in the left panel of Fig. 8. As already pointed out by Herbst et al. (2019a), substantial differences between the GCR-and GLE-induced ionization profiles can be observed (their Fig. 6). As pointed out, differences of up to 11 orders of magnitudes could be observed at high altitudes. The left panel of Fig. 8 shows that the differences within the CR-induced atmospheric radiation dose profiles at these altitudes are not as dramatic. Here, increases of up to six orders of magnitude can be found. Of particular interest for this study is the influence of strong SEP events on the Venusian HZ. Here, differences caused by the varying shapes of the GLE spectra can be observed. While   GLE05 (23 February 1956), the strongest GLE event directly measured at Earth so far, would have led to drastic increases in the absorbed radiation dose of microbial cells of two orders of magnitude within the HZ. We note, however, that strong SEP events are often accompanied by so-called Forbush decreases (FDs), which are sudden decreases of the GCR background intensity caused by interplanetary shocks or magnetized ejecta of coronal mass ejections (see, e.g., Cane 2000). In addition to this, the GCR-induced radiation dosage would be reduced, and thus the total absorbed radiation dose a microbial cell would suffer within the HZ during strong SEP events.

Total absorbed dose rate induced by historic extreme SEP events
However, from historical records like the white-light flare observation on 1 September 1859, for instance, which was accompanied by what is known as the Carrington event, or from the records of the so-called cosmogenic radionuclides 10 Be, 14 C, and 36 Cl, produced by interactions of CRs with the terrestrial atmosphere, we know that much stronger SEP events have occurred over the past 10,000 years. In particular, three exceptionally strong increases have been discussed in the literature: the AD774-775 event (see, e.g., Miyake et al. 2012), the AD993-994 event (see, e.g., Mekhaldi et al. 2015), and the BC660 event (see O'Hare et al. 2019). With a production increase of more than 15 % above the GCR-induced background, the AD774-775 event is so far the strongest event ever recorded. Model studies by Kovaltsov & Usoskin (2014) and Herbst et al. (2015) showed that the event was most likely more than one order of magnitude stronger than GLE05.
However, no information about the spectral shape of these historical extreme events is available. Therefore, assumptions about the spectral shape and its temporal evolution have to be made. A common approach is to rescale the spectra of known recent GLE events. For example, Dartnell et al. (2015) scaled the the event spectra of GLE05 (23 February 1956), GLE24 (4 August 1972), and GLE45 (24 October 1989) to match the estimated E>30 MeV integral fluence of the Carrington event (F 30 = 2×10 10 protons cm −2 Cliver & Dietrich 2013) and the estimated E> 200 MeV fluence of the AD774-775 event (F 200 = 8×10 9 protons cm −2 ) by Kovaltsov & Usoskin (2014). The duration of the events also has to be estimated in order to account for average differential spectra. Following the approach by Dartnell et al. (2015), an event duration of 20 hours was assumed for both historical events. Using these spectra (see, e.g., Fig. 3 in Dartnell et al. 2015), we here calculate the altitude-dependent absorbed dose rates microbial cells would suffer during such extreme solar events.
The results of this study are shown in the right panel of Fig. 8. Here the total absorbed dose rates of the Carrington and the AD774-775 event based on the different energy spectra are displayed. With increases of up to five orders of magnitude above the GCR-induced background, these extreme events would have a tremendous impact on the absorbed dose rate values within the HZ. The AD 774-775 event rescaled ton the basis of the hard spectrum of GLE05 would have had the most substantial influence throughout the entire HZ, while the Carrington event rescaled based on the soft spectrum of GLE24 would have had almost no impact on the absorbed dose rates within the HZ.
According to Dartnell et al. (2015), at the top of the HZ, the radiation dose of the GCR-induced absorbed dose is about 9.5·10 −4 mGy/h, while the GLE05-based Carrington-and AD775 events show values of 1.6 mGy/h and 4.9 mGy/h, respectively. A comparison with the values a water-based microbial cell would suffer at the lower and the upper boundary of the potential Venusian HZ as given in Table 3 shows that the radiation hazard for any putative Venusian cloud-based life is about five times higher than previously reported. This difference most likely arises because a Venusian-air-based phantom was used and the quality factors of different radiation types were therefore not accounted for.
In order to directly compare the Venusian radiation hazard on water-based life, we can compare our results to the absorbed dose rates measured at Earth. Here, the GCR-induced background leads to absorbed dose rates between 5×10 −4 mGy/h at the ground (∼ 1033 g/cm 2 , see, e.g., Banjac et al. 2019b) and 0.3 mGy/h at flight altitudes where the dose rate maximum is located (∼ 120 g/cm 2 see, e.g., Herbst et al. 2019b). Thus, the GCR-induced Venusian absorbed dose rates at the top of the HZ, for example, are well below the magnitudes that can be found at flight altitudes at Earth. As has been pointed out by Dartnell et al. (2015), even for the worst-case scenario of the GLE05-based AD774/775 spectrum, the induced absorbed radiation dose rates of about 33.9 mGy/h at the top of the habitable zone are in addition well below those of the impact of GLE05 on the terrestrial absorbed dose rate values (41.7 mGy/h, see Herbst et al. 2019b). However, as has been pointed out by Dartnell et al. (2015), even radiation-sensitive microorganisms, for example, Shewanella oneidensis, can survive much higher radiation dosages than those induced during this worst-case scenario. According to Ghosal et al. (2005), 90% of Shewanella oneidensis cells do not survive 24 hours of being exposed to a radiation of 70 Gy. The cells can therefore survive absorbed radiation doses below 2.9 Gy/h.
Consequently, all of the so-far known extreme SEP events (AD774-775, AD 992-993, and BC660) would not have had a hazardous effect on possible microorganisms within the Venusian habitable zone.

Summary and conclusions
Here we provided updated altitude-dependent Venusian radiation dose profiles based on the newly developed full 3D Monte Carlo simulation code AtRIS (Banjac et al. 2019b). We used the most recent numerical physics models (Geant4 10.5, FTFP_BERT_HP) as well as the primary GCR spectra of hydrogen (Z=1) to nickel (Z=28) from the CREME2009 model (see, e.g., Tylka et al. 1997) with energies between 1 MeV and 10 TeV. To study the influence of strong solar events on the Venusian radiation field, we also studied the effects of 67 space-age GLE events, and two extreme historical events.
We showed that the contribution of GCR protons and alpha particles to the induced Venusian radiation field, in particular in the HZ, are in the order of up to 75% and 20%, respectively. The contribution of Z>2 particles is less than 10 % and can therefore be neglected.
For the first time, we investigated the Venusian absorbed dose rates in different phantoms. Here we also took phantoms into account that are built out of water, which is the key ingredient for life as we know it from Earth. An ICRUsphere-like liquid-water phantom, mimicking the human body, and a microbial cell of 50 µm, were also developed and implemented in AtRIS. We showed that the computations by Dartnell et al. (2015) lead to an underestimation of the absorbed dose rates that water-based organisms would suffer within the Venusian atmosphere by up to 45% for a 50µm microbial cell, and up to 60% compared for a water-based ICRU-sphere phantom. However, as discussed below, compared to the limitations of lethal radiation levels, these differences may only play a minor role.
As recently shown by Herbst et al. (2019a), modeling the influence of strong SEP events on the Venusian atmosphere can be restricted to studying the influence of solar protons. Using the spectra of GLE events measured at Earth (provided by, e.g., Raukunen et al. 2018), we investigated the impact of the strongest SEP events measured throughout the space era (1942 to 2012) on the radiation field and thus also any putative organisms in the Venusian atmosphere. We showed that the strongest GLE event measured so far, the event of 23 February 1956, with a rather hard spectrum, would have led to a drastic increase in the absorbed dose rates of up to two orders of magnitude within the potential upper and middle habitable zone. From the cosmogenic radionuclide records of 10 Be, 14 C, and 36 Cl, however, we know that much stronger events have occurred in the past. The dramatic production rate increases of more than 15% above the GCR-induced background around AD774-775, AD993-994, and BC660 are of particular interest in this context (Miyake et al. 2012;Mekhaldi et al. 2015;O'Hare et al. 2019, respectively). As discussed in Nordheim et al. (2015) and Herbst et al. (2019a), these historical events would have caused dramatic changes in the Venusian ionization rates, of about up to six orders of magnitude. A comparison to Earth-bound measurements shows that the induced Venusian radiation dose values are much lower than the induced values within the Earth's atmosphere. We showed that all known extreme SEP events would not have had a hazardous effect on potential microorganisms, and even radiation-sensitive microorganisms like Shewanella oneidensis would have survived. This further supports the previous findings by Dartnell et al. (2015).
However, most recently, Notsu et al. (2019) found that old Sun-like stars can produce superflares with energies above 5×10 34 erg (5×10 −1 W/m 2 ) once every 2000 to 3000 years. Thus, much stronger SEP events than the imprints that have so far been found in the cosmogenic radionuclide records of 10 Be, 14 C, and 36 Cl around AD774-5, AD992-3, and 660 BC (see, e.g., Miyake et al. 2012;Mekhaldi et al. 2015;O'Hare et al. 2019, respectively) may have occurred on the Sun. According to Herbst et al. (2019c), such energetic flares could lead to particle events that are more intense by than two orders of magnitude than GLE05, which would drastically increase the absorbed radiation dose values within microbial cells in the HZ of the Venusian atmosphere. Even more intense flares with energies up to 10 36 erg (10 W/m 2 ) have been detected in young G-type stars like our Sun (see Notsu et al. 2019). This means that if they were produced by our Sun earlier in its history, such intense flares may have produced particle events five orders of magnitude stronger than GLE05 (see Herbst et al. 2019c). These superflare events could have had drastic impacts on any potential airborne life in the Venusian clouds.