Let $\mathcal{M}$ be a smooth manifold of positive dimension $n$ equipped with a smooth density $dμ_{\mathcal{M}}$. Let $A$ be a polyhomogeneous elliptic pseudo-differential operator of positive order $m$ on $\mathcal{M}$ which is symmetric for the $L^2$ scalar product defined by $dμ_{\mathcal{M}}$. For each $L>0$, the space $U_L=\bigoplus_{λ\leq L}Ker(A-λId)$ is a finite dimensional subspace of $C^\infty(\mathcal{M})$. Let $Π_L$ be the spectral projector onto $U_L$. Given $s\in\mathbb{R}$, we compute the asymptotics of the integral kernel $K_L$ of $Π_LA^{-s}$ in the cases where $n>ms$ and $n=ms$ respectively. Next, assuming that $\mathcal{M}$ is closed, let $(e_n)_{n\in\mathbb{N}}$ and $(λ_n)_{n\in\mathbb{N}}$ be the sequence of $L^2$ normalized eigenfunctions and eigenvalues of $A$ where the latter sequence organized in increasing order. Let $(ξ_n)_{n\in\mathbb{N}}$ be a sequence of independent centered gaussians of variance $1$. We fix a parameter $s\in\mathbb{R}$ such that $n\geq ms$ and consider the family $(φ_L)_{L>0}$ of smooth random fields on $\mathcal{M}$ defined by \[φ_L=\sum_{0<λ_j\leq L}λ_j^{-\frac{s}{2}}ξ_je_j\] for each $L>0$. It turns out that the covariance function of $φ_L$ is $K_L$. Using this information, we apply the derived asymptotics to study the zero set of $φ_L$. If $n>ms$ then the number of components of the zero set of $φ_L$ concentrates around $aL^{\frac{n}{m}}$ for some positive constant $a$. On the other hand, if $n=ms$, each Betti number of the zero set has an expectation bounded by $C\ln\Big(L^{\frac{1}{m}}\Big)^{-\frac{1}{2}}L^{\frac{n}{m}}$ where $C$ is an explicit constant. When $\mathcal{M}$ is a closed surface with a Riemmanian metric, $A$ is the Laplacian and $dμ_\mathcal{M}$ is the Riemmanian volume, $C$ equals $\frac{1}{4π^2}\sqrt{\frac{3}{2}}Vol(\mathcal{M})$.