Profound understanding of behaviors of organic matter from sources to multistage rivers assists watershed management for improving water quality of river networks in rural areas. Ninety-one water samples were collected from the three orders of receiving rivers in a typical combined polluted subcatchment (diffuse agricultural pollutants and domestic sewage) located in China. Then, the fluorescent dissolved organic matter (FDOM) information for these samples was determined by the excitation-emission matrix coupled with parallel factor analysis (EEM-PARAFAC). Consequently, two typical humic-like (C1 and C2) and other two protein-like (C3 and C4) components were separated. Their fluorescence peaks were located at λ ex/em = 255(360)/455, <250(320)/395, 275/335, and <250/305 nm, which resembled the traditional peaks of A + C, A + M, T, and B, respectively. In addition, C1 and C2 accounted for the dominant contributions to FDOM (>60 %). Principal component analysis (PCA) further demonstrated that, except for the autochthonous produced C4, the allochthonous components (C1 and C2) had the same terrestrial origins, but C3 might possess the separate anthropogenic and biological sources. Moreover, the spatial heterogeneity of contamination levels was noticeable in multistage rivers, and the allochthonous FDOM was gradually homogenized along the migration directions. Interestingly, the average content of the first three PARAFAC components in secondary tributaries and source pollutants had significantly higher levels than that in subsequent receiving rivers, thus suggesting that the supervision and remediation for secondary tributaries would play a prominent role in watershed management works.