Radioactivity of soil in Croatia II: 137Cs, 40K, and absorbed dose rate.

We took samples of uncultivated soil from the surface layer (0-10 cm) at 138 sites from all over Croatia and measured their radionuclide activity concentrations with high-resolution gamma-ray spectrometry. This second part of our report brings the results on 40K and 137Cs to complement those on the 232Th and 238U decay chains addressed in the first part. Together they give the most complete picture of radioactivity of Croatian soil so far. Activity concentrations of 40K were the highest in the Pannonian region, and there was an opposite trend for 137Cs. We found that the concentrations of 137Cs tended to increase with altitude, annual precipitation, and vegetation density. The concentration ratio of 137Cs and K in soil, which indicates the potential for 137Cs entering food chains via uptake by plants, was the lowest in agriculturally important areas in the east of the Pannonian region. In addition, we used the obtained results on activity concentrations to calculate the related absorbed dose rate as a measure of external exposure to ionising radiation from soil. The sum of the absorbed dose rates for naturally occurring radionuclides and 137Cs showed that external exposure was generally the highest in the Dinaric region and Istrian Peninsula.

Exposure to ionising radiation in the environment stems from both naturally occurring and anthropogenic radionuclides. Since life on Earth has evolved in the absence of anthropogenic radionuclides, their presence in any quantity is threatening. Nuclear fission is the main source of anthropogenic radionuclides, of which 137Cs is considered particularly harmful to living organisms due to its long half-life T1/2 of 30.1 years, large production rate, and chemical properties similar to those of other alkali metals. In fact, activity concentrations of 137Cs in the environment are – together with those of 222Rn – in the focus of radioecology (1). The ecological impact of 137Cs is determined not only by its radioactivity (it emits both beta and gamma radiation) but also by the physical and chemical properties that govern its propagation through the environment and biota. The Cs+ ionic state facilitates potassium and sodium (biogenic elements) substitution with 137Cs in living organisms, which implies that 137Cs concentrations in a given medium should be compared with the concentrations of K at least, which is possible by measuring radioactive 40K.

Under non-accidental conditions, most of external exposure to ionising radiation originates from soil. In spite of a vast amount of related studies, radioactive properties of soil continue to attract attention and are investigated with regard to agronomy (2, 3), human activities (4, 5, 6, 7, 8, 9), and the ecology of intact nature (10, 11). In the first part of the report (12) on the radioactivity of uncultivated soil in Croatia, we presented the activity concentrations (A) of the most significant members of the 238U and 232Th naturally occurring radioactive chains. Here we complement these results with those for 137Cs and 40K in the same samples and also present calculated absorbed dose rate for external exposure to ionising radiation from soil. Overall, our results for naturally occurring radionuclides and 137Cs give the most complete picture of soil radioactivity in Croatia so far. As in (12), we use maps to present our results. Detailed numerical data can be found elsewhere (13) and are also available from the authors upon request.

Geomorphological, biogeographical, and climatological characteristics of Croatia

Figure 1 shows a widely used division of Croatia with respect to geomorphological, biogeographical, and climatological characteristics. A more detailed description of the main regions and subregions has been given in the first part of this report (12), and here we repeat the main points only.

Figure 1

Main geomorphological, biogeographical, and climatological regions of Croatia. Region I belongs to the Pannonian Plain, with subregion Ia comprising hilly areas and subregion Ib being a flat lowland. Region II belongs to the Dinaric Alps. In subregion IIa, the climate is cold continental and vegetation is subalpine, whereas in subregion IIb, karst prevails, and the climate is Mediterranean. Region III is Mediterranean in both climate and vegetation. In subregion IIIa the influence of regions I and II is stronger than in subregion IIIb. Sizeable areas of dense forests are indicated by letter F

Region I is part of the Pannonian Plain. Its western subregion (Ia) comprises hilly and flat terrains, whereas subregion Ib in the east is a flat lowland. The vegetation, soil types, and climate are typical of Central Europe, and the annual precipitation decreases from the west (~1200 mm) to the east (~500 mm) of the region. Region II is part of the Dinaric Alps, underlain by limestone and dolomite bedrock. In subregion IIa the vegetation is subalpine and the climate is cold continental, with annual precipitation mostly exceeding 1500 mm and being as high as 3500 mm in densely forested mountains. Subregion IIb is transitional and influenced by the vicinity of the Mediterranean. Annual precipitation is 1000–1500 mm. Region III comprises the eastern Adriatic coast, accompanying islands, and parts of the hinterland. The bedrock is again limestone and dolomite, whereas the climate and vegetation are Mediterranean. The annual precipitation ranges from ~700 mm in the southeast to ~1500 mm in the highest northeast of the Istrian Peninsula. Subregion IIIa is affected by the vicinity of subregion IIa, whereas subregion IIIb is barely influenced by the continental climate.

Regions I and II contain sizeable areas of dense forests, and these are indicated by letter F in Figure 1.

Materials and methods

Our experimental procedure has been explained in more detail in the first part of the report (12). Briefly, we sampled the surface layer of soil (0–10 cm) at 138 locations throughout Croatia in 2015 and 2016 following a recommended procedure for sampling and preparations for measurements (14). Radionuclide activity concentrations were determined by means of high-resolution gamma-ray spectrometry using a setup based on a high-purity germanium coaxial detector. We also applied corrections for self-attenuation (15) and coincidence summing effects (16). Detection limits were typically 2 Bq/kg for 40K and 0.3 Bq/kg for 137Cs.

Results and discussion

40K

Figure 2 shows the distribution of the measured A of 40K (T1/2 =1.3 billion years) in the surface layer of soil in Croatia. The range was 62–769 Bq/kg, and the average was Ā=423 Bq/kg. Globally, the reported median is 400 Bq/kg and the range 140–850 Bq/kg, and our results do not depart from those obtained elsewhere (including the rest of Europe) (17). Apart from some local deviations, there was a trend of A being above average in region I and below average in regions II and III. The A of 40K was generally lower in the regions underlain by karst, where limestone and dolomite bedrock dominate (18). This behaviour was opposite to that found for the 238U and 232Th decay chains (12).

Figure 2

Distribution of the A of 40K in Croatian soil

40K accounts for 0.012 % of the total K found in nature, and its distribution is actually the same as that of K. Since K is a biogenic element, its concentration in soil depends not only on underlying geology but also on the exchange of matter between soil and living organisms. While the A of 40K in the environment is comparatively large, which applies to our soil as well, this is usually not considered to be a threat, because 40K is a primordial, naturally occurring radionuclide that has been present throughout the evolution of life on Earth. Generally, 40K contributes significantly to absorbed dose rate for external exposure to ionising radiation from soil, but in our case, no health hazard can be associated with measured activity concentrations.

137Cs

137Cs, in contrast, represents a serious radioecological threat. The origin of most of the present-day 137Cs in Croatia is radionuclide release during the Chernobyl disaster in 1986. Much smaller contributions come from nuclear weapons tests and the Fukushima accident (9, 19). The distribution of the A of 137Cs in Croatian soil (in 2015–16) is shown in Figure 3. The measured values ranged from 0 (below the detection limit) to 205 Bq/kg, with Ā=25 Bq/kg.

Figure 3

Distribution of the A of 137Cs in Croatian soil

Generally, the A of 137Cs in our samples was lower than that of 232Th, 238U, 226Ra, 210Pb, and 40K. It did not depend on the distance from Chernobyl, as the lowest values (~4 Bq/kg in average) were measured in subregion Ib, closest to Chernobyl.

In order to understand results in Figure 2, one has to address several factors which might have affected the deposition of (initially airborne) 137Cs onto the ground and its fixation in soil. These factors are expected to be primarily related to atmospheric phenomena and geomorphological properties, including physical and chemical properties of soil (2). 137Cs is still abundant in the upper atmosphere and is deposited on the ground mainly via precipitation. When we compare results in Figure 2 with altitude and annual precipitation (20) at sampling locations, we notice that higher altitudes and more precipitation have A > Ā. The effect of altitude might be a consequence of the pinning of clouds by high mountains, so that 137Cs could be deposited onto the ground directly from the moist air, without actual precipitation. However, there are other variables involved, such as specific soil properties that favour 137Cs fixation on the surface layer, as discussed for chernozem in subregion Ib (2).

Due to similar electronic structures K and Cs have similar chemical properties, and at insufficient concentrations of K in soil, plants tend to take up chemically similar elements, including 137Cs (21). This is, in fact, why fertilisers rich in K have been used to reduce 137Cs uptake in heavily polluted areas (22). For the same reason the ratio R of atomic concentrations of 137Cs and total K is used as a measure of 137Cs potential to enter food chains through plant uptake. This ratio can be calculated as

where ACs and AK are activity concentrations of 137Cs and 40K, respectively. Figure 4 shows that R is lower in region I than regions II and III, which is a consequence of a relatively high AK and low ACs. It is of particular importance that R is the lowest in subregion Ib, where the majority of intense agricultural production is taking place in Croatia.

Figure 4

Ratio of the concentrations of 137Cs and K in Croatian soil

Absorbed dose rate

Internal exposure to ionising radiation due to soil radioactivity depends on numerous factors and has to be estimated on the basis of specific conditions. On the other hand, there is a well-established method of assessing external exposure from measured activity concentrations of radionuclides in soil and their gamma emissions. For naturally occurring radionuclides, the absorbed dose rate DN at 1 m above the ground can be calculated using the following equation;

where DN is expressed in nGy/h, and activity concentrations AK, AU, and ATh of 40K, 238U, and 232Th, respectively, are expressed in Bq/kg. Therefore, our results are sufficient for calculating the distribution of DN for Croatia, as shown in Figure 5. DN ranged from 23 to 131 nGy/h. Its average D̅ N was 63 nGy/h, which moderately exceeds the global median value of 51 nGy/h (17).

Figure 5

Absorbed dose rate due to naturally occurring radionuclides in Croatian soil

In the case of 137Cs, there is no universal equation for calculating the related absorbed dose rate DCs from ACs. Hence, we have to rely on estimates based on other results and models. There are arguments that 1 Bq/m2 due to 137Cs in soil results in about 0.001 nGy/h at 1 m above the ground (24, 25). Other studies indicate that only the topmost ~10 cm layer of soil contributes to DCs, since ACs significantly decreases with depth (26, 27) and self-attenuation effects (15). Total 137Cs activity in the upper 10 cm of a 1 m2 surface layer of soil therefore equals 0.1ρACs numerically, where ρ is soil mass density in kg/m3. This implies that

expressed in nGy/h. By taking into account measured ACs and ρ, we obtained results shown in Figure 6. DCs ranged from 0 to 24 nGy/h, and the average D̅ Cs was 2.7 nGy/h, which is considerably lower than DN. Therefore, while 137Cs in Croatia may still pose a certain threat regarding internal exposure due to its competition with K in food chains, its contribution to external exposure is comparatively small.

Figure 6

Absorbed dose rate due to 137Cs in Croatian soil

Figure 7 shows the distribution of the total absorbed dose rate (Dtot = DN +DCs) due to soil radioactivity in Croatia. It ranged from 24 to 136 nGy/h, and its average D̅ tot was 66 nGy/h. In most of region II, Dtot was above average, which was not surprising, because only activity concentrations of 40K were not elevated there. Dtot in the northwest of subregion IIIa, i.e., on the Istrian Peninsula, were also above average at most of the sampling spots (along the western coast in particular).

Figure 7

Total absorbed dose rate due to radionuclides in Croatian soil

Conclusions

In this second part of our report on uncultivated soil radioactivity in Croatia we focus on 137Cs, 40K, and absorbed dose rate. Since 137Cs and K have similar chemical properties, they compete in uptake by plants, which implies that their relative concentrations in soil are of interest for the propagation of 137Cs through food chains. Measured activity concentrations of 40K were the highest in the Pannonian part of Croatia, and 137Cs showed the opposite trend. Their ratio was the smallest in the Pannonian region, especially in its eastern part, where the majority of Croatian intense agricultural production has been taking place. Activity concentrations of 137Cs tended to be above average at high altitudes, in areas with dense vegetation and annual precipitation above average. By combining these results with those for the 232Th and 238U decay chains presented in the first part, we determined that the highest dose rates were in the Dinaric region and on the Istrian Peninsula. Since this study relied on long-lived radionuclides, our findings are likely to remain current for a prolonged period of time and may also serve for reference in case of accidental radioactivity bursts in the environment of Croatia.