Tm Water Chemistry and Characterization
In the trial to characterize Tm concentrations in RM only (H. azteca not exposed), the Tm-Rec concentrations were generally less than the planned nominal concentrations (Table 1). As nominal concentrations increase the proportion of recovered Tm was reduced (except at the highest concentration, Table 1). Measured Tm-T and Tm-D concentrations generally showed increasing amounts of precipitation as the nominal and Tm-Rec concentrations increased. Immediately after solutions were created (0 h time) there were differences between Tm-T and Tm-D concentrations but by 24 h the solution Tm concentrations appeared to stabilize and there were no differences between Tm-T and Tm-D concentrations over the next 96 h (Table 1). Following equilibration there was a clear pattern of increased precipitation of Tm at higher concentrations with an apparent solubility threshold of approximate 3-4 µM Tm-Rec. At the highest of the prepared concentrations 84% (24 h) and 87% (120 h) of the Tm that had been added to the beaker (as measured by Tm-Rec) was in an insoluble form.
Overall, the bench test results illustrated that there was likely precipitation occurring in the stock solution used to prepare exposure solutions (nominal vs Tm-Rec concentration differences in Table 1). This was anticipated because the stock solution was adjusted to pH 7.3. From experience, the use of an acidified stock solution results in recovered concentrations being much closer to nominal however it also produces significant differences in solution pH across the concentration range and the process of solution-by-solution pH adjustment can be exceptionally time-consuming. By using a stock solution with a pH that matched that of the RM we reduced the need to adjust the pH of each test solution but the trade-off was variation in Tm concentrations. Of key importance is measurement of Tm-T and Tm-D at the beginning and the end of the test. Data from the bench test indicated that 24 h was sufficient for test solution equilibration as there were few differences between 120 h and 24 h concentrations (Table 1). There were indications of significant precipitation and at the highest concentrations with only 13-17% of the Tm was in the dissolved form (Table 1).
Precipitation is recognized as a complicating factor in REE toxicity tests (Gonzalez et al. 2014). Borgmann et al. (2005) reported a nominal LC50 of 721 ug Tm/L with a measured LC50 of 0.01 ug Tm/L. In aquatic toxicity tests La readily precipitated out of solution and measured values were always less than 30% of the nominal concentrations (Barry and Meehan 2000). Ce, Gd and Nd precipitated from solution as concentration increased and also over time (up to 72 h, Blinova et al. 2019). Vukov et al. (2016) found that precipitation of Dy correlated with increased exposure concentration. It was reported that at high exposure concentrations of Dy, dissolved concentrations were less than 34% of total concentrations (Vukov et al. 2016). Precipitation can likely be accounted for by the pH and the carbonate content of RM and the formation of insoluble salts (Jiang and Ji 2012; Janssen and Verweij 2003). The study of Gonzalez et al. (2015) also reported the formation of “insoluble species” in REE test solutions and used geochemical modelling to show that hydroxide and carbonate precipitation was expected.
Tm acute toxicity in RM
H. azteca mortality increased with increasing concentrations of Tm (Table 2). Based on the results from the bench test, we expected measured Tm concentrations at test beginning and at test end to be relatively similar however, this was not always the case (Table 2). An equilibration time of 24 h may not have been sufficient and it is unknown if the addition of organisms influenced the geochemistry of test solutions. Calculations of the standard acute toxicity endpoint (96 h LC50) were done with measured concentrations and this was possible for both Tm-T or Tm-D and either at the beginning of the test or at the end (Table 2). In general, the measured concentrations were lowest at the end of the tests and therefore, as a conservative approach, calculations were based on samples collected at 96 h.
In unaltered RM the 96 h LC50 value for measured dissolved was 573 µg/L (3.4 µM) with 95% confidence interval from 482 to 663 µg/L (2.9-3.9 µM; Fig 1, light gray bar at 0.5 mM Ca). This concentration is much higher than the 7 d LC50 value of 0.01 µg/L previously reported (Borgmann et al. 2005) for measured Tm-D in soft water. Compared to our study, Borgmann et al. (2005) conducted static tests for a longer duration (4 vs 7 d) in softer water (hardness value 60 vs 12 mg CaCO3/L) and with feeding part way through (our tests were without food). While exposure duration, geochemistry and the provision of food are recognized as potential influences on the bioavailability and toxicity of metals, it is unknown if these can explain the dramatically different results (57,000 fold difference in LC50 values). Interestingly, the nominal 7 d LC50 values reported by Borgmann et al. (2005) were 721 µg/L in soft water and 739 µg/L in hard water and this compares well to our estimate of 1062 µg/L for the 96 h LC50 on a nominal basis (using the Tm-D to nominal Tm relationship in the test to estimate). The very low Tm-D value previously reported by Borgmann et al. (2005) remains unexplained and it is worth noting that the authors of the work do not offer any discussion of it what-so-ever. Their analysis and discussion of Tm toxicity used the nominal based toxicity endpoints exclusively.
Influence of Cations on Tm toxicity
The results for effects of Ca on Tm toxicity (Fig. 1) were somewhat difficult to interpret. On the basis of Tm-T concentrations the lowest LC50 was at the 0.5 mM Ca treatment and toxicity was significantly reduced at both lower and higher concentrations (0.25 and 1.5 respectively). However, for Tm-D there was no trend evident for the low Ca exposure because the LC50 values at 0.25 and 0.5 Ca were not significantly different. At the higher Ca treatment Tm-D showed higher toxicity (Fig 1). There was no protective effect with increasing Na concentrations (Fig. 2). Similarly, we did not see a protective effect with Mg as there were no significant differences in LC50 values across the range of Mg tested (Fig. 3).
We hypothesized that increases in Ca would have a protective effect to Tm toxicity but there was no consistent trend across the range of Ca tested (Fig 1). Other studies have shown that Ca provides significant protection against the toxicity of REEs. For example, Vukov et al. (2016) showed a 1.8-fold decrease in Dy toxicity to H. azteca over a 3-fold increase in Ca concentration. Barry and Meehan (2000) showed that La toxicity to Daphnia carinata was reduced as hardness increased. Cardon et al. 2019 and Ma et al. 2016, with Y and Ce (respectively) similarly demonstrated reduced toxicity to D. magna with increased hardness. In these studies the changes in toxicity are linked to hardness and so cannot be exclusively attributed to Ca. As discussed above, changes in exposure hardness did not result in changes of Tm toxicity (Borgmann et al 2005) and this is consistent with our results but only nominal LC50 values are available. It may be that Tm uptake and toxicity is not influenced by competitive interaction with Ca and therefore it is unlike other REEs (e.g. Dy, La, Ce and Y) that are. There was no protective effect with increasing Na concentrations (Fig. 2) and this result is different to the Vukov et al., (2016) study where a 3 fold increase of Na significantly decreased Dy toxicity by a factor of 1.4 times. However, these results were based on total Dy concentrations and the study states that LC50s for dissolved Dy concentrations were much less clear (Vukov et al. 2016). Out results with Mg align with those of Vukov et al. (2016) where Mg additions did not show a protective effect on Dy toxicity.
We had hypothesized that Tm toxicity would be influenced by cationic competition, particularly Ca2+. Previous studies on the toxicity of inorganic forms of metals attribute the toxicity reduction achieved by cations to direct competition at the site of uptake of essential ions such as Ca2+ and Na+ (Niyogi and Wood 2004). Toxic metals have similar characteristics (e.g. ionic radius and charge) compared to essential ions and via ionic mimicry (Bridges and Zalups 2005) divalent metal cations such as Cd2+, Zn2+ and Pb2+ inhibit Ca2+ uptake while monovalent metals such as Ag+ disrupting Na+ uptake (Niyogi and Wood 2004). From this perspective the mechanism by which a trivalent REE free ion would interact (compete) with an essential divalent (Ca or Mg) or monovalent (Na) cation is not clear. However, studies have highlighted the similar properties of trivalent REEs, particularly in relation to Ca2+ (Evans 1983) and it is well known that La3+ is an effective analogue. There is some evidence of direct competition between Ca2+ and trivalent REEs in algae. Ca2+ competitively inhibit La3+ and Ce3+ uptake and protect against toxicity in Chlorella fusca (Aharchaou et al. 2020). A 1000 fold increase in Ca concentration (1 µM - 1 mM) resulted in a 2 fold (La exposure) and 3 fold (Ce exposure) increase in cell density over 120 h of exposure (Aharchaou et al. 2020). Tests with C reinhardtii demonstrated competition with reduced Nd3+ (Yang and Wilkinson 2018) and Sm3+ (Tan et al. 2017) uptake as either Ca2+ or Mg2+ increased. In our experiments neither Ca, Mg nor Na influenced Tm toxicity to Hyalella. Understanding the mechanisms of uptake of Tm in aquatic invertebrates would also be valuable in understanding the potential for toxicity mitigation in natural waters.
Influence of DOM on Tm toxicity
In solutions with added Luther Marsh DOM toxicity was significantly reduced above concentrations of 3 mg DOC/L (Fig 4). The addition of DOM also altered the relative concentrations of Tm-T and Tm-D and it appeared that the precipitation threshold may have been increased with elevated DOC content (Table 2). DOM has been shown to reduce the toxicity of numerous metals in a concentration dependent manner (Wood et al. 2011). It is a complex heterogeneous molecule with a variety of negatively charged moieties that are capable of interacting with cationic metals. Complexation of the free ion form of the metal reduces the availability for uptake thereby reducing toxicity. This was evident in our study for Tm (Fig 4) and we assumed that mitigation of toxicity was due to reduced Tm3+ concentrations. However, this is only an assumption as the bioavailable forms of Tm associated with toxicity and the mechanism of uptake at the biotic surface are unknown.
While we did not measure free ion concentrations of Tm3+ in solution, we did use the geochemical equilibrium modeling software WHAM (Windermere Humic Aqueous Model, Ver 7.02; Tipping et al. 2011) to estimate Tm3+. Using the Tm-D and measured DOC concentrations at the end of the test as model inputs as described by Stockdale et al (2015), WHAM predicted virtually complete complexation of Tm. The predicted Tm3+ concentrations in test solutions with added DOM were at least 140 fold lower (highest Tm-D with lowest DOC) and ranged up to 2.5x106 fold lower (lowest Tm-D with highest DOC) than the corresponding Tm3+ concentrations in solutions with no added DOM. Clearly the predicted Tm3+ estimates were not linked to the acute toxicity of Tm. One possible conclusion is that Tm3+ is not associated with toxicity in Hyalella and that other (or additional) geochemical forms are. It is also possible that WHAM is predicting a much higher level of complexation of Tm3+ than is actually occurring in our test solutions. Either way, DOM significantly reduces Tm toxicity and an improved understanding of the geochemical speciation of Tm in relation to acute toxicity is required.
There are relatively few studies on the effects of REEs on aquatic biota and even fewer investigating the potential influence of DOM on toxicity. DOM has been highlighted as an important factor to include in water quality derivations for La (Hermann et al 2016). Vukov et al (2016) used Suwannee River DOM to show 3-4 fold reductions of Dy toxicity to Hyalella at a DOC concentration of 13 mg/L. The biouptake of Sm3+ (as measured directly by ion exchange technique) into the unicellular green algae Chlamydomonas reinhardtii was significantly reduced by DOM in a DOC concentration dependent manner (Rowell et al 2018). In that study four different sources of DOM were tested, including Luther Marsh DOM, and even very small additions of 0.5 mg DOC/L dramatically reduced uptake by 10 fold (Rowell et al 2018). In tests with the synthetic organic ligands malic acid, diglycolic acid and citric acid, Sm uptake to C. reinhardtii was reduced but the possibility of complexed Sm also being taken up could not be ruled out (Tan et al 2016). Similar reduced uptake results into Chlorella vulgarize were shown for La, Gd and Y using the organic ligands citrate, nitriloacetic acid and ethylenediamine tetraacetic acid (Sun et al 1997). In natural waters the important role that DOM has in complexing REEs is well recognized (Moermond et al 2001; Tang and Johannesson 2003, Tipping and Filella 2020) and it is generally assumed that complexation will reduce toxicity (Herrmann et al 2016). Given the ubiquitous nature of DOM in natural waters and the significant reduction in toxicity we observed, further study on DOM-Tm interactions would contribute to the development of water quality thresholds for assessing the environmental risk of this REE.