[step:Verify that the scalar Dolbeault operator squares to zero in coordinates]We recall the local scalar computation. Choose a holomorphic coordinate chart $(V,\varphi)$ with $V \subset U$ and holomorphic coordinates $z_1,\dots,z_n$ on $V$. For each $i \in \{1,\dots,r\}$, write the scalar form $\alpha_i|_V$ as
\begin{align*}
\alpha_i|_V = \sum_{|I|=p,\, |J|=q} a_{i,I,J}\, dz_I \wedge d\bar z_J,
\end{align*}
where $I$ and $J$ range over strictly increasing multi-indices, each coefficient is a smooth map $a_{i,I,J}: V \to \mathbb{C}$, and $dz_I \wedge d\bar z_J$ denotes the corresponding coordinate basis form.
For a fixed coefficient $a_{i,I,J}$, the scalar Dolbeault operator differentiates only in the anti-holomorphic coordinate directions. Thus $\bar\partial^2$ contributes second derivatives of the form
\begin{align*}
\frac{\partial^2 a_{i,I,J}}{\partial \bar z_k \partial \bar z_\ell}\, d\bar z_\ell \wedge d\bar z_k \wedge dz_I \wedge d\bar z_J
\end{align*}
summed over $k,\ell \in \{1,\dots,n\}$, up to the fixed sign convention determined by the placement of the $d\bar z_k$ factors. The terms with $k=\ell$ vanish because $d\bar z_k \wedge d\bar z_k=0$. For $k \neq \ell$, the term indexed by $(k,\ell)$ cancels the term indexed by $(\ell,k)$ because smoothness of $a_{i,I,J}$ gives equality of mixed partial derivatives, while
\begin{align*}
d\bar z_\ell \wedge d\bar z_k = - d\bar z_k \wedge d\bar z_\ell.
\end{align*}
Therefore $\bar\partial^2 \alpha_i = 0$ on $V$. Since these coordinate neighbourhoods cover $U$, $\bar\partial^2 \alpha_i = 0$ on $U$.[/step]