I took a basic multivariable calculus class as an undergrad where I saw Green's theorem, Stokes theorem, and Divergence theorem without proof. I know how to use and make sense of them in dimension $n \leq 3$.
I have recently started learning PDEs from Evans and take the 'Guass-Green' theorem for granted since I haven't learnt about differential forms/manifolds yet.
$$\int_{U}{\frac{\partial u}{\partial x_i}}dx=\int_{\partial U}u\nu^{i}dS\;\;\;\;(i=1,\ldots,n),$$
I am trying to understand the proof of MVP for a harmonic function $u:U \to \mathbb{R}$, where $U \subseteq \mathbb{R}^n$ and $u \in C^2(\overline{U})$ which says $$u(x)=\def\avint{\mathop{\,\rlap{-}\!\!\int}\nolimits} \avint_{\partial B(x,r)} u(y) d\sigma $$ for all balls $\overline{B(x,r)} \subseteq U$ where $d\sigma$ is the surface element.
The proof begins by fixing $x$ and defining
$$ \phi(r):= \def\avint{\mathop{\,\rlap{-}\!\!\int}\nolimits} \avint_{\partial B(x,r)} u(y) d\sigma$$ and showing $\phi'(r)=0$ to get
$$\phi(r)=\phi(0)=u(x)=\def\avint{\mathop{\,\rlap{-}\!\!\int}\nolimits} \avint_{\partial B(x,r)} u(y) d\sigma$$.
My question is what is the expression for $d \sigma$? How does it depend on $y$ in the above integrals?
I know when using hyperspeherical coordinates, we get $$dy=dy_1dy_2 \ldots dy_n=r^{n-1}\sin^{n-2}(\theta_1)\sin^{n-3}(\theta_2) \ldots \sin^{n-2}(\theta_{n-2})dr d\theta_1d\theta_2 \ldots d\theta_{n-1}=r^{n-1}dr d\sigma$$
How can I write $d \sigma$ explicitly in terms of $y$ and make sure I am using change of variables correctly?