Watershed features shape spatial patterns of fish tissue mercury in a boreal river network

Freshwater mercury (Hg) contamination is a widespread environmental concern but how proximate sources and downstream transport shape Hg spatial patterns in riverine food webs is poorly understood. We measured total Hg (THg) in slimy sculpin (Cottus cognatus) across the Kuskokwim River, a large boreal river in western Alaska and home to subsistence fishing communities which rely on fish for primary nutrition. We used spatial stream network models (SSNMs) to quantify watershed and instream conditions influencing sculpin THg. Spatial covariates for local watershed geology and slope accounted for 55 % of observed variation in sculpin THg and evidence for downstream transport of Hg in sculpins was weak. Empirical semivariograms indicated these spatial covariates accounted for most spatial autocorrelation in observed THg. Watershed geology and slope explained up to 70 % of sculpin THg variation when SSNMs accounted for instream spatial dependence. Our results provide network-wide predictions for fish tissue THg based largely on publicly available geospatial data and open-source software for SSNMs, and demonstrate how these emerging models can be used to understand contaminant behavior in spatially complex aquatic ecosystems.