Wavelet analysis is used to identify and describe variability in annual Canadian streamflows and to gain insights into the dynamical link between the streamflows and the dominant modes of climate variability in the Northern Hemisphere. Results from applying continuous wavelet transform to mean annual streamflows from 79 rivers selected from the Canadian Reference Hydrometric Basin Network (RHBN) reveal striking climate-related features before the 1950s and after the 1970s. The span of available observations, 1911--1999, allows for depicting variance for periods up to 12 years. Scale-averaged wavelet power spectra are used to simultaneously assess the interannual and spatial variability in 79 annual streamflow time series. The most striking feature, in the 2--3 year band and in the 3--6 year band (the 6--12 year band is dominated by white noise (since 1950) and is not considered further) is a net distinction between the timing and intensity of the interannual variability in western, central, and eastern streamflows, which is shown to be linked to the regional climatology. It is found that for the 2--3 year band, the Canadian streamflows are depicted mainly by the Pacific North America (PNA) during 1950--1999, and the Northern Hemisphere Annular Mode (NAM) only prior to 1950, and the North Atlantic Oscillation (NAO) after 1970. Similarly, in the 3--6 year band, the streamflows are depicted mostly by the NAO, the sea surface temperature anomalies over the Ni¿o-3 region (ENSO3) and the PNA during the period 1950--1999, and the NAM prior to 1950. Furthermore, strong local correlations between teleconnection patterns and western, central, and eastern streamflows are also revealed in both the 2--3 and 3--6 year bands with striking changes around 1950 and 1970. The correlation analysis in the 2--3 year and 3--6 year bands revealed the presence of two change points in the west and east streamflows occurring around 1950 and 1970.