Monte Carlo calculation of the radiation field at aircraft altitudes

Energy spectra of secondary cosmic rays are calculated for aircraft altitudes and a discrete set of solar modulation parameters and rigidity cut-off values covering all possible conditions. The calculations are based on the Monte Carlo code FLUKA and on the most recent information on the interstellar cosmic ray flux including a detailed model of solar modulation. Results are compared to a large variety of experimental data obtained on the ground and aboard aircraft and balloons, such as neutron, proton, and muon spectra and yields of charged particles. Furthermore, particle fluence is converted into ambient dose equivalent and effective dose and the dependence of these quantities on height above sea level, solar modulation, and geographical location is studied. Finally, calculated dose equivalent is compared to results of comprehensive measurements performed aboard aircraft.     

Averaged particle dose conversion coefficients in air crew dosimetry

The MCNPX Monte Carlo code was used to calculate energy-dependent fluence-to-effective dose conversion coefficients for neutrons, protons, electrons, photons, charged pions and muons. The FLUKA Monte Carlo code was used to calculate the spectral particle fluences of secondary cosmic rays for different altitudes, and for different combinations of solar modulation and vertical cut-off rigidity parameters. The energy-averaged fluence-to-dose conversion coefficients were obtained by folding the particle fluence spectra with the conversion coefficients for effective dose and ambient dose equivalent. They show a slight dependence on altitude, solar activity and location in the geomagnetic field.     

The influence of the type of filling gas on the response of ionisation chambers to a mixed high-energy radiation field

Radiation protection dosimetry in radiation fields behind the shielding of high-energy accelerators such as CERN is a challenging task and the quantitative understanding of the detector response used for dosimetry is essential. Measurements with ionisation chambers are a standard method to determine absorbed dose (in the detector material). For applications in mixed radiation fields, ionisation chambers are often also calibrated in terms of ambient dose equivalent at conventional reference radiation fields. The response of a given ionisation chamber to the various particle types of a complex high-energy radiation field in terms of ambient dose equivalent depends of course on the materials used for the construction and the chamber gas used. This paper will present results of computational studies simulating the exposure of high-pressure ionisation chambers filled with different types of gases to the radiation field at CERN's CERN-EU high-energy reference field facility. At this facility complex high-energy radiation fields, similar to those produced by cosmic rays at flight altitudes, are produced. The particle fluence and spectra calculated with FLUKA Monte Carlo simulations have been benchmarked in several measurements. The results can be used to optimise the response of ionisation chambers for the measurement of ambient dose equivalent in high-energy mixed radiation fields.     

Calculation of the radiation environment caused by galactic cosmic rays for determining air crew exposure

The spectra of secondary particles resulting from interactions of primary galactic cosmic rays with the nuclei in the atmosphere have been calculated using the Monte Carlo transport code FLUKA. The simulations have been carried out at solar minimum and solar maximum activity, for several values of the vertical geomagnetic cut-off. The effective dose rate and the ambient dose equivalent rate as a function of geomagnetic cut-off and altitude have been obtained using appropriate sets of conversion coefficients. The calculated results are discussed and compared with experimental data and other calculations. A simple method is proposed to calculate the radiation exposure at aircraft altitudes.     

FLUKA simulation of TEPC response to cosmic radiation

The aircrew exposure to cosmic radiation can be assessed by calculation with codes validated by measurements. However, the relationship between doses in the free atmosphere, as calculated by the codes and from results of measurements performed within the aircraft, is still unclear. The response of a tissue-equivalent proportional counter (TEPC) has already been simulated successfully by the Monte Carlo transport code FLUKA. Absorbed dose rate and ambient dose equivalent rate distributions as functions of lineal energy have been simulated for several reference sources and mixed radiation fields. The agreement between simulation and measurements has been well demonstrated. In order to evaluate the influence of aircraft structures on aircrew exposure assessment, the response of TEPC in the free atmosphere and on-board is now simulated. The calculated results are discussed and compared with other calculations and measurements.     

Galactic and solar radiation exposure to aircrew during a solar cycle

An on-going investigation using a tissue-equivalent proportional counter (TEPC) has been carried out to measure the ambient dose equivalent rate of the cosmic radiation exposure of aircrew during a solar cycle. A semi-empirical model has been derived from these data to allow for the interpolation of the dose rate for any global position. The model has been extended to an altitude of up to 32 km with further measurements made on board aircraft and several balloon flights. The effects of changing solar modulation during the solar cycle are characterised by correlating the dose rate data to different solar potential models. Through integration of the dose-rate function over a great circle flight path or between given waypoints, a Predictive Code for Aircrew Radiation Exposure (PCAIRE) has been further developed for estimation of the route dose from galactic cosmic radiation exposure. This estimate is provided in units of ambient dose equivalent as well as effective dose, based on E/H x (10) scaling functions as determined from transport code calculations with LUIN and FLUKA. This experimentally based treatment has also been compared with the CARI-6 and EPCARD codes that are derived solely from theoretical transport calculations. Using TEPC measurements taken aboard the International Space Station, ground based neutron monitoring, GOES satellite data and transport code analysis, an empirical model has been further proposed for estimation of aircrew exposure during solar particle events. This model has been compared to results obtained during recent solar flare events.     

Dosimetric considerations on TEPC fluka-simulation and measurements

The response of a tissue equivalent proportional counter (TEPC) has been simulated with the Monte Carlo transport code FLUKA. The absorbed dose distribution of lineal energy y has been determined for several monoenergetic photon and neutron sources. The agreement between the calculated results and the measurements carried out with different well-known sources is well demonstrated. Work is in progress in order to evaluate the response of the instrument in the cosmic ray environment.     

Simulation studies on a prototype ionisation chamber for measurement of personal dose equivalent, Hp(10)

A prototype ionisation chamber for direct measurement of the personal dose equivalent, Hp(10), similar to the one developed by the Physikalisch-Technische Bundesantalt (PTB), was designed and constructed by the Metrological Laboratory of Ionizing Radiation (LMRI) of Nuclear and Technological Institute (ITN). Tests already performed have shown that the behaviour of this chamber is very similar to the PTB chamber, mainly the energy dependence for the X-ray radiation qualities of the ISO 4037-1 narrow series N-30, N-40, N-60, N-80, N-100 and N-120 and also for gamma radiation of 137Cs and 60Co. However, the results obtained also show a dependence on the energy and angles of incident radiation and a low magnitude of the electrical response of the ionisation chamber. In order to optimise the performance of the chamber, the LMRI initiated numerical simulation of this ionisation chamber by Monte Carlo method using the MCNPX code.     

A mathematical model of aircraft for evaluating the effects of shielding structure on aircrew exposure

To investigate the influence of the aircraft structures and contents on the exposure of aircrew to the galactic component of cosmic rays, a mathematical model of an aeroplane has been developed. The irradiation of the mathematical model in the cosmic ray environment has been simulated using the Monte Carlo transport code FLUKA. Effective dose andambient dose-equivalent rates have been determined inside the aircraft at several locations along the fuselage at a typicaI civil aviation altitude. A significant effect of the shielding of aircraft structures has been observed on the ambient dose-equivalent rates, while the impact on the effective dose rates seems to be minor. Care should be taken in positioning the detectors onboard when the measurements are aimed at validating the codes.     

Study of the ratio of non-neutron to neutron dose components of cosmic radiation at typical commercial flight altitudes

CIEMAT, in close co-operation with Iberia Airlines, carried out an extensive programme of in-flight measurements, covering both hemispheres, during the years 2001 and 2002. Although the instrumentation onboard included different active devices, the results presented here were obtained from a polyethylene/tungsten-moderated rem meter (SWENDI2; Eberline) and an ionisation chamber (RSS-131; Reuter-Stokes) used for measuring the ambient dose equivalent due to the neutron and the non-neutron components of cosmic radiation, respectively. This paper presents a study of each of the dose components mentioned as a function of the vertical cut-off rigidity and the flight altitude. The ratio between the two components is also presented to determine the variations in cosmic radiation composition as a function of the aforementioned parameters. The experimental results have also been compared with those predicted by the code EPCARD3.2 for the non-neutron and the neutron components of the ambient dose equivalent.     

Field calibration studies for ionisation chambers in mixed high-energy radiation fields

The monitoring of ambient doses at work places around high-energy accelerators is a challenging task due the complexity of the mixed stray radiation fields encountered. At CERN, mainly Centronics IG5 high-pressure ionisation chambers are used to monitor radiation exposure in mixed fields. The monitors are calibrated in the operational quantity ambient dose equivalent H*(10) using standard, source-generated photon- and neutron fields. However, the relationship between ionisation chamber reading and ambient dose equivalent in a mixed high-energy radiation field can only be assessed if the spectral response to every component and the field composition is known. Therefore, comprehensive studies were performed at the CERN-EU high-energy reference field facility where the spectral fluence for each particle type has been assessed with Monte Carlo simulations. Moreover, studies have been performed in an accessible controlled radiation area in the vicinity of a beam loss point of CERN's proton synchrotron. The comparison of measurements and calculations has shown reasonable agreement for most exposure conditions. The results indicate that conventionally calibrated ionisation chambers can give satisfactory response in terms of ambient dose equivalent in stray radiation fields at high-energy accelerators in many cases. These studies are one step towards establishing a method of 'field calibration' of radiation protection instruments in which Monte Carlo simulations will be used to establish a correct correlation between the response of specific detectors to a given high-energy radiation field.     

New calculations of the atmospheric cosmic radiation field--results for neutron spectra

The propagation of primary cosmic rays through the Earth's atmosphere and the energy spectra of the resulting secondary particles have been calculated using the Monte Carlo transport code FLUKA with several novel auxiliary methods. Solar-modulated primary cosmic ray spectra were determined through an analysis of simultaneous proton and helium measurements made on spacecraft or high-altitude balloon flights. Primary protons and helium ions are generated within the rigidity range of 0.5 GV-20 TV, uniform in cos2theta. For a given location, primaries above the effective angle-dependent geomagnetic cut-off rigidity, and re-entrant albedo protons, are transported through the atmosphere. Helium ions are initially transported using a separate transport code called HEAVY to simulate fragmentation. HEAVY interfaces with FLUKA to provide interaction starting points for each nucleon originating from a helium nucleus. Calculated cosmic ray neutron spectra and consequent dosimetric quantities for locations with a wide range of altitude (atmospheric depth) and geomagnetic cut-off are presented and compared with measurements made on a high-altitude aeroplane. Helium ion propagation using HEAVY and inclusion of re-entrant albedo protons with the incident primary spectra significantly improved the agreement of the calculated cosmic ray neutron spectra with measured spectra. These cosmic ray propagation calculations provide the basis for a new atmospheric ionising radiation (AIR) model for air-crew dosimetry, calculation of effects on microelectronics, production of cosmogenic radionuclides and other uses.     

The results of cosmic radiation in-flight TEPC measurements during the CAATER flight campaign and comparison with simulation

The European-Commission-supported project DOSMAX (Dosimetry of Aircrew Exposure to Radiation During Solar Maximum) was aimed at measuring aircrew exposure to cosmic radiation on-board the aircraft during solar maximum. During a dedicated international comparison mission (Co-ordinated Access to Aircraft for Transnational Environmental Research; CAATER) different measurement techniques have been compared by six European institutes (Results of the CAATER Mission, DOSMAX Meeting, Dublin, June 2004). In this paper, we present the tissue-equivalent proportional counter (TEPC) measurements carried out by ARC Seibersdorf research (ARCS), Austria, and Institut de Radioprotection et de S?reté Nucléaire (IRSN), France, together with a comparison with simulation results under the same conditions. The whole flight campaign consists of four different in-flight investigations performed at two different geographical positions at 12.2 km (FL 400) and 9.8 km (FL 320). One location was chosen above Rome (42 degrees North, 12 degrees East), Italy, for high cut-off rigidity (6.4 GV) and the second above Aalborg (57 degrees North, 10 degrees East), Denmark, for low cut-off rigidity (1.8 GV). The TEPC measurements are presented in terms of absorbed dose and ambient dose equivalent as well as microdosimetric spectra as a function of lineal energy. For the same conditions of the CAATER flights the response of the TEPC has also been simulated by using the Monte Carlo Transport Code FLUKA (version 2003). The results from simulations are compared with measurements and they show a reasonable agreement.     

A method applicable to effective dose rate estimates for aircrew dosimetry

The inclusion of cosmic radiation as occupational exposure under ICRP Publication 60 and the European Union Council Directive 96/29/Euratom has highlighted the need to estimate the exposure of aircrew. According to a report of the Group of Experts established under the terms of Article 31 of the European Treaty, the individual estimates of dose for flights below 15 km may be done using an appropriate computer program. In order to calculate the radiation exposure at aircraft altitudes, calculations have been performed by means of the Monte Carlo transport code FLUKA. On the basis of the calculated results, a simple method is proposed for the individual evaluation of effective dose rate due to the galactic component of cosmic radiation as a function of latitude and altitude.     

Modelling of aircrew radiation exposure from galactic cosmic rays and solar particle events

Correlations have been developed for implementation into the semi-empirical Predictive Code for Aircrew Radiation Exposure (PCAIRE) to account for effects of extremum conditions of solar modulation and low altitude based on transport code calculations. An improved solar modulation model, as proposed by NASA, has been further adopted to interpolate between the bounding correlations for solar modulation. The conversion ratio of effective dose to ambient dose equivalent, as applied to the PCAIRE calculation (based on measurements) for the legal regulation of aircrew exposure, was re-evaluated in this work to take into consideration new ICRP-92 radiation-weighting factors and different possible irradiation geometries of the source cosmic-radiation field. A computational analysis with Monte Carlo N-Particle eXtended Code was further used to estimate additional aircrew exposure that may result from sporadic solar energetic particle events considering real-time monitoring by the Geosynchronous Operational Environmental Satellite. These predictions were compared with the ambient dose equivalent rates measured on-board an aircraft and to count rate data observed at various ground-level neutron monitors.     

Monte Carlo calculation of the angular distribution of cosmic rays at flight altitudes

The angular distribution of the secondary radiation produced by the galactic component of cosmic rays has been determined by simulating the penetration of the primary spectra in the Earth's atmosphere. The simulations have been carried out with the latest version of the FLUKA code. Particles have been scored at various altitudes according to their angle of incidence for some significant values of vertical cut-off rigidity and solar modulation parameter. The calculated results at typical cruise altitudes for a civil aircraft are presented. The data at 10.7 km have been fitted with simple mathematical equations. It has been demonstrated that the major contribution to the doses at aviation altitudes arises from downward-directed particles. The isotropic irradiation usually assumed for the evaluation of aircrew exposure could be a very poor approximation.     

Variations observed in environmental radiation at ground level

To investigate and monitor environmental radiation at ground level, Physikalisch-Technische Bundesanstalt (PTB) has installed several dosemeters and particle detectors at the new Ambient Radiation Dosimetry Site. The separation of the total ambient dose equivalent rate H*10(env) of environmental radiation into the different contributions is achieved by comparing the data of different detectors: the muon detector MUDOS, a modified neutron dosemeter, proportional counters and ionisation chambers. The response of the latter two dosemeter systems to cosmic radiation was determined at the Cosmic Radiation Dosimetry Site on a lake near PTB. Besides the increase of the ambient dose equivalent rate during rainfall, variations owing to air pressure, solar activity and temperature changes in the upper atmosphere are observed. Without rain and solar effects, smooth variations of the cosmic component at ground level of +/-6.9 nSv h(-1) should be treated as naturally occurring variations during an entire year.     

Aircrew exposure from cosmic radiation on commercial airline routes

As a result of the recent recommendations of the ICRP 60, and in anticipation of possible regulation on occupational exposure of Canadian-based aircrew, an extensive study was carried out by the Royal Military College of Canada over a one-year period to measure the cosmic radiation at commercial jet altitudes. A tissue-equivalent proportional counter was used to measure the ambient total dose equivalent rate on 62 flight routes, resulting in over 20,000 data points at one-minute intervals at various altitudes and geomagnetic latitudes (i.e. which span the full cut-off rigidity of the Earth's magnetic field). These data were then compared to similar experimental work at the Physikalisch Technische Bundesanstalt, using a different suite of equipment, to measure separately the low and high linear energy transfer components of the mixed radiation field, and to predictions with the LUIN transport code. All experimental and theoretical results were in excellent agreement. From these data, a semiempirical model was developed to allow for the interpolation of the dose rate for any global position, altitude and date (i.e. heliocentric potential). Through integration of the dose rate function over a great circle flight path, a computer code was developed to provide an estimate of the total dose equivalent on any route worldwide at any period in the solar cycle.     

Shielding variation effects for 250 MeV protons on tissue targets

This paper provides results of computer simulation studies with the goal to analyse issues regarding radiation protection for personnel, patients and third persons involved in hadron therapy treatment. The treatment room and the patient are modelled by simple cylindrical geometries at incident proton energies of 250 MeV. Monte Carlo simulations of the energy and angular dependence of proton, neutron and photon radiation fields and resulting ambient dose equivalent distributions outside the shielding walls are performed. In order to investigate systematic uncertainties due to the shielding materials and inherent to the computer models, various concrete compositions, densities and water contents are modelled, and the influence of simulation parameters on the results obtained is determined. Generally, good agreement is found between results provided by MCNPX and FLUKA computer codes. Variations in neutron ambient dose attenuation from -50 to +/-30% are found due to varying concrete composition. Changes in the water content of the concrete in the order of 8% may cause variations up to 20%.