Abstract

The occurrence of toxic metal(loid)s in surface freshwater is a global concern due to its impacts on human and ecosystem health. Conceptual and quantitative metal(loid) models are needed to assess the impact of metal(loid)s in watersheds affected by acid rock drainage. Few case studies have focused on arid and semiarid headwaters, with scarce hydrological and hydrochemical information. This work reports the use of WASP8 (US EPA) to model Al, Fe, As, Cu, and SO42? concentrations in the Upper Elqui River watershed in north–central Chile. Calibrated model performance for total concentrations was “good” (25.9, RRMSE; 0.7, R2-d) to “very good” (0.8–0.9, R2-d). The dissolved concentrations ranged between “acceptable” (56.3, RRMSE), “good” (28.6, RRMSE; 0.7 d), and “very good” (0.9, R2-d). While the model validation achieved mainly “very good” (0.8–0.9, R2-d) predictions for total concentrations, the predicted dissolved concentrations were less accurate for all indicators. Sensitivity analysis showed that the partition coefficient is a sensitive constant for estimating dissolved concentrations, and that integrating sorption and sediment interaction reduces the model error. This work highlights the need for detailed and site-specific information on the reactive and hydrodynamic properties of suspended solids, which directly impact the partition coefficient, sedimentation, and resuspension velocity calibration. © 2025 by the authors.