Below is an example demonstrating the differences in mean curvatures and one other quantity, and below is my question. Please comment if you need me to explain in more detail or clarify notation, or if there are any mistakes.
Example 1 Consider the manifold with coordinates given by $$(x,y,x+y, x^2).$$ At each point, there are two unit normal vectors $$n_1 = \frac{\sqrt{3}}{3} [1, 1, -1, 0]^T$$ and $$n_2 = \frac{\sqrt{3}}{3\sqrt{8x^2+3}} [4x,-2x,2x,-3].$$ The vector of mean curvatures in the directions of each of these normals is given by $$c=\left[0, -\frac{2\sqrt{3}}{(8x^2+3)^{3/2}}\right]^T.$$ Further, if we let $q$ be the vector with components $q^r = \operatorname{Tr}{N (Dn_r) N^T}$ where $Dn_r$ is the Jacobian derivative of the vector $n_r$, $r=1,2$ and $N$ is the matrix with $r$-th row equal to $n_r$, then we get that $q=\vec{0}$.
Brownian motion on this manifold has drift $\mu =N^T(q+c)$, i.e. $$\mu = \frac{1}{(8x^2+3)^2}[-8x, 4x, -4x, 6]^T.$$
Example 2 Now, permute the last two coordinates so that instead we have $$(x,y,x^2,x+y).$$ Now, we can obtain unit normal vectors $$n_1 = \frac{1}{\sqrt{4x^2+1}} [2x,0,-1,0]^T,$$ $$n_2 = \frac{1}{\sqrt{(4x^2+1)(8x^2+3)}}\left[1, 4x^2+1, 2x, -(4x^2+1)\right]^T,$$ mean curvature vector $c$, $$c=\left[-\frac{1}{(4x^2+1)^{3/2}}, -\frac{2x(16x^2+5)}{(4x^2+1)^{3/2}(8x^2+3)^{3/2}}\right]^T$$ and now a non-zero $q$, as defined above, $$q=\frac{1}{\sqrt{(4x^2+1)(8x^2+3)}}\left[\frac{1}{(4x^2+1)\sqrt{8x^2+3}}, -\frac{2x}{4x^2+1}\right]^T$$ and drift of a Brownian motion $$\mu = \frac{1}{(8x^2+3)^2}[-8x, 4x, -4x, 6]^T,$$ again from $\mu=N^T(q+c)$.
My questions:
Why are the mean curvatures different for these two permutations of coordinates of, essentially, the same manifold? (From a probabilist perspective, I get the same SDE for Brownian motion on both manifolds--this might be a blasphemous viewpoint to differential geometers, so please excuse my ignorance!)
What is the geometric interpretation of $q^r = \operatorname{Tr}(N (D n_r)N^T)$ and why does it change depending on the permutation of coordinates (or rather on the order of orthonormalizing the rows of $Df$ where $f(x,y,z,w)=0$ defines the manifold, via component definitions e.g. $f^1(x,y,z,w)=x+y-z$ and and $f^2(x,y,z,w)=x^2-w$) for the first presentation.
Question 1: The mean curvature measures the average deviation of the surface from its tangent planes, and it depends on the normal vectors and their derivatives. In Example 1, the normal vectors have a specific form that is different from Example 2, and this difference in normal vectors leads to different mean curvature vectors. Essentially, the way the surface deviates from its tangent planes is different for the two examples due to the different normal vectors.
Provided your probabilistic perspective, the SDE for Brownian motion only depends on the Riemannian metric of the manifold, which is the same for both examples. However, the drift term of the SDE, which depends on the mean curvature, is different for the two examples due to the difference in normal vectors and mean curvature vectors.