The relationship between 3d Chern-Simons theory on the product of the disk and the real line ($D\times \mathbb{R}$) and the chiral WZW model on $S^1\times \mathbb{R}$ was shown in Elitzur et al Nucl.Phys. B326 (1989) 108 (the main details can be found in the top answer to this question, whose notation we follow).
The essential point is that the Chern-Simons gauge field $\tilde{a}$ on the disk can be shown to obey the flatness constraint $\tilde{f}=d\tilde{a}+\tilde{a}\wedge\tilde{a}=0$, which is solved by $$\tilde{a}=-\tilde{d}UU^{-1}. \tag{1}\label{1}$$ Substituting this into the action leads to the chiral WZW model on $S^1\times \mathbb{R}$.
For the quantum theory, one also ought to show that the path integral measure does not have a nontrivial Jacobian, i.e., $$\int \mathcal{D}\tilde{a}~\delta{(\tilde{f})}=\int \mathcal{D}U, \tag{2}\label{2}$$ where $\int \mathcal{D}U$ comes from the Haar measure. My question is, how do we show this?
My attempt: In 1, one way to show this was explained on page 111. They claim that for the change of variables $$\tilde{a}=-\tilde{d}UU^{-1}+\epsilon \tag{3}\label{3}$$ (where $\epsilon$ is a small variation transverse to the space of flat gauge fields), the path integral Jacobian is proportional to $|\textrm{det}(\partial_z+\partial_zU U^{-1})|^2$, and they further claim that this cancels the factor obtained in converting $\delta(\tilde{f})$ to $\delta(\epsilon)$. Here, $z$ is a complex coordinate on $S^1\times \mathbb{R}$.
Here is my attempt to show this. Firstly, by varying \eqref{3}, I obtain $$\delta \tilde{a}=-\tilde{d}(\delta U) U^{-1}+\tilde{d}UU^{-1}\delta U U^{-1}+\delta\epsilon .$$ If the Jacobian is understood to be $|\textrm{det}\frac{\delta\tilde{a}}{\delta U}|$ this seems to imply that $$\int \mathcal{D}\tilde{a}_z\int \mathcal{D}\tilde{a}_{\bar{z}}=|\textrm{det}(\partial_z-\partial_zU U^{-1})U^{-1}|^2\int\mathcal{D}U,$$ which is not exactly what we want.
Next, to show that $$\delta(\tilde{f})\propto \frac{\delta(\epsilon)}{|\textrm{det}(\partial_z+\partial_zU U^{-1})|^2},$$ it seems that we should use a formula of the form $\delta(f(x)) = \delta(x-x_0)/|f'(x_0)|$, together with the explicit form of $|\textrm{det}\frac{\delta \tilde{f}}{\delta \epsilon_z}|$. I was able to show that $$\tilde{f}=\partial_z\epsilon_{\bar{z}}-[\partial_z U U^{-1},\epsilon_{\bar{z}}]-\partial_{\bar{z}}\epsilon_{{z}}+[\partial_{\bar{z}} U U^{-1},\epsilon_{{z}}],$$ but I am not sure how to generalize $\delta(f(x)) = \delta(x-x_0)/|f'(x_0)|$ appropriately.