Now that we can parameterize surfaces and we can calculate their surface areas, we are able to define surface integrals. First, let’s look at the surface integral of a scalar-valued function. Informally, the surface integral of a scalar-valued function is an analog of a scalar line integral in one higher dimension. The domain of integration of a scalar line integral is a parameterized curve (a one-dimensional object); the domain of integration of a scalar surface integral is a parameterized surface (a two-dimensional object). Therefore, the definition of a surface integral follows the definition of a line integral quite closely. For scalar line integrals, we chopped the domain curve into tiny pieces, chose a point in each piece, computed the function at that point, and took a limit of the corresponding Riemann sum. For scalar surface integrals, we chop the domain region (no longer a curve) into tiny pieces and proceed in the same fashion.

Let [latex]S[/latex] be a piecewise smooth surface with parameterization [latex]{\bf{r}}(u,v)=\langle{x}(u,v),y(u,v)z(u,v)\rangle[/latex] with parameter domain [latex]D[/latex] and let [latex]f(x, y, z)[/latex] be a function with a domain that contains [latex]S[/latex]. For now, assume the parameter domain [latex]D[/latex] is a rectangle, but we can extend the basic logic of how we proceed to any parameter domain (the choice of a rectangle is simply to make the notation more manageable). Divide rectangle [latex]D[/latex] into subrectangles [latex]D_{ij}[/latex] with horizontal width [latex]\Delta{u}[/latex] and vertical length [latex]\Delta{v}[/latex]. Suppose that [latex]i[/latex] ranges from [latex]1[/latex] to [latex]m[/latex] and [latex]j[/latex] ranges from [latex]1[/latex] to [latex]n[/latex] so that [latex]D[/latex] is subdivided into [latex]mn[/latex] rectangles. This division of [latex]D[/latex] into subrectangles gives a corresponding division of [latex]S[/latex] into pieces [latex]S_{ij}[/latex]. Choose point [latex]P_{ij}[/latex] in each piece [latex]S_{ij}[/latex], evaluate [latex]P_{ij}[/latex] at [latex]f[/latex], and multiply by area [latex]S_{ij}[/latex] to form the Riemann sum

[latex]\large{\displaystyle\sum_{i=1}^m\displaystyle\sum_{j=1}^nf(P_{ij})\Delta{S_{ij}}}[/latex].

To define a surface integral of a scalar-valued function, we let the areas of the pieces of [latex]S[/latex] shrink to zero by taking a limit.

Again, notice the similarities between this definition and the definition of a scalar line integral. In the definition of a line integral we chop a curve into pieces, evaluate a function at a point in each piece, and let the length of the pieces shrink to zero by taking the limit of the corresponding Riemann sum. In the definition of a surface integral, we chop a surface into pieces, evaluate a function at a point in each piece, and let the area of the pieces shrink to zero by taking the limit of the corresponding Riemann sum. Thus, a surface integral is similar to a line integral but in one higher dimension.

The definition of a scalar line integral can be extended to parameter domains that are not rectangles by using the same logic used earlier. The basic idea is to chop the parameter domain into small pieces, choose a sample point in each piece, and so on. The exact shape of each piece in the sample domain becomes irrelevant as the areas of the pieces shrink to zero.

Scalar surface integrals are difficult to compute from the definition, just as scalar line integrals are. To develop a method that makes surface integrals easier to compute, we approximate surface areas [latex]\Delta{S_{ij}}[/latex] with small pieces of a tangent plane, just as we did in the previous subsection. Recall the definition of vectors [latex]{\bf{t}}_u[/latex] and [latex]{\bf{t}}_v[/latex]:

[latex]\large{{\bf{t}}_u=\left\langle\frac{\partial{x}}{\partial{u}},\frac{\partial{y}}{\partial{u}},\frac{\partial{z}}{\partial{u}}\right\rangle\text{ and }{\bf{t}}_v=\left\langle\frac{\partial{x}}{\partial{v}},\frac{\partial{y}}{\partial{v}},\frac{\partial{z}}{\partial{v}}\right\rangle}[/latex].

From the material we have already studied, we know that

[latex]\large{\Delta{S}_{ij}\approx||{\bf{t}}_u(P_{ij})\times{\bf{t}}_v(P_{ij})||\Delta{u}\Delta{v}}[/latex].

Therefore,

[latex]\large{\displaystyle\iint_Sf(x,y,z)dS\approx\displaystyle\lim_{m,n\to\infty}\displaystyle\sum_{i=1}^m\displaystyle\sum_{j=1}^nf(P_{ij})||{\bf{t}}_u(P_{ij})\times{\bf{t}}_v(P_{ij})||\Delta{u}\Delta{v}}[/latex].

This approximation becomes arbitrarily close to [latex]\displaystyle\lim_{m,n\to\infty}\displaystyle\sum_{i=1}^m\displaystyle\sum_{j=1}^nf(P_{ij})\Delta{S}_{ij}[/latex] as we increase the number of pieces [latex]S_{ij}[/latex] by letting [latex]m[/latex] and [latex]n[/latex] go to infinity. Therefore, we have the following equation to calculate scalar surface integrals:

[latex]\large{\displaystyle\iint_Sf(x,y,z)dS=\displaystyle\iint_Df({\bf{r}}(u,v)||{\bf{t}}_u\times{\bf{t}}_v||}dA[/latex].

The surface integral of a scalar-valued function allows us to calculate a surface integral by transforming it into a double integral. This equation for surface integrals is analogous to the surface integral equation for line integrals:

[latex]\large{\displaystyle\iint_Cf(x,y,z)ds=\displaystyle\int_a^bf({\bf{r}}(t))||{\bf{r}}^\prime(t)||dt}[/latex].

In this case, vector [latex]{\bf{t}}_u\times{\bf{t}}_v[/latex] is perpendicular to the surface, whereas vector [latex]{\bf{r}}^\prime[/latex] is tangent to the curve.

Example: calculating the surface integral of a cylinder

Calculate surface integral [latex]\displaystyle\iint_S(x+y^2)dS[/latex], where [latex]S[/latex] is cylinder [latex]x^2+y^2=4[/latex], [latex]0\leq{z}\leq3[/latex] (Figure 1).

Figure 1. Integrating function [latex]f(x,y,z)=x+y^2[/latex] over a cylinder.

Show Solution

To calculate the surface integral, we first need a parameterization of the cylinder. Following Example “Parameterizing a Cylinder”, a parameterization is

[latex]{\bf{r}}(u,v)=\langle2\cos{u},2\sin{u},v\rangle, \ 0\leq{u}\leq2\pi, \ 0\leq{v}\leq3[/latex].

The tangent vectors are [latex]{\bf{t}}_u=\langle\sin{u},\cos{u},0\rangle[/latex] and [latex]{\bf{t}}_v=\langle0,0,1\rangle[/latex]. Then,

[latex]{\bf{t}}_u\times{\bf{t}}_v=\begin{vmatrix}{\bf{i}}&{\bf{j}}&{\bf{k}} \\ -\sin{u}&\cos{u}&0 \\ 0&0&1\end{vmatrix}=\langle\cos{u},\sin{u},0\rangle[/latex]

and [latex]||{\bf{t}}_u\times{\bf{t}}_v||=\sqrt{\cos^2u+\sin^2u}=1[/latex]. By the surface integral equation,

[latex]\begin{aligned} \displaystyle\iint_Sf(x,y,z)dS&=\displaystyle\iint_Df({\bf{r}}(u,v))||{\bf{t}}_u\times{\bf{t}}_v|| \ dA \\ &=\displaystyle\int_0^3\displaystyle\int_0^{2\pi}(2\cos{u}+4\sin^2u)dydv \\ &=\displaystyle\int_0^3[2\sin{u}+2u-\sin(2u)]_0^{2\pi}dv=\displaystyle\int_0^34\pi{d}v=12\pi \end{aligned}[/latex].

Example: calculating the surface integral of a piece of a sphere

Calculate surface integral [latex]\displaystyle\iint_Sf(x,y,z)dS[/latex], where [latex]f(x, y, z)=z^{2}[/latex] and [latex]S[/latex] is the surface that consists of the piece of sphere [latex]x^{2}+y^{2}+z^{2}=4[/latex] that lies on or above plane [latex]z=1[/latex] and the disk that is enclosed by intersection plane [latex]z=1[/latex] and the given sphere (Figure 2).

Figure 2. Calculating a surface integral over surface [latex]S[/latex].

Show Solution

Solution

Notice that [latex]S[/latex] is not smooth but is piecewise smooth; [latex]S[/latex] can be written as the union of its base [latex]S_1[/latex] and its spherical top [latex]S_2[/latex], and both [latex]S_1[/latex] and [latex]S_2[/latex] are smooth. Therefore, to calculate [latex]\displaystyle\iint_Sz^2dS[/latex], we write this integral as [latex]\displaystyle\iint_Sz^2d{S_1}+\displaystyle\iint_{S_2}z^2dS[/latex] and we calculate integrals [latex]\displaystyle\iint_{S_1}z^2dS[/latex] and [latex]\displaystyle\iint_{S_2}z^2dS[/latex].

First, we calculate [latex]\displaystyle\iint_{S_1}z^2dS[/latex]. To calculate this integral we need a parameterization of [latex]S_1[/latex]. This surface is a disk in plane [latex]z=1[/latex] centered at [latex](0, 0, 1)[/latex]. To parameterize this disk, we need to know its radius. Since the disk is formed where plane [latex]z=1[/latex] intersects sphere [latex]x^{2}+y^{2}+z^{2}=4[/latex], we can substitute [latex]z=1[/latex] into equation [latex]x^{2}+y^{2}+z^{2}=4[/latex]:

[latex]x^2+y^2+1=4\Rightarrow{x}^2+y^2=3[/latex].

Therefore, the radius of the disk is [latex]\sqrt3[/latex] and a parameterization of [latex]S_1[/latex] is [latex]{\bf{r}}(u,v)=\langle{u}\cos{v},u\sin{v},1\rangle[/latex], [latex]0\leq{u}\leq\sqrt3[/latex],[latex]0\leq{v}\leq2\pi[/latex]. The tangent vectors are [latex]{\bf{t}}_u=\langle\cos{v},\sin{v},0\rangle[/latex] and [latex]{\bf{t}}_v=\langle-u\sin{v},u\cos{v},0\rangle[/latex], and thus

[latex]{\bf{t}}_u\times{\bf{t}}_v=\begin{vmatrix}{\bf{i}}&{\bf{j}}&{\bf{k}} \\ \cos{v}&\sin{v}&0 \\-u\sin{v}&u\cos{v}&0\end{vmatrix}=\langle0,0,u\cos^2v+u\sin^2v\rangle=\langle0,0,u\rangle[/latex].

The magnitude of this vector is [latex]u[/latex]. Therefore,

[latex]\begin{aligned} \displaystyle\iint_{S_1}z^2dS&=\displaystyle\int_0^{\sqrt3}\displaystyle\int_0^{2\pi}f({\bf{r}}(u,v))||{\bf{t}}_u\times{\bf{t}}_v|| \ dvdu \\ &=\displaystyle\int_0^{\sqrt3}\displaystyle\int_0^{2\pi}u \ dvdu \\ &=2\pi\displaystyle\int_0^{\sqrt3} \ udu \\ &=3\pi \end{aligned}[/latex].

Now we calculate [latex]\displaystyle\iint_{S_2} \ dS[/latex]. To calculate this integral, we need a parameterization of [latex]S_2[/latex]. The parameterization of full sphere [latex]x^{2}+y^{2}+z^{2}=4[/latex] is

[latex]{\bf{r}}(\phi,\theta)=\langle2\cos\theta\sin\phi,2\sin\theta\sin\phi,2\cos\phi\rangle, \ 0\leq\theta\leq2\pi, \ 0\leq\phi\leq\pi[/latex].

Since we are only taking the piece of the sphere on or above plane [latex]z=1[/latex], we have to restrict the domain of [latex]\phi[/latex]. To see how far this angle sweeps, notice that the angle can be located in a right triangle, as shown in Figure 3 (the [latex]\sqrt3[/latex] comes from the fact that the base of [latex]S[/latex] is a disk with radius [latex]\sqrt3[/latex]). Therefore, the tangent of [latex]\phi[/latex] is [latex]\sqrt3[/latex], which implies that [latex]\phi[/latex] is [latex]\pi/6[/latex]. We now have a parameterization of [latex]S_2[/latex]:

[latex]{\bf{r}}(\phi,\theta)=\langle2\cos\theta\sin\phi,2\sin\theta\sin\phi,2\cos\phi\rangle, \ 0\leq\theta\leq2\pi, \ 0\leq\phi\leq\pi/3[/latex].

Figure 3. The maximum value of [latex]\phi[/latex] has a tangent value of [latex]{\sqrt3}[/latex]

The tangent vectors are

[latex]{\bf{t}}_\phi=\langle2\cos\theta\cos\phi,2\sin\theta\cos\phi,-2\sin\phi\rangle\text{ and }{\bf{t}}_\theta=\langle-2\sin\theta\sin\phi,u\cos\theta\sin\phi,0\rangle[/latex],

and thus

[latex]\begin{aligned} {\bf{t}}_\phi\times{\bf{t}}_\theta&=\begin{vmatrix} {\bf{i}}&{\bf{j}}&{\bf{k}} \\ 2\cos\theta\cos\phi&2\sin\theta\cos\phi&-2\sin\phi \\ -2\sin\theta\sin\phi&2\cos\theta\sin\phi&0 \end{vmatrix} \\ &=\langle4\cos\theta\sin^2\phi,4\sin\theta\sin^2\phi,4\cos^2\theta\sin\phi+4\sin^2\theta\cos\phi\sin\phi\rangle \\ &=\langle4\cos\theta\sin^2\phi,4\sin\theta\sin^2\phi,4\cos\phi\sin\phi\rangle \end{aligned}[/latex].

The magnitude of this vector is

[latex]\begin{aligned} ||{\bf{t}}_\phi\times{\bf{t}}_\theta||&=\sqrt{16\cos^2\theta\sin^4\phi+16\sin^2\theta\sin^4\phi+16\cos^2\phi\sin^2\phi} \\ &=4\sqrt{\sin^4\phi+\cos^2\phi\sin^2\phi} \end{aligned}[/latex].

Therefore,

[latex]\begin{aligned} \displaystyle\iint_{S_2}z^2dS&=\displaystyle\int_0^{\pi/3}\displaystyle\int_0^{2\pi}f({\bf{r}}(\phi,\theta))||{\bf{t}}_\phi\times{\bf{t}}_\theta|| \ d\theta{d}\phi \\ &=\displaystyle\int_0^{\pi/3}\displaystyle\int_0^{2\pi}16\cos^2\phi\sqrt{\sin^4\phi+\cos^2\phi\sin^2\phi}d\theta{d}\phi \\ &=32\pi\displaystyle\int_0^{\pi/3}\cos^2\phi\sqrt{\sin^4\phi+\cos^2\phi\sin^2\phi}d\phi \\ &=32\pi\displaystyle\int_0^{\pi/3}\cos^2\phi\sin\phi\sqrt{\sin^2\phi+\cos^2\phi}d\phi \\ &=32\pi\displaystyle\int_0^{\pi/3}\cos^2\phi\sin\phi{d}\phi \\ &=32\pi\left[-\frac{\cos^3\phi}3\right]_0^{\pi/6}=32\pi\left[\frac13-\frac{\sqrt3}8\right]=\frac{28\pi}3 \end{aligned}[/latex].

Since

[latex]\displaystyle\iint_Sz^2 \ dS=\displaystyle\iint_{S_1}z^2 \ dS+\displaystyle\iint_{S_2}z^2 \ dS=3\pi\frac{28\pi}3=\frac{37\pi}3[/latex].

Analysis

In this example we broke a surface integral over a piecewise surface into the addition of surface integrals over smooth subsurfaces. There were only two smooth subsurfaces in this example, but this technique extends to finitely many smooth subsurfaces.

Watch the following video to see the worked solution to the above Try It

Scalar surface integrals have several real-world applications. Recall that scalar line integrals can be used to compute the mass of a wire given its density function. In a similar fashion, we can use scalar surface integrals to compute the mass of a sheet given its density function. If a thin sheet of metal has the shape of surface [latex]S[/latex] and the density of the sheet at point [latex](x, y, z)[/latex] is [latex]\rho(x,y,z)[/latex], then mass [latex]m[/latex] of the sheet is [latex]m=\displaystyle\iint_S\rho(x,y,z)dS[/latex].