Let $\phi:\mathbb{R}^2\to\mathbb{R}$ such that $\phi^{-1}(0)\neq\emptyset$ and $\phi\in C^1(W)$ where $W$ is a compact neighborhood of $\phi^{-1}(0)$, with $\nabla\phi\neq 0$ in $W$. So there is some $t_0>0$ such that for all $|t|<t_0$ we have that the set $\{\phi=t\}\subset W$ (even compact inclusion can be assumed). How can we prove that for any $f\in L^1(\mathbb{R}^2)$ the following relation holds:
$$\lim\limits_{t\to 0} \int_{\{\phi=t\}} f(x)\ d\mathcal{H}^1=\int_{\{\phi=0\}} f(x)\ d\mathcal{H}^1 $$
?
I found this result proved in the appendix of the article of G.H.Cottet, E. Maitre - A level set method for fluid-structure interaction with immersed surface (that can be downloaded from here for free https://www.researchgate.net/publication/229059739_A_level_set_method_for_fluid-structure_interactions_with_immersed_surfaces). But the proof uses a $C^2$ regularity for $\phi$. The idea of the proof used there (to generate a family of diffeomorphisms that push $\phi=t$ to $\phi=0$ using orthogonal trajectories) cannot be extended to $C^1$ regularity. The same idea is posted also on MS:https://math.stackexchange.com/questions/505332/parametrization-of-level-sets-of-a-smooth-function
However this result seems to be a consequence of the inverse function theorem. I made a similar post here: Continuity of Hausdorff measure on level sets. The solution posted there has a problem: the inverse function theorem cannot be applied in all points of $\phi=0$ (only on those $x\in\phi^{-1}(0)$ in which $\dfrac{\nabla\phi(x)}{|\nabla\phi (x)|}\neq\pm (0,1)$).
Maybe the following question can be useful Continuity of $r\mapsto\int_{\Sigma\cap B_r(x)}f^2d\mu$.