In the context of $L^2$ space, it is usually stated that any square-integrable function can be expanded as a linear combination of Spherical Harmonics: $$ f(\theta,\varphi)=\sum_{\ell=0}^\infty \sum_{m=-\ell}^\ell f_\ell^m \, Y_\ell^m(\theta,\varphi)\tag 2 $$ where $Y_\ell^m( \theta , \varphi )$ are the Laplace spherical harmonics.
The context here is important because this equality holds only in the sense of the $L^2$-norm:
$$\lim_{N\to\infty} \int_0^{2\pi}\int_0^\pi \left|f(\theta,\varphi)-\sum_{\ell=0}^N \sum_{m=- \ell}^\ell f_\ell^m Y_\ell^m(\theta,\varphi)\right|^2\sin\theta\, d\theta \,d\varphi = 0.$$
Do we also have pointwise convergence almost everywhere?