A well known change of variable theorem of integral calculus says the following.
Let $X$ and $Y$ be two open sets in $\mathbb{R}^2$ and $\phi: X \rightarrow Y$ be a diffeomorphism. Then \begin{eqnarray}\tag{1} \int\limits_{Y} f(y_1,y_2) dy_1 dy_2=\int\limits_{X} f(\phi(x_1,x_2)) \left|J_{\phi}(x_1,x_2)\right| dx_1 dx_2, \end{eqnarray} where $J_{\phi}$ denotes the Jacobian matrix of $\phi.$
What is the corresponding change of variable formula if the integrand involves derivatives. More precisely, for $f \in C^1$ what is the change of variable formula if the integrand is \begin{eqnarray} \int\limits_{Y} \left(\partial_1f(y_1,y_2)+\partial_2f(y_1,y_2)\right) dy_1 dy_2=\int\limits_{X} \cdot\cdot\cdot\cdot dx_1 dx_2. \end{eqnarray} How to prove the change of variable formula for such integrand using the change of formula (1) ?
One element of confusion is coming from a somewhat imprecise notation. Writing $\partial_1 f$ is asking for trouble when you invoke a change of variables. If instead you write $\frac{\partial f}{\partial y_1}$ or $\partial_{y_1} f$ this removes (in principle) any ambiguity. For example if $f(y_1,y_2)=y_1^2 y_2$ then (omitting domains) $$\int \partial_{y_1} f (y_1,y_2) dy_1\, dy_2 =\int 2 y_1 y_2 \, dy_1\, dy_2 $$ when changing variables gives $$ \int \partial_{y_1}f\circ \phi(x_1,x_2) \; J(\phi)\; dx_1 dx_2 = \int 2\phi_1(x_1,x_2)\phi_2(x_1,x_2)\; J(\phi)\;dx_1 dx_2$$