Zonal averaging of partial derivatives

Given a quantity $Q(x,y,t)$, the zonal average operator $[Q] = \frac{1}{2\pi}\int_0^{2\pi} Q\:\mathrm{d} \lambda$, and zonal anomaly $Q^\star$ such that $Q = [Q] + Q^\star$, my text book says that zonally averaging

$$\frac{\partial Q}{\partial t} + \frac{\partial}{\partial x}(uQ) + \frac{\partial}{\partial y}(vQ) = S$$ gives $$\frac{\partial [Q]}{\partial t} + [v] \frac{\partial [Q]}{\partial y} + \frac{\partial}{\partial y}[v^\star Q^\star] = [S]$$ I'm unable to see how to arrive at this answer. What is the method or intuition for applying the zonal averaging?

I figured it out. It's clear that the zonally averaged $\partial/\partial x$ is zero. The terms involving $\partial/\partial y$ come from expanding $vQ = ([v] + v^\star)([Q] + Q^\star)$, then noticing that $[v^\star [Q]] = [[v] Q^\star] = 0$. I also had to convince myself that $[\partial Q/\partial t] = \partial [Q] / \partial t$. I'd still be interested in a more rigorous answer, however.