Experimental and theoretical studies of the statistical properties of surface plasmon polaritons (SPOs) are described. Both classical and non-classical properties of surface plasmons are analysed. The temporal statistical behaviour at low excitation level, as measured by detecting the SPO emitted photon statistics as expressed by the correlation function and the temporal photon count distribution, show that the SPOs preserve the photon statistics of the laser. In the spatial distribution of the plasmon field as measured by an STM, squeezing, i.e. non-classical properties, were found. Independent simple model calculations confirmed the existence of both enhanced EM fields of surface plasmons and their squeezed character.