Issue 
J. Space Weather Space Clim.
Volume 5, 2015



Article Number  A10  
Number of page(s)  16  
DOI  https://doi.org/10.1051/swsc/2015011  
Published online  27 May 2015 
Research Article
Numerical model for computation of effective and ambient dose equivalent at flight altitudes
Application for dose assessment during GLEs
^{1}
ReSoLVE Center of Excellence, University of Oulu, Finland
^{2}
Sodankylä Geophysical Observatory (Oulu Unit), University of Oulu, Finland
^{*} Corresponding author: alex_mishev@yahoo.com
Received:
8
December
2014
Accepted:
20
April
2015
A numerical model for assessment of the effective dose and ambient dose equivalent produced by secondary cosmic ray particles of galactic and solar origin at commercial aircraft altitudes is presented. The model represents a full chain analysis based on groundbased measurements of cosmic rays, from particle spectral and angular characteristics to dose estimation. The model is based on newly numerically computed yield functions and realistic propagation of cosmic ray in the Earth magnetosphere. The yield functions are computed using a straightforward full Monte Carlo simulation of the atmospheric cascade induced by primary protons and αparticles and subsequent conversion of secondary particle fluence (neutrons, protons, gammas, electrons, positrons, muons and charged pions) to effective dose or the ambient dose equivalent. The ambient dose equivalent is compared with reference data at various conditions such as rigidity cutoff and level of solar activity. The method is applied for computation of the effective dose rate at flight altitude during the ground level enhancement of 13 December 2006. The solar proton spectra are derived using neutron monitor data. The computation of the effective dose rate during the event explicitly considers the derived anisotropy i.e. the pitch angle distribution as well as the propagation of the solar protons in the magnetosphere of the Earth.
Key words: Atmospheric cascade simulation / Yield function / Effective dose rate / Ambient dose equivalent / Ground level enhancement
© A. Mishev and I. Usoskin, Published by EDP Sciences 2015
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Our planet is constantly hit by energetic particles, known as cosmic rays (CRs). Primary CR particles are mostly protons and αparticles with a small addition of heavier nuclei. They penetrate deep into the atmosphere, producing large amount of secondary particles via a complicated nuclearelectromagneticmuon atmospheric cascade. In such a cascade only a fraction of the initial primary particle energy reaches the ground as secondaries. Most of the primary particle’s energy is released in the atmosphere by ionization and excitation of the air (Dorman 2004; Bazilevskaya et al. 2008; Usoskin et al. 2009; Vainio et al. 2009). There is a general agreement that the majority of CRs originate from the Galaxy, called galactic cosmic rays (GCRs). Their intensity varies by 10–20% modulated by the solar wind i.e. it depends on the level of the solar activity. Therefore it follows the 11year solar cycle and responds to solarwind variations and transient phenomena (Forbush 1958).
The Earth is also hit sporadically by solar energetic particles (SEPs), which may sometimes produce an atmospheric cascade in a similar way. Solar energetic particles are accelerated during explosive energy releases on the Sun (see Reames 1999; Cliver et al. 2004; Aschwanden 2012, and references therein). While the energy of SEPs is usually of the order of a few tens of MeV, in some cases it can reach a few GeV. According to a recent study, there is an evidence of possible proton acceleration up to energies in excess of 20 GeV (Bostanjyan et al. 2007). Therefore, they penetrate deep into the atmosphere or even reach the ground, leading to the socalled ground level enhancements (GLEs). Such events occur roughly once a year with a higher probability during a solar maximum (Shea & Smart 1990). Cosmic ray particles significantly affect the radiation environment and accordingly exposure at commercial flight altitudes, specifically during GLEs (O’Brien et al. 1997; Shea & Smart 2000, 2012; Spurny et al. 2002). In addition, recently an assessment of thunderstorm radiation environment at aviation flight altitudes has been carried out, which could be significant under some circumstances (Dwyer et al. 2010; Drozdov et al. 2013).
Since aircrews are subject to increased exposure comparing to the sea level due to the increased intensity of secondary CR, the radiation protection of aircrews has been widely discussed. As a result, new regulations appeared, namely Publication 60 of the International Commission on Radiological Protection (ICRP), where the exposure of flying personnel to cosmic radiation is recommended to be regarded as occupational (ICRP 1991). Accordingly in EU, the Euratom Directive 96/12 suggests measures to assess the individual doses of cockpit and cabin crew (EURATOM 1996).
A standard way to assess the exposure is computation based on CR measurements. In this paper, we describe a numerical model to estimate effective and ambient dose equivalent at flight altitudes, based on newly computed yield functions and reconstruction procedure of SEPs based on ground measurements. We focus on dose assessment during GLEs.
2. Numerical model for computation of effective dose rate and ambient dose equivalent at flight altitude
Accessment of primary CRs to the Earth atmosphere is influenced by the spatialtemporal variability of complex magnetospheric and interplanetary conditions. Therefore, the task to estimate the aircrew exposure due to CR of galactic and solar origin is not trivial, since it depends on geographic position, altitude and solar activity (Spurny et al. 1996, 2002; Vainio et al. 2009). An assessment of the radiation dose hazard due to CR involves several consecutive steps. In the first step fluxes of GCR and SEPs are determined outside the magnetosphere. Then their propagation through the atmosphere and magnetosphere of the Earth is modelled, resulting in a realistic evaluation of the secondary CR flux. In the work presented here we apply an appropriate magnetospheric model (see below) and new numerically computed yield functions for conversion of the secondary particle fluence to dose (effective/ambient dose equivalent). We also present a convenient procedure for estimation of SEP spectral and angular characteristics based on a neutron monitor data analysis.
2.1 Definition of the dose rate
The dose rate can be computed as a function of the geomagnetic rigidity cutoff and altitude using a full Monte Carlo simulation of the atmospheric cascade (Ferrari et al. 2001). At present, several models have been proposed, aiming to estimate the dose rate (effective and/or ambient dose equivalent) at flight altitudes due to primary CR radiation (Schraube et al. 2000; Ferrari et al. 2001; Roesler et al. 2002; Lewis et al. 2005; Matthiä et al. 2008; Sato et al. 2008; Mishev & Hristova 2012; Mertens et al. 2013; Al Anid et al. 2014; Mishev 2014; Mishev et al. 2015).
According to the definition, the absorbed dose is the energy deposited in a medium by ionizing radiation per unit mass. It is measured as joules per kilogram, i.e. according to SI the unit is Gy. In order to assess the biological risk due to radiation exposure it is convenient to use the effective dose. Since the effective dose is not a measurable quantity, ICRP suggests for operational purpose the ambient dose equivalent (ICRP 2007) denoted as H*(d). It represents the dose equivalent that would be produced by the corresponding expanded and aligned field at a depth d in an International Commission on Radiation Units and Measurements (ICRU) sphere (a sphere with diameter of 30 cm made of tissue equivalent material with a density of 1 g cm^{3} and a mass composition of 76.2% oxygen, 11.1% carbon, 10.1% hydrogen and 2.6% nitrogen) on the radius vector opposing the direction of the aligned field. The unit for both effective dose and ambient dose equivalent is Sv. According to recent studies the ambient dose equivalent at a depth of d = 10 mm H*(10) is recommended as a reasonable proxy for the effective dose (Pelliccioni 2000). In general it slightly overestimates the effective dose. However, it is not a conservative estimate for cosmic radiation exposure at aviation altitudes according to Mertens et al. (2013) and the work presented here, specifically in high rigidity regions. Nevertheless, it is regarded as an acceptable approximation for effective dose at aircraft altitudes (Meier et al. 2009; Mertens et al. 2013).
2.2. Formalism
The effective dose rate at a given atmospheric depth h induced by a primary CR particle can be determined as:(1)where T′ is the energy of the primary CR particle arriving from zenith angle θ and azimuth angle φ, J _{ i } (T′) is the differential energy spectrum of the primary CR at the top of the atmosphere for i component (proton and/or αparticle), λ _{ m } is the geomagnetic latitude, Ω(θ, φ) is a solid angle and Y _{ i } is the effective dose yield function.
Accordingly, the effective dose yield function Y _{ i } is defined as:(2)where C _{ j }(T ^{*}) is the fluence to effective dose conversion coefficient for a secondary particle of type j (neutron, proton, γ, e ^{ − }, e ^{ + }, μ ^{ − }, μ ^{ + }, π ^{−}, π ^{+}) with energy T ^{*}, F _{ i,j } (h, T′, T ^{*}, θ, φ) is the secondary particle fluence of type j, produced by a primary particle of type i (proton and/or αparticle) with a given primary energy T′. The conversion coefficients C _{ j }(T ^{*}) (see Pelliccioni 2000; PetoussiHenss et al. 2010, and references therein) for a secondary particle of type j are obtained on the basis of extensive Monte Carlo simulations with various codes including a precise crosscheck, namely: FLUKA (Fasso et al. 2005; Battistoni et al. 2007), MCNPx (Briesmeister 1997; Waters et al. 2007), PHITS (Iwase et al. 2002), GEANT4 (Agostinelli et al. 2003) and EGSnrc (Kawrakow 2001) performed by the DOCAL task Group assuming a reference computational phantom (ICRP 2009). The secondary particle fluence at a given altitude above the sea level (a.s.l.) is obtained here with extensive simulations of the atmospheric cascade with GEANT4 (Agostinelli et al. 2003) based PLANETOCOSMICS code (Desorgher et al. 2005) assuming a realistic atmo spheric model NRLMSISE2000 (Picone et al. 2002). Here we assume an isotropic secondary particle flux, which is in good agreement with a realistic case for aircrew dosimetry (Sato et al. 2011). The atmospheric cascade simulations have been performed in a wide energy range, namely between 0.5 GeV/nucleon and 1 TeV/nucleon for primary protons and αparticle, assuming isotropic particle incidence. The effective dose yield function considers explicitly the complexity of the atmospheric cascade development, since it brings information of particle fluence and spectrum at a given altitude in the atmosphere, considering explicitly the secondary particle attenuation.
Accordingly, the ambient dose equivalent H*(10) i.e. the dose equivalent produced at a depth of 10 mm in a ICRU sphere, at given atmospheric altitude h induced by a primary cosmic ray particle can be determined as:(3)where is the ambient dose equivalent yield function determined in a way similar to Eq. (2). Here we adopt the data sets for an isotropic particle fluence from Appendix 2 of Pelliccioni (2000) i.e. the fluence to ambient dose equivalent conversion coefficients for the described above secondaries. Accordingly the fluence to ambient dose equivalent conversion coefficients is obtained using Monte Carlo simulations with the FLUKA code (Fasso et al. 2005; Battistoni et al. 2007).
Yield functions for the effective dose for primary protons and αparticles at altitude of 35 kft (typical commercial flights ≈ 10.5 km) and altitude of 50 kft (supersonic flights ≈ 5 km) are presented in Figure 1 with details given in Table A.1 of Appendix A. Accordingly, Table A.2 gives details for the ambient dose equivalent yield function. For a realistic interpretation of SEP with oblique particle incidence we present also effective dose yield function for primary protons with vertical and 15°, 30° and θ = 45° incidence (Clem 1997; Cramp et al. 1997; Vashenyuk et al. 2006). Accordingly, the details are given in Table A.3.
Fig. 1.
Effective dose yield function as a function of the energy per nucleon for primary CR protons and αparticles at two different altitudes a.s.l., namely 35 kft (typical commercial flight ≈ 10.5 km) and 50 kft (supersonic flight ≈ 15 km). 
As mentioned above, it is necessary to have a precise model for the CR spectrum. Here we use a force field model for the GCR (see below). Since the flux and spectra of SEPs vary considerably between solar eruptive events, we describe below a method for determination of their spectral and angular characteristics based on neutron monitor (NM) data.
2.3. GCR spectrum
A detailed description of force field model for GCR spectrum is given elsewhere (Usoskin et al. 2005). The GCR flux is affected by the interplanetary magnetic field and solar wind. It results in modulation of GCR in the vicinity of the Earth. The modulation varies with solar activity (Gleeson & Axford 1968; CaballeroLopez & Moraal 2004). The only explicit parameter of this model is the modulation potential ϕ given in units of MV. The value of Zeϕ (Z is the charge number) corresponds to the average energy loss of cosmic rays inside the heliosphere. Here we consider the nucleonic ratio of heavier particles including αparticles to protons in the interstellar medium as 0.3 (Gaisser & Stanev 2010) similarly to Kovaltsov et al. (2012). The model with the corresponding parametrization of LIS provides very good fitting of the measured spectra (for details see Usoskin et al. 2005, and references therein). For the computations, a realistic mass composition of the primary GCR (Gaisser & Stanev 2010) is assumed, where the contribution of heavier species is scaled to that for αparticles (Usoskin & Kovaltsov 2006; Mishev & Velinov 2011; Kovaltsov et al. 2012).
2.4. Model of the Earth magnetosphere
Both GCR and SEPs are shielded by the geomagnetic field. The shielding is most effective near the geomagnetic equator and absent close to the geomagnetic poles. The capacity of the shielding is approximately quantified by the effective vertical rigidity cutoff R _{ c }, which varies with the geographical location (Cooke et al. 1991).
Since the precise particle propagation in the magnetosphere is important for both GLE analysis and access of SEPs, which typically possess an essential anisotropic part, to a given location, a realistic geomagnetospheric model is necessary (Smart et al. 2000; Desorgher et al. 2009; Mishev & Usoskin 2013). Several models, tools, algorithms to calculate geomagnetic cutoff rigidity, particle trajectories and the asymptotic viewing cones were proposed during the years (McCracken et al. 1962, 1968; Shea et al. 1965; Smart et al. 2000). As the internal field we consider the IGRF geomagnetic model, which is a Gauss spherical harmonic model of the geomagnetic field, based on magnetic field measurements from geomagnetic stations, magnetometers and satellites (Langel 1987). As the external field model we use a semiempirical Tsyganenko (1989) model which is based on a large number of satellite observations (Tsyganenko 1989). The model takes into account contributions from external magnetospheric sources: ring current, magnetotail current system, magnetopause currents and a largescale system of fieldaligned currents. The model takes into consideration seasonal and diurnal changes of the magnetospheric field as well as the geomagnetic activity level K _{ p }. It provides seven different states of the magnetosphere corresponding to different levels of geomagnetic activity. It is driven only by the geomagnetic activity index K _{ p } and provides perfect balance between simplicity (Nevalainen et al. 2013) and realism (Kudela & Usoskin 2004; Kudela et al. 2008). All computations of particle’s transport in the geomagnetic field are performed with the MAGNETOCOSMICS code (Desorgher et al. 2005).
It is known that during major SEP events the configuration of the magnetosphere could be changed dramatically. Accordingly a magnetospheric model has been recently proposed (Tsyganenko & Sitnov 2005). The influence of a stormtime geomagnetic field model (Tsyganenko & Sitnov 2005), specifically on computation of neutron monitor asymptotic cones and the subsequent analysis of GLEs, is beyond the topic of this work and should be a subject of future study.
2.5. Analysis of GLEs using data from the global neutron monitor network
Details for the analysis of GLE events based on data from several NM stations are given in (Shea & Smart 1982; Cramp et al. 1997; Bombardieri et al. 2006; Vashenyuk et al. 2006, 2008). Detailed description of the method used in this study is given by Mishev et al. (2014). The relative count rate increase of a given NM to a SEP event can be expressed as:(4)where J _{sep}(P) is the primary solar particles rigidity spectrum in the direction of a maximal flux, J _{GCR}(P, t) is the rigidity spectrum of GCR at given time t, Y _{NM} (P) is the NM yield function (Mishev et al. 2013), G(α(P)) is the pitch angle distribution of SEPs, N is the NM count rate due to GCR and ΔN(P _{cut}) is the count rate increase due to solar particles. In Eq. (4), P _{cut} is the minimum rigidity cutoff of the station and P _{max} is the maximum rigidity considered in the model, assumed to be 20 GV. The fractional increase of the count rate of a NM station represents the ratio between the NM count rate due to SEP and that due to the GCR averaged over 2 h before the event’s onset. For the GCR we apply the described above force field model and reconstructed solar modulation parameter according to Usoskin et al. (2011). A modified power law rigidity spectrum of SEP is assumed in the analysis similarly to Cramp et al. (1997); Bombardieri et al. (2006); Vashenyuk et al. (2008):(5)where J _{sep} is the particle flux arriving from the Sun along the pitch angle distribution (PAD) axis of symmetry defined by geographic coordinate angles Ψ and Λ, γ is the powerlaw spectral exponent at rigidity P = 1 GV and δγ is the rate of the spectrum steepening. The pitch angle distribution is assumed as a superposition of two Gaussians:(6)where α is the pitch angle, σ _{1} and σ _{2} are parameters corresponding to the width of the pitch angle distribution and B is a parameter. At B → 0, PAD becomes a Gaussian similar to Vashenyuk et al. (2006). The pitch angle is defined as the angle between the asymptotic direction and the axis of anisotropy. Therefore according to Eqs. (4)–(6), nine model parameters have to be determined: J _{0}, γ, δγ, Ψ and Λ, α′, σ _{1}, σ _{2} and B. This PAD allows to consider a bidirectional anisotropy similarly to Cramp et al. (1997), results from bouncing of particles inside a looped magnetic field structure (magnetic bottle). In the study presented here we obtained the best fit with a simple Gaussian for the PAD (see Table B.1).
3 Application of the model for computation of the dose rate at various conditions
Since the flux of GCR is isotropic, it is necessary to compute the rigidity cutoff at a given condition and to apply the described above model to assess the exposure due to GCR. During the GLE events, the exposure is a superposition of CR particles of galactic and solar origin. Since SEPs possess an essential anisotropic part, it is necessary to compute the asymptotic cones in the region of interest, in our case in a grid of 5° × 5° in order to consider explicitly the anisotropy using the derived PAD.
3.1. Exposure due to GCR and comparison of the model with reference data
In order to compare the model with the reference data (Menzel 2010) we perform computations in various conditions. The ambient dose equivalent is computed at altitude of 35 kft a.s.l. using the corresponding yield function. The estimated ambient dose equivalent as a function of the rigidity cutoff is presented in Figures 2a–2c. A good agreement is achieved, specifically in mid and high rigidity cutoff regions.
Fig. 2.
Ambient dose equivalent H*(10) at the altitude of 35 kft a.s.l. as a function of the rigidity cutoff computed at various conditions, compared with reference data Menzel (2010). (a) Ambient dose equivalent H*(10) for January 1998. Modulation potential 427 MV, (b) ambient dose equivalent H*(10) for January 2000. Modulation potential 752 MV, (c) ambient dose equivalent H*(10) for January 2002. Modulation potential 977 MV. 
Fig. 3.
Effective dose rate at the altitude of 35 kft a.s.l. as a function of the rigidity cutoff computed at various conditions. (a) January 1998, modulation potential 427 MV, (b) January 2000, modulation potential 752 MV, (c) January 2002, modulation potential 977 MV. 
Fig. 4.
The effective dose rate (μSv h^{−1}) at the altitude of 35 kft a.s.l. during the GLE 70 on 13 December 2006 in a region with R _{ c } ≤ 1 GV. (a) initial phase of the event, (b) main phase of the event, (c) late phase of the event. 
Fig. 5.
The relative difference in % of the effective dose rate at the altitude of 35 kft a.s.l. during the GLE 70 on 13 December 2006 assuming different spectral and angular characteristics of SEPs. (a) initial phase of the event, (b) main phase of the event. 
In general, our model computations demonstrate good agreement with reference data (see Figs. 2 and A.1). The modelled ambient dose equivalent tends to slightly overestimate the dose in the low rigidity cutoff region and underestimate it in high rigidity cutoff regions, similarly to Mertens et al. (2013). The slight underestimation of the reference data at the high level of the rigidity cutoff could be explained by the contribution of thunderstorm radiation, widely discussed in Dwyer et al. (2010) and Drozdov et al. (2013). The detailed study of the contribution of thunderstorm radiation within our model to exposure is a topic of a further work. The observed difference for 1998 and a marginal difference for 2000 in the region of low rigidity cutoff are mostly due to modulation effects and some specifics (cascade maximum) of atmospheric cascade development of heavier nuclei within their scaling to αparticles (Mishev & Velinov 2014). Accordingly, Figures 3a–3c show the computed effective dose rate as a function of the rigidity cutoff.
Some details are given in Appendix A. As one can see in Figure A.2, there is a general agreement between the model computations and reference data. The difference between the effective dose E and the ambient dose equivalent H*(10) is shown in Figure A.3. One can see that the H*(10) is a conservative approach for the exposure in a low rigidity cutoff regions, while the effective dose is a conservative approach above some 5 GV in accord with a recent study (Mertens et al. 2013). Moreover, our model computations demonstrate better agreement with the reference data (Figs. A1 and A4) comparing to some recent works Mertens et al. (2013).
Our model also demonstrates good agreement with a recently proposed model for numerical computation of the radiation exposure from GCRs at aviation altitudes (Matthiä et al. 2014). Both models demonstrate good agreement with experimental/reference data, although the calculations slightly overestimate the measurements in the polar region, which are mostly due to modulation/GCR model effects.
3.2. Effective dose rate computed during GLE 70 on 13 December 2006
The event of 13 December 2006 occurred during the decline phase of solar cycle 23. The event was characterized by a large anisotropy in its initial phase. It has been well studied in the sense of exposure, specifically at flight altitudes (Beck 2009; Bütikofer et al. 2009; Matthiä et al. 2009). We derive the spectral and angular characteristics of SEPs using NM data as described above. The computed asymptotic cones of the NMs used for the analysis, the direction of heliospheric magnetic field (IMF) derived from ACE satellite measurements at 03:00 UT, the derived apparent source position and lines of equal pitch angles relative to the derived anisotropy axis for 30°, 60°, 150° and 120° are shown in Figure B.1 of Appendix B. Accordingly, the estimated apparent source position throughout the event is shown in Figure B.2. The spectral and angular characteristics are presented in Figure B.3. Accordingly, Table B.1 summarizes the details.
Subsequently, the effective dose rate at altitude of 35 kft a.s.l. was computed on the basis of the described above model and the computed yield function (Appendix A). Since the initial phase of the event was very anisotropic, taking into account the falling SEP spectrum, we estimated the effective dose rate in a region with R _{ c } ≤ 1 GV, where the expected exposure is maximal. The computed effective dose rate at the flight altitude during the initial phase of GLE 70 on 13 December 2006 was about 40–50 μSv h^{−1} (Fig. 4a), which is in a very good agreement with previous estimations (Matthiä et al. 2009). The computed effective dose rate during the main and late phases of the event is considerably lower (Figs. 4b and 4c).
It is apparent that the contribution of SEP to the exposure due to CR is important, specifically during the initial phase of GLE 70. During the main and late phases of GLE 70 the contribution of SEPs to the exposure is slightly greater compared to the average due to GCR.
As demonstrated recently the computed effective dose rate at flight altitude during GLEs is model dependent, specifically by the derived SEP spectra, angular distribution and apparent source position (Bütikofer & Flückiger 2013). In this connection, we compared the estimated effective dose rate at altitude of 35 kft a.s.l. during GLE 70 assuming various SEP spectra (Vashenyuk et al. 2008). Figure 5a shows the relative difference of the computed effective dose rate during the initial phase of GLE 70 assuming our derived spectra (Table B.1) and spectra by Vashenyuk et al. (2008). Similarly, the relative difference for the main phase is presented in Figure 5b.
One can see that the relative difference does not exceed ~20%, which is in very good agreement with the results of Bütikofer & Flückiger (2013). During the initial phase of the event this difference is mostly due to the difference in the apparent source position determination, rather than the derived SEP spectra. The dramatic change of the observed difference in antisun direction is due to the point like PAD of SEPs during the initial phase of the event (Fig. B.3 and Table B.1). During the main and late phases of the event the derived by us SEP spectra lead to a more conservative approach of the estimated effective dose (Fig. 5b).
4 Discussion
According to the common definition, space weather refers to the dynamic, variable conditions on the Sun, solar wind and Earth’s magnetosphere and ionosphere, that can compromise the performance and reliability of spacecraft and groundbased systems and can endanger human life or health (Baker 1998; Lilensten & Belehaki 2009; Lilensten & Bornarel 2009). An important feature of space weather is the assessment of aircrew exposure due to CR, specifically during GLEs. At present there is no evidence of any immediate threat on human health from SEPs, and the biological risk of radiation doses accumulated by aircrew is still a matter of scientific debate (Sigurdson & Ron 2004; Ballarini et al. 2007; Hammer et al. 2009; Pukkala et al. 2012; Dos Santos Silva et al. 2013). In addition, a recent survey of the exposure due to CRs demonstrated that none exceeded the recommended level of 6 mSv/y and the majority of aircrew received around 3 mSv/y (Bennett et al. 2013a, 2013b). However, according to several studies discovering chromosome aberrations in aircrew members, the topic deserves further investigations (Bolzan et al. 2008; Yong et al. 2009) as well as additional verification (Wolf et al. 1999). Therefore, the development of new models as well as the improvement and validation of existing models for dose assessment, specifically during GLEs, is of a big importance.
Since the difference between effective dose and ambient dose equivalent is slightly above the intrinsic atmospheric cascade simulations, specifically in the lower energy range, the application of both of them for dose assessment at flight altitudes is quite reasonable. The ambient dose equivalent looks more conservative in a low rigidity cutoff, while effective dose at high rigidity cutoff. Therefore the application of ambient dose equivalent is more conservative in the case of exposure assessment during SEP events, because the expected radiation environment change is more important in polar and subpolar regions. Although, the application of effective dose responds to regulations and it is fully applicable as it is carried out in the presented here work.
5 Conclusions
In this paper we describe a fully operational model for exposure assessment of aircrew based on newly computed yield functions for effective and/or ambient dose equivalent. The model represents a full chain analysis of CR measurements with NMs, from spectral and angular characteristics of SEPs evaluation to dose assessment.

The model is validated and compared with reference data.

New yield functions for effective and/or ambient dose equivalent are presented.

The effective dose rate during GLE 70 on 13 December is computed as an example of application of the model.
Acknowledgments
This work was supported by the Center of Excellence ReSoLVE (Project No. 272157) and the University of Oulu Finland. The authors are grateful to colleagues from Institute for Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences and International Atomic Energy Agency in Vienna for technical support and discussions. The editor and authors thank T. Karapetyan and A. Drozdov for their assistance in evaluating this paper.
Appendix A.
Effective dose and ambient dose equivalent tables
Table A.1 presents the effective dose yield function Y [Sv cm^{2} sr nucleon^{−1}] for primary protons and αparticles with isotropic incidence. The effective dose yield function is computed at the altitudes of 35 kft a.s.l. (typical commercial flights ≈ 10.5 km a.s.l.) and 50 kft a.s.l. (supersonic flights ≈ 15 km a.s.l.). Accordingly, the ambient dose equivalent H*(10) yield function Y* [Sv cm^{2} sr nucleon^{−1}] for primary protons and αparticles at the altitude of 35 kft a.s.l. is given in Table A.2.
Fig. A.1.
Relative difference between computed ambient dose equivalent H*(10) at the altitude of 35 kft a.s.l. and the reference data (Menzel 2010). (a) January 1998, modulation potential 427 MV, (b) January 2000, modulation potential 752 MV, (c) January 2002, modulation potential 977 MV. The modulation potential is according to Usoskin et al. (2011). 
Fig. A.2.
Computed effective dose E, ambient dose equivalent H*(10) and reference data (Menzel 2010) at the altitude of 35 kft a.s.l. (a) January 1998, (b) January 2000, (c) January 2002. 
Fig. A.3.
Relative difference between computed effective dose E and ambient dose equivalent H*(10) at the altitude of 35 kft a.s.l. (a) January 1998, (b) January 2000, (c) January 2002. 
Fig. A.4.
Relative difference between computed effective dose E at the altitude of 35 kft a.s.l. and reference data (Menzel 2010). (a) January 1998, (b) January 2000, (c) January 2002. 
Fig. A.5.
Computed effective dose rate E at the altitudes of 50 kft a.s.l. and 35 kft a.s.l. (a) January 1998, (b) January 2000, (c) January 2002. 
Effective dose yield function Y [Sv cm^{2} sr nucleon^{−1}] for primary proton and αparticles at the altitudes of 35 kft and 50 kft a.s.l.
Ambient dose equivalent H*(10) yield function Y* [Sv cm^{2} sr nucleon^{−1}] for primary protons and αparticles at the altitude of 35 kft a.s.l.
Effective dose yield function Y [Sv cm^{2} sr nucleon^{−1}] for primary protons with various incidence at the altitude of 35 kft a.s.l.
Appendix B.
Spectral and angular characteristics of SEPs during GLE 70 on 13 of December 2006
Here we use the data from several NM stations for the analysis of the GLE 70 on 13 December 2006. They are given in Table B.2. Most of the data are retrieved from NM database (Mavromichalaki et al. 2011).
Fig. B.1.
Calculated NM asymptotic directions during GLE 70 on 13 December 2006 at 03:00 UT. The cross represents the direction of heliospheric magnetic field (IMF) derived from ACE satellite measurements at 03:00 UT. The small circle represents the derived apparent source position. The lines of equal pitch angles relative to the derived anisotropy axis are plotted for 30°, 60°, 150° and 120°. The asymptotic directions of polar NMs are plotted with solid lines, while midlatitude NMs are plotted with dashed lines. 
Fig. B.2.
Derived apparent source position (red line) and heliospheric magnetic field (crosses) from ACE satellite measurements throughout GLE 70 on 13 December 2006. 
Fig. B.3.
The derived spectral and angular characteristics of SEPs for GLE 70 on 13 December 2006 for 03:00 UT throughout 06:00 UT as denoted in the legend. 
Derived spectral and angular characteristics for GLE 70 on 13 December 2006.
Neutron monitors with corresponding geomagnetic rigidity cutoffs used in the analysis.
References
 Agostinelli, S., J. Allison, and K. Amako. Geant4 – a simulation toolkit. Nucl. Instrum. Methods Phys. Res.: Sect. A, 506 (3), 250–303, 2003, DOI: 10.1016/S01689002(03)013688. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Al Anid, H., B. Lewis, L. Bennett, M. Takada, and M. Duldig. Aircrew radiation dose estimates during recent solar particle events and the effect of particle anisotropy. Radiat. Prot. Dosim., 158 (3), 355–367, 2014, DOI: 10.1093/rpd/nct234. [CrossRef] [Google Scholar]
 Aschwanden, M. GeV particle acceleration in solar flares and ground level enhancement (GLE) events. Space Sci. Rev., 171 (14), 3–21, 2012, DOI: 10.1007/s112140119865x. [NASA ADS] [CrossRef] [Google Scholar]
 Baker, D. What is space weather? Adv. Space Res., 22 (1), 7–16, 1998. [Google Scholar]
 Ballarini, F., D. Alloni, A. Facoetti, A. Mairani, R. Nano, and A. Ottolenghi. Radiation risk estimation: modelling approaches for “targeted” and “nontargeted” effects. Adv. Space Res., 40 (9), 1392–1400, 2007, DOI: 10.1016/j.asr.2007.04.021. [CrossRef] [Google Scholar]
 Battistoni, G., S. Muraro, P. Sala, F. Cerutti, A. Ferrari, S. Roesler, A. Fasso, and J. Ranft. The FLUKA code: description and benchmarking. In: Proceedings of the Hadronic Shower Simulation Workshop 2006, 6–8 September 2006, M., Albrow and R. Raja Editors, AIP Conference, 896, 31–49, 2007. [Google Scholar]
 Bazilevskaya, G.A., I.G. Usoskin, E. Flückiger, R. Harrison, L. Desorgher, et al. Cosmic ray induced ion production in the atmosphere. Space Sci. Rev., 137, 149–173, 2008, DOI: 10.1007/s112140089339y. [Google Scholar]
 Beck, P. Overview of research on aircraft crew dosimetry during the last solar cycle. Radiat. Prot. Dosim., 136 (4), 244–250, 2009, DOI: 10.1093/rpd/ncp158. [CrossRef] [Google Scholar]
 Bennett, L., B. Lewis, B. Bennett, M. McCall, M. Bean, L. Dor, and I. Getley. Cosmic radiation exposure survey of an air force transport squadron. Radiat. Meas., 48 (1), 35–42, 2013a, DOI: 10.1016/j.radmeas.2012.10.012. [CrossRef] [Google Scholar]
 Bennett, L., B. Lewis, B. Bennett, M. McCall, M. Bean, L. Dor, and I. Getley. A survey of the cosmic radiation exposure of Air Canada pilots during maximum galactic radiation conditions in 2009. Radiat. Meas., 49 (1), 103–108, 2013b, DOI: 10.1016/j.radmeas.2012.12.004. [CrossRef] [Google Scholar]
 Bolzan, A., M. Bianchi, E. Gimenez, M. Flaque, and V. Ciancio. Analysis of spontaneous and bleomycininduced chromosome damage in peripheral lymphocytes of longhaul aircrew members from Argentina. Mutation Research – Fundamental and Molecular Mechanisms of Mutagenesis, 639 (1–2), 64–79, 2008, DOI: 10.1016/j.mrfmmm.2007.11.003. [CrossRef] [Google Scholar]
 Bombardieri, D., M. Duldig, K. Michael, and J. Humble. Relativistic proton production during the 2000 July 14 solar event: the case for multiple source mechanisms. Astrophys. J., 644 (1), 565–574, 2006, DOI: 10.1086/501519. [Google Scholar]
 Bostanjyan, N., A. Chilingarian, V. Eganov, and G. Karapetyan. On the production of highest energy solar protons at 20 January 2005. Adv. Space Res., 39 (9), 1456–1459, 2007, DOI: 10.1016/j.asr.2007.03.024. [CrossRef] [Google Scholar]
 Briesmeister, J. MCNP A general Monte Carlo Transport code (version 4B). Tech. Rep. LA 12625M, Los Alamos National Laboratory, 1997. [Google Scholar]
 Bütikofer, R., and E. Flückiger. Differences in published characteristics of GLE60 and their consequences on computed radiation dose rates along selected flight paths. J. Phys: Conf. Ser., 409 (1), 012166, 2013. [Google Scholar]
 Bütikofer, R., E. Flückiger, L. Desorgher, M. Moser, and B. Pirard. The solar cosmic ray groundlevel enhancements on 20 January 2005 and 13 December 2006. Adv. Space Res., 43 (4), 499–503, 2009, DOI: 10.1016/j.asr.2008.08.001. [Google Scholar]
 CaballeroLopez, R., and H. Moraal. Limitations of the force field equation to describe cosmic ray modulation. J. Geophys. Res., 109, A01101, 2004, DOI: 10.1029/2003JA010098. [Google Scholar]
 Clem, J. Contribution of Obliquely incident particles to neutron monitor counting rate. J. Geophys. Res., 102, 919, 1997. [Google Scholar]
 Cliver, E., S. Kahler, and D. Reames. Coronal shocks and solar energetic proton events. Astrophys. J., 605, 902–910, 2004, DOI: 10.1086/382651. [Google Scholar]
 Cooke, D., J. Humble, M. Shea, D. Smart, N. Lund, I. Rasmussen, B. Byrnak, P. Goret, and N. Petrou. On cosmicray cutoff terminology. IL Nuovo Cimento C, 14 (3), 213–234, 1991. [Google Scholar]
 Cramp, J., M. Duldig, E. Flückiger, J. Humble, M. Shea, and D. Smart. The October 22, 1989, solar cosmic ray enhancement: an analysis of the anisotropy spectral characteristics. J. Geophys. Res., 102 (A11), 24237–24248, 1997. [Google Scholar]
 Desorgher, L., E. Flückiger, M. Gurtner, M. Moser, and R. Bütikofer. A Geant 4 code for computing the interaction of cosmic rays with the Earth’s atmosphere. Int. J. Mod. Phys. A, 20 (A11), 6802, 2005, DOI: 10.1142/S0217751X05030132. [NASA ADS] [CrossRef] [Google Scholar]
 Desorgher, L., K. Kudela, E. Flückiger, R. Bütikofer, M. Storini, and V. Kalegaev. Comparison of Earth’s magnetospheric magnetic field models in the context of cosmic ray physics. Acta Geophys., 57 (1), 75–87, 2009, DOI: 10.2478/s1160000800653. [CrossRef] [Google Scholar]
 Dorman, L. Cosmic Rays in the Earth’s Atmosphere Underground, Kluwer Academic Publishers, Dordrecht, ISBN 1402020716, 2004. [Google Scholar]
 Dos Santos Silva, I., B. De Stavola, C. Pizzi, A. Evans, and S. Evans. Cancer incidence in professional flight crew and air traffic control officers: disentangling the effect of occupational versus lifestyle exposures. International Journal of Cancer, 132 (2), 374–384, 2013, DOI: 10.1002/ijc.27612. [CrossRef] [Google Scholar]
 Drozdov, A., A. Grigoriev, and Y. Malyshkin. Assessment of thunderstorm neutron radiation environment at altitudes of aviation flights. J. Geophys. Res. [Space Phys.], 118 (2), 947–955, 2013, DOI: 10.1029/2012JA018302. [CrossRef] [Google Scholar]
 Dwyer, J.R., D.M. Smith, M.A. Uman, Z. Saleh, B. Grefenstette, B. Hazelton, and H.K. Rassoul. Estimation of the fluence of highenergy electron bursts produced by thunderclouds and the resulting radiation doses received in aircraft. J. Geophys. Res. [Atmos.], 115 (D9), 2010, DOI: 10.1029/2009JD012039. [Google Scholar]
 EURATOM. Council Directive 96/29/EURATOM of 13 May 1996 laying down basic safety standards for protection of the health of workers and the general public against the dangers arising from ionising radiation. Official Journal of the European Communities, 39 (L159), 1996. [Google Scholar]
 Fasso, A., A. Ferrari, J. Ranft, and P. Sala. FLUKA: a multiparticle transport code. SLACR773 200510, CERN, CERN, Geneva, 2005. [Google Scholar]
 Ferrari, A., M. Pelliccioni, and T. Rancati. Calculation of the radiation environment caused by galactic cosmic rays for determining air crew exposure. Radiat. Prot. Dosim., 93 (2), 101–114, 2001. [Google Scholar]
 Forbush, S. Cosmicray intensity variations during two solar cycles. J. Geophys. Res., 63 (4), 651–669, 1958. [CrossRef] [Google Scholar]
 Gaisser, T.K., and T. Stanev. Cosmic rays. In: K. N. et al., ed., Review of Particle Physics, 269–275, J. Phys. G, 37, 269–275, 2010. [Google Scholar]
 Gleeson, L., and W. Axford. Solar modulation of galactic cosmic rays. Astrophys. J., 154, 1011–1026, 1968. [Google Scholar]
 Hammer, G., M. Blettner, and H. Zeeb. Epidemiological studies of cancer in aircrew. Radiat. Prot. Dosim., 136 (4), 232–239, 2009, DOI: 10.1093/rpd/ncp125. [CrossRef] [Google Scholar]
 ICRP. ICRP Publication 60: 1990 recommendations of the international commission on radiological protection. Ann. ICRP, 21 (1–3), 1991. [Google Scholar]
 ICRP. ICRP Publication 103: the 2007 recommendations of the international commission on radiological protection. Ann. ICRP, 37 (2–4), 2007. [Google Scholar]
 ICRP. ICRP Publication 110: adult reference computational phantoms. Ann. ICRP, 39 (2), 2009. [Google Scholar]
 Iwase, H., K. Niita, and T. Nakamura. Development of generalpurpose particle and heavy ion transport Monte Carlo code. J. Nucl. Sci. Technol., 39 (11), 1142–1151, 2002. [Google Scholar]
 Kawrakow, I. Electron impact ionization cross sections for EGSnrc. Med. Phys., 29, 1230, 2001. [Google Scholar]
 Kovaltsov, G., A. Mishev, and I. Usoskinc. A new model of cosmogenic production of radiocarbon 14C in the atmosphere. Earth Planet. Sci. Lett., 337, 114–120, 2012, DOI: 10.1016/j.epsl.2012.05.036. [Google Scholar]
 Kudela, K., R. Bučik, and P. Bobik. On transmissivity of low energy cosmic rays in disturbed magnetosphere. Adv. Space Res., 42 (7), 1300–1306, 2008, DOI: 10.1016/j.asr.2007.09.033. [Google Scholar]
 Kudela, K., and I. Usoskin. On magnetospheric transmissivity of cosmic rays. Czech. J. Phys., 54 (2), 239–254, 2004, DOI: 10.1023/B:CJOP.0000014405.61950.e5. [Google Scholar]
 Langel, R. Main Field in Geomagnetism. In: Geomagnetism, chap. 1, 249–512, J.A. Jacobs Academic Press, London, 1987. [Google Scholar]
 Lewis, B., L. Bennett, A. Green, A. Butler, M. Desormeaux, F. Kitching, M. McCall, B. Ellaschuk, and M. Pierre. Aircrew dosimetry using the Predictive Code for Aircrew Radiation Exposure (PCAIRE). Radiat. Prot. Dosim., 116 (1–4), 320–326, 2005, DOI: 10.1093/rpd/nci024. [Google Scholar]
 Lilensten, J., and A. Belehaki. Developing the scientific basis for monitoring, modelling and predicting space weather. Acta Geophys., 57 (1), 1–14, 2009, DOI: 10.2478/s1160000800813. [CrossRef] [Google Scholar]
 Lilensten, L., and J. Bornarel. Space Weather, Environment and Societies, Springer, Dordrecht, ISBN 9781402043321, 2009. [Google Scholar]
 Matthiä, D., B. Heber, G. Reitz, L. Sihver, T. Berger, and M. Meier. The ground level event 70 on December 13th, 2006 and related effective doses at aviation altitudes. Radiat. Prot. Dosim., 136 (4), 304–310, 2009, DOI: 10.1093/rpd/ncp141. [Google Scholar]
 Matthiä, D., M. Meier, and G. Reitz. Numerical calculation of the radiation exposure from galactic cosmic rays at aviation altitudes with the PANDOCA core model. Space Weather, 12 (3), 161–171, 2014, DOI: 10.1002/2013SW001022. [CrossRef] [Google Scholar]
 Matthiä, D., L. Sihver, and M. Meier. MonteCarlo calculations of particle fluences and neutron effective dose rates in the atmosphere. Radiat. Prot. Dosim., 131 (2), 222–228, 2008, DOI: 10.1093/rpd/ncn130. [Google Scholar]
 Mavromichalaki, H., A. Papaioannou, C. Plainaki, C. Sarlanis, G. Souvatzoglou, et al. Applications and usage of the realtime Neutron Monitor Database. Adv. Space Res., 47, 2210–2222, 2011, DOI: 10.1016/j.asr.2010.02.019. [Google Scholar]
 McCracken, K., V. Rao, B. Fowler, M. Shea, and D. Smart. Cosmic ray tables (asymptotic directions, etc.). In: Annals of the IQSY, chap. 1, 198–214, MIT Press, Cambridge, MA, USA, 1968. [Google Scholar]
 McCracken, K., V. Rao, and M. Shea. The trajectories of cosmic rays in a high degree simulation of the geomagnetic field. Technical Report 77, Massachusetts Institute of Technology, Cambridge, MA, USA, 1962. [Google Scholar]
 Meier, M., M. Hubiak, D. Matthiä, M. Wirtz, and G. Reitz. Dosimetry at aviation altitudes (2006–2008). Radiat. Prot. Dosim., 136 (4), 1–35, 2009, DOI: 10.1093/rpd/ncp142. [Google Scholar]
 Menzel, H. The international commission on radiation units and measurements. Journal of the ICRU, 10 (2), 1–35, 2010. [CrossRef] [Google Scholar]
 Mertens, C., M. Meier, S. Brown, R. Norman, and X. Xu. NAIRAS aircraft radiation model development, dose climatology, and initial validation. Space Weather, 11 (10), 603–635, 2013, DOI: 10.1002/swe.20100. [CrossRef] [Google Scholar]
 Mishev, A. Computation of radiation environment during ground level enhancements 65, 69 and 70 at equatorial region and flight altitudes. Adv. Space Res., 54 (3), 528–535, 2014, DOI: 10.1016/j.asr.2013.10.010. [CrossRef] [Google Scholar]
 Mishev, A., F. Adibpour, I. Usoskin, and E. Felsberger. Computation of dose rate at flight altitudes during ground level enhancements no. 69, 70 and 71. Adv. Space Res., 55 (1), 354–362, 2015, DOI: 10.1016/j.asr.2014.06.020. [Google Scholar]
 Mishev, A., and E. Hristova. Recent gamma background measurements at high mountain altitude. J. Environ. Radioact., 113, 77–82, 2012, DOI: 10.1016/j.jenvrad.2012.04.017. [CrossRef] [Google Scholar]
 Mishev, A., L. Kocharov, and I. Usoskin. Analysis of the ground level enhancement on 17 May 2012 using data from the global neutron monitor network. J. Geophys. Res., 119, 670–679, 2014, DOI: 10.1002/2013JA019253. [Google Scholar]
 Mishev, A., and I. Usoskin. Computations of cosmic ray propagation in the Earth’s atmosphere, towards a GLE analysis. J. Phys: Conf. Ser., 409, 012152, 2013. [Google Scholar]
 Mishev, A., I. Usoskin, and G. Kovaltsov. Neutron monitor yield function: new improved computations. J. Geophys. Res., 118, 2783–2788, 2013, DOI: 10.1002/jgra.50325. [Google Scholar]
 Mishev, A., and P. Velinov. Influence of hadron and atmospheric models on computation of cosmic ray ionization in the atmosphereExtension to heavy nuclei. J. Atmos. Sol. Terr. Phys., 120, 111–120, 2014, DOI: 10.1016/j.jastp. 2014.09.007. [CrossRef] [Google Scholar]
 Mishev, A., and P.I. Velinov. Normalized ionization yield function for various nuclei obtained with full Monte Carlo simulations. Adv. Space Res., 48 (1), 19–24, 2011, DOI: 10.1016/j.asr.2011.02.008. [Google Scholar]
 Nevalainen, J., I. Usoskin, and A. Mishev. Eccentric dipole approximation of the geomagnetic field: application to cosmic ray computations. Adv. Space Res., 52 (1), 22–29, 2013, DOI: 10.1016/j.asr.2013.02.020. [Google Scholar]
 O’Brien, K., W. Friedberg, H. Sauer, and D. Smart. Atmospheric cosmic rays and solar energetic particles at aircraft altitudes. Environ. Int., 22 (Suppl. 1), S9–S44, 1997. [CrossRef] [Google Scholar]
 Pelliccioni, M. Overview of fluencetoeffective dose and fluencetoambient dose equivalent conversion coefficients for high energy radiation calculated using the FLUKA Code. Radiat. Prot. Dosim., 88 (4), 279–297, 2000. [Google Scholar]
 PetoussiHenss, N., W. Bolch, K. Eckerman, A. Endo, N. Hertel, J. Hunt, M. Pelliccioni, H. Schlattl, and M. Zankl. Conversion coefficients for radiological protection quantities for external radiation exposures. Ann. ICRP, 40 (2–5), 1–257, 2010. [Google Scholar]
 Picone, J., A. Hedin, D. Drob, and A. Aikin. NRLMSISE00 empirical model of the atmosphere: statistical comparisons and scientific issues. J. Geophys. Res., 107 (A12), 2002. [Google Scholar]
 Pukkala, E., M. Helminen, T. Haldorsen, N. Hammar, K. Kojo, A. Linnersj, V. Rafnsson, H. Tulinius, U. Tveten, and A. Auvinen. Cancer incidence among Nordic airline cabin crew. International Journal of Cancer, 131 (12), 2886–2897, 2012, DOI: 10.1002/ijc.27551. [CrossRef] [Google Scholar]
 Reames, D. Particle acceleration at the Sun and in the heliosphere. Space Sci. Rev., 90 (3–4), 13–491, 1999. [Google Scholar]
 Roesler, S., W. Heinrich, and H. Schraube. Monte Carlo calculation of the radiation field at aircraft altitudes. Radiat. Prot. Dosim., 98 (4), 367–388, 2002. [Google Scholar]
 Sato, T., A. Endo, M. Zankl, N. PetoussiHenss, H. Yasuda, and K. Niita. Fluencetodose conversion coefficients for aircrew dosimetry based on the new ICRP Recommendations. Progress in Nuclear Science and Technology, 1 , 134–137, 2011. [Google Scholar]
 Sato, T., H. Yasuda, K. Niita, A. Endo, and L. Sihver. Development of PARMA: PHITSbased analytical radiation model in the atmosphere. Radiat. Res., 170, 244–259, 2008, DOI: 10.1667/RR1094.1. [Google Scholar]
 Schraube, H., G. Leuthold, W. Heinrich, S. Roesler, and D. Combecher. European program package for the calculation of aviation route doses, version 3.0. Tech. Rep. D85758, National Research Center for Environment and Health Institute of Radiation Protection, Neuherberg, Germany, 2000. [Google Scholar]
 Shea, M., and D. Smart. Possible evidence for a rigiditydependent release of relativistic protons from the solar corona. Space Sci. Rev., 32, 251–271, 1982. [Google Scholar]
 Shea, M., and D. Smart. A summary of major solar proton events. Sol. Phys., 127, 297–320, 1990. [Google Scholar]
 Shea, M., and D. Smart. Cosmic ray implications for human health. Space Sci. Rev., 93 (1–2), 187–205, 2000, DOI: 10.1023/A:1026544528473. [Google Scholar]
 Shea, M., and D. Smart. Space weather and the groundlevel solar proton events of the 23rd solar cycle. Space Sci. Rev., 171, 161–188, 2012, DOI: 10.1007/s112140129923z. [CrossRef] [Google Scholar]
 Shea, M., D. Smart, and K. McCracken. A study of vertical cut off rigidities using sixth degree simulations of the geomagnetic field. J. Geophys. Res., 70, 4117–4130, 1965. [CrossRef] [Google Scholar]
 Sigurdson, A., and E. Ron. Cosmic radiation exposure and cancer risk among flight crew. Cancer Investigation, 22 (5), 743–761, 2004, DOI: 10.1081/CNV200032767. [CrossRef] [Google Scholar]
 Smart, D., M. Shea, and E. Flückiger. Magnetospheric models and trajectory computations. Space Sci. Rev., 93 (1), 305–333, 2000. [NASA ADS] [CrossRef] [Google Scholar]
 Spurny, F., I. Votockova, and J. BottollierDepois. Geographical influence on the radiation exposure of an aircrew on board a subsonic aircraft. Radioprotection, 31 (2), 275–280, 1996. [Google Scholar]
 Spurny, F., T. Dachev, and K. Kudela. Increase of onboard aircraft exposure level during a solar flare. Nuclear Energy Safety, 10 (48), 396–400, 2002. [Google Scholar]
 Tsyganenko, N. A magnetospheric magnetic field model with a warped tail current sheet. Planet. Space Sci., 37 (1), 5–20, 1989. [Google Scholar]
 Tsyganenko, N., and M. Sitnov. Modeling the dynamics of the inner magnetosphere during strong geomagnetic storms. J. Geophys. Res. [Space Phys.], 110 (A3), 2005, DOI: 10.1029/2004JA010798. [Google Scholar]
 Usoskin, I., K. AlankoHuotari, G. Kovaltsov, and K. Mursula. Heliospheric modulation of cosmic rays: monthly reconstruction for 1951–2004. J. Geophys. Res., 110 (A12108), 2005, DOI: 10.1029/2005JA01125. [Google Scholar]
 Usoskin, I., G. Bazilevskaya, and G.A. Kovaltsov. Solar modulation parameter for cosmic rays since 1936 reconstructed from groundbased neutron monitors and ionization chambers. J. Geophys. Res., 116 (A02), 104, 2011, DOI: 10.1029/2010JA016105. [Google Scholar]
 Usoskin, I., and G. Kovaltsov. Cosmic ray induced ionization in the atmosphere: full modeling and practical applications. J. Geophys. Res., 111 (D21206), 2006, DOI: 10.1029/2006JD007150. [Google Scholar]
 Usoskin, I.G., L. Desorgher, P. Velinov, M. Storini, E. Flückiger, R. Bütikofer, and G. Kovaltsov. Ionization of the Earth’s atmosphere by solar and galactic cosmic rays. Acta Geophys., 57 (1), 88–101, 2009, DOI: 10.2478/s1160000800199. [NASA ADS] [CrossRef] [Google Scholar]
 Vainio, R., L. Desorgher, D. Heynderickx, M. Storini, E. Flückiger, et al. Dynamics of the Earth’s particle radiation environment. Space Sci. Rev., 147 (3–4), 187–231, 2009, DOI: 10.1007/s1121400994967. [Google Scholar]
 Vashenyuk, E., Y. Balabin, B. Gvozdevsky, and L. Schur. Characteristics of relativistic solar cosmic rays during the event of December 13, 2006. Geomag. Aeron., 48 (2), 149–153, 2008, DOI: 10.1007/s1147800820036. [Google Scholar]
 Vashenyuk, E., Y. Balabin, J. PerezPeraza, A. GallegosCruz, and L. Miroshnichenko. Some features of the sources of relativistic particles at the Sun in the solar cycles 21–23. Adv. Space Res., 38 (3), 411–417, 2006, DOI: 10.1016/j.asr.2005.05.012. [Google Scholar]
 Waters, L., G. McKinney, J. Durkee, M. Fensin, J. Hendricks, M. James, R. Johns, and D. Pelowitz. The MCNPX Monte Carlo radiation transport code. AIP Conference Proceedings, 896 (1), 81–90, 2007, DOI: 10.1063/1.2720459. [CrossRef] [Google Scholar]
 Wolf, G., G. Obe, and L. Bergau. Cytogenetic investigations in flight personnel. Radiat. Prot. Dosim., 86 (4), 275–278, 1999. [CrossRef] [Google Scholar]
 Yong, L., A. Sigurdson, E. Ward, M. Waters, E. Whelan, M. Petersen, P. Bhatti, M. Ramsey, E. Ron, and J. Tucker. Increased frequency of chromosome translocations in airline pilots with longterm flying experience. Occupational and Environmental Medicine, 66 (1), 56–62, 2009, DOI: 10.1136/oem.2008.038901. [CrossRef] [Google Scholar]
Cite this article as: Mishev A & Usoskin I. Numerical model for computation of effective and ambient dose equivalent at flight altitudes. J. Space Weather Space Clim., 5, A10, 2015, DOI: 10.1051/swsc/2015011.
All Tables
Effective dose yield function Y [Sv cm^{2} sr nucleon^{−1}] for primary proton and αparticles at the altitudes of 35 kft and 50 kft a.s.l.
Ambient dose equivalent H*(10) yield function Y* [Sv cm^{2} sr nucleon^{−1}] for primary protons and αparticles at the altitude of 35 kft a.s.l.
Effective dose yield function Y [Sv cm^{2} sr nucleon^{−1}] for primary protons with various incidence at the altitude of 35 kft a.s.l.
Neutron monitors with corresponding geomagnetic rigidity cutoffs used in the analysis.
All Figures
Fig. 1.
Effective dose yield function as a function of the energy per nucleon for primary CR protons and αparticles at two different altitudes a.s.l., namely 35 kft (typical commercial flight ≈ 10.5 km) and 50 kft (supersonic flight ≈ 15 km). 

In the text 
Fig. 2.
Ambient dose equivalent H*(10) at the altitude of 35 kft a.s.l. as a function of the rigidity cutoff computed at various conditions, compared with reference data Menzel (2010). (a) Ambient dose equivalent H*(10) for January 1998. Modulation potential 427 MV, (b) ambient dose equivalent H*(10) for January 2000. Modulation potential 752 MV, (c) ambient dose equivalent H*(10) for January 2002. Modulation potential 977 MV. 

In the text 
Fig. 3.
Effective dose rate at the altitude of 35 kft a.s.l. as a function of the rigidity cutoff computed at various conditions. (a) January 1998, modulation potential 427 MV, (b) January 2000, modulation potential 752 MV, (c) January 2002, modulation potential 977 MV. 

In the text 
Fig. 4.
The effective dose rate (μSv h^{−1}) at the altitude of 35 kft a.s.l. during the GLE 70 on 13 December 2006 in a region with R _{ c } ≤ 1 GV. (a) initial phase of the event, (b) main phase of the event, (c) late phase of the event. 

In the text 
Fig. 5.
The relative difference in % of the effective dose rate at the altitude of 35 kft a.s.l. during the GLE 70 on 13 December 2006 assuming different spectral and angular characteristics of SEPs. (a) initial phase of the event, (b) main phase of the event. 

In the text 
Fig. A.1.
Relative difference between computed ambient dose equivalent H*(10) at the altitude of 35 kft a.s.l. and the reference data (Menzel 2010). (a) January 1998, modulation potential 427 MV, (b) January 2000, modulation potential 752 MV, (c) January 2002, modulation potential 977 MV. The modulation potential is according to Usoskin et al. (2011). 

In the text 
Fig. A.2.
Computed effective dose E, ambient dose equivalent H*(10) and reference data (Menzel 2010) at the altitude of 35 kft a.s.l. (a) January 1998, (b) January 2000, (c) January 2002. 

In the text 
Fig. A.3.
Relative difference between computed effective dose E and ambient dose equivalent H*(10) at the altitude of 35 kft a.s.l. (a) January 1998, (b) January 2000, (c) January 2002. 

In the text 
Fig. A.4.
Relative difference between computed effective dose E at the altitude of 35 kft a.s.l. and reference data (Menzel 2010). (a) January 1998, (b) January 2000, (c) January 2002. 

In the text 
Fig. A.5.
Computed effective dose rate E at the altitudes of 50 kft a.s.l. and 35 kft a.s.l. (a) January 1998, (b) January 2000, (c) January 2002. 

In the text 
Fig. B.1.
Calculated NM asymptotic directions during GLE 70 on 13 December 2006 at 03:00 UT. The cross represents the direction of heliospheric magnetic field (IMF) derived from ACE satellite measurements at 03:00 UT. The small circle represents the derived apparent source position. The lines of equal pitch angles relative to the derived anisotropy axis are plotted for 30°, 60°, 150° and 120°. The asymptotic directions of polar NMs are plotted with solid lines, while midlatitude NMs are plotted with dashed lines. 

In the text 
Fig. B.2.
Derived apparent source position (red line) and heliospheric magnetic field (crosses) from ACE satellite measurements throughout GLE 70 on 13 December 2006. 

In the text 
Fig. B.3.
The derived spectral and angular characteristics of SEPs for GLE 70 on 13 December 2006 for 03:00 UT throughout 06:00 UT as denoted in the legend. 

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