Any homogeneous polynomial P (x, y, z) of degree d, being restricted to a unit sphere S 2 , admits essentially a unique representation of the form, where L kj 's are linear forms in x, y and z and λ k 's are real numbers. The coefficients of these linear forms, viewed as 3D vectors, are called multipole vectors of P . In this paper we consider similar multipole representations of polynomial and analytic functions on other quadratic surfaces Q(x, y, z) = c, real and complex. Over the complex numbers, the above representation is not unique, although the ambiguity is essentially finite. We investigate the combinatorics that depicts this ambiguity. We link these results with some classical theorems of harmonic analysis, theorems that describe decompositions of functions into sums of spherical harmonics. We extend these classical theorems (which rely on our understanding of the Laplace operator ∆ S 2 ) to more general differential operators ∆ Q that are constructed with the help of the quadratic form Q(x, y, z). Then we introduce modular spaces of multipoles. We study their intricate geometry and topology using methods of algebraic geometry and singularity theory. The multipole spaces are ramified over vector or projective spaces, and the compliments to the ramification sets give rise to a rich family of K(π, 1)-spaces, where π runs over a variety of modified braid groups.