Institute of Physics Institute of Botany Lithuania The
Institute of Physics, Institute of Botany Lithuania The LIETDOS-BIO assessment approach to Environment protection from ionizing radiation T. Nedveckaite, D. Marciulioniene Institute of Physics, Savanoriu 231, Lt-2058, Vilnius, Lithuania, E-mail: tatjana@cablenet. lt; Institute of Botany, Žaliųjų ežerų 49, LT-2021 Vilnius, Lithuania, E-mail: dmarciulione@yahoo. com
Ignalina NPP with decommissioning, spent fuel storage and waste disposal problems; Maisiagala RADON type radioactive waste repository amongst others with earlier (1963 -1973) military radioactive waste from Lithuania, Kaliningrad and Belarus. Estonia Russian Federation Latvia Lithuania Kaliningrad INPP Maisiagala repository Poland Belarus
Ignalina NPP, Druksiai Lake – cooling pond, spent fuel storage RSR – rain sewage release, Intake - cooling water intake, Release - cooling water release, PFo. HE - purification facility of the household effluents of INPP and Visaginas, HSS – household sewerage release after biological treatment, SFSF - Spent Fuel Storage Facility.
The sampling points in Ignalina NPP cooling pond-Druksiai Lake
The spatial pattern of activated corrosion products 54 Mn and 90 Co in bottom sediments of Druksiai Lake. The highest activity concentrations correspond to sampling points: 2 -sewerage outflow and 4 -cooling water outflow. Other monitoring radioniclides and heavy metals: 54 Mn, 60 Co, 137 Cs, 90 Sr, 210 Pb, 210 Po, 238 U, 226 Ra, 232 Th Pb, Cd, Cr-VI, Ni, Cu, Zn-II, Zn
The projected Ignalina NPP decommissioning low-level repository sites
MAISIAGALA radioactive waste repository
Tritium in underground water in close vicinity of the repository reservoir (January-March 2007) Borehole No.
LIETDOS and LIETDOS-BIO – a software to calculate the radiological dose for human and non-human biota due to the presence in LITHUANIA radiation safety problems
LIETDOS –segmental atmospheric diffusion-convection model-code tested using data on I-131 released from Hanford (BIOMASS)
LIETDOS-human exposure model used for the regulating document in Lithuania LAND 41 -2001 : “Limitation of Radioactive Discharges from nuclear facilities and Order of Issuance of Permits for Discharges” • Nedveckaite T. , Motiejunas S. , Kucinskas V. , Mazeika J. , Filistovic V. , Jusciene D. , Maceika E. , Morkeliunas L. , Hamby D. M. Environmental releases of radioactivity and the incidence of thyroid disease at the Ignalina Nuclear Power Plant. Health Physics, Vol. 79, No. 6, 666 -684 (2000)
LIETDOS-BIO model a tool for calculating radiation doses to Freshwater and Terrestrial Biota with input file that is subsequently read by MCNPX General/site-specific parameters Source (nuclide) Ecosystem Preparing MCNPX input file for DCC calculations Crystal Ball MCNPX code Internal radionuclide concentration DCC calculations Internal dose (rate) Organisms/reference organisms Dynamic/steady-state transfer modelling Radionuclide concentration in medias External dose (rate) Weighting Total dose (rate)
LIETDOS-BIO code • LIETDOS-BIO code is designed to be consistent with MCNPX code. MCNPX is commonly used general purpose radiation transport code that tracks all particles at all energies for internal and external exposure of any geometry organisms calculation. • LIETDOS-BIO code is designed to be consistent with Crystal Ball software (Crystal Ball is a software for uncertainty analyses) for statistical calculation approaches. EXPERT CONTROLLED INFORMATION This is a limited distribution unclassified documents with code and may not be given to other persons. LIETDOS-BIO code Libraries and Data Bases: Nuclides Library - Electron, Positron, Beta, Alpha and Gamma Energy's data (ICRP 38) Reference Organisms Library: Terrestrial and Freshwater Ecosystems Concentration ratios (CR) – 3 Data Bases Partitioning coefficients (Kd) Library and etc.
LIETDOS-BIO (3 CR Data Bases) • • • Site-specific radionuclide CR distribution-geometric mean values : Data Bases No. 1. Generic radionuclide CR values (FASSET, ERICA, BWG Ex 2 CR values): Data Bases No. 2. Stable nuclide (heavy metals) CR values (ERICA, sitespecific). As demonstrate presented picture after the Chernobyl accident stable CR value may be used after few years: Data Bases No. 3.
FRESHWATER ECOSYSTEM 90 Sr CR for different type of Macrophytes – site-specific values Based on these data statistical CR approach has been used
The distribution of site specific Sr-90 CR values (macrophytes) involved in DB Nr 1
Co-60 CR (FW) site-specific distribution (macrophytes): LIETDOS-BIO site-specific DB (example) Statsistics: Value Trials 20000 Mean 8. 61 Geometric Mean 3. 87 Median Mode Standard Deviation Variance Skewness Kurtosis Coeff. of Variability 3. 565 --17. 114 292. 875 7. 81 115. 80 1. 99 Range Minimum 0. 008 Range Maximum 544. 488 Range Width 544. 480 Mean Std. Error 0. 121
TERRESTRIAL ECOSYSTEM Plume of the Chernobyl accident formation by meteorological conditions for instantaneous releases on dates and times (GMT) indicated
Aerial Gamma - Spectrometry Test Results. Spots Denote Areas Where the Density of Cs-137 is at least 3700 Bq/square meter (the first contour), with a Contour-spacing of 1100 Bq/m 2
Soil Cs-137 activity distribution Bq/kg, 1986 -1996 Lithuania Statistic Value Trials 10 000 Mean 12. 4 Median Mode o 12. 28 --- Standard Deviation 1. 57 Variance 2. 47 Skewness 0. 44 Kurtosis 3. 4 Coeff. of Variability 0. 13 Range Minimum 7. 28 Range Maximum 19. 83 Range Width 12. 56 Mean Std. Error 0. 02
The distribution of site specific Cs-137 CR values involved in DB Nr 1 (Wild boar) Trials 10 000 Mean 2. 41 Geom. Mean 2. 37 Median 2. 36 Mode --- Stand. Deviation 0. 46 Variance 0. 21 Skewness 0. 61 Kurtosis 3. 72 Range Minimum 1. 26 Range Maximum 4. 98
Internal Dose rate calculation • Biota activity concentration and internal dose rate was calculated as follows:
Crystal Ball statistical technique was used with 10 000 number of trials and the Latin Hypercube sampling method. An example: Forecast: Co-60 Internal dose rate (macrophyte) Trials 20000 Mean 2. 83 E-04 Geometric Mean 2. 07 E-04 Median 2. 07 E-04 Mode --- Standard Deviation 2. 63 E-04 Variance 6. 90 E-08 Skewness Kurtosis Coeff. of Variability 3. 47 28. 18 0. 93 Range Minimum 9. 15 E-06 Range Maximum 5. 61 E-03 Range Width 5. 60 E-03
External dose rate calculation • External dose rate from sediment (FW) : • External dose rate from sediment (DW):
An example of model for MCNPX code DCC calculation fish (ellipsoid on sediment) in the middle of water cylinder R = 100 cm Hw = 100 cm Water Hc=(10 +2 l); Dc= (5 +2 l) cm Hf = 10; Df = 5 cm Wc = (60 + 2 l) cm Fish 5 cm 10 cm Wf = 60 cm Sediment
An example of model for MCNPX code external exposure DCC calculation: fish (ellipsoid) in the middle of water cylinder R = 100 cm Hc=(12 +2 l); Dc= (4 +2 l) cm 50 cm Hf = 12; Df = 4 cm Wc = (20 + 2 l) cm Hw = 100 cm Fish Wf = 20 cm Water 10 cm Sediment
An example of model for MCNPX code DCC calculation rooted submerged hydrophytes = 1 cm R = 100 cm Hw = 10 0 cm Plant Water 10 cm Sediment
An example of model for MCNPX code DCC calculation unrooted submerged hydrophytes h 0 = 0. 1 cm h = 0. 2 cm R = 100 cm = 10 cm Hw = 10 0 cm Plant Water 10 cm Sediment
Crystal Ball statistical technique was used with 10 000 number of trials and the Latin Hypercube sampling method. An example: Forecast: Co-60 External dose rate (macrophyte) Trials 20000 Mean 1. 38 E-02 Geometric Mean 9. 72 E-03 Median 9. 73 E-03 Mode --- Standard Deviation 1. 39 E-02 Variance 1. 93 E-04 Skewness Kurtosis Coeff. of Variability 4. 23 42. 67 1. 01 Range Minimum 3. 57 E-04 Range Maximum 3. 39 E-01 Range Width 3. 38 E-01 Mean Std. Error 9. 83 E-05
LIETDOS-W subroutine • Taking into account tritium activity concentrations in the Ignalina Lake water and groundwater as well as groundwater activity in the vicinity of Maisiogala repository the LIETDOS -W code was developed to predict freshwater BIOTA exposure. • LIETDOS-W code was tested during Perch Lake scenario (EMRAS working group on Modeling of tritium transfer to biota and man).
LIETDOS-W subroutine Steady-state tritium transfer modelling Freshwater HTO Biota OBT Biota HTO Sediment Layer HTO LIETDOS-BIO
LIETDOS-W: HTO and OBT concentrations in clams (EMRAS TWG Pearch Lake scenario: L – LIETDOS-W data) HTO concentrations in clams. The model predictions are shown as solid diamonds with the vertical lines representing 95% confidence intervals as estimated by the modelers. The solid horizontal line is the observation with the 95% confidence interval indicated by the dashed lines. OBT concentrations in clams. The model predictions are shown as solid diamonds with the vertical lines representing 95% confidence intervals as estimated by the modelers. The solid horizontal line is the observation with the 95% confidence interval indicated by the dashed lines.
The main LIETDOS-BIO features • As the organism for which DCC calculated is defined by user, there is possibility for optional using references as well as other organisms or criteria. • The application of stochastic and probabilistic approach is one of possibility to improve the correctness of biota exposure predictions. • Radiation Protection and radiological dose assessment tools are well established for humans through the ICRP and national recommendations. Up to now Environment protection from ionizing radiation in Lithuania is accepted as much as it is interconnected with human exposure. THANK YOU FOR ATTENTION!
- Slides: 33