Monofractal and multifractal analysis were combined to quantitatively investigate soil particle size distribution (PSD) in hilly and gully areas of the Loess Plateau, and the effect of land uses on PSD parameters was analyzed. The results indicated that the PSD showed heterogeneous distribution, had multifractal feature, and multifractal spectrum functions f(α)-α were asymmetrical convex curve. There was no significant relationship between capacity dimension (D0) and soil particle contents. Clay contents were positively correlated with the entropy dimension (D1), correlation dimension (D2) and the width of multifractal spectra Δα(P<0.01). The clay contents, silt contents, sand contents and dimension of silt domain (Dsilt) values were significantly influenced by land uses (P<0.01). The multifractal analysis provided an accurate method to characterize the PSD.