Towards a better understanding of trace element transport in permafrost-affected Earth surface environments, we sampled ∼ 60 large and small rivers (< 100 to ≤ 150 000 km2 watershed area) of Western Siberia Lowland (WSL) during spring flood and summer and winter base-flow across a 1500 km latitudinal gradient covering continuous, discontinuous, sporadic and permafrost-free zones. Analysis of ∼ 40 major and trace elements in dissolved (< 0.45 μm) fraction allowed establishing main environmental factors controlling the transport of metals and trace elements in rivers of this environmentally important region. No statistically significant effect of the basin size on most TE concentration was evidenced. Three category of trace elements were distinguished according to their concentration – latitude pattern: (i) increasing northward in spring and winter (Fe, Al, Ga (only winter), Ti (only winter), REEs, Pb, Zr, Hf, Th (only winter)), linked to leaching from peat and/or redox processes and transport in the form of Fe-rich colloids, (ii) decreasing northward during all seasons (Sr, Mo, U, As, Sb) marking the underground water influence of river feeding and (iii) elements without distinct trend from S to N whose variations within each latitude range were higher than the difference between latitudinal ranges (B, Li, Ti (except summer), Cr, V, Mn, Zn, Cd, Cs, Hf, Th). In addition to these general features, specific, northward increase during spring period was mostly pronounced for Fe, Mn, Co, Zn and Ba and may stem from a combination of enhanced leaching from the topsoil and vegetation and bottom waters of the lakes (spring overturn). A spring time northward decrease was observed for Ni, Cu, Zr, Rb. The southward increase in summer was strongly visible for Fe, Ni, Ba, Rb and V, probably due to peat/moss release (Ni, Ba, Rb) or groundwater feeding (Fe, V). The Principal Component Analysis demonstrated two main factors potentially controlling the ensemble of TE concentration variation. The first factor, responsible for 16–20 % of overall variation, included trivalent and tetravalent hydrolysates, Cr, V, and DOC and presumably reflected the presence of organo-mineral colloids, as also confirmed by previous studies in Siberian rivers. The 2nd factor (8–14 % variation) was linked to the latitude of the watershed and acted on elements affected by the groundwater feeding (DIC, Sr, Mo, As, Sb, U), whose concentration decreased significantly northward during all seasons. Overall, the rank of environmental factors on TE concentration in western Siberian rivers was latitude (3 permafrost zones) > season > watershed size. The effect of the latitude was minimal in spring for most TE but highly visible for Sr, Mo, Sb and U. The main factors controlling the shift of river feeding from surface and subsurface flow to deep underground flow in the permafrost-bearing zone were the depth of the active (unfrozen) seasonal layer and its position in organic or mineral horizons of the soil profile. In the permafrost-free zone, the relative role of carbonate mineral-bearing base rock feeding vs. bog water feeding determined the pattern of trace element concentration and fluxes in rivers of various size as a function of season. Comparison of obtained TE fluxes in WSL rivers with those of other subarctic rivers demonstrated reasonable agreement for most trace elements; the lithology of base rocks was the major factor controlling the magnitude of TE fluxes. The climate change in western Siberia and permafrost boundary migration will affect essentially the elements controlled by underground water feeding (DIC, alkaline-earth elements (Ca, Sr), oxyanions (Mo, Sb, As) and U). The thickening of the active layer may increase the export of trivalent and tetravalent hydrolysates in the form of organo-ferric colloids. Plant litter-originated divalent metals present as organic complexes may be retained via adsorption on mineral horizon. However, due to various counterbalanced processes controlling element source and sinks in plants – peat – mineral soil – river systems, the overall impact of the permafrost thaw on TE export from the land to the ocean may be smaller than that foreseen by merely active layer thickening and permafrost boundary shift.