Application of the Weibull distribution to describe the vertical distribution of cesium-137 on a slope under permanent pasture in Luxembourg

The application of solutions of the Advection-Dispersion-Equation (ADE) for soil profiles is sometimes questionable. An alternative method, based on the Weibull distribution, has been developed, which can approximate the vertical distribution of radiocesium (137Cs) and allows extrapolation to account for the complete inventory. The structure of the equation allows simple parameters describing the soil depth profile to be derived. Reliable estimates of the total 137Cs inventory can help to explain the lateral distribution. This equation was used to analyse the 137Cs-activities of Chernobyl fall-out measured on a slope under permanent pasture in Luxembourg.     

A simplified 137Cs transport model for estimating erosion rates in undisturbed soil

(137)Cs is an artificial radionuclide with a half-life of 30.12 years which released into the environment as a result of atmospheric testing of thermo-nuclear weapons primarily during the period of 1950s-1970s with the maximum rate of (137)Cs fallout from atmosphere in 1963. (137)Cs fallout is strongly and rapidly adsorbed by fine particles in the surface horizons of the soil, when it falls down on the ground mostly with precipitation. Its subsequent redistribution is associated with movements of the soil or sediment particles. The (137)Cs nuclide tracing technique has been used for assessment of soil losses for both undisturbed and cultivated soils. For undisturbed soils, a simple profile-shape model was developed in 1990 to describe the (137)Cs depth distribution in profile, where the maximum (137)Cs occurs in the surface horizon and it exponentially decreases with depth. The model implied that the total (137)Cs fallout amount deposited on the earth surface in 1963 and the (137)Cs profile shape has not changed with time. The model has been widely used for assessment of soil losses on undisturbed land. However, temporal variations of (137)Cs depth distribution in undisturbed soils after its deposition on the ground due to downward transport processes are not considered in the previous simple profile-shape model. Thus, the soil losses are overestimated by the model. On the base of the erosion assessment model developed by Walling, D.E., He, Q. [1999. Improved models for estimating soil erosion rates from cesium-137 measurements. Journal of Environmental Quality 28, 611-622], we discuss the (137)Cs transport process in the eroded soil profile and make some simplification to the model, develop a method to estimate the soil erosion rate more expediently. To compare the soil erosion rates calculated by the simple profile-shape model and the simple transport model, the soil losses related to different (137)Cs loss proportions of the reference inventory at the Kaixian site of the Three Gorge Region, China are estimated by the two models. The over-estimation of the soil loss by using the previous simple profile-shape model obviously increases with the time period from the sampling year to the year of 1963 and (137)Cs loss proportion of the reference inventory. As to 20-80% of (137)Cs loss proportions of the reference inventory at the Kaixian site in 2004, the annual soil loss depths estimated by the new simplified transport process model are only 57.90-56.24% of the values estimated by the previous model.     

A blind test of the MOIRA lake model for radiocesium for Lake Uruskul, Russia, contaminated by fallout from the Kyshtym accident in 1957

This paper presents results of a model-test carried out within the framework of the COMETES project (EU). The tested model is a new lake model for radiocesium to be used within the MOIRA decision support system (DSS; MOIRA and COMETES are acronyms for EU-projects). This model has previously been validated against independent data from many lakes covering a wide domain of lake characteristics and been demonstrated to yield excellent predictive power (see H?kanson, Modelling Radiocesium in Lakes and Coastal Areas. Kluwer, Dordrecht, 2000, 215 pp). However, the model has not been tested before for cases other than those related to the Chernobyl fallout in 1986, nor for lakes from this part of the world (Southern Urals) and nor for situations with such heavy fallout as this. The aims of this work were: (1) to carry out a blind test of the model for the case of continental Lake Uruskul, heavily contaminated with 90Sr and 137Cs as a result of the Kyshtym radiation accident (29 September 1957) in the Southern Urals, Russia, and (2) if these tests gave satisfactory results to reconstruct the radiocesium dynamics for fish, water and sediments in the lake. Can the model provide meaningful predictions in a situation such as this? The answer is yes, although there are reservations due to the scarcity of reliable empirical data. From the modelling calculations, it may be noted that the maximum levels of 137Cs in fish (here 400 g ww goldfish), water and sediments were about 100,000 Bq/kg ww, 600 Bq/l and 30,000 Bq/kg dw, respectively. The values in fish are comparable to or higher than the levels in fish in the cooling pond of the Chernobyl NPP. The model also predicts an interesting seasonal pattern in 137Cs levels in sediments. There is also a characteristic "three phase" development for the 137Cs levels in fish: first an initial stage when the 137Cs concentrations in fish approach a maximum value, then a phase with relatively short ecological half-lives followed by a final phase with long ecological half-lives more or less corresponding to the physical decay of radiocesium.     

Uncertainty analysis of the external gamma-dose rate due to the variability of the vertical distribution of 137Cs in the soil

The external gamma-dose rate at 1 m height above a flat area due to the presence of fallout radiocesium in the soil is frequently calculated from the observed depth profile of the 137Cs activity as well as the soil mass per unit area. At a given site, these depth profiles may, however, vary considerably, thus introducing an uncertainty to the external gamma-dose calculated in this way. To assess this source of uncertainty for a typical grassland site, the activity of Chernobyl-derived 137Cs and the wet bulk density in the three upper soil layers at 100 plots in a 100 m x 100 m pasture were determined. Analysis of these data shows that the frequency distribution of the dose rates calculated from the corresponding depth profiles of all plots is similar to a log-normal distribution (mean 25 nGy h-1, median 22 nGy h-1, standard deviation 11 nGy h-1; range 1.6-56 nGy h-1). The various sources which contribute to the uncertainty of the dose rate are quantified. The semi-variogram indicates that any spatial dependence of the dose rates occurs on this pasture only over distances that are smaller than the shortest sampling interval (here about 10 m). It is estimated which errors have to be expected for the median dose rate when the depth profiles of 137Cs and of the wet bulk density are determined only for a small number of plots. It is preferable to calculate the mean dose rate as a mean from the n individual dose rates rather than from an averaged 137Cs depth profile of the n plots.     

Vertical distribution and inventories of 137Cs in the Syrian soils of the eastern Mediterranean region

Vertical distribution and inventories of (137)Cs have been determined using radiocesium distributions in presumably undistributed soil profiles, collected from 36 sites distributed all over Syria (eastern Mediterranean region). Vertical distributions of (137)Cs in the collected profiles were found to be strongly correlated with soil type and five groups were identified. Based on these profiles, total (137)Cs inventory (bomb test and Chernobyl) varied between 320 Bq m(-2) and 9,647 Bq m(-2). Geographical mapping of (137)Cs inventories showed that the highest values were found in the coastal, middle and north-east regions of Syria indicating that Chernobyl atmospheric contribution to the total (137)Cs deposition in the region is predominant. In contrast, the lowest values were found in the south-east region (Syrian Badia), where a relatively uniform distribution was observed, which may only be attributed to the past global nuclear bomb test. The measured inventories were also compared with a mathematical model for estimating bomb derived (137)Cs reference inventories.     

Inventory and vertical migration of 90Sr fallout and 137Cs/90Sr ratio in Spanish mainland soils

In this paper the inventory of ⁹⁰Sr in 34 points distributed along the Spanish peninsular territory is presented. Obtained values range between 173 Bq/m² and 2047 Bq/m². From these data set and those ¹³⁷Cs data obtained in a previous work the ¹³⁷Cs/⁹⁰Sr activity ratio has been established, laying this value between 0.9 and 3.6. Also the migration depth of both radionuclides has been analysed obtaining for ¹³⁷Cs an average value 57% lower than that obtained for ⁹⁰Sr.Additionally, this paper presents the results obtained in 11 sampling points in which the activity vertical profile has been measured. These profiles have been analysed to state the behaviour of strontium in soils and after, by using a convective-diffusive model, the parameters of the model which governs the vertical migration of ⁹⁰Sr in the soil, v (apparent convection velocity) and D (apparent diffusion coefficient) have been evaluated. Mean values obtained are 0.20 cm/year and 3.67 cm²/year, respectively.     

Determining floodplain sedimentation rates using Cs in a low fallout environment dominated by channel- and cultivation-derived sediment inputs, central Queensland, Australia

Fallout ¹³⁷Cs has been widely used to determine floodplain sedimentation rates in temperate environments, particularly in the northern hemisphere. Its application in low fallout, tropical environments in the southern hemisphere has been limited. In this study we assess the utility of ¹³⁷Cs for determining rates of floodplain sedimentation in a dry-tropical catchment in central Queensland, Australia. Floodplain and reference site cores were analysed in two centimetre increments, depth profiles were produced and total ¹³⁷Cs inventories calculated from the detailed profile data. Information on the rates of ¹³⁷Cs migration through local soils was obtained from the reference site soil cores. This data was used in an advection–diffusion model to account of ¹³⁷Cs mobility in floodplain sediment cores. This allowed sedimentation rates to be determined without the first year of detection for ¹³⁷Cs being known and without having to assume that ¹³⁷Cs remains immobile following deposition. Caesium-137 depth profiles in this environment are demonstrated to be an effective way of determining floodplain sedimentation rates. The total ¹³⁷Cs inventory approach was found to be less successful, with only one of the three sites analysed being in unequivocal agreement with the depth profile results. The input of sediment from catchment sources that have little, or no, ¹³⁷Cs attached results in true depositional sites having total inventories that are not significantly different from those of undisturbed reference sites.     

Radiocesium accumulation in mosses from highlands of Serbia and Montenegro: chemical and physiological aspects

The aim of this work was (i) to determine the activity levels of 137Cs in mosses from highland ecosystems of Serbia and Montenegro, (ii) to find out if radiocesium is associated with essential biomacromolecules, and (iii) to investigate 137Cs distribution among intracellular compartments. It was found that biomolecules of mosses do not bind significant amounts of radiocesium (2.3-3.3% of the absorbed 137Cs), a behavior that was independent of the moss species. Cellular fractionation of mosses showed that membranes are the primary 137Cs-binding sites at the cellular level. They contained 26.1-43.1% of the initial radiocesium activity. It seems that 137Cs-binding molecules in different mosses are of similar chemical nature, and their distribution between various cellular compartments is not species specific.     

Modelling 137Cs uptake in plants from undisturbed soil monoliths

A model predicting 137Cs uptake in plants was applied on data from artificially contaminated lysimeters. The lysimeter data involve three different crops (beans, ryegrass and lettuce) grown on five different soils between 3 and 5 years after contamination and where soil solution composition was monitored. The mechanistic model predicts plant uptake of 137Cs from soil solution composition. Predicted K concentrations in the rhizosphere were up to 50-fold below that in the bulk soil solution whereas corresponding 137Cs concentration gradients were always less pronounced. Predictions of crop 137Cs content based on rhizosphere soil solution compositions were generally closer to observations than those based on bulk soil solution composition. The model explained 17% (beans) to 91% (lettuce) of the variation in 137Cs activity concentrations in the plants. The model failed to predict the 137Cs activity concentration in ryegrass where uptake of the 5-year-old 137Cs from 3 soils was about 40-fold larger than predicted. The model generally underpredicted crop 137Cs concentrations at soil solution K concentration below about 1.0 mM. It is concluded that 137Cs uptake can be predicted from the soil solution composition at adequate K nutrition but that significant uncertainties remain when soil solution K is below 1 mM.     

Inventory and vertical migration of ¹³⁷Cs in Spanish mainland soils

In this study the total activity of ¹³⁷Cs deposited per unit area over the Spanish peninsular territory was analysed using a 150 × 150 km² mesh grid, with samples taken from 29 points. The deposited activities ranged between 251 and 6074 Bq/m². A linear relationship was obtained between these values and the mean annual rainfall at each sampling point which allowed a map to be drawn, using GIS software, which shows the distribution of total deposited ¹³⁷Cs activity across the Spanish mainland. At twelve of these sampling points the vertical migration profile of ¹³⁷Cs was obtained. These profiles are separated into two groups with different behaviour, one of which includes clay and loam soils and the other containing sandy soils. For both groups of profiles the parameters of the convective-diffusive model, which describes the vertical migration of ¹³⁷Cs in the soil, v (apparent convection velocity) and D (apparent diffusion coefficient) were calculated.     

A new generic sub-model for radionuclide fixation in large catchments from continuous and single-pulse fallouts, as used in a river model

This paper presents a new general sub-model for fixation in catchment areas to be used within the framework of a river model for substances such as radionuclides and metals from continuous and single-pulse fallouts. The model has been critically tested using data from 27 European river sites covering a very wide geographical area and contaminated by radiocesium and radiostrontium from the Chernobyl accident and from the nuclear weapons tests (NWT fallout). This modelling approach gives radionuclide concentrations in water (total, dissolved and particulate phases) at defined sites on a monthly basis. The overall river model is based on processes in the upstream river stretch and in the catchment area. The catchment area is differentiated into inflow (approximately dry land) areas and outflow (approximately wetland) areas. The model has a general structure, which can be used for all radionuclides or substances. It is simple to apply in practice since all driving variables may be readily accessed from maps and standard monitoring programs. The driving variables are: latitude, altitude, catchment area, mean annual precipitation and fallout. Note that for large catchments, this model does not require data on the characteristic soil type or the percentage of outflow areas (wet lands) in the catchment, as in most previous models, since in practice it is very difficult to obtain reliable data on characteristic soil type or percentage of outflow areas, especially in large and topographically complex catchments. Modelled values have been compared to empirical data from rivers sites covering a wide domain (catchment areas from 3000 to 3,000,000 km2, precipitation from 400 to 1700 mm/year; fallouts from 1600 to 280,000 Bq/m2; altitudes from 0 to 1000 m.a.s.l. and latitudes from 41 degrees to 72 degrees N). The river model with its sub-model for fixation predicts close to the uncertainty factors given by the empirical data, which have been shown to be about a factor of 1.6 for 137Cs and a factor of 2.2 for 90Sr in river water. The obtained characteristic uncertainty factors for 137Cs from the Chernobyl fallout is 2.4, for 137Cs from the NWT fallout it is 1.3 and for the 90Sr results from the NWT fallout it is 3 using the new model.     

Radiocesium storage in soil microbial biomass of undisturbed alpine meadow soils and its relation to 137Cs soil-plant transfer

This study focuses on radiocesium storage in soil microbial biomass of undisturbed alpine meadow sites and its relation to the soil-to-plant transfer. Soil and plant samples were taken in August 1999 from an altitude transect (800-1600m.a.s.l.) at Gastein valley, Austria. Soil samples were subdivided into 3-cm layers for analyses of total, K(2)SO(4)-extractable and microbially stored (137)Cs. Microbial biomass was measured by the fumigation extraction method, and fungal biomass was quantified using ergosterol as biomarker molecule. In general, the quantity of (137)Cs stored in the living soil microbial biomass was relatively small. At the high-altitude meadows, showing high amounts of fungal biomass, microbially stored (137)Cs amounted to 0.64+/-0.14kBqm(-2) which corresponds to about 1.2-2.7% of the total (137)Cs soil inventory. At lower altitudes, microbial (137)Cs content was distinctly smaller and in most cases not measurable at all using the fumigation extraction method. However, a positive correlation between the observed soil-to-plant aggregated transfer factor, microbially stored (137)Cs and fungal biomass was found, which indicates a possible role of fungal biomass in the storage and turnover of (137)Cs in soils and in the (137)Cs uptake by plants.     

Analysis of Factors Determining Accumulation of 137Cs by Woody Plants

The accumulation of ¹³⁷Cs by coniferous (Pinus sylvestris) and deciduous (Betula pendula) trees was studied in the area contaminated with radioactive products after the Chernobyl accident. The availability of ¹³⁷Cs to woody plants is determined by its vertical distribution in the soil profile, the proportion of its exchangeable form, and the distribution of tree roots in the soil. The index of biological availability of ¹³⁷C in the soil is proposed, and its linear dependence on the coefficients of ¹³⁷Cs transfer to wood is demonstrated.     

Radiocesium contamination behavior and its effect on potassium absorption in tropical or subtropical plants

The accumulation and long-term decline of radiocesium contamination in tropical plant species was studied through measurements of gamma-ray spectra from pomegranate (Punica granatum) and chili pepper (Capsicum fructescens) trees. The plants were originally grown at a (137)Cs contaminated site (where a radiological accident occurred in the city of Goiania, Brazil, in 1987), and transplanted to uncontaminated soil, so that the main source of contamination of the new leaves and fruits would be the fraction of the available radiocesium in the body of the plants. Measurements of (137)Cs and (40)K concentrations along the roots, main trunk, twigs, leaves and fruits before and after the transplant process of both plant species indicated a direct competition between Cs and K ions, suggesting that these elements could have a common accumulation mechanism. Cesium transfer factors from soil to pomegranate, green and red chili pepper fruits were evaluated as 0.4 +/- 0.1, 0.06 +/- 0.01 and 0.05 +/- 0.01, respectively. Biological half-life values due to (137)Cs translocation from the tree reservoir (BHL(T)) were calculated as 0.30 years for pomegranate, 0.12 years and 0.07 years for red and green peppers, respectively.     

Soil organic horizons as a major source for radiocesium biorecycling in forest ecosystems

Here we review some of the main processes and key parameters affecting the mobility of radiocesium in soils of semi-natural areas. We further illustrate them in a collection of soil surface horizons which largely differ in their organic matter contents. In soils, specific retention of radiocesium occurs in a very small number of sorbing sites, which are the frayed edge sites (FES) born out of weathered micaceous minerals. The FES abundance directly governs the mobility of trace Cs in the rhizosphere and thus its transfer from soil to plant. Here, we show that the accumulation of organic matter in topsoils can exert a dilution of FES-bearing minerals in the thick humus of some forest soils. Consequently, such accumulation significantly contributes to increasing 137Cs soil-to-plant transfer. Potassium depletion and extensive exploration of the organic horizons by plant roots can further enhance the contamination hazard. As humus thickness depends on both ecological conditions and forest management. our observations support the following ideas: (1) forest ecosystems can be classified according to their sensitivity to radiocesium bio-recycling, (2) specific forest management could be searched to decrease such bio-recycling.     

Seasonality of Cs in roe deer from Austria and Germany

Empirical data on the ¹³⁷Cs activity concentration in meat of roe deer (Capreolus capreolus) roaming in 3 spruce forest areas and one peat bog area are presented and compared. They cover time series of nearly 20 years after a spike contamination in 1986 originating from Chernobyl. A model is presented which considers three soil compartments to describe the change of the availability of ¹³⁷Cs with time.     

Cs and K in fruiting bodies of different fungal species collected in a single forest in southern Poland

Fruiting bodies of fungi belonging to more than 70 species were collected within a few thousand square meter area of one forest during 2006 and 2007. The soil profile was collected to check the cumulative deposition of ¹³⁷Cs, which was relatively high, equal to 64 ± 2 kBq/m² (calculated for October 2006). The majority of this activity was in the first 6 cm. Fruitbodies were analyzed for radiocesium and ⁴⁰K by means of gamma-spectrometry. The highest ¹³⁷Cs activity was 54.1 ± 0.7 kBq/kg (dry weight) for a sample of Lactarius helvus collected in 2006. The results for 2006 were higher than those for 2007. In a few cases the traces of short-lived (T_1/2 = 2.06 a) ¹³⁴Cs were still found in samples. The importance of mycorrhizal fungi for radiocesium accumulation is confirmed. The differences in activity among the species are discussed in relation to observations and predictions from previous studies, where the change in relative accumulation between fruiting bodies of different species was at least partially explained by the differences in the depth of the mycelium localization in a litter/soil system. It is concluded that in some cases, such as Boletus edulis and Xerocomus badius, this prediction is fulfilled and therefore this explanation confirmed.     

Numerical simulation of 137Cs and (239,240)Pu concentrations by an ocean general circulation model

We simulated the spatial distributions and the temporal variations of 137Cs and (239,240)Pu concentrations in the ocean by using the ocean general circulation model which was developed by National Center of Atmospheric Research. These nuclides are introduced into seawaters from global fallout due to atmospheric nuclear weapons tests. The distribution of radioactive deposition on the world ocean is estimated from global precipitation data and observed values of annual deposition of radionuclides at the Meteorological Research Institute in Japan and several observed points in New Zealand. Radionuclides from global fallout have been transported by advection, diffusion and scavenging, and this concentration reduces by radioactive decay in the ocean. We verified the results of the model calculations by comparing simulated values of 137Cs and (239,240)Pu in seawater with the observed values included in the Historical Artificial Radionuclides in the HAM database, which has been constructed by the Meteorological Research Institute. The vertical distributions of the calculated 137Cs concentrations were in good agreement and are in good agreement with the observed profiles in the 1960s up to 250 m, in the 1970s up to 500 m, in the 1980s up to 750 m and in the 1990s up to 750 m. However, the calculated 137Cs concentrations were underestimated compared with the observed 137Cs at the deeper layer. This may suggest other transport processes of 137Cs to deep waters. The horizontal distributions of 137Cs concentrations in surface water could be simulated. A numerical tracer release experiment was performed to explain the horizontal distribution pattern. A maximum (239,240)Pu concentration layer occurs at an intermediate depth for both observed and calculated values, which is formed by particle scavenging. The horizontal distributions of the calculated (239,240)Pu concentrations in surface water could be simulated by considering the scavenging effect.     

Global fallout (137)Cs accumulation and vertical migration in selected soils from South Patagonia

The spatial distribution and vertical migration of global fallout (137)Cs were studied in soils from South Patagonia at the austral region of South America in semi-natural and natural environments located between 50-54 degrees S and 68-74 degrees W. The (137)Cs areal activity density varied from 222 to 858 Bq m(-2), and was found to be significantly positively correlated (p<0.001) with the mean annual precipitation rate. The fraction of the total activity density observed in steppe grass varied from <0.03% to 0.12% (median <0.07%) and is considerably lower than the results obtained at the South Shetland Islands (median 8%) and in other temperate environments in south-central Chile (median 0.2%). The median of the convection velocity v(s) of (137)Cs in the soil in such polar isotundra climate has been determined to be 0.056 cm y(-1). This value is higher than v(s) determined under polar climate (-0.012 cm y(-1)) and is near to the upper limit of v(s)-values determined in temperate environments from Chile (0.019 cm y(-1)). The median value of the diffusion coefficient D(s) (0.048 cm(2) y(-1)) is similar to D(s) observed in an Antarctic region (0.043 cm(2) y(-1)) and lower than D(s) in temperate regions of Chile (1.24 cm(2) y(-1)). About 35 years after the highest depositions, (137)Cs had penetrated to a depth of 6-14 cm in the Patagonian soils and can be expected to remain in the rooting zone of grass for many decades. Nevertheless, because of its low transfer to steppe grass observed at this region, the radioecological sensitivity of this ecosystem with respect to fallout radiocesium seems to be lower than in other polar regions.     

Validation of I ecological transfer models and thyroid dose assessments using Chernobyl fallout data from the Plavsk district, Russia

Within the project “Environmental Modelling for Radiation Safety” (EMRAS) organized by the IAEA in 2003 experimental data of ¹³¹I measurements following the Chernobyl accident in the Plavsk district of Tula region, Russia were used to validate the calculations of some radioecological transfer models. Nine models participated in the inter-comparison. Levels of ¹³⁷Cs soil contamination in all the settlements and ¹³¹I/¹³⁷Cs isotopic ratios in the depositions in some locations were used as the main input information. 370 measurements of ¹³¹I content in thyroid of townspeople and villagers, and 90 measurements of ¹³¹I concentration in milk were used for validation of the model predictions.A remarkable improvement in models performance comparing with previous inter-comparison exercise was demonstrated. Predictions of the various models were within a factor of three relative to the observations, discrepancies between the estimates of average doses to thyroid produced by most participant not exceeded a factor of ten.