Estimation of apparent rate coefficients for radionuclides interacting with marine sediments from Novaya Zemlya

Kinetic reactive transport explains distinct subsurface deposition patterns of pollutants in sediments. The case of the Sellafield-derived 236U, 137Cs and 239,240Pu in the Esk Estuary, UK

The kinetics of the uptake of pollutants by solids in sediments interacts with transitional eddy diffusivity in the pore fluid, leading to different depth-distribution patterns. This work aims to gain insights into the still poorly understood behaviour in the marine environment of the anthropogenic ²³⁶U, a recently postulated tracer of water masses. It is hypothesized that the transition from mobile U(VI) to highly particle-reactive U(IV) in the anoxic zone of the sediment produces a subsurface deposition pattern. A novel numerical model for kinetic reactive transport in sediments, which merges diagenetic processes for transport and box models for the uptake, is used for concept demonstration. It is applied to synthetic environments with high eddy diffusivity to obtain the singular depth-distribution patterns of pulsed inputs of tracers that mimic the anthropogenic ^239,240Pu, ¹³⁷Cs, and ²³⁶U. While the first is retained in the upper cm, the second shows an exponential penetration pattern over few cm, and ²³⁶U is deposited with a Gaussian-like pattern centred below the transition to the anoxic zone. These patterns are then merged into a diagenetic model to compute the depth distribution at decadal or centennial scales of dissolved and particle-bound inputs of these radiotracers. It is successfully applied to a real case using literature data for a sediment core from the Esk Estuary, UK, affected by radioactive releases from the Sellafield nuclear reprocessing plant. This work provides insight into until now poorly understood field data and provides a novel view of broad implications in the study of the behaviour of pollutants in surficial aquatic sediments.

Solution to the particle concentration effect on determining Kd value of radionuclides

The particle concentration effect on Kd values of radionuclides has been observed but the underlying mechanism remains controversial. The hope is to use the relationship between particle concentration, adsorption-desorption isotherms and reversibility, in combination with surface component activity of model (SCA model), to solve this issue. 137Cs, 60Co, 90Sr were used as tracers, batch experiments were conducted in freshwater-sediment and seawater-sediment. The experiment of each radionuclide was designed with five different particle concentrations Cp, and for each Cp there were seven different initial concentrations C0. After adsorption experiments, four consecutive desorption experiments were carried out. At the fourth desorption experiment, radionuclide concentrations in the supernatant and sediment were measured. The results showed that adsorption and single desorption data of 137Cs, 60Co, 90Sr might be described by linear isotherms. 137Cs was reversible in the seawater-sediment, so hysteresis angles of the five-particle concentration were approximately 0°, all adsorption and desorption data could be classified into one line. In the remaining systems, besides the adsorption and single desorption isotherms moved upward with the decrease of particle concentration, hysteresis angles and irreversibility also increased, thus, the particle concentration effect was obvious. The reversible and resistant component concentrations calculated by adsorption, single desorption and consecutive desorption isotherm were linear functions of equilibrium concentration Ce1, respectively. Data from adsorption and desorption experiments with particle concentration effect could be classified into the same line using the Freundlich-SCA model. The results of this study indicate that the particle concentration effect is related to reversibility. When adsorption isotherm and single desorption isotherm are both linear, consecutive desorption isotherm, reversible and resistant component concentrations approach linearity too. After the Freundlich-SCA model eliminated the particle concentration effect on adsorption and desorption data, the data can be used to predict the adsorption, single desorption isotherm and Kd value at any particle concentration.

Modelling the kinetic reactive transport of pollutants at the sediment-water interface. Applications with atmospheric fallout radionuclides

Understanding the behaviour of particulate matter and chemicals at the sediment-water interface (SWI) is of interest in environmental studies and risk assessments. These processes are still poorly understood, and this work aims to gain relevant insights by using a kinetic reactive transport model. It merges early diagenetic processes and box models for the uptake kinetics. Numerical solutions have been found for synthetic scenarios and for studying real cases from the literature (²¹⁰Pb and Chernobyl fallout radionuclides in Lake Sniardwy, Poland, and ⁷Be in sediments from Tema Harbour, Ghana). The study identifies a series of factors that dynamically interact to govern the final fate of tracers in the SWI region, leading to a wide diversity of behaviours. When a term of eddy diffusivity is included in the upper regions of the pore fluid, which seems feasible for some energetic scenarios, it is possible to explain the observed large penetration depths for Cs and Be, while high particle-reactive elements are retained in thinner sediment layers. Desorption from the sediment occurs through the pore fluid as diffusive fluxes. Transient depth profiles of tracer concentrations can last from months up to a year, and they can show subsurface maxima at positions unrelated with the accretion rate. In the application cases, the model explained a wide set of observational data that was beyond the capabilities of other approaches involving physical mixing of solids and equilibrium k_d. This modelling study could provide useful guidance for future research works.

Coastal transport of river-discharged radionuclides: Impact of speciation and transformation processes in numerical model simulations

Following a potential nuclear accident, river run-off may potentially become a significant source of radionuclide contamination to the coastal marine environment. In the present work, code for radionuclide speciation and dynamic transfer of radionuclides between the different species was implemented in a Lagrangian marine dispersion model. A case study was performed where the model system utilized ocean circulation fields at relatively high spatial (160 m × 160 m in horizontal direction) and temporal resolution (1 hour), considering a hypothetical accident scenario including river discharges of ¹³⁷Cs to the marine environment. Results from a number of simulations were compared to identify how factors associated with radionuclide speciation and transfer between the model compartments could affect the predicted radiocesium activity concentrations. The results showed that by including dynamic transfer of radionuclides between the model compartments, the total activity concentrations at far-field sites could vary with more than two orders of magnitude, demonstrating that this model configuration enables prediction of potential local hot-spots. However, the total activity concentration near the river outlets was less affected (< factor 10). The radionuclide speciation in the river discharges and the parameterization of ¹³⁷Cs particle affinity greatly affected the specie distribution (> factor 10³ increase in concentration of particle-associated ¹³⁷Cs) as well as the settling of radionuclides towards the seabed (up to factor 10² increase in ¹³⁷Cs sediment concentrations). These factors were therefore identified as important contributors to the overall uncertainty.

Bioavailability of sediment-associated and low-molecular-mass species of radionuclides/trace metals to the mussel Mytilus edulis

Sediments can act as a sink for contaminants in effluents from industrial and nuclear installations or when released from dumped waste. However, contaminated sediments may also act as a potential source of radionuclides and trace metals to the water phase due to remobilisation of metals as dissolved species and resuspension of particles. The marine mussel Mytilus edulis is a filter-feeding organism that via the gills is subjected to contaminants in dissolved form and from contaminants associated to suspended particles via the digestive system. In this paper the bioavailability of sediment-associated and seawater diluted Cs, Co, Cd and Zn radioactive tracers to the filtering bivalve M. edulis has been examined. The mussels were exposed to tracers diluted in ultrafiltered (<10kDa) seawater (Low Molecular Mass form) or to tracers associated with sediment particles from the Stepovogo Fjord at Novaya Zemlya in short-term uptake experiments, followed by 1-month depuration experiments in flow-through tanks. A toxicokinetic model was fitted to the uptake and depuration data, and the obtained parameters were used to simulate the significance of the two uptake pathways at different suspended sediment loads and sediment-seawater distribution coefficients. The results of the model simulations, assuming steady state conditions, suggest that resuspended particles from contaminated sediments can be a highly significant pathway for mussels in the order (109)Cd everse congruent(65)Zn<(134)Cs<(60)Co. The significance increases with higher suspended sediment load and with higher K(d). Furthermore, the experimental depuration data suggest that Cs is retained longer and Co, Cd and Zn shorter by the mussels when associated with ingested sediments, than if the metals are taken up from the low molecular mass (LMM) phase.

Kinetic box models for the uptake of radionuclides and heavy metals by suspended particulate matter: equivalence between models and its implications

In recent years an increasing experimental effort has been paid to the study of the sorption process of radionuclides and heavy metals by particulate matter in aquatic environments. This has led to the development of different kinetic box models. Most of them are variations of two basic approaches: one containing several (up to three) parallel reactions while the other involves consecutive reactions. All the reactions are reversible (irreversibility is contained as a particular case) with concentration independent coefficients. The present work provides analytical solutions and demonstrates that both approaches are mathematically equivalent. That is, both models produce the same analytical solution for the uptake curve (time course of the concentrations in the dissolved phase), which is illustrated using literature data. This result unifies the description of the observed behaviour, but it brings up the question of the physical meaning of the involved coefficients. Finally, the mathematical relationship developed here serves to discuss some limitations found in recent attempts in literature devoted to distinguish the actual uptake mechanism.

Trends of radionuclide sorption by estuarine sediments. Experimental studies using 133Ba as a tracer

Sediments play an important (but still poorly understood) role in the dispersion and final fate of radionuclides and other hazardous materials in aquatic systems. Adopting an experimental point of view, the present work deals with the transfer of a radioactive tracer from water to sediments. Sediments and waters were sampled in the Odiel and Tinto estuaries (South-West Spain) with anthropogenic-enhanced 226Ra concentrations. 133Ba was used as a tracer since it is a gamma emitter and a good analogue of the environmental behaviour of 226Ra. Laboratory experiments have been carried out to quantify the uptake kinetics of 133Ba by sediments in aqueous suspensions and by sediment cores under a water column at rest. Depth distributions of 133Ba in sediments have been studied with different contact times and using sediment samples with different grain sizes. The results reveal an important and fast uptake by suspended sediments (up to 40% in 10 min with a 20 g l(-1) suspension) and sediment cores (up to 70% in a few minutes). The kinetics of the uptake by suspended sediments could be reasonably described by a model of two parallel and reversible reactions followed by a weakly-reversible reaction. The total uptake and the rate of reaction decreased with the increment of grain size. Furthermore, uptake following diffusion through the interstitial water was hardly detectable beyond 1 cm depth. For the case with more experimental results, the depth distribution could be reasonably described by the analytical solution based on the use of an effective diffusion coefficient (4.2 x 10(-12) m2 s(-1)) and the measured intrinsic distribution coefficient (kd = 63 +/- 6 l kg(-1)).

Experimental and modelling study on the uptake and desorption kinetics of 133Ba by suspended estuarine sediments from southern Spain

Dispersion of pollutants in aquatic environments depends on their uptake by suspended solids. This work deals with the uptake kinetics of 133Ba (gamma-emitter and a good analogue of 226Ra) by suspended estuarine sediments (which can be resuspended into the water column under certain conditions). This study presents a wide set of tracing experiments, including second tracing, decantation and desorption processes. The purpose is to characterize 133Ba uptake by sediments and to investigate the use and limitations of box models in order to describe the uptake kinetics. Water and sediment samples were collected in the Huelva estuary (Spain), where environmental 226Ra concentrations have been increased by two phosphate fertilizer industries. Samples were characterized by granulometric, organic carbon content, cation exchange capacity and XRF-EP analyses. Results revealed three-step kinetics, with characteristic times of minutes, hours and days. These results enabled the selection and calibration of a suitable box model and facilitated the testing of its use as a fully predictive tool.

Testing the behaviour of different kinetic models for uptake/release of radionuclides between water and sediments when implemented in a marine dispersion model

Three kinetic models for adsorption/release of (137)Cs between water and sediments have been tested when they are included in a previously validated dispersion model of the English Channel. Radionuclides are released to the Channel from La Hague nuclear fuel reprocessing plant (France). The kinetic models are a 1-step model consisting of a single reversible reaction, a 2-step model consisting of two consecutive reversible reactions and an irreversible model consisting of three parallel reactions: two reversible and one irreversible. The models have been tested under three typical situations that correspond to the source terms that can generally be found: instantaneous release, continuous release and redissolution of radionuclides from contaminated sediments. Differences between the models become more evident when contact times between water and sediments are larger (continuous release) and in the case of redissolution from sediments. Time scales for the redissolution process are rather different between the three models. The 1-step model produces a redissolution that is too fast when compared with experimental evidence. The irreversible model requires that saturation effects of the irreversible phase are included. Probably, the 2-step model represents the best compromise between ease and level of detail of the description of sorption/release processes.

Fixation of Cs to marine sediments estimated by a stochastic modelling approach

Dumping of nuclear waste in the Kara Sea represents a potential source of radioactive contamination to the Arctic Seas in the future. The mobility of 137Cs ions leached from the waste will depend on the interactions with sediment particles. Whether sediments will act as a continuous permanent sink for released 137Cs, or contaminated sediments will serve as a diffuse source of 137Cs in the future, depends on the interaction kinetics and binding mechanisms involved. The main purpose of this paper is to study the performance of different stochastic models using kinetic information to estimate the time needed for Cs ions to become irreversibly fixed within the sediments. The kinetic information was obtained from 134Cs tracer sorption and desorption (sequential extractions) experiments, conducted over time, using sediments from the Stepovogo Fjord waste dumping site, on the east coast of Novaya Zemlya. Results show that 134Cs ions interact rapidly with the surfaces of the Stepovogo sediment, with an estimated distribution coefficient Kd(eq) of 300 ml/g (or 13m2/g), and the 134Cs ions are increasingly irreversibly fixed to the sediment over time. For the first time, stochastic theory has been utilised for sediment-seawater systems to estimate the mean residence times (MRTs) of Cs ions in operationally defined sediment phases described by compartment models. In the present work, two different stochastic models (i) a Markov process model (MP) being analogous to deterministic compartment models, and (ii) a semi-Markov process model (SMP) which should be physically more relevant for inhomogeneous systems, have been compared. As similar results were obtained using the two models, the less complicated MP model was utilised to predict the time needed for an average Cs ion to become irreversibly fixed in the Stepovogo sediments. According the model, approximately 1100 days of contact time between Cs ions and sediments is needed before 50% of the 134Cs ion becomes fixed in the irreversible sediment phase. while about 12.5 years are needed before 99.7% of the Cs ions are fixed. Thus, according to the model estimates the contact time between 137Cs ions leached from dumped waste and the Stepovogo Fjord sediment should be about 3 years before the sediment will act as an efficient permanent sink. Until then a significant fraction of 137Cs should be considered mobile. The stochastic modelling approach provides useful tools when assessing sediment-seawater interactions over time, and should be easily applicable to all sediment-seawater systems including a sink term.