Plutonium and thorium isotopes in the bottom sediments of some Mazurian Lakes (Poland)

Presented are results on the Pu and Th isotopes activity concentration found in the upper part of bottom sediments collected from a deep part of 29 lakes in N-E Poland by diving in 2000. Analyses of Pu isotopic ratios allowed for the discussion of Pu origin. Maximum percentage of 239+240Pu activity from Chernobyl fallout was 24%. Surface deposition of 239+240Pu was calculated. No relationship was found between the plutonium and main chemical matrix components of sample nor with the trophic status of the lake. Pu activities were weakly correlated with measured previously 137Cs activities.


Introduction
Study on the plutonium concentration in forest litter [1][2][3] conducted on the area of North-Eastern Poland, especially in Augustów Primeval Forest, showed the significant deposition of Chernobyl plutonium on some areas there. The activities were up to 30 Bq/m 2 for alpha emitting isotopes: 238?239?240 Pu [1] and up to 1 kBq/m 2 (in 1986) for 241 Pu [2]. Deposition of Pu isotopes from Chernobyl most likely took place in a form of small ''fuel like hot particle'' and therefore by nature was rather non homogeneous. Earlier deposition of 239?240 Pu alpha emitters from global fallout [4] was in average equal to about 60 Bq/m 2 , and remains of that for relatively short living (T 1/2 = 14.4 a) beta emitter, the 241 Pu was equal to about 100 Bq/m 2 in year 2001. Global deposition of 238 Pu from SNAP 9A satellite accidental re-entry in 1964 was small [4], for our latitude it was equal to about 2.5 Bq/m 2 .
One of the aims of present work was to get some data on Pu activity in examined upper part of bottom sediments for which radiocaesium data were obtained [5]. It seemed for authors interesting to get general recognition on level of Pu present in Mazurian lakes deep sediments for any further analyses in sediments of shallow waters. Another aim was to compare the data obtained from analyses conducted in the forest environment with the results from lake sediments. It seemed interesting for us also to complete the results of Pu and Th isotopes to those for radiocaesium [5] obtained before for the same bottom sediment samples.
Samples of bottom sediments were taken from 29 lakes in North-Eastern Poland. The 17 lakes were eutrophic [6] (Bro_ zane, Dębniak, Guzki, Jałówek, Jędzelewo, Krzywe, Lepaki, Oleckie, Pobojno, Rajgrodzkie, Rogale, Sawinda Wielka, Tajty, Tobołowo, Toczyłowo, Wi_ zajny, Zawadzkie), 9 were mesotrophic lakes [6] (Dobskie, Gaładuś, Garbaś, Jałowo, Płaskie, Serwy, Sumowo, Sunnowo, Szelment Du_ zy) and the 3 were oligotrophic lakes [6] (Aszarynis, Białe Filipowskie, Białe Wigierskie). In one case (Sumowo) two samples were taken in different parts of lake. All lakes are lying in the Eastern part of Mazurian lakes area. A map of North-Eastern Poland showing the lakes locations is presented in Fig. 1. Some data characterizing the lakes are presented in Table 1. The lakes have different hydrological status [6]. Six of them are endoreic (landlocked) lakes, 12 are transit lakes (i.e. there is at least one river coming in and one going out the lake) and 11 are outflow lakes (i.e. there is no river coming in). More details on the lakes themselves might be found in the monographic publication on Polish lakes [6]. To make it more easy the code numbers given to the lakes in that book are also presented in Table 1.
Sediments were taken in central part of each lake, offsite of littoral zone, where on bottom no vegetation was observed [5]. Diving were conducted till the depth enabling visibility sufficient for sample collection. For oligotrophic lakes it was nearly 20 m. In other cases it was less, not more than 10 m. Sediments were collected and taken to the surface as monolithic cores in plastic tubes of 10 cm diameter and 30 cm height [5]. In given site sampling repetitions nearby (from few to dozens of meters) were done. Cylinders were pressed into the bottom sediment, blocked from both ends and then taken out to the water surface. In majority of cases collected material had a form of a sludge. Samples were air-dried in laboratory for several weeks until they become a dry. Shrank sediment samples were sliced into 2 cm height layers, homogenised and taken for gamma spectrometric measurements [5]. The contraction factor (ratio of original total sample height to height of dried material) varied from 2 to 2.8. For 19 lakes only top 2 cm were taken. After performing analyses for all 29 lakes the 10 lakes for which additional profile taken close to first one was available (Białe Wigierskie, Gaładuś, Garbaś, Guzki, Płaskie, Rogale Du_ ze, Sawinda Wielka, Sumowo, Toczyłowo, Zawadzkie) we did repeated analyses for plutonium using more than one layer. Those samples were also sliced in 2 cm thick layers. Typically, three slices were obtained, however in two cases collected material allowed to use only two samples and in one case four layers were obtained. This repeated analyses were performed for determination of alpha emitting plutonium isotopes only. All presented data refers to dried material.
For quality assurance a reference material IAEA 300 Baltic Sea sediment [7] was analysed as well.

Procedures and measurements
Applied procedure for Pu and Th determination was an our implementation [1,3] of the radiochemical method [8] used by J.J. LaRosa and his co-workers for the analyses of Chernobyl Project samples in the IAEA Laboratories at Seibersdorf. The procedure started with incineration in 600°C lasting for at least for 8 h. Ten gram ash was taken from each sample for further radiochemical procedure. The plutonium tracers ( 236 Pu of 9.3 mBq activity at 1.03.2001 in first set and 3.4 mBq of 242 Pu in case of repetited samples) was added to each sample and subsequent wet mineralization in hot concentrated acids (HF, HNO 3 , HCl ? H 3 BO 3 ) was conducted. The oxidation adjustment procedure [1,8] was performed using hydrazine and NaNO 2 . Thorium and plutonium was separated by anionexchange chromatography using DOWEX 1 9 8, 200 mesh, Clform (Sigma Ltd.). Thorium was stripped off the column with 10 M HCl. Plutonium was later striped off the column with 0.1 M HCl-0.1 M HF solution. Thin alpha spectrometric sources were prepared form both Th and Pu stripes using NdF 3 co-precipitation method [8,9] using 50 nm pore diameter filters.
The alpha spectrometric measurements were performed using Silena AlphaQuattro alpha-spectrometer with silicon detectors (Canberra PIPS, 450 mm 2 ). Alpha-sources were measured in about 2 mm source to detector distance. Obtained spectra were analysed using home made ALF software [1]. Mean value of FWHM was equal to 68 keV with standard deviation of 16 keV.
In the aim of the determination of 241 Pu activity Pu alpha-sources were then dissolved in small amounts of hot concentrated HNO 3 with traces of H 3 BO 3 added, evaporated to dryness, re-dissolved in 1 M HNO 3 , mixed with liquid scintillation cocktail Wallac HiSafe3 and measured for the 241 Pu using Wallac 1414-003 Guardian liquid scintillation spectrometer. More details on procedure was presented previously [2]. The spectrometer can distinguish [10] scintillations caused by alpha particles from those caused by beta particles and therefore a simultaneous separate acquisition of alpha and beta spectra is possible. The total alpha emitting Pu isotopes activity present in source, which was already known here as the result of the alpha spectrometric measurement, was used to control the recovery. The tritium protocol [2,11] together with builtup library data was applied for all quenching and efficiency problems. Details of applied method of 241 Pu activity calculation were described already elsewhere [2]. Measurement time was set for 10 5 s.

Calculation of Chernobyl component
The analysis of plutonium activity concentration ratios ( 238 Pu to 239,240 Pu) allowed us to distinguish the Chernobyl and global fallout components of plutonium activities. Assumed 238 Pu to 239,240 Pu activity ratio in global fallout was equal to 0.03, and that for Chernobyl fallout was taken equal to 0.50 (decay corrected for year 2000), from 0.55 in fresh fallout, the percentage F o of Chernobyl plutonium can be calculated [1] as follow: where A 238 -activity concentration of 238 Pu, A 239 -activity concentration of 239,240 Pu. Global fallout component is equal to 100%-F o . The analysis of another plutonium activity concentration ratios, 241 Pu to 239,240 Pu also allows to distinguish the Chernobyl and global fallout components of plutonium activities. The 241 Pu to 239,240 Pu activity ratio in global fallout [4] was equal to about 12 in 1963, and that for Chernobyl fallout was equal 86 [2] in 1986. Corrected for the decay, at the year 2001 those ratio should be equal to 1.7 and 41.7, for global and Chernobyl fallouts, respectively. Thus, the percentage F 1 of Chernobyl plutonium based on 241 Pu activity can be calculated as follow: where A 241 -activity concentration of 241 Pu, A 239 -activity concentration of 239,240 Pu.
Global fallout component is equal to 100%-F 1 .

Calculation of surface deposition
The surface activities A surf in sediments were calculated from the activity concentrations A conc for radiocesium and then plutonium accordingly to the following formula: where q-density for dry material; V-volume of dry mass taken for analyses; S-surface area of a sampling cylinder. Jałówek  [12], where the maximum of Pu deposition was at about 3 cm, and whole deposit was in fist 6 cm of sediment. The variety of plutonium deposition observed in examined surface layer samples of bottom sediments reflects most likely the differences in the sediment accumulation rate, which depends on several parameters [13]. Results are lower than those for Obleije and _ Zuwintas lakes in Southern Lithuania [17] where maximum activities of 7500 ± 400 mBq/kg and 1800 ± 500 mBq/kg were found for 239?240 Pu and 238 Pu, respectively.
Regarding 241 Pu measurements only in one case, for Zawadzkie Lake, the 241 Pu was measurable by LSC without doubts, and the activity was equal to 50.5 ± 3.2 Bq/kg. Another result above the detection limit is for Garbaś Lake, where it was 5.3 ± 4.7 Bq/kg with the minimum detectable activity of 2.2 Bq/kg. In all other cases the activities were below the detection limits, which varied from 2 to 6 Bq/kg.
Since our initial study was limited to the top layer ony because of the amount of 236 Pu tracer present in laboratory after about three years we returned to the analyses using new tracer, the 242 Pu by NIST. We suspected that subsequent amount of plutonium can be present in deeper layers. The check this a repetition run of Pu measurements was performed for 10 lake sediments, for which we had spare profiles. Results are presented in Table 3. Again, the last row contains results for RM IAEA-300 obtained during this run. Our result is 3400 ± The minimum value for 239?240 Pu was found in depth 6-8 cm in Rogale Du_ ze lake (9 ± 7 mBq/kg) whereas maximum was found for top 2 cm of the same lake (3670 ± 350 mBq/kg) In general those results did not confirmed results obtained for the same lakes for samples collected in the same time, when we analysed only uppermost 2 cm. It means, that very large variety of Pu content can be observed in deep sediments at distance of few meters only. Perhaps the sediments were so much not settled that the diving itself makes disturbance to it. The results confirmed that subsequent plutonium activity is present in deeper layers. In all cases but the lake Garbaś, where majority of Pu was found in top 2 cm. However in this lake the activity concentration in top leayer of second set is more than five times lower that it was found in first run. Table 4 shows the activity ratio of 238 Pu to 239?240 Pu and 241 Pu to 239?240 Pu as well as the estimated percentages F o (see Eq. 1) and F 1 (see Eq. 2) of Chernobyl-origin Pu in samples obtained from first analyses. Maximum 238 Pu to 239?240 Pu activity ratio does not exceed 15% (Zawadzkie Lake), what yields in F o = (24 ± 2) %. Such percentage of Pu Chernobyl component in this sample was confirmed by the result from 241 Pu analyses, where F 1 = (29 ± 3) %. For 3 other samples Chernobyl component was about 15%, for two another it was between 5 and 10%. For those of remains samples, for which it was possible to calculate F o value, it was below 5%. General not high content of Chernobyl fallout was rather surprising for us and perhaps unexpected, especially if one would compare present results to those for the forest litter samples, in which up to 100% of Chernobyl fraction was observed in surface layer of forest litter collected in North-Eastern Poland [1]. Such sort of discrepancy between our results in two environments suggest, that hot particles which contained Pu from Chernobyl were somehow eliminated from the examined part of sediment, or rather even never were there. One of possible way of elimination is a relatively high velocity of migration down the profile, similar to that found for forest soil from Poland [3,14]. In the case of forest soil it was explained by a greater diameter of particles which makes them more mobile with water flow downward. Another, and even more possible in our opinion explanation of lower than expected levels of Chernobyl plutonium is the difference in the aerodynamic properties of tree canopies and lake which might influence the intensity of hot particle fallout on different surfaces.
Comparison of total and Chernobyl-origin inventory found for first 2 cm layer during first analyses and for deeper profile is presented in Table 5. This comparison shows clearly lack of reproducibility of results obtained. There is some similarity of results only in case of three lakes (Sumowo, Sawinda Wielka, Płaskie). In five cases the inventory of Pu found in top 2 cm is larger than then found in whole deeper profile. This shows problem related to inhomogeneity of sediments which has form of not well settled sludge. In particular, relatively high content of Chernobyl origin Pu in lake Zawadzkie was not confirmed.
It seems that study should be repeated with less disturbing core sampling procedure (not man-diving but core sampler) and with at least five repetitions of sampling from one location.
Results for thorium activities are presented in Table 6. They were calculated under assumption of complete, equal to 100%, recoveries of thorium. They vary a lot, ranging (for 232 Th) from 0.09 ± 0.01 Bq/kg (Białe Wigierskie Lake) to almost 11 Bq/kg for a couple of lakes like Aszarynis, Szelment, Sumowo, Sunnowo. Despite a big difference in activities between the samples, isotopes from  Fig. 2). The Persons correlation factor r = 0.947 with significance level p \ 0.0001. Thorium isotopes are in almost constant proportion, close to one to one for all samples but few, namely those of Dobskie and Pobojno Lakes, which are significantly above the correlation line and Aszarynis and Dębniak Lakes, which are significantly below the correlation line. This results suggests almost uniform geological structure of examined bottom sediments. Two samples from one lake (Sumowo) lake had identical thorium isotopes activity and almost identical 239?240 Pu activity, but they differ each to other for 238 Pu activity. Since the main source of this isotope is Chernobyl, it might be caused by scattered, particle nature of Chernobyl Pu fallout. There is a complete lack of correlation between 232 Th and 239?240 Pu (see Fig. 3), what might suggests, that the origin of radionuclides (which is different for those nuclides) but not the accumulation properties of sediments (expected similar for both radionuclides) is the main factor which governs the observed variety of activity levels.
A significant correlation was found for surface activity between the activity of radiocaesium determined previously [5] and 239?240 Pu (Fig. 4). The Persons correlation factor r was equal to 0.582 with the level of significance p = 0.001. Less significant (r = 0.337, p = 0.077) correlation calculated for activity concentrations of those   (Fig. 5). Low correlation characterized by R = 0.350 and p = 0.068 was found between the pH in sediment measured in water and the 239?240 Pu in sediment (Fig. 5). Perhaps it is caused by a higher mobility of humic acid in more basic environment. The humic acid is well known complexing agent for plutonium [15,16]. No correlation was found between the Pu content and basic matrix content of samples like: SiO 2 , Ca, Mg, Na, K, C, P, Fe.
In principle, like for radiocaesium, differentiation of plutonium deposition in the lake bottom sedimentations may be due to the morphometric parameters of the lakes and the type of the drainage area. A differentiation has been made on the basis of the trophic levels of lakes. It was found previously for radiocesium [5]

Conclusions
In average about 30% or 40% of expected amount of global fallout plutonium was observed in examined 2 cm slice of top layer of bottom sediment from 29 Mazurian lakes. Chernobyl plutonium was not so much abundant as was expected from previous investigations conducted on this area using forest litter samples, reaching maximum content of 24%. Moreover, later performed analyses of 10 sediments from lakes, where deeper layers were also available did not confirmed initial analyses revealing even lower Chernobyl content and indicating problems with homogenity of Pu content in bottom sediments. This was result mainly of very wet, sludge form. This difference between forest litter and lake sediment seems to be the most interesting observation of present work. One of possible explanation is the difference in the aerodynamic properties of forest and lake which might influence the intensity of hot particle fallout.
Unlikely the radiocaesium [5], the plutonium content in top layer sediment of lakes seems to have in presented here cases no clear relationship with the trophic status of lake. However, a weak but significant correlation between plutonium and cesium surface activity was noticed. Hardly no correlation was observed between plutonium activity and any of known for given lake parameter, among them such like main content of matrix. Only a weak correlation was found between Pu content and pH of sediment determined in water.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creative commons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.