Coupling river concentration simulations with a toxicokinetic model effectively predicts the internal concentrations of wastewater-derived micropollutants in field gammarids.
Although the exposure assessment of wastewater-derived micropollutants via chemical, bioanalytical, and modeling methods in environmental compartments are becoming more frequent, the whole-body burden (i.e. internal concentrations) in non-target organisms resulting from micropollutant exposure is rarely assessed. An understanding of the internal concentration fluctuation is especially important when exploring the mechanistic linkage between exposure and effects. In this study, we coupled a simple river model with a first-order toxicokinetic (TK) model to predict the concentrations of wastewater-derived micropollutants in freshwater invertebrates (Gammarus spp.). We applied Monte Carlo simulations and conducted laboratory experiments to account for the uncertain input data and the lack of uptake/depuration rate constants required for the TK model. The internal concentrations in field gammarids were predicted well and the estimates varied only by a factor of 0.1-1.9. Fast equilibrium may also be assumed such that bioconcentration factors (BCFs) is used together with the daily river dilution patterns to predict internal concentrations. While this assumption is suitable for compounds observed in our experiment to reach steady-state within 48 h in gammarids, it overpredicted the concentrations of substances that become in steady-state after longer periods. Nevertheless, this approach provides conservative estimates, and simplifies the coupling of models as BCFs are slightly more accessible than the rate constants. However, if one is interested in a more detailed exposure information (e.g., peak concentration and the whole-body burden recovery after a spill ), then the non-steady state formulation should be employed.