Speciation of rare earth elements in natural terrestrial waters: assessing the role of dissolved organic matter from the modeling approach

Abstract

Humic Ion-Binding Model V, which focuses on metal complexation with humic and fulvic acids, was modified to assess the role of dissolved natural organic matter in the speciation of rare earth elements (REEs) in natural terrestrial waters. Intrinsic equilibrium constants for cation-proton exchange with humic substances (i.e., p K MHA for type A sites, consisting mainly of carboxylic acids), required by the model for each REE, were initially estimated using linear free-energy relationships between the first hydrolysis constants and stability constants for REE metal complexation with lactic and acetic acid. p K MHA values were further refined by comparison of calculated Model V “fits” to published data sets describing complexation of Eu, Tb, and Dy with humic substances. A subroutine that allows for the simultaneous evaluation of REE complexation with inorganic ligands (e.g., Cl −, F −, OH −, SO 4 2−, CO 3 2−, PO 4 3−), incorporating recently determined stability constants for REE complexes with these ligands, was also linked to Model V. Humic Ion-Binding Model V’s ability to predict REE speciation with natural organic matter in natural waters was evaluated by comparing model results to “speciation” data determined previously with ultrafiltration techniques (i.e., organic acid–rich waters of the Nsimi-Zoetele catchment, Cameroon; dilute, circumneutral-pH waters of the Tamagawa River, Japan, and the Kalix River, northern Sweden). The model predictions compare well with the ultrafiltration studies, especially for the heavy REEs in circumneutral-pH river waters. Subsequent application of the model to world average river water predicts that organic matter complexes are the dominant form of dissolved REEs in bulk river waters draining the continents. Holding major solute, minor solute, and REE concentrations of world average river water constant while varying pH, the model suggests that organic matter complexes would dominate La, Eu, and Lu speciation within the pH ranges of 5.4 to 7.9, 4.8 to 7.3, and 4.9 to 6.9, respectively. For acidic waters, the model predicts that the free metal ion (Ln 3+) and sulfate complexes (LnSO 4 +) dominate, whereas in alkaline waters, carbonate complexes (LnCO 3 + + Ln[CO 3] 2 −) are predicted to out-compete humic substances for dissolved REEs. Application of the modified Model V to a “model” groundwater suggests that natural organic matter complexes of REEs are insignificant. However, groundwaters with higher dissolved organic carbon concentrations than the “model” groundwater (i.e., >0.7 mg/L) would exhibit greater fractions of each REE complexed with organic matter. Sensitively analysis indicates that increasing ionic strength can weaken humate-REE interactions, and increasing the concentration of competitive cations such as Fe(III) and Al can lead to a decrease in the amount of REEs bound to dissolved organic matter.