Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

In spherical polars, $$x=r\cos(\phi)\sin(\theta)$$ $$y=r\sin(\phi)\sin(\theta)$$ $$z=r\cos(\theta)$$ I want to work out an integral over the surface of a sphere - ie $r$ constant. I'm able to derive through scale factors, ie $\delta(s)^2=h_1^2\delta(\theta)^2+h_2^2\delta(\phi)^2$ (note $\delta(r)=0$), that: $$h_1=r\sin(\theta),h_2=r$$ $$dA=h_1h_2=r^2\sin(\theta)$$

I'm just wondering is there an "easier" way to do this (eg. Jacobian determinant when I'm varying all 3 variables). I know you can supposedly visualize a change of area on the surface of the sphere, but I'm not particularly good at doing that sadly.

share|cite|improve this question
You can try having a look here, perhaps you'll find something useful: – Giuseppe Negro Apr 14 '12 at 14:36
up vote 23 down vote accepted

I've come across the picture you're looking for in physics textbooks before (say, in classical mechanics). A bit of googling and I found this one for you!area element on sphere

Alternatively, we can use the first fundamental form to determine the surface area element. Recall that this is the metric tensor, whose components are obtained by taking the inner product of two tangent vectors on your space, i.e. $g_{i j}= X_i \cdot X_j$ for tangent vectors $X_i, X_j$. We make the following identification for the components of the metric tensor, $$ (g_{i j}) = \left(\begin{array}{cc} E & F \\ F & G \end{array} \right), $$ so that $E = <X_u, X_u>, F=<X_u,X_v>,$ and $G=<X_v,X_v>.$

We can then make use of Lagrange's Identity, which tells us that the squared area of a parallelogram in space is equal to the sum of the squares of its projections onto the Cartesian plane: $$|X_u \times X_v|^2 = |X_u|^2 |X_v|^2 - (X_u \cdot X_v)^2.$$
Here's a picture in the case of the sphere:

Projecting surface element onto the Cartesian plane

This means that our area element is given by $$ dA = | X_u \times X_v | du dv = \sqrt{|X_u|^2 |X_v|^2 - (X_u \cdot X_v)^2} du dv = \sqrt{EG - F^2} du dv. $$

So let's finish your sphere example. We'll find our tangent vectors via the usual parametrization which you gave, namely, $X(\phi,\theta) = (r \cos(\phi)\sin(\theta),r \sin(\phi)\sin(\theta),r \cos(\theta)),$ so that our tangent vectors are simply $$ X_{\phi} = (-r\sin(\phi)\cos(\theta),r\cos(\phi)\sin(\theta),0), \\ X_{\theta} = (-r\cos(\phi)\sin(\theta),r\sin(\phi)\cos(\theta),-r\sin(\theta)) $$ Computing the elements of the first fundamental form, we find that $$ E = r^2 \sin^2(\theta), \hspace{3mm} F=0, \hspace{3mm} G= r^2. $$ Thus, we have $$ dA = \sqrt{r^4 \sin^2(\theta)}d\theta d\phi = r^2\sin(\theta) d\theta d\phi $$

share|cite|improve this answer
Yea I saw that too, I'm just wondering if there's some other way similar to using Jacobian (if someday I'm asked to find it in a self-invented set of coordinates where I can't picture it) – Rcwt Apr 14 '12 at 15:10
@R.C. I've edited my response for you. Perhaps this is what you were looking for... ? – tentaclenorm Apr 26 '12 at 4:58

When you have a parametric representatuion of a surface $$S:\quad (u,v)\ \mapsto\ {\bf x}(u,v)$$ then an infinitesimal rectangle $[u, u+du]\times [v,v+dv]$ in the parameter plane is mapped onto an infinitesimal parallelogram $dP$ having a vertex at ${\bf x}(u,v)$ and being spanned by the two vectors ${\bf x}_u(u,v)\, du$ and ${\bf x}_v(u,v)\,dv$. The area of this parallelogram is $${\rm d}\omega:=|{\bf x}_u(u,v)\times{\bf x}_v(u,v)|\ {\rm d}(u,v)\ .$$ If you are given a "surface density“ ${\bf x}\mapsto \rho({\bf x})$ $\ ({\bf x}\in S)$ then the integral $I(S)$ of this density over $S$ is then given by $$I(S)=\int_B \rho\bigl({\bf x}(u,v)\bigr)\ {\rm d}\omega = \int_B \rho\bigl({\bf x}(u,v)\bigr)\ |{\bf x}_u(u,v)\times{\bf x}_v(u,v)|\ {\rm d}(u,v)\ ,$$ where $B$ is the parameter domain corresponding to the exact piece $S$ of surface.

Now this is the general setup. When your surface is a piece of a sphere of radius $r$ then the parametric representation you have given applies, and if you just want to compute the euclidean area of $S$ then $\rho({\bf x})\equiv1$.

Therefore in your situation it remains to compute the vector product ${\bf x}_\phi\times {\bf x}_\theta$ as a function of $\phi$ and $\theta$, resp., the absolute value of this product, and then you have to integrate over the desired parameter domain $B$.

You have explicitly asked for an explanation in terms of "Jacobians". The vector product $\times$ is the appropriate surrogate of that in the present circumstances, but in the simple case of a sphere it is pretty obvious that ${\rm d}\omega=r^2\sin\theta\,{\rm d}(\theta,\phi)$.

share|cite|improve this answer

There is yet another way to look at it using the notion of the solid angle. Then the area element has a particularly simple form: $$dA=rd\Omega$$

share|cite|improve this answer
Shouldn't be $r^2$? – Physicist137 Jul 31 '15 at 13:33

Because of the striking fact that

(a) The area of [a slice of the spherical surface between two parallel planes (within the poles)] is proportional to its width.

. . . here's a rarely (if ever) mentioned way to integrate over a spherical surface. We assume the radius = 1.

(b) Note that every point on the sphere is uniquely determined by its z-coordinate and its counterclockwise angle phi, 0 <= phi < 2pi, from the half-plane y = 0, x >= 0.

From (a) and (b) it follows that an element of area on the unit sphere centered at the origin in 3-space is just dphi dz.

Then the integral of a function f(phi,z) over the spherical surface is just

  Integral_{-1 <= z <= 1, 0 <= phi <= 2pi} f(phi,z) dphi dz.
share|cite|improve this answer

The answers above are all too formal, to my mind. The straightforward way to do this is just the Jacobian. The Jacobian is the determinant of the matrix of first partial derivatives. The first row is $\partial r/\partial x$, $\partial r/\partial y$, etc, the second the same but with $r$ replaced with $\theta$ and then the third row replaced with $\phi$. So to compute each partial you hold the other variables constant and just differentiate with respect to the variable in the denominator, e.g. $r=\sqrt{x^2+y^2+z^2}$

so $\partial r/\partial x = x/r $. You then just take the determinant of this 3-by-3 matrix, which can be done by cofactor expansion for instance.

share|cite|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.