This study aims at evaluating the hydrological balance of large watersheds of the Canadian Shield in the James Bay area in Northwestern Quebec, Canada. The focus is set on six rivers of the Canadian Shield altogether draining more than 185,000 km² of the Boreal Shield, Taiga Shield and Hudson Plains ecozones of Canada. River discharge measurements, geochemical data (δ²H, δ¹⁸O and electrical conductivity [EC] of water), remote sensing, and GIS models are used jointly to calculate water balances. The approach allows for partitioning the influence of rainwater, snowmelt, surface runoff, evaporation, transpiration, and groundwater discharge to the hydrological balances of watersheds. On an annual basis, the results suggest that runoff from rainwater (30–61 % of total precipitation) and snowmelt (18–40 % of total precipitation) are the main contributions to river discharge, while the contribution of groundwater discharge to rivers represents < 12 % of the total precipitation. Over the study area, this contribution represents 2–5 km³ of water. The stable isotope mass balances allow for estimating watershed-scale evaporation over inflow ratios ranging between 2 and 10 % and suggest that transpiration has an isotopic composition close to summer rainwater. The hydrological balances further suggest that the total pool of water stored in the active portion of watersheds represents 10–20 % of the total annual precipitation, while the exports of groundwater beyond the limits of surface watersheds are negligible. The seasonal trends in the hydrological balances of monitored watersheds were further documented to provide insights into the sensitivity of watersheds as they face climate change. The observations are used to propose recommendations for monitoring of rivers in the Canadian Shield and to identify future research needs.