How to get the equations of motion with the Kalb-Ramond 3-form field strength by varying the bosonic string action?

Question

I have the bosonic string action:

$S= - {1 \over 4 \pi \alpha'} \int_\Sigma d^2\sigma \sqrt{-g} g^{ab} \partial_a X^\mu \partial_b X^\nu G_{\mu\nu}(X)\, +\epsilon^{ab} B_{\mu\nu}(X) \partial_a X^\mu \partial_b X^\nu$

where $G_{\mu\nu}$is the metric for the background spacetime, $g_{ab}$is the worldsheet metric, $B_{\mu\nu}$is the Kalb-Ramond 2-form field and $\epsilon^{ab}$is the totally antisymmetric tensor density.

I'm supposed to get the equations of motion:

$$\partial_a \partial^a X^\mu + \Gamma^\mu_{\nu\rho}\partial_a X^\nu \partial^a X^\rho - \frac{1}{2} H^\mu_{\nu\rho} \partial_a X^\nu \partial_b X^\rho \epsilon^{ab}= 0$$

where $\Gamma^\mu_{\nu\rho}$are the connection components for the background spacetime and $H^{\mu}_{\mu\rho} = G^{\mu\xi}H_{\xi\nu\rho}$are the components of the Kalb-Ramond field strength $H=dB$. When I vary the second term of the action I get a term of the form  $\nabla_a B_{\mu\nu} \partial_b X^\nu \epsilon^{ab} = \nabla_\rho B_{\mu\nu} \partial_a X^\rho \partial_bX^\nu \epsilon^{ab}$, but I don't know how to get the other terms of $H_{\mu\nu\rho} = \frac{1}{2} \nabla_{[\mu} B_{\nu\rho]}$.

What I do instead is assume the worldsheet is the boundary of some 3-dimensional region. Then the second term of the action turns into an integral of $H$ over this region and the variation gives the desired equations of motions.

This seems a bit hacky, as I'm not sure I can assume the worldsheed is a boundary even for a closed string, and obviously not for the open string. I'd also like to analyze the two terms of the action in parallel. Is there a way to get the equations of motion in this form by direct variation of the action?