A special case is this:

Let $u\in\mathcal{E}'(\mathbb{R}^{n+m})$. Then $supp(u)\subseteq\mathbb{R}^n\times\{0\}$ iff $u$ is a linear combination of the form $\sum_{\substack{\alpha\in\{0\}\times\mathbb{N}^m \\ |\alpha|\leq k}} c_\alpha \partial^\alpha \iota(v_\alpha)$ for certain $v_\alpha\in\mathcal{E}'(\mathbb{R}^n)$. Here $\iota(v)$ denotes the distribution on $\mathbb{R}^{n+m}$ defined by $\langle \iota(v),\phi\rangle:=\langle v,\phi_{|\mathbb{R}^n}\rangle$.

The proof idea is to use a taylor expansion of $\phi$ orthogonal to $\mathbb{R}^n$, to show that $u$ vanishes on those functions whose orthogonal derivates up to the $k$th order vanish on $\mathbb{R}^n\times\{0\}$ and deduce the result.

(I suspect that a similar thing could be true for $u\in\mathcal{D}'$. In that case the order of $u$ is only locally bounded, to one would expect that that linear combination has to be replaced by some infinite (but locally finite) linear combination of derivatives)

By localizing and some application of change of coordinates I think this should extend to the following general result: If $u\in\mathcal{E}'(M,E)$ is a distributional section of a vector bundle $E$ on the riemannian manifold $M$ (in your case $E$ is the bundle of $n$-forms) and $S\subseteq M$ an embedded submanifold, then $supp(u)\subseteq S$ iff it is of the form $\sum_{a=0}^k \mathcal{L}_{X_{i1}} \cdots \mathcal{L}_{X_{ij_a}}\iota(v_a)$ for certain $v_a\in\mathcal{E}'(S,E_{|S})$ and certain vector fields $X_{ij} \in \Gamma(NS)$ where $NS\subseteq TM$ is the normal bundle of $S$ and $\iota:\mathcal{E}'(S,E_{|S})\to\mathcal{E}'(M,E)$ is the analogue embedding. Again, it might be the case that some version of this also holds for $\mathcal{D}'$ if one replace the finite sum by an locally finite series.