The paper presents the results of measurements the specific concentration of lead-210 for six peat profiles representing four peat bogs from two regions of Poland and the problem of creating age-depth models. For the construction of age-depth models, it is proposed to use mathematical functions, the best fit to the measured activity. The F-statistics were used as a measure of the match quality. The obtained models are visualized in two ways – showing the age calculated on the basis of direct measurements of activity and indicating points that are the results of the used approximation. Such visualization is important to clearly distinguish the places of the age-depth model that result from the measurements of activity from those places that are the result of the approximation used.

This paper proposes and tests the application of activity modelling for the cores shorter than the range of the lead method. The paper also outlines the limitations and potential dangers related to the interpretation of core dating results i) obtained by using the activity approximation resulting in the smoothing of the age-depth profile, and ii) for the cores of a length smaller than the depth of the presence of the unsupported lead.

Additionally, the proposed models were compared with the models obtained by using SIT, Mod-Age, OxCal, Clam, and Bacon programs.

#### Keywords

- Pb-210 chronologies
- young peat profiles
- Pb dating
- CRS
- age-depth model

The dating method using the radioactive ^{210}Pb is widely applied for dating geological sediments in the environmental research. It is one of the fundamental and leading methods allowing the determination of the chronology of samples whose ages do not exceed approximately 200 years (the range of the method is determined by the half-life of ^{210}Pb – T_{1/2} = 22.6 years).

Historically, the investigations of young sediments for which the concentration of the radioactive isotope of lead was measured, were conducted mainly with the goal of reconstructing the changes of the rate of sedimentation of ice, marine, lacustrine, and fluvial sediments. The lead method was successfully applied for water basins of varying area, depth, and rate of sedimentation (e.g. Luque and Julià, 2002; Xiang

However, the lead method finds a much wider range of applications, for example for the dating of peat sediments. The last decade has seen a dynamic development of research of peat bogs. They are declining reservoirs that had recorded the history of industrialization which is one of the symptoms of anthropopression (Turetsky ^{210}Pb in soils, air samples, ice cores as independent reservoirs or as complementary measurements in the mentioned sedimentary environments and peat bogs are an essential aspect of geochronometric and climatological investigations (Ebaid and Khater, 2006; Olid

The goal of the presented work was to determine precise age-depth models (e.g. with 1 cm resolution) for the cores where the samples for activity measurements were collected with a lower resolution. The common approach is to divide the core into 1 cm slices and to determine the specific activity of ^{210}Pb in each of them (Putyrskaya

Currently, the typical approach in ^{210}Pb dating method is to calculate the age of the selected sediment layers based on the activity measurements for these layers. Appropriate models are used for this – e. g. CRS – Constant Rate of Supply (Appleby and Oldfield, 1978; Liu

Increasingly, in the next step, the age for the intermediate layers is modelled using the Bayes or Monte Carlo method (e.g. Blaauw, 2010; Blaauw and Christen, 2011; Hercman and Pawlak, 2012, Hercman

This work offers a slightly different approach. Namely, it introduces matching the mathematical function that best reflects the values of activity measurements (see Section 3). Based on the selected function, it is possible to calculate the activity for all intermediate layers with the necessary resolution. This approach is less statistically expanded than the ones mentioned above. On the other hand, using ready-made computer programs to create age-depth models may create the threat of too much automation in the application of the the ^{210}Pb dating method to create such models. In the author's opinion, this stage of dating requires special care, regardless of whether we approximate the activity or the age. The performed analyses allowed determining an age model for non-standard profiles that are difficult to interpret. For this purpose, a classical CRS dating model was used. The paper also discusses the possibilities and limitations connected to the interpretations of the ages obtained for the profiles where sampling encompassed only the layers younger than the range of the method.

The investigated peat bogs, Bagno Bruch and Bagno Mikołeska (^{°}C. In winter, the average temperature is about –1^{°}C, and in summer it is +18^{°}C approximately (Fiałkiewicz-Kozieł

The cores from these peat bogs were collected in 2016, in the central parts of the investigated areas. For comparison purposes, two cores were acquired at a distance of 2 m from each other, from each peat bog – the cores

The investigated peat bogs (^{°}C and +17.9^{°}C, respectively (Kaszewski, 2002). The average annual precipitation in the park is 551 mm (Kaszewski, 2002; Mieczan, 2007). The peat bogs are supplied with water mainly by precipitation. However, groundwater supply cannot be excluded, especially in the case of Bagno Bubnów peat bog, drained by an unnamed watercourse. The cores from these peat bogs were collected in 2014.

^{°}21’17”; E: 23^{°}16’26”) has a surface area of 2309 ha. It is a calcareous fen. Due to their rarity in Poland, the most interesting communities are connected to the calcium carbonate present in the soil. These are the species such as

The profile

^{°}27’01”; E: 23^{°}13’23”), which has a surface area of 213 ha, is a continental raised bog. The species such as

There are two basic methods used to determine the specific activity of ^{210}Pb – namely alpha spectrometry and gamma spectrometry of the radiation emitted by the sample (e.g. Sikorski and Bluszcz, 2003, 2008; Le Roux and Marshall, 2010).

The most frequently applied approach is an indirect measurement of the ^{210}Pb decay product ^{210}Po (T_{1/2}=138 days) using alpha spectrometry, also applied in the presented work.

The activity was measured using a quad alpha spectrometer (Canberra model 7404) with Passivated Implanted Planar Silicon (PIPS) detector. The sensitive area of the detector is 300 mm^{2}, and its energy resolution is ^{210}Po (E = 5.308 MeV) present in the sample and ^{208}Po (E = 5.105 MeV) added as a standard. Because the half-lives of both polonium isotopes are relatively short (T_{208} = 2.898 years, T_{210} = 138.4 days), the results of the measurements are corrected into ^{208}Po radioactive decay from the moment of its calibration, and ^{210}Po disappearance from the moment of polonium extraction from the sediment.

All measurement results are presented according to ISO/IEC Guide 98–3 (2008), so the uncertainties are given in brackets. The decimal separator was not used when giving the uncertainties. Therefore, one needs to keep in mind that the uncertainty is always related to the last two decimal places of the measurement result.

The samples were mineralised and chemically processed in order to transfer polonium into the solution. Polonium was subsequently deposited through spontaneous deposition on silver. In the first step of preparation, ^{208}Po standard (manufactured by AERE Harwell, UK), of known specific activity, was added to 1 g of the dry sample. Next, the sample was treated with hot, concentrated HNO_{3} (90^{°}C, 65%), and placed in a sealed container. In order to separate the desired fraction, the mixture was centrifuged for _{3} was eliminated from the sample and replaced by an HCl solution. The remainder of the last rinse was dissolved in 40 ml 0.5 M HCl with an addition of 0.5 g of hydroxylamine (H_{3}NO^{⋅}HCl) and 0.5 g sodium citrate (C_{6}H_{5}O_{7}Na_{3}^{⋅}2H_{2}O) in order to achieve the reduction of inhibitors. Polonium was deposited on a silver disc (Flynn, 1968; Sikorski and Bluszcz, 2003; De Vleeschouwer ^{°}C for

The obtained radioactive sources underwent alpha spectrometric measurements. This methodology was used for all measured peat samples.

For the ^{210}Pb age calculation, the CRS model was used (Appleby and Oldfield, 1978; Sikorski and Bluszcz, 2003). This model assumes that the ^{210}Pb flux remained constant through time, but the peat accumulation rate may have changed.

The results usually are not given as a function of depth, but rather the depth mass, in order to eliminate the possible effect of sediment compaction or the changes in its density. Such approach was also used in the presented work.

In the simplest case, the activity versus mass depth should be described by the exponential function. Unfortunately, exponential fit to the plots of the specific unsupported ^{210}Pb activity vs. mass depth for examined peat profiles did not provide satisfactory results (F-statistics starting from the value = 30 and worse for some other functions).

In the current study, it was proposed to approximate the activity for each depth by matching the mathematical function that best reflects the values of activity measurements. The quality of fitting was tested by the observed value of the F-statistics. In the next step, the age of each layer was calculated using the CRS dating model.

Matching the best mathematical function was determined for each of the cores independently, by the application of the measured activities. The selection of the matching function is crucial not only for the modelling of ^{210}Pb concentration but also for the obtained age-depth model.

The program ^{210}Pb and their uncertainties), the program matches these functions for the input data for each sediment core, by the application of the weighted least-squares procedure. The results of fitting are presented as a statistically ranked list of candidate equations. The range of the search was limited to the models consisting of at most four parameters (

For the available cores the activity was measured in every second analysed layer. As a result, 70 activity values were obtained for peat bogs Bagno Bruch (BB) and Bagno Mikołeska (BM), and 24 for Bagno Bubnów and Durne Bagno (TW).

^{210}Pb

The specific activity of the supported ^{210}Pb for each of the BB and BM cores was determined for the measurement from the deep layers (for the depth of 50 cm), and from the measurements of the deepest layers, where the activity reached the value of the deep layer. These values were averaged using a weighted mean. For sites BB3, BB4, BM3, and BM4 these activities A_{a} with their uncertainties u(A_{a}) were determined to be 12.2(81); 15.3(83); 20.8(95); 23.4(96) Bq/kg respectively. The results obtained for the cores within each peat bog do not differ significantly – they agree within measurement uncertainties. The values of the Student’s t-test:

where:

t_{S}

A_{a1}, A_{a2} ^{210}Pb (Bq/kg)

u(A_{a1}), u(A_{a2}) ^{α} = 0.05, whose value, in this case, is 2 (uncertainties of A_{a1}and A_{a2} are known from the law of decay hence we may assume that the testing variable has the normal distribution).

Therefore, the results of the t-test allow the calculation of the weighted average for each of the investigated peat bogs. For the cores taken in close vicinity, the following values were obtained: 13.7(82) Bq/kg for BB and, 22.1(95) Bq/kg for BM.

For TW1 and TW2 the activity of supported ^{210}Pb was determined for measurements performed for deep (45 cm) layers of the peat profiles. They are equal to 32(11) Bq/kg and 33(11) Bq/kg, respectively. Although an agreement within individual uncertainties is seen, the results were not averaged because they were taken from two different peat bogs and slightly different sedimentary environments (see Section 2).

It is interesting to observe that the best extrapolation results were obtained for the functions described by the same formula for all cores:

where:

^{2})

^{210}Pb (Bq/kg)

The only difference is the value of

The basic criterion for determining the match quality was the value of the F-statistics.

The values of F-statistic for the investigated peat profiles are as follows:

Bagno Bruch: F = 192 (BB3) and F = 552 (BB4),

Bagno Mikołeska: F = 251 (BM3) and F = 195 (BM4),

Bagno Bubnów F = 571 (TW1),

Durne Bagno F = 169 (TW2).

The Fisher-Snedecor probability distribution was used to calculate the likelihood of an equal or higher F value occurring by chance, and as a result the following probabilities were obtained:

Bagno Bruch: P = 7.7·10^{–10} (BB3) and P = 1.5·10^{–12} (BB4),

Bagno Mikołeska: P = 7.0·10^{–10} (BM3) and P = 3.3·10^{–9} (BM4),

Bagno Bubnów P = 2.3·10^{–9} (TW1),

Durne Bagno P = 2.9·10^{–7} (TW2).

In this case, the null hypothesis is that the observed relationship occurred by chance. These probability values are extremely small, hence we may reject the null hypothesis and assume that the model fits properly with the experimental results.

For mass depths _{i}_{i}

Subsequently, uncertainty of these activities _{i})

where

The presented results of specific activity measurements originating from profiles BB3, BB4, BM3, BM4, TW1, TW2 are shown in

Taking into account the fact that

Parameters a, b, c of the functions obtained in TableCurve program.

Core | |||
---|---|---|---|

BB3 | 6.093(45) | –0.391(52) | 2.5 |

BB3 shortened | 6.148(52) | –0.486(77) | 2.5 |

BB4 | 6.175(26) | –0.499(38) | 2.5 |

BB4 shortened | 6.220(26) | –0.618(45) | 2.5 |

BM3 | 6.063(39) | –0.420(61) | 2.5 |

BM3 shortened | 6.062(49) | –0.417(80) | 2.5 |

BM4 | 6.022(36) | –0.198(35) | 3.0 |

BM4 shortened | 6.018(40) | –0.187(41) | 3.0 |

The above-described procedure of measurements, activity approximation and age determination was applied to the sites located in close vicinity – peat bogs BB and BM.

The used methodology of sediment measurements performed for every second layer and approximation of activity (or age) for the missing layers along the profile has its limitations. However, one has to bear in mind that such approach may lead to missing some changes occurring in the profile, blurring of information, and smoothing of the age-depth model. In such cases, extreme events are not visible. It is particularly dangerous for the profiles

with low sedimentation rate. On the other hand, because the number of dated profiles can be doubled in this way, this approach seems to be justified.

The ages of the layers sampled from the same depths were first compared within each peat bog, and, subsequently, between both peat bogs, BB and BM. It was observed that the ages for the pairs of cores from the same peat bog (BB3 and BB4 as well as BM3 and BM4) are in agreement, whereas for layers located at 46 depths the agreement is within a single uncertainty, and for the remaining 12 – within two uncertainties. A similar agreement is observed between the peat bogs.

The analysed peat bogs are close to each other, and hence the vegetation conditions of the peat are similar, and a similar rate of accumulation would be expected. The agreement of the results obtained for two independent peat bogs suggests a similar deposition history for both localities.

In the cases where the oldest layers of the investigated profile are younger than the range of the method, this fact has to be taken into account, by appropriately modifying the age model. Such situation is often seen in young, mostly artificial (anthropogenic) reservoirs, when their age is known. In the CRS model, the activity for the whole thickness of the sediment has to be assumed in a way to yield the age of the reservoir for the sample found at its bottom (Sikorski, 2003, Sikorski and Bluszcz, 2003). However, the determination of the activity of the supported ^{210}Pb in such situation is very difficult.

Similar problem is encountered when a core that is too short has been collected. In this situation, an approximation of the activity for deeper samples is necessary. For these profiles, the approximation method described in the Methods section was used. The values obtained in such a way can be then used in the age model. Such problem was encountered in the profiles TW1 and TW2 from the peat bogs from Poleski National Park, where the measurements of activity and concentration of supported ^{210}Pb revealed that radioactive lead reaches deeper than 20 cm, a depth to which the core was taken. For obvious reasons, the ages obtained for layers deeper than 20 cm cannot be a result of direct alpha spectrometric measurements. Modelling has to be applied, and the results cannot be treated with equal validity as in the first case.

It has to be stressed that these are only estimates and should be treated with caution. In environmental research, particular attention and restraint have to be exercised in the interpretation of the results obtained for each layer. Such dating should be treated roughly, as one way of providing the information about the range of the lead method, or determining the average rate of peat accumulations for the whole core.

Additional tests were conducted for the short cores, in order to compare the measured and estimated level of supported ^{210}Pb. The measurement results for BB and BM cores were used for this purpose, but they were modified in such a way that short cores were created from them, by the removal of the results for the layers deeper than 20 cm (1.7 g/cm^{2}) from the set of measurement data. For these new data sets, the activities were extrapolated independently for each profile, in ^{210}Pb, determined in such a way for BB3, BB4, BM3, and BM4 stands, are presented in

Comparison of measured and approximated values of supported ^{210}Pb.

Core | Depth (cm) | Measured values of supported ^{210}Pb (Bq/kg) |
Approximated values of supported ^{210}Pb (Bq/kg) |
t-test | p-value |
---|---|---|---|---|---|

BB3 | 29.5 | 12.2(81) | 14.7(91) | 0.21 | 0.84 |

BB4 | 29.5 | 15.3(83) | 13.8(90) | 0.12 | 0.91 |

BM3 | 31.5 | 20.8(95) | 22.9(99) | 0.15 | 0.88 |

BM4 | 29.5 | 23.4(96) | 23.1(98) | 0.02 | 0.98 |

The

Two methods of determining the age were compared in the presented study. In the first, classical approach, the age is calculated using CRS model directly for the measured activity. In the second approach, the approximated values of the activity were used as an input for the CRS model.

The ^{c}

The insufficient number of dating results may result in the need to use the methods based on Bayesian approach or Monte Carlo simulations to determine the age of non-measurable samples. The approximation of the activity enables the use of the lead method as a robust and comprehensive technique of dating.

In practice, the age values obtained for the youngest sediments by the lead method are often combined with the results of the dating of older profile layers (e.g. using the ^{14}C method) to create comprehensive age-depth models. In this situation, the results of ^{210}Pb dating are modelled for the second time. In the first step ^{210}Pb ages for each layer were modelled based on activity measurements and in the next step these dates are used together with ^{14}C dates to create new age-depth model – again for the same layers.

The graphs of the age-depth models calculated on the basis of the CRS model used in lead dating are shown in two variants – with the indication of points only for the depths from which actual samples were taken (

The dating results obtained by two methods are consistent within the limits of the measurement uncertainties. However it can be observed that ages calculated by using the chosen mathematical function describing the changes of activity are younger than the ones calculated by using measured activities for BB and BM cores before 1900 AD. This observation is in line with the observations of other researchers using the ^{210}Pb lead method (e.g. Gunten

Average growth rates of the investigated peat bogs.

Average growth rates (mm/y) | ||||
---|---|---|---|---|

Core | by 1967 | from 1967 to present times | ||

weighted average | weighted average | |||

BB3 | 1.03(12) | 1.10(13) | 3.41(19) | 3.21(18) |

BB4 | 1.21(15) | 3.04(18) | ||

BM3 | 1.15(13) | 0.99(12) | 3.26(18) | 3.40(18) |

BM4 | 0.87(11) | 3.54(18) | ||

TW1 | 1.34(16) | |||

TW2 | 1.41(17) |

Finally, the obtained age-depth model for each core was compared with the models generated by the following software: SIT (Carroll and Lerche, 2003), OxCal (P_Sequence procedure, Bronk Ramsey, 2008), Clam (Blaauw, 2010), Bacon (Blaauw and Christen, 2011), and Mod-Age (Hercman and Pawlak, 2012).

By the application of OxCal, Clam and Bacon software, the age-depth model was determined on the basis of the age of the sample and its location (depth) in the sediment profile. The approach presented in this paper (activity extrapolation), SIT and Mode-Age programs enable the direct determination of the chronology on the basis of the activity of unsupported ^{210}Pb and the location of the sample within the profile. The comparison is presented in

The values with normal distribution of uncertainty were applied as the input data in the used models, with the exception of Mod-Age software that accepts the dates with other distribution. It enabled age-depth model construction using the real lead age distributions, for each depth taking account its uncertainties. The actual age distribution is obtained as a result of dating by using Monte-Carlo simulations of activity assuming its normal distribution (Hercman and Pawlak, 2012).

All the models used in the presented research yielded comparable results (

The chronologies obtained by independent methods for the older samples are generally comparable within the limits of the single measurement uncertainty with the results obtained from the application of activity approximation, with one exception. The deviation from this conformity is the dependence obtained by the application of SIT model. In this case, the dates are significantly rejuvenated, even up to 100 years. Only for the SIT model time marker was introduced: the year when the zinc and lead plant was built – see the previous chapter (1967 AD for the following depths: BB3: 16 cm, BB4: 15 cm, BM3: 16 cm, BM4: 18 cm). The data obtained without this limitation were even younger, by additional 50 years, and unrealistic in relation to the other models. The SIT model, in the presented research, demonstrated the tendency of linearization of the generated age profile for the discussed measurement sites. It causes difficulties in the application of this model, especially for the sites where the deposition rate changes over time. The limitations of SIT model were discussed by other authors (e.g. Abril, 2015), and the discrepancies between the results of SIT and CRS models for older samples were described by e. g. Tylmann

The general tendency to be younger of the results obtained by the application of Bacon software can also be observed. It is also the result of the linearization tendency of the age-depth model – especially for TW cores and the older layers of the other cores (

The age characteristics obtained for Mod-Age model are moved in the direction of older dates, in comparison to the ages obtained after the activity approximation. The most probable reason for this divergence is the fact that this model takes into account the skewness of ^{210}Pb age distribution.

The paper presents the results of the analysis performed for six peat profiles representing two independent research sites. These sites have different geographical location and thus require independent interpretation of the results.

The conducted research can be summarised as follows:

A frequent practice of the users of lead dates is the sole presentation in publications of the final dating results, although such a form is insufficient. The dates should be presented without the points obtained by approximation, and when such points are included it should be clearly indicated how many samples were measured and at which depths they were taken.

In the case of multi-proxy analyses, the best practice is the measurement of activity for each of the analysed layers along the core. Performing measurements for every second layer is acceptable. Naturally, such simplification generates a smoothing of the age-depth profile in the final result, and thus extreme events can be missed.

The lead dating method, taking advantage of the activity approximation, similar to the methods using age approximation, is a tool that enables obtaining an age for any layer of the profile that lies within the age range.

The employment of the activity approximation method enables the estimation of activity, even the supported activity for the cores shorter than the range of the method. The tests conducted on the artificially shortened cores confirmed the usefulness of this method. However, the awareness of the associated limitations is necessary. Such ages are only estimates and should be treated with caution. For example, as the information about the range of the lead method or average deposition rate for the whole core, especially for monotonic profiles with a constant sedimentation rate.

It is worth noticing that for six analysed cores, the function that best fits the measurement activities has the form: ^{c}

A comparison of the activity approximation to creating the age-depth model with the models obtained using the Monte Carlo simulation and Bayesian theorem (using Mod-Age, OxCal, Clam, and Bacon software), showed their compatibility within the limits of a single uncertainty. One should note that among the tested software, Mod-Age program is the only one that offers the possibility of taking into account the skewness of age distribution, introduced into the model.

The results of modelling with the application of SIT software are in accordance with the results obtained from other programs for the upper layers of the researched cores only. For the lower layers, the differences become more pronounced with the depth. It suggests that SIT software applies too much linearization.

#### Parameters a, b, c of the functions obtained in TableCurve program.

Core | |||
---|---|---|---|

BB3 | 6.093(45) | –0.391(52) | 2.5 |

BB3 shortened | 6.148(52) | –0.486(77) | 2.5 |

BB4 | 6.175(26) | –0.499(38) | 2.5 |

BB4 shortened | 6.220(26) | –0.618(45) | 2.5 |

BM3 | 6.063(39) | –0.420(61) | 2.5 |

BM3 shortened | 6.062(49) | –0.417(80) | 2.5 |

BM4 | 6.022(36) | –0.198(35) | 3.0 |

BM4 shortened | 6.018(40) | –0.187(41) | 3.0 |

#### Average growth rates of the investigated peat bogs.

Average growth rates (mm/y) | ||||
---|---|---|---|---|

Core | by 1967 | from 1967 to present times | ||

weighted average | weighted average | |||

BB3 | 1.03(12) | 1.10(13) | 3.41(19) | 3.21(18) |

BB4 | 1.21(15) | 3.04(18) | ||

BM3 | 1.15(13) | 0.99(12) | 3.26(18) | 3.40(18) |

BM4 | 0.87(11) | 3.54(18) | ||

TW1 | 1.34(16) | |||

TW2 | 1.41(17) |

#### Comparison of measured and approximated values of supported 210Pb.

Core | Depth (cm) | Measured values of supported ^{210}Pb (Bq/kg) |
Approximated values of supported ^{210}Pb (Bq/kg) |
t-test | p-value |
---|---|---|---|---|---|

BB3 | 29.5 | 12.2(81) | 14.7(91) | 0.21 | 0.84 |

BB4 | 29.5 | 15.3(83) | 13.8(90) | 0.12 | 0.91 |

BM3 | 31.5 | 20.8(95) | 22.9(99) | 0.15 | 0.88 |

BM4 | 29.5 | 23.4(96) | 23.1(98) | 0.02 | 0.98 |

^{210}Pb-based chronology in recent-sediment cores? ^{210}Pb-based chronology in recent-sediment cores?

^{210}Pb and ^{137}Cs radiometric techniques. ^{210}Pb and ^{137}Cs radiometric techniques

^{210}Pb to the sediment. ^{210}Pb to the sediment

^{210}Pb chronologies for two shallow Argentine lakes in Southern Pampas. ^{210}Pb chronologies for two shallow Argentine lakes in Southern Pampas

^{210}Pb in environmental samples. ^{210}Pb in environmental samples

^{210}Pb method. ^{210}Pb method

^{137}Cs and ^{210}Pb radiochronology to estimate recent changes in sedimentation rates in Kushiro Mire, Northern Japan, resulting from land use change and river channelization. ^{137}Cs and ^{210}Pb radiochronology to estimate recent changes in sedimentation rates in Kushiro Mire, Northern Japan, resulting from land use change and river channelization

^{210}Pbdating of peat cores. ^{210}Pbdating of peat cores

^{210}Pb determined by gamma-spectrometry: progress and problems. ^{210}Pb determined by gamma-spectrometry: progress and problems

^{210}Pb chronologies and accumulation rates of sediments and peat bogs. ^{210}Pb chronologies and accumulation rates of sediments and peat bogs

^{210}Pb measurements). PhD thesis, Silesian University of Technology Gliwice: 130pp.

^{210}Pb method to date sediment of human–made lake Kozłowa Góra. ^{210}Pb method to date sediment of human–made lake Kozłowa Góra

^{210}Pb method to model sedimentation in artificial retention reservoir. ^{210}Pb method to model sedimentation in artificial retention reservoir

^{210}Pb activities in the lake deposits. ^{210}Pb activities in the lake deposits

^{210}Pb dating results with varve chronology and independent chronostratigraphic markers: Problems and implications. ^{210}Pb dating results with varve chronology and independent chronostratigraphic markers: Problems and implications

^{137}Cs and ^{210}Pb measurements. ^{137}Cs and ^{210}Pb measurements

^{210}Pb geochronology. ^{210}Pb geochronology