Holomorphic dependence for the Beltrami equation in Sobolev spaces

We prove that, given a path of Beltrami differentials on $\mathbb C$ that live in and vary holomorphically in the Sobolev space $W^{l,\infty}_{loc}(Ω)$ of an open subset $Ω\subset \mathbb C$, the canonical solutions to the Beltrami equation vary holomorphically in $W^{l+1,p}_{loc}(Ω)$ for admissible $p > 2$. This extends a foundational result of Ahlfors and Bers (the case $l = 0$). As an application, we deduce that Bers metrics on surfaces depend holomorphically on their input data.