Background information about polonium-210 – considerations on its biokinetics and internal dosimetry

The mysterious death of the former secret service agent Alexander Litvinenko on 23 November 2006 suddenly attracted the attention of the public to the rather unknown radionuclide 210Po. In this respect, this paper presents some general background information about this radionuclide as well as some details about its biokinetic behaviour and internal dosimetry. The results are discussed in the context of the Litvinenko case with regard to possible health effects on people from Germany being potentially contaminated. However, no evidence for contamination with 210Po could be found for any people tested at the GSF-National Research Center for Environment and Health (since 01.01.2008 Helmholtz Zentrum München).

Differences in exhaled gas profiles between patients with type 2 diabetes and healthy controls

AIMS

Recent advances in analytical technology allow the detection of several hundred volatile organic compounds (VOCs) in human exhaled air, many of which reflect unidentified endogenous pathways. This study was performed to determine whether a breath gas analysis using proton transfer reaction-mass spectrometry (PTR-MS) could serve as a noninvasive method to distinguish between patients with type 2 diabetes mellitus and healthy controls.

METHODS

Breath and room air samples were measured from 21 patients with insulin-treated type 2 diabetes and 26 healthy controls. VOCs in the mass range of 20-200 atomic mass units were analyzed using PTR-MS.

RESULTS

We identified eight masses characteristic of endogenous VOCs that showed significant differences in the gas profiles of patients with type 2 diabetes and healthy control subjects. Using these VOCs for linear discriminant analysis, the sensitivity and specificity were found to be 90% and 92%, respectively.

CONCLUSIONS

These results suggest that it is possible to separate patients with diabetes mellitus type 2 from healthy controls by multivariate analysis of exhaled endogenous VOCs. This is a first step towards the development of a noninvasive test using breath gas of at-risk persons and making it an attractive option for large-scale testing of at-risk populations. However, the establishment of exhaled volatiles as metabolic markers requires additional confirmatory investigations.

Measurements of secondary neutrons from cosmic radiation with a Bonner sphere spectrometer at 79 degrees N

Air crew members and airline passengers are continuously exposed to cosmic radiation during their flights. Particles ejected by the sun during so-called solar particle events (SPEs) in periods of high solar activity can contribute to this exposure. In rare cases the dose from a single SPE might even exceed the annual dose limit of 1 mSv above which dose monitoring of air crews is legally required in Germany. Measurements performed by means of neutron monitors have already shown that the relative intensity of secondary neutrons from cosmic radiation is enhanced during an SPE, particularly at regions close to the magnetic poles of the Earth where shielding of the cosmic radiation by the geomagnetic field is low. Here we describe a Bonner sphere spectrometer installed at the Koldewey station at 79 degrees N, i.e. about 1,000 km from the geographic North pole, which is designed to provide first experimental data on the time-dependent energy spectrum of neutrons produced in the atmosphere during an SPE. This will be important to calculate doses from these neutrons to air crew members. The system is described in detail and first results are shown that were obtained during quiet periods of sun activity.

Promotion of initiated cells by radiation-induced cell inactivation

Cells on the way to carcinogenesis can have a growth advantage relative to normal cells. It has been hypothesized that a radiation-induced growth advantage of these initiated cells might be induced by an increased cell replacement probability of initiated cells after inactivation of neighboring cells by radiation. Here Monte Carlo simulations extend this hypothesis for larger clones: The effective clonal expansion rate decreases with clone size. This effect is stronger for the two-dimensional than for the three-dimensional situation. The clones are irregular, far from a circular shape. An exposure-rate dependence of the effective clonal expansion rate could come in part from a minimal recovery time of the initiated cells for symmetric cell division.

Comparisons of calculations with PARTRAC and NOREC: transport of electrons in liquid water

Monte Carlo computer models that simulate the detailed, event-by-event transport of electrons in liquid water are valuable for the interpretation and understanding of findings in radiation chemistry and radiation biology. Because of the paucity of experimental data, such efforts must rely on theoretical principles and considerable judgment in their development. Experimental verification of numerical input is possible to only a limited extent. Indirect support for model validity can be gained from a comparison of details between two independently developed computer codes as well as the observable results calculated with them. In this study, we compare the transport properties of electrons in liquid water using two such models, PARTRAC and NOREC. Both use interaction cross sections based on plane-wave Born approximations and a numerical parameterization of the complex dielectric response function for the liquid. The models are described and compared, and their similarities and differences are highlighted. Recent developments in the field are discussed and taken into account. The calculated stopping powers, W values, and slab penetration characteristics are in good agreement with one another and with other independent sources.

Measurement techniques for tracer kinetic studies with stable isotopes of zirconium

Biokinetic models are used in radiation protection to assess internal radiation doses. Experiments with stable isotopes as tracers can be performed to obtain characteristic parameters of these models. Two methods for the measurement of zirconium isotopes in human biological samples are presented--thermal ionisation mass spectrometry (TIMS) and proton nuclear activation analysis (PNA). Descriptions include sample preparation, operating conditions, relative uncertainties and method detection limits as well as important properties of both methods.

Measurements of daily urinary uranium excretion in German peacekeeping personnel and residents of the Kosovo region to assess potential intakes of depleted uranium (DU)

Following the end of the Kosovo conflict, in June 1999, a study was instigated to evaluate whether there was a cause for concern of health risk from depleted uranium (DU) to German peacekeeping personnel serving in the Balkans. In addition, the investigations were extended to residents of Kosovo and southern Serbia, who lived in areas where DU ammunitions were deployed. In order to assess a possible DU intake, both the urinary uranium excretion of volunteer residents and water samples were collected and analysed using inductively coupled plasma-mass spectrometry (ICP-MS). More than 1300 urine samples from peacekeeping personnel and unexposed controls of different genders and age were analysed to determine uranium excretion parameters. The urine measurements for 113 unexposed subjects revealed a daily uranium excretion rate with a geometric mean of 13.9 ng/d (geometric standard deviation (GSD)=2.17). The analysis of 1228 urine samples from the peacekeeping personnel resulted in a geometric mean of 12.8 ng/d (GSD=2.60). It follows that both unexposed controls and peacekeeping personnel excreted similar amounts of uranium. Inter-subject variation in uranium excretion was high and no significant age-specific differences were found. The second part of the study monitored 24 h urine samples provided by selected residents of Kosovo and adjacent regions of Serbia compared to controls from Munich, Germany. Total uranium and isotope ratios were measured in order to determine DU content. (235)U/(238)U ratios were within +/-0.3% of the natural value, and (236)U/(238)U was less than 2 x 10(-7), indicating no significant DU in any of the urine samples provided, despite total uranium excretion being relatively high in some cases. Measurements of ground and tap water samples from regions where DU munitions were deployed did not show any contamination with DU, except in one sample. It is concluded that both peacekeeping personnel and residents serving or living in the Balkans, respectively, were not exposed to significant amounts of DU.

Daily uranium excretion in German peacekeeping personnel serving on the Balkans compared to ICRP model prediction

An investigation was performed to assess a possible health risk of depleted uranium (DU) for residents and German peacekeeping personnel serving on the Balkans. In order to evaluate a possible DU intake, the urinary uranium excretions of volunteers were collected and analysed using inductively coupled plasma mass spectrometry (ICP-MS). In total, more than 1300 urine samples from soldiers, civil servants and unexposed controls of different genders and ages were analysed to determine uranium excretion parameters. All participating volunteers, aged 3-92 y, were grouped according to their gender and age for evaluation. The results of the investigation revealed no significant difference between the unexposed controls and the peacekeeping personnel. In addition, the geometric means of the daily urinary excretion in peacekeeping personnel, ranging from 3 to 23 ng d(-1) for different age groups, fall toward the lower end of renal uranium excretion values published for unexposed populations in literature. The measured data were compared with the International Commission on Radiological Protection prediction for the intake of natural uranium by unexposed members of the public. The two data sets are in good agreement, indicating that no relevant intake of additional uranium, either natural or DU, has appeared for German peacekeeping personnel serving on the Balkans.

Comparisons of 239Pu inhalation doses calculated with ICRP 67 and proposed systemic models

The International Commission on Radiological Protection (ICRP) has issued an age-specific systemic biokinetic model for plutonium (Pu), which was later modified to give better agreement with measured urinary excretion data. Recently, the current ICRP systemic Pu model was improved by Leggett et al. based on recently developed data. Incorporation of 239Pu in the human body may result in significant internal radiation exposure. In the present work, the retentions in organs and tissues, the equivalent dose and effective dose from 239Pu for workers and members of the public were estimated and compared under the current ICRP and the proposed models. 239Pu contents in liver and in other soft tissue calculated with the proposed model are higher than predicted by the ICRP model, whereas bone content is lower than predicted by the ICRP model. Based on the proposed model, the inhalation equivalent dose coefficient in some organs, e.g. liver and kidneys, is increased, but there is no significant change in the effective inhalation dose coefficients of 239Pu for workers and members of the public.

Long-term measurements of cosmic ray neutrons by means of a Bonner spectrometer at mountain altitudes - first results

A Bonner multi-sphere spectrometer has been installed in 2005 at the Environmental Research Station 'Schneefernerhaus' (2660 m above sea level) on the Zugspitze mountain, Germany, to measure the energy spectrum of cosmic-ray neutrons at high altitudes continuously. The system can be used to investigate small temporal variations in the cosmic radiation intensity. For example, measurements were done during periods of 2 Forbush decreases of the cosmic radiation intensity in July and September 2005, respectively. The results were compared with those obtained by using neutron monitors, and neutron fluence spectra measured during these events are presented and discussed.

Lung cancer risk in mice: analysis of fractionation effects and neutron RBE with a biologically motivated model

Data from Argonne National Laboratory on lung cancer in 15,975 mice with acute and fractionated exposures to gamma rays and neutrons are analyzed with a biologically motivated model with two rate-limiting steps and clonal expansion. Fractionation effects and effects of radiation quality can be explained well by the estimated kinetic parameters. Both an initiating and a promoting action of neutrons and gamma rays are suggested. While for gamma rays the initiating event is described well with a linear dose-rate dependence, for neutrons a nonlinear term is needed, with less effectiveness at higher dose rates. For the initiating event, the neutron RBE compared to gamma rays is about 10 when the dose rate during each fraction is low. For higher dose rates this RBE decreases strongly. The estimated lifetime relative risk for radiation-induced lung cancers from 1 Gy of acute gamma-ray exposure at an age of 110 days is 1.27 for male mice and 1.53 for female mice. For doses less than 1 Gy, the effectiveness of fractionated exposure to gamma rays compared to acute exposure is between 0.4 and 0.7 in both sexes. For lifetime relative risk, the RBE from acute neutrons at low doses is estimated at about 10 relative to acute gamma-ray exposure. It decreases strongly with dose. For fractionated neutrons, it is lower, down to about 4 for male mice.

Human exposure to space radiation: role of primary and secondary particles

Human exposure to space radiation implies two kinds of risk, both stochastic and deterministic. Shielding optimisation therefore represents a crucial goal for long-term missions, especially in deep space. In this context, the use of radiation transport codes coupled with anthropomorphic phantoms allows to simulate typical radiation exposures for astronauts behind different shielding, and to calculate doses to different organs. In this work, the FLUKA Monte Carlo code and two phantoms, a mathematical model and a voxel model, were used, taking the Galactic Cosmic Rays (GCR) spectra from the model of Badhwar and O'Neill. The time integral spectral proton fluence of the August 1972 Solar Particle Event (SPE) was represented by an exponential function. For each aluminium shield thickness, besides total doses the contributions from primary and secondary particles for different organs and tissues were calculated separately. More specifically, organ-averaged absorbed doses, dose equivalents and a form of 'biological dose', defined on the basis of initial (clustered) DNA damage, were calculated. As expected, the SPE doses dramatically decreased with increasing shielding, and doses in internal organs were lower than in skin. The contribution of secondary particles to SPE doses was almost negligible; however it is of note that, at high shielding (10 g cm(-2)), most of the secondaries are neutrons. GCR organ doses remained roughly constant with increasing Al shielding. In contrast to SPE results, for the case of cosmic rays, secondary particles accounted for a significant fraction of the total dose.

Role of DNA/chromatin organisation and scavenging capacity in USX- and proton- induced DNA damage

DNA higher-order structures and (non-histonic) *;OH radical scavengers have well known protective effects in the induction of single- and double-strand breaks by ionising radiation. In a previous work, such protective roles have been quantified for gamma radiation (Valota et al., Int. J. Radiat. Biol. 79, 2003). As a starting base for the simulations, we used the PARTRAC Monte Carlo code, developed within a collaboration involving the University of Pavia and the GSF institute. The code can reproduce the track structure of photons, electrons, protons and heavier ions in liquid water, and it can simulate the DNA content of a human cell at different organisation levels, based on an atom-by-atom approach. In this work we extended the calculations to Ultra-Soft X rays (USX) and protons, separately analysing the effects of different radiation types on various DNA structures (i.e. linear DNA, SV40 'minichromosomes' and compact chromatin) as a function of the *OH scavenging capacity (SC). Both for USX and protons, the calculated damage yields decreased by increasing the SC for the three considered target types. Such decrease can be ascribed to the competition between the reactions *OH-DNA and *OH-scavenger, which becomes more and more likely by increasing the SC. Furthermore, linear DNA was found to be more radiosensitive than SV40 'minichromosomes', which in turn were more radiosensitive than compact chromatin, which is protected by histones. Comparisons with experimental data by Fulford et al. (Int. J. Radiat. Biol. 77, 2001) relative to USX irradiation showed very good agreement. The dependence of the modulating role played by DNA organisation and scavenging capacity on radiation quality is presented and discussed.

Bone cancer risk in mice exposed to 224Ra: protraction effects from promotion

This paper analyzes data for the osteosarcoma incidence in life-time experiments of (224)Ra injected mice with respect to the importance of initiating and promoting action of ionizing high LET-radiation. This was done with the biologically motivated two step clonal expansion (TSCE) model of tumor induction. Experimentally derived osteosarcoma incidence in 1,194 mice following exposure to (224)Ra with different total radiation doses and different fractionation patterns were analyzed together with incidence data from 1,710 unirradiated control animals. Effects of radiation on the initiating event and on the clonal expansion rate, i.e. on promotion were found to be necessary to explain the observed patterns with this model. The data show a distinct inverse protraction effect at high doses, whereas at lower doses this effect becomes insignificant. Such a behavior is well reproduced in the proposed model: At dose rates above 6 mGy/day a longer exposure produces higher ERR per dose, while for lower rates the reverse is the case. The TSCE model permits the deduction of several kinetic parameters of a postulated two-step bone tumorigenesis process. Mean exposure rates of 0.13 mGy/day are found to double the baseline initiation rate. At rates above 100 mGy/day, the initiation rate decreases. The clonal expansion rate is doubled at 8 mGy/day, and it levels out at rates beyond 100 mGy/day.

Solubility of airborne radioactive fuel particles from the Chernobyl reactor and implication to dose

Airborne particles of nuclear fuel from the Chernobyl reactor that had been collected on air filters and stored, were characterised using in vitro dissolution tests to assess effective doses after their inhalation. As solvent, the Gamble biological fluid was used to simulate lung fluid. The solubility of the measured radionuclides decreased in the order (137)Cs>(90)Sr>>(241)Am>or=(239+240)Pu in the simulated lung fluid. The dissolution rate constant of e.g. (239+249)Pu ranged from 0.72 to 5.4 x 10(-6) g x cm(-2) d(-1) and decreased (for all nuclides) with increasing particle size as predicted from theoretical considerations. Considering the inhalation dose, decreasing dose with size and increasing doses with lower solubility may partly counterbalance each other for (137)Cs and (90)Sr. On the other hand, for (239)Pu and (241)Am larger particles and associated lower solubility both change the resulting dose in the same direction towards lower values. The comparison of the experimentally determined dose coefficients with ICRP values indicates that nuclear fuel particles closely resemble type M material characteristics for (137)Cs and (90)Sr and type S material characteristics for (239)Pu and (241)Am.

Soil-to-grain transfer of fallout 137Cs for 28 winter wheat cultivars as observed in field experiments

In order to find wheat cultivars with a minimum soil-to-grain transfer of fallout (137)Cs, 28 winter wheat cultivars were investigated at 3 different sites with different soil types in Bavaria, Germany. Each cultivar was grown on an area of 10 m(2) and harvested in August 1999. The soil-to-grain concentration ratios (C(r)) of (137)Cs varied by a factor of up to 3 from cultivar to cultivar at a given site and from site-to-site for a given cultivar. The mean C(r) values at the three sites, 4.2 x 10(-4), 4.9 x 10(-4) and 7.5 x 10(-4), differed significantly. The fact that no cultivar showed similar C(r) values at the three sites indicates a strong influence of the soil on C(r). The cultivars Flair, Kornett and Previa showed a minor uptake of (137)Cs compared with the mean of all cultivars at each site. Unlike (137)Cs, the (40)K concentrations in the wheat grains varied only within a small range (122-190 Bq kg(-1)) at each site, which is due to the potassium regulation by the plants. For both radionuclides, the differences between the root uptake characteristics of the cultivars may not only be explained by an inter-cultivar variability due to genetic differences between the cultivars, but also by an intra-cultivar variability due to different soil conditions.

Role of shielding in modulating the effects of solar particle events: Monte Carlo calculation of absorbed dose and DNA complex lesions in different organs

Distributions of absorbed dose and DNA clustered damage yields in various organs and tissues following the October 1989 solar particle event (SPE) were calculated by coupling the FLUKA Monte Carlo transport code with two anthropomorphic phantoms (a mathematical model and a voxel model), with the main aim of quantifying the role of the shielding features in modulating organ doses. The phantoms, which were assumed to be in deep space, were inserted into a shielding box of variable thickness and material and were irradiated with the proton spectra of the October 1989 event. Average numbers of DNA lesions per cell in different organs were calculated by adopting a technique already tested in previous works, consisting of integrating into "condensed-history" Monte Carlo transport codes--such as FLUKA--yields of radiobiological damage, either calculated with "event-by-event" track structure simulations, or taken from experimental works available in the literature. More specifically, the yields of "Complex Lesions" (or "CL", defined and calculated as a clustered DNA damage in a previous work) per unit dose and DNA mass (CL Gy-1 Da-1) due to the various beam components, including those derived from nuclear interactions with the shielding and the human body, were integrated in FLUKA. This provided spatial distributions of CL/cell yields in different organs, as well as distributions of absorbed doses. The contributions of primary protons and secondary hadrons were calculated separately, and the simulations were repeated for values of Al shielding thickness ranging between 1 and 20 g/cm2. Slight differences were found between the two phantom types. Skin and eye lenses were found to receive larger doses with respect to internal organs; however, shielding was more effective for skin and lenses. Secondary particles arising from nuclear interactions were found to have a minor role, although their relative contribution was found to be larger for the Complex Lesions than for the absorbed dose, due to their higher LET and thus higher biological effectiveness.

Interpretation by modelling of observations in radon radiation carcinogenesis

Biophysical models for radon-related induction of lung cancer are developed with the aim of reducing the uncertainties in current risk estimates at low doses by a better understanding of the relevant mechanisms. These models can make use of the full dosimetric information when extracting information on, say, age-at-exposure, protraction or fractionation effects. It is found that irradiation by radon and its progeny does act on the initiating event of carcinogenesis (e.g. mutation), but its dominating effect is via promoting the division of already initiated cells. Data show that the concept of a unit of exposure giving, in an additive way, a unit of lung cancer risk is too limited, while relatively simple mechanistic assumptions described in this article do yield an adequate description of observations. Exposures in epidemiological data sets are measured with error. For various error models it has been shown that likelihood-based techniques of correction work reliably; likewise for biologically based cancer models. When several parameters are allowed to be exposure dependent, for example, initiation and promotion, then their relative importance is influenced.

Modelling study on the protective role of OH radical scavengers and DNA higher-order structures in induction of single- and double-strand break by gamma-radiation

PURPOSE

To quantify the protective effects of (non-histonic) OH-radical scavengers and DNA higher-order structures in induction of single- (ssbs) and double-strand breaks (dsbs) by gamma-rays.

MATERIALS AND METHODS

Spatial distributions of energy depositions by gamma-rays in liquid water were modelled with the track structure modules of the biophysical simulation code PARTRAC. Such distributions were superimposed on different DNA structure models (e.g. linear DNA, SV40 'minichromosomes' and compact chromatin), and direct energy depositions in the sugar-phosphate were considered as potential (direct) ssbs. The diffusion and interaction of the main chemical species produced in liquid water radiolysis were explicitly simulated, and reactions of *OH with the sugar-phosphate were considered as potential (indirect) ssbs. Two ssb on opposite DNA strands within 10 base pairs were considered as one dsb. Yields of ssb and dsb Gy(-1) Dalton(-1) in different DNA target structures were calculated as a function of the *OH mean lifetime, whose inverse value was taken as representative of the scavenging capacity of the DNA environment.

RESULTS AND CONCLUSIONS

A further validation of the models implemented in the PARTRAC code has been provided, thus allowing a better understanding of the mechanisms underlying DNA damage. More specifically, the protection due to *OH scavengers was separately quantified with respect to that due to histones and chromatin folding, which could be 'switched off' in the simulations. As expected, for a given value of the environment scavenging capacity, linear DNA was more susceptible to strand breakage than SV40 minichromosomes, which in turn showed higher damage yields with respect to cellular DNA due to the larger accessibility offered to *OH. Furthermore, by increasing the scavenging capacity, the break yields decreased in all structures and tended to coincide with direct damage yields. Very good agreement was found with available experimental data. Comparisons with data on 'nucleoid' DNA (i.e. unfolded and histone-depleted DNA) also suggested that the experimental procedures used to obtain such structures might lower the environment scavenging capacity owing to the loss of cellular scavengers.

Simulation of 125I induced DNA strand breaks in a CAP-DNA complex

The E. coli catabolite gene activator protein (CAP)-DNA complex with 125I located at the position of the H5 atom of the cytosine near the centre was incorporated into the PARTRAC track structure code. DNA strand breaks due to irradiation were calculated by track structure and radical attack simulations; strand breaks due to neutralisation of the highly charged 125Te ion were derived from a semi-empirical distribution. According to the calculations, the neutralisation effect dominates the strand breakage frequency at 2 bases away from the 125I decay site on both strands. The first breakage distribution counted from a 32P labelled end on the strand with 125I agreed well with experimental data, but on the opposite strand, the calculated distribution is more concentrated around the decay site and its yield is about 20% larger than the measured data.

Simulation of DNA damage after proton and low LET irradiation

A module for proton track structure simulation in liquid water was implemented in the biophysical model PARTRAC. Simulated tracks of energy deposition events from the radiation under investigation were superimposed on a higher-order DNA target model describing the whole genome inside a human cell. The parameters controlling DNA damage from direct and indirect effects were adapted to agree with yields and pathway contributions derived from gamma ray irradiation experiments. Single and double strand break (DSB) induction was simulated for irradiations by protons, photons and electrons over a wide range of initial energies. The relative biological effectiveness for DSB induction after proton irradiation was found to rise from 1.2 at 5 keV.micron-1 to about 2.5 at 70 keV.micron-1. About half of this growth resulted from an increased production of DSB clusters associated with small (< 10 kbp) fragments.

Monte Carlo simulation of DNA damage by low LET radiation using inhomogeneous higher order DNA targets

To test possible effects of the heterogeneous nature of the cell nucleus on simulation results of radiation-induced DNA damage, inhomogeneous targets have been implemented in the biophysical code PARTRAC. The geometry of the DNA and the histones was defined by spheres around the constituent atoms. Electron cross sections in liquid water were scaled according to the mass density of the different materials, whereas photon cross sections were derived from the sum of the cross sections for the constituent atoms. In the case of higher energy electrons the simulations show an increase of energy deposition in the DNA proportional to its high mass density. For photons with energies in the range of the carbon and the oxygen K-shell (0.28-0.53 keV), cross sections of DNA are larger than those of water, leading to an increased yield of strand breaks per average absorbed dose in the cell nucleus.

The two-stage clonal expansion model as an example of a biologically based model of radiation-induced cancer

A model with two stages and clonal expansion (TSCE) is reviewed as a prototype for biologically based models of cancer development. Applications of the TSCE model to data sets for animals and humans for particle radiation (alpha particles) are presented. The results suggest that the radiation not only influences the initiating mutation, but may also act as a promoter. A possible mechanism for the promoting action is described. The consequences of these results for the shapes of the radiation dose-response curves at low doses and dose rates are discussed.

Track structures and dose distributions from decays of (131)I and (125)I in and around water spheres simulating micrometastases of differentiated thyroid cancer

The disintegration of the radionuclides (131)I and (125)I and the subsequent charged-particle tracks left behind in water (as a model substance for a biological cell) are simulated by the Monte Carlo track structure simulation code PARTRAC, using new inelastic electron scattering cross sections for condensed water. Every photon and electron emitted was followed in detail, event by event, down to 10 eV. From the spatial information on the track structures, absorbed dose distributions per (131)I and (125)I decay were calculated in and around water spheres simulating micrometastases as well as in the tissue surrounding such metastases. These radionuclides were assumed to be distributed uniformly inside spheres of different diameters (0.01, 0.03, 0.1, 0.3, 1.0 and 3.0 mm). The respective electron degradation spectra, the nearest-neighbor distance distributions between inelastic events, and the distance distributions for all activations for both iodine radionuclides were calculated. The absorbed fractions of the initial electron energies, absorbed doses and energy depositions, and single-event distributions, F(1)(epsilon), inside the six water spheres described above and in the surrounding tissue were also calculated. The absorbed doses per decay inside the six water spheres, i.e., the calculated S values (listed from 0.01 to 3.0 mm), were 6.8 x 10(-4), 7.2 x 10(-5), 5.5 x 10(-6), 4.9 x 10(-7), 3.1 x 10(-8) and 1.8 x 10(-9) Gy Bq(-1) s(-1) for (131)I, and 3.4 x 10(-3), 1.7 x 10(-4), 5.1 x 10(-6), 2.0 x 10(-7), 5.6 x 10(-9) and 2.2 x 10(-10) Gy Bq(-1) s(-1) for (125)I. It is concluded that, in the treatment of thyroid cancer, the geometrical track structure properties of (125)I might be superior to those of (131)I in micrometastases with diameters less than 0.1 mm; however, in this medical context, many other factors also have to be considered.

Radiation-induced cell inactivation can increase the cancer risk

Radiation can inactivate cells that are replaced by dividing neighboring cells. If cells on the way to malignancy can fill the deficit faster than healthy cells, their number increases. A major part of the radon-induced lung cancers in the Colorado miners can be explained by a moderate increase in the replacement probability.

Simulation of exon deletion mutations induced by low-LET radiation at the HPRT locus

Friedland, W., Li, W. B., Jacob, P. and Paretzke, H. G. Simulation of Exon Deletion Mutations Induced by Low-LET Radiation at the HPRT Locus. Radiat. Res. 155, 703-715 (2001). The induction of HPRT mutants with exon deletions after irradiation with photons was simulated using the biophysical radiation track structure model PARTRAC. The exon-intron structure of the human HPRT gene was incorporated into the chromatin fiber model in PARTRAC. After gamma and X irradiation, simulated double-stranded DNA fragments that overlapped with exons were assumed to result in exon deletion mutations with a probability that depended on the genomic or the geometric distance between the breakpoints. The consequences of different assumptions about this probability of deletion formation were evaluated on the basis of the resulting fractions of total, terminal and intragenic deletions. Agreement with corresponding measurements was obtained assuming a constant probability of deletion formation for fragments smaller than about 0.1 Mbp, and a probability of deletion formation decreasing with increasing geometric or genomic distance between the end points for larger fragments. For these two assumptions, yields of mutants with exon deletions, size distributions of deletions, patterns of deleted exons, and patterns of deleted STS marker sites surrounding the gene were calculated and compared with experimental data. The yields, size distributions and exon deletion patterns were grossly consistent, whereas larger deviations were found for the STS marker deletion patterns in this comparison.

Age- and sex-specific relative thyroid radiation exposure to 131I in Ukraine after the Chernobyl accident

The age- and sex-dependence of the 131I induced count rates is determined for the measurements performed in Ukraine after the Chernobyl accident on the thyroids of over 60,000 persons. For this, the individual measurements are scaled in such a way that the mean values over age and sex on one side and the mean values over measurement series on the other side are normalized to one. The resulting distribution of all scaled measurements is roughly log-normal. Half of them lie within a factor 1.6 of the median. 131I induced count rates have a minimum at birth year 1986, about half the value of adults. The maximum count rates with about 30% above adults are reached for males around age 16 y. The count rates are up to about 40% (at age 14-17 y) higher for males than for females. The results are within statistical uncertainties independent of the geographical area and the urban or rural nature of the settlements. Starting from the relative count rates, the age- and sex-dependence is calculated for the thyroid activities 1 mo after the accident for the integrated activities and for the doses. The dose of young children is a factor of about 6.5 higher than that of adults. Uncertainties are estimated throughout.

Analysis of lung tumour risk in radon-exposed rats: an intercomparison of multi-step modelling

Three carcinogenesis modelling groups have both jointly and separately applied a multi-step carcinogenesis model with clonal expansion to one data set of lung tumours in rats exposed to radon (CEA, France). This study was designed to investigate the differences in modelling approach and fitting procedures used by the three groups in detail, and to explore possible discrepancies in the results. Using the same model assumptions and a (linear) radiation dependence on the first model step only, the three groups arrived at identical best fits, proving that the mathematical formalisms and fitting procedures do not lead to different results. However, when each group was allowed to find its own preferred fit for this data set, all three found a significantly better, but different fit to the data. All solutions indicated radiation to be an initiating agent and found additional radiation action necessary. The character of this additional radiation dependence, however, could not be unambiguously pinpointed. Tumour incidence data were described equally well when radiation dependence was taken into account in clonal expansion ("promotion") or in the second mutational step ("transformation"); extension to three model stages also resulted in an adequate description. The study showed that, although the three groups used one carcinogenesis model in principle, different model assumptions and/or different methods of finding the "best fit" could result in different descriptions of experimental data. This implies that on statistical grounds, different interpretations can be given for the action that radiation had in this data set. Different data, i.e. other data sets with age-dependent tumour data and/or information from cellular radiobiology experiments, are needed to specifically pin down the radiation dependence in the multi-step carcinogenesis process.

Stochastic aspects and uncertainties in the prechemical and chemical stages of electron tracks in liquid water: a quantitative analysis based on Monte Carlo simulations

A new physical module for the biophysical simulation code PARTRAC has recently been developed, based on newly derived electron inelastic-scattering cross-sections in liquid water. In the present work, two modules of PARTRAC describing the production, diffusion and interaction of chemical species were developed with the specific purpose of quantifying the role of the uncertainties in the parameters controlling the early stages of liquid water radiolysis. A set of values for such parameters was identified, and time-dependent yields and frequency distributions of chemical species produced by electrons of different energies were calculated. The calculated yields were in good agreement with available data and simulations, thus confirming the reliability of the code. As the primary-electron energy decreases down to 1 keV, the *OH decay kinetics were found to get faster, reflecting variations in the spatial distribution of the initial energy depositions. In agreement with analogous works, an opposite trend was found for energies of a few hundred eV, due to the very small number of species involved. The spreading effects shown at long times by *OH frequency distributions following 1 keV irradiation were found to be essentially due to stochastic aspects of the chemical stage, whereas for 1 MeV tracks the physical and pre-chemical stages also were found to play a significant role. Relevant differences in the calculated e(aq) -yields were found by coupling the physics of PARTRAC with descriptions of the pre-chemical and chemical stages adopted in different models. This indicates a strict interrelation of the various stages, and thus a strong dependence of the parameter values on the assumptions made for the preceding and subsequent stages of the process. Although equally acceptable results can be obtained starting from different assumptions, it is necessary to keep control of such uncertainties, since they can significantly influence the modeling of radical attack on DNA and, more generally, radiobiological damage estimation. This study confirms the need for new, independently derived data on specific steps of water radiolysis, to be included in comprehensive biophysical simulation codes.

Biologically based analysis of the data for the Colorado uranium miners cohort: age, dose and dose-rate effects

This study is a comprehensive analysis of the latest follow-up of the Colorado uranium miners cohort using the two-stage clonal expansion model with particular emphasis on effects related to age and exposure. The model provides a framework in which the hazard function for lung cancer mortality incorporates detailed information on exposure to radon and radon progeny from hard rock and uranium mining together with information on cigarette smoking. Even though the effect of smoking on lung cancer risk is explicitly modeled, a significant birth cohort effect is found which shows a linear increase in the baseline lung cancer risk with birth year of the miners in the cohort. The analysis based on the two-stage clonal expansion model suggests that exposure to radon affects both the rate of initiation of intermediate cells in the pathway to cancer and the rate of proliferation of intermediate cells. However, in contrast to the promotional effect of radon, which is highly significant, the effect of radon on the rate of initiation is found to be not significant. The model is also used to study the inverse dose-rate effect. This effect is evident for radon exposures typical for mines but is predicted to be attenuated, and for longer exposures even reversed, for the more protracted and lower radon exposures in homes. The model also predicts the drop in risk with time after exposure ceases. For residential exposures, lung cancer risks are compared with the estimates from the BEIR VI report. While the risk estimates are in agreement with those derived from residential studies, they are about two- to fourfold lower than those reported in the BEIR VI report.

Childhood exposure due to the Chernobyl accident and thyroid cancer risk in contaminated areas of Belarus and Russia

The thyroid dose due to 131I releases during the Chernobyl accident was reconstructed for children and adolescents in two cities and 2122 settlements in Belarus, and in one city and 607 settlements in the Bryansk district of the Russian Federation. In this area, which covers the two high contamination spots in the two countries following the accident, data on thyroid cancer incidence during the period 1991-1995 were analysed in the light of possible increased thyroid surveillance. Two methods of risk analysis were applied: Poisson regression with results for the single settlements and Monte Carlo (MC) calculations for results in larger areas or sub-populations. Best estimates of both methods agreed well. Poisson regression estimates of 95% confidence intervals (CIs) were considerably smaller than the MC results, which allow for extra-Poisson uncertainties due to reconstructed doses and the background thyroid cancer incidence. The excess absolute risk per unit thyroid dose (EARPD) for the birth cohort 1971-1985 by the MC analysis was 2.1 (95% CI 1.0-4.5) cases per 10(4) person-year Gy. The point estimate is lower by a factor of two than that observed in a pooled study of thyroid cancer risk after external exposures. The excess relative risk per unit thyroid dose was 23 (95% CI 8.6-82) Gy(-1). No significant differences between countries or cities and rural areas were found. In the lowest dose group of the settlements with an average thyroid dose of 0.05 Gy the risk was statistically significantly elevated. Dependencies of risks on age-at-exposure and on gender are consistent with findings after external exposures.

Simulation of DNA fragment distributions after irradiation with photons

The Monte Carlo track structure code PARTRAC has been further improved by implementing electron scattering cross-sections for liquid water and by explicitly modelling the interaction of water radicals with DNA. The model of the genome inside a human cell nucleus in its interphase is based on the atomic coordinates of the DNA double helix with an additional volume for the water shell. The DNA helix is wound around histone complexes, and these nucleosomes are folded into chromatin fibres and further to fibre loops, which are interconnected to build chromosomes with a territorial organisation. Simulations have been performed for the irradiation of human fibroblast cells with carbon K and aluminium K ultrasoft x-rays, 220 kVp x-rays and 60Co gamma-rays. The ratio single-strand breaks to double-strand breaks (ssb/dsb) for both types of ultrasoft x-rays is lower than for gamma-rays by a factor of 2. The contributions of direct and indirect effects to strand break induction are almost independent of photon energy. Strand break patterns from indirect effects reflect differences in the susceptibility of the DNA helix to OH* attack inside the chromatin fibre. Distributions of small DNA fragments (<3 kbp) are determined by the chromatin fibre structure irrespective of whether direct or indirect effects are causing the breaks. In the calculated fragment size distributions for larger DNA fragments (>30 kbp), a substantial deviation from random breakage is found only for carbon K irradiation, and is attributed to its inhomogeneous dose distribution inside the cell nucleus. For the other radiation qualities, the results for larger fragments can be approximated by random breakage distributions calculated for a yield of dsb which is about 10% lower than the average for the whole genome. The excess of DNA fragments detected experimentally in the 8-300 kbp region after x-ray irradiation is not seen in our simulation results.

Time trends of thyroid cancer incidence in Belarus after the Chernobyl accident

The rates of childhood thyroid cancer incidence observed in Belarus during the period 1986 to 1995 are described as a function of time after exposure, age at exposure, and sex. Conclusions are drawn for the excess absolute risk function. After a minimum latent period of about 3 years after exposure, this risk function has a linear increase with time for at least 6 years. After correction for the dependence of average doses on age, the radiation-induced absolute thyroid risk in Gomel is about a factor of 3 higher for children up to age 10 at exposure compared to older ones; this may be due in part to different case-collection quality. In addition, in the group up to 10 years at exposure, the thyroid of girls is more sensitive to radiation by a factor of about 1.5 than the thyroid of boys on an absolute scale. Risk estimates from external exposure are consistent with risk estimates from Gomel assuming that the increase in excess cases reaches a plateau soon.

Two-step model for the risk of fatal and incidental lung tumors in rats exposed to radon

Data from 4276 rats with radon exposures up to 10,000 WLM at rates up to 1000 WL are analyzed with a two-step clonal expansion model. The age dependences of the hazard for the risks for fatal and for incidental tumors are very different. Therefore, two different parameterizations of the model are used in the two cases. In both cases radiation acts only on the initiating mutation and the clonal expansion, but not on the second mutation. Average exposure rates of 5 WL for fatal tumors and 0.5 WL for incidental tumors double the rate of spontaneous mutations. While the fatal tumors show a linear increase in the effective clonal expansion rate up to about 100 WL average exposure rate and a saturation at higher exposure rates, the incidental tumors follow a step-like behavior of this parameter. It is proposed that only the fatal lung tumors among the rats be used for generalizations to models for lung cancer in humans. The fitted model for fatal tumors shows an inverse dose-rate effect at average exposure rates above 20 WL. However, below 10 WL the lung cancer risk per unit exposure decreases with increasing duration of exposure. Between 10 and 20 WL, the difference in ERR/WLM between acute and protracted exposure is small.

Measurement of resuspended aerosol in the Chernobyl area. Part III. Size distribution and dry deposition velocity of radioactive particles during anthropogenic enhanced resuspension

During anthropogenic activities, such as agricultural soil management and traffic on unpaved roads, size distribution measurements were performed of atmospheric particulate radionuclides at a site in the Chernobyl 30-km exclusion zone. Analysis of cascade impactor measurements showed an increase of the total atmospheric radioactivity. In the cases of harrowing by a tractor and traffic on unpaved roads, a common shape of the size distribution was found with two maxima, the first in the 2-4 microm range, the second in the 12-20 microm range. The size distributions were compared to measurements during wind-driven resuspension. Particle number concentration measurements with an Aerodynamic Particle Sizer showed a dynamic dependence of the particle concentration in different size ranges on anthropogenic action. The increase of the mean concentration was for the large particles more than one order of magnitude higher than for fine particles during anthropogenic enhanced resuspension. From the measurement of the mass concentration, the radioactive loading could be estimated. An enrichment of radionuclides on resuspended particles (compared to soil particles) was found, with the highest enrichment for large particles. Micrometeorological considerations showed that large particles may frequently be subject to medium range transport. The dry deposition velocity was measured; the mean value of 0.026 m s(-1) +/- 0.016 m s(-1) is typical for 6-9 microm diameter particles.

Monte Carlo simulation of the production of short DNA fragments by low-linear energy transfer radiation using higher-order DNA models

A realistic DNA target model has been developed and implemented in the biophysical simulation code PARTRAC. It describes five levels of the B-DNA structure (nucleotides, DNA helices, nucleosomes, chromatin fiber structure and chromatin fiber loops) on an atomic level for the whole genome inside a mammalian cell nucleus. The model is capable of describing regular solenoidal, crossed-linker or zigzag structures as well as repeating stochastic arrangements of nucleosomes in the chromatin fiber. Electron tracks resulting from monoenergetic electrons with energies up to 100 keV and from 220 kVp X rays, starting at random positions in the cell, were superimposed on four DNA target models with different chromatin fiber structures. The yields of SSBs, DSBs and short single- and double-stranded DNA fragments were determined from spatial coincidences with strand atoms. Two parameters of the model-the energy necessary to create an SSB and the distance between two breaks that would be scored as a DSB-were adapted to equate simulated and measured strand break yields after X irradiation of human fibroblast cells. The integral fractions of short single- and double-stranded fragments were rather similar for all condensed chromatin fiber structures; they agreed with experimental data for DNA fragments below 2 kbp. The simulated fragment size distributions in the range from 0.1 to 1.5 kbp reflected the fiber structure irrespective of strandedness or electron energy. The distributions using a stochastic arrangement of nucleosomes in the chromatin fiber were found to be in better accordance with experimental data than those obtained with regular fiber structures.

Measurement of resuspended aerosol in the Chernobyl area. Part II. Size distribution of radioactive particles

Size distribution measurements of particulate radionuclides were performed at two sites in the Chernobyl 30-km exclusion zone using several cascade impactors. The results obtained in the period September 1986 till June 1993 were discussed with regard to the general assumption of a log-normal activity size distribution in inhalation dose assessment. At Zapolie (a site 14 km from the Chernobyl reactor) a bimodal distribution was observed in 91% of all measured distributions. In most cases the medians were about 4 microns and in the range 20-30 microns. According to soil granulometric data this finding was explained by superimposing two processes: local resuspension and advective transport of radioactive aerosol from highly contaminated territories. The mean air concentration showed an increasing proportion of inhalable particles over the years since the accident. In 1993 the inhalable fraction was about 48% of the total concentration. At Pripyat, a site situated within a highly contaminated area, unimodal types of size distributions were predominant with the median diameters in the range 5-10 microns for 137Cs. For the three nuclides 137Cs, 144Ce and 106Ru, very similar types of distribution were observed. Apparently, the radioactive aerosol was of fuel origin. During a forest fire at a distance of 17 km, the majority of the radioactivity was associated with submicrometer particles with median diameters in the range 0.28-0.50 micron.

PARATI--a dynamic model for radiological assessments in urban areas. Part III: Parameter uncertainty analysis

The uncertainties in the exposure predictions after contamination of an urban area due to the variabilities in environmental transfer parameters and in dose conversion factors have been estimated. This was done using the "Latin Hypercube" sampling scheme and the computer codes PRISM and PARATI. For the scenario 'urban contamination by 137Cs' and the population groups considered, the main sources contributing to the uncertainty in the resulting exposures are the limited knowledge of the initial deposition and retention, the weathering processes, the actual urban environments, and the characteristics and habits of the population. The effect of the parameter uncertainties on the variability of the dose is almost constant over time.

Health impacts of large releases of radionuclides. Late somatic health effects

This chapter reviews the risks of radiation-induced cancer for the dose range likely to occur after releases of radionuclides into the environment. Epidemiological evidence from exposed workers and the atomic bomb survivors of Hiroshima and Nagasaki is surveyed. Influences on such risk functions of individual related quantities (e.g. age, sex, nationality, time since exposure and organs exposed) and of radiation modality-related quantities (e.g. dose, dose rate and radiation quality) are also briefly discussed.

Measurement of resuspended aerosol in the Chernobyl area. I. Discussion of instrumentation and estimation of measurement uncertainty

Results of measurements of the resuspended radioactive aerosols in the Chernobyl area are presented which were obtained soon after the Chernobyl reactor accident and in a European project in 1992-1993. The measurements were carried out with the intention of obtaining a data base for dose assessment of resuspended radioactive particles. Potential significant dose contributions may result from inhalation and secondary contamination due to resuspended radionuclides. In this first article of a series of three papers, the instrumentation and the measurement uncertainties are discussed. An effort was made to sample quantitatively giant aerosol particles (particles larger than 10 microns aerodynamic diameter) as well. The comparison of the samplers shows, in general, an agreement of concentration measurements of 137Cs and 7Be within a factor of two. One sampler was identified with larger discrepancies and needs additional investigation of its sampling characteristics; for another device, the recalibration of the analysing system is recommended. Ordinary integrating samplers have a loss of about 30% in 137Cs activity compared to an isokinetic sampler collecting giant particles as well. The mean ratio of 137Cs activity concentration between an instrument sampling only particles larger than 10 microns and an ordinary integrating sampler is 0.39 +/- 0.15 during anthropogenic-enhanced resuspension. These findings demonstrate the significant contribution of giant particles to resuspended airborne radioactivity. The results of this study concerning integral measurements during wind-driven resuspension proved to be in good agreement with previously published data on resuspension.