Let $(\Sigma,g)$ be an $n$-dimensional Riemannian manifold. For a differential form $\alpha$, given $d^{*}\alpha=0$, where $d^*$ is the codifferential with respect to $g$, can we rewrite the equation $d^*d\alpha=0$ as a divergence-form strongly elliptic system of equations using local coordinates (just like the case when $\alpha$ is a function) ?

My purpose is to improve as much as possible the regularity of $\alpha$ given that $\alpha\in L^2$ in case $g\in W^{1,n}$ is not continuous. We have $\Delta \alpha=0$, but in local coordinates this is an elliptic system of non-divergence form, which is not convenient to apply regularity results as in the function case. That's why I want to see if it's possible to insert the condition $d^*\alpha=0$ to $d^*d\alpha=0$ to get a divergence-form system.