Fingerprint technique is effective in relative sediment contribution estimation. As recently developed methods, Bayesian model and analytical solutions for single factor received more and more attentions. However, compared with most often used multivariate mixing models, their performance and stability in calculating the proportion of sediment sources kept unknow. In order to make clear the contribution of topsoil and subsoil in a typical small watershed in the black soil region of Northeast China, where distributed vast farmlands experiencing serious soil erosion and a number of gullies developing very fast, these three methods were introduced to provide estimation. In this studied watershed with 27.60 km² area, totally 69 samples from the sediment sources area and 30 samples from the sedimentation area were collected. Sediment sources covered 45 topsoil samples in farmland, woodland and grassland, and 24 subsoil samples from gullies. Based on the analysis for 33 properties in these samples, the optimal composite fingerprints including P, Ce, Ga, Rb and ¹³⁷Cs were screened by non-parametric test and multiple discriminant analysis. Taken these 5 fingerprints as group I, radionuclides ¹³⁷Cs and ²¹⁰Pb_ex as group II, Walling-Collins model, the representative of multivariate mixing models, and Bayesian model were used to calculate the sediment contribution for the two potential sources, respectively. Taken each property in the optimal composite fingerprints as group III, analytical solutions for single factor were also used to provide such estimation. The results showed that the contribution ratios of sediment sources based on different fingerprint group and various methods were similar, for example, the ratio provided by group II (topsoil 47.5% and subsoil 52.5%) kept consistent with that calculated by group I (topsoil 44.6% and subsoil 55.4%), while the multivariate mixing model was adopted, the contribution estimated by group II (about half to half) was slightly different with that based on group I (topsoil 58.8% and subsoil 41.2%) in Bayesian model. However, the results provided by these two models were relative great-up to 14.2%, while taken group I as tracers. In group III, P, Ce, Ga, and ¹³⁷Cs could differentiate the sources while applying for the single factor in analytical solutions, the results that half contribution from topsoil and half from gully, were not completely same, but very closed. The differences in sediment contribution might be caused by different principles of models. There were still some spaces to make improvement for different models in order to obtain more reliable results. Attention should also be paid to the gully, as it has caused severe soil erosion and contributed about 50% sediments with less than 1% area ratio in the whole watershed. To reduce sediment derived from the watershed, it is necessary to strengthen prevention and treatment for the gullies to the future land management and soil erosion controls.