Europium sorption to sediments in the presence of natural organic matter: A laboratory and modeling study

Cellulosic materials, such as wood, paper products and cardboard that have been co-disposed with low-level nuclear waste have been shown to produce leachate with natural organic matter (NOM) concentrations of hundreds of mg/L C and, as such, have the potential to influence the fate and transport of radionuclides in the subsurface environment. The objective of this study was to examine the influence of NOM on the sorption of Eu (an analogue for trivalent radionuclides) to two coastal plain sediments from the US Department of Energy’s Savannah River Site. Particular attention was directed at quantifying Eu interactions with NOM sorbed to sediments (NOM_sed) in laboratory experiments and developing conditional stability constants for that interaction using the thermodynamic equilibrium speciation model MINTEQA2. Europium sorption to the two sediments systematically increased as pH increased from 3.9 to 6.7. With increasing additions of NOM to the aqueous phase from 0 to 222 mg/L C, Eu sorption initially increased to a maximum at 10 mg/L C NOM_aq and then decreased with increasing NOM_aq concentrations. Increases in Eu sorption at low NOM additions was attributed to the sorption of NOM to the sediment surface increasing the number of sorption sites on the low cation-exchange capacity sediments and/or increasing the association constant (log K) for the Eu-sediment surface reaction. Decreases in Eu sorption at higher NOM levels was attributed to Eu_aq complexation to NOM_aq being more favored than Eu sorption to the solid phase. A component additivity model was developed to describe the Eu–NOM-sediment system by the additive effects of the three binary system models: Eu–NOM, Eu-sediment and NOM-sediment. The model generally captured the data trends in the ternary system. Conditional stability constants developed from the experimental data for the complexation of Eu to NOM_sed were as much as four orders of magnitude greater than Eu complexation with NOM_aq, presumably due to the NOM_sed deriving additional negative (attractive) charge from the sediment surface. At high initial NOM_aq levels, >99 mg/L C, the model captured the trend of reduced Eu sorption but tended to over-estimate Eu sorption. The additivity approach of combining binary models to form a ternary model was only successful when the unique complexation properties of the NOM_sed were properly calculated.