We study the problem of the influence of the pairing interaction on two- and four-particle surface correlations in deformed nuclei. Taking the overlap between the two- or four-particle wave function and a di-nucleon or alpha-particle wave function, respectively, one obtains cluster probability distributions, as functions of the center-of-mass coordinate of the considered particles. For particles moving in pure Nilsson orbits the probability is localized in the intrinsic frame in different region of the nuclear surface, according to the K-quantum numbers of the considered orbitals. The inclusion of the pairing interaction leads to a probability distributed over the entire surface but at the same time to values of the total spectroscopic factors orders of magnitude larger than those associated with pure Nilsson orbits.