Surface Integrals of Scalar Fields - Calculus 3

SECTION 1

Historical Context & Motivation

The idea of integrating a function over a surface grew naturally from the broader program of extending calculus to higher dimensions. By the late eighteenth century, mathematicians and physicists routinely computed line integrals to find work done along curves, but many physical problems—heat conduction through a curved shell, gravitational attraction of a thin dome, total charge on a bent conductor—demanded integration over two-dimensional surfaces embedded in three-dimensional space. The conceptual leap was significant: whereas a line integral reduces to a single parameter, a surface integral requires parameterizing a patch of surface with two independent variables and accounting for the way the surface stretches or compresses relative to the parameter plane.

1760 • Euler's Work on Surfaces

Leonhard Euler developed systematic methods for describing curved surfaces analytically, introducing parametric representations that would later underpin surface integration.

1813 • Gauss & Surface Area

Carl Friedrich Gauss formalized the computation of surface area elements through differential geometry, establishing the area magnification factor |

r_u

r_v

| central to modern surface integrals.

1828 • Green's Theorem Extended

George Green published results linking volume integrals to surface integrals, motivating precise definitions of integration over closed surfaces and connecting scalar surface integrals to potential theory.

1850s • Riemann & Integration Theory

Bernhard Riemann placed surface integrals on rigorous footing by extending his theory of integration to higher-dimensional manifolds, treating the surface element dS as a genuine measure.

1860s • Maxwell's Electrodynamics

James Clerk Maxwell's equations required surface integrals of both scalar and vector fields over arbitrary closed surfaces, cementing these tools as essential in mathematical physics.

The central question this lesson addresses is: given a scalar-valued function f(x, y, z) and a curved surface S in ℝ³, how do we define and compute the integral of f over S? This is the natural generalization of integrating a function over a flat region in the plane, but the curvature and orientation of the surface introduce a geometric correction factor—the surface area element dS—that makes the computation both richer and more demanding.

SECTION 2

Core Principles & Definitions

A surface integral of a scalar field generalizes the familiar double integral ∬_D f(x, y) dA over a flat region D to an integral ∬_S f(x, y, z) dS over a curved surface S. The key ideas revolve around parameterizing the surface, computing the infinitesimal area element, and accumulating the contributions of f across the entire surface. Before diving into formulas, it is important to understand the foundational concepts that make this construction well-defined.

Parametric Surface

A surface S is described by a vector-valued function r(u, v) = ⟨x(u,v), y(u,v), z(u,v)⟩ where (u, v) ranges over a domain D in the parameter plane. This maps each point of D to a point on S.

Tangent Vectors & Normal

The partial derivatives

r_u

and

r_v

are tangent vectors to S. Their cross product

r_u

r_v

is a normal vector to the surface at each point, and its magnitude measures local area stretching.

Surface Area Element dS

The infinitesimal area element is dS = |

r_u

r_v

| du dv. This factor corrects for the distortion introduced by the parameterization, converting parameter-plane area into true surface area.

Scalar Surface Integral

The integral ∬_S f dS = ∬_D f(r(u,v)) |

r_u

r_v

| du dv sums the values of f weighted by the area element over the surface, yielding a single scalar quantity.

Independence of Parameterization

The value of ∬_S f dS is intrinsic to the surface and the function f; it does not depend on which parameterization of S is chosen, provided the parameterization is smooth and regular.

✦

KEY TAKEAWAY

Think of a surface integral as weighing a thin, curved sheet. If f(x, y, z) represents the density at each point, the surface integral gives the total mass. Just as a cartographer uses a scale factor to translate areas on a flat map to areas on Earth's curved surface, the factor |r_u × r_v| converts parameter-plane area into true surface area, ensuring each infinitesimal patch is weighted correctly.

SECTION 3

Visual Explanation

From Flat Region to Curved Surface

The diagram below illustrates the fundamental geometry behind a surface integral. On the left, you see the parameter domain D in the uv-plane, subdivided into small rectangles of area Δu Δv. The mapping r(u, v) sends each small rectangle to a curved patch on the surface S (shown on the right). Because the surface may be stretched or compressed, each flat rectangle of area Δu Δv maps to a curvilinear parallelogram whose true area is approximately |r_u × r_v| Δu Δv. The cross product of the two tangent vectors r_u and r_v produces the normal vector at each point, whose magnitude encodes the area magnification.

The parameter domain D (left) is mapped onto the surface S (right) by the vector function r(u, v). The highlighted rectangle Δu Δv in the parameter plane maps to a curved patch on S. The tangent vectors r_u (pink) and r_v (amber) span the tangent plane, while their cross product r_u × r_v (green) is the outward normal whose magnitude gives the area magnification factor.

SECTION 4

Mathematical Framework

We now formalize the construction sketched in the previous sections. Let S be a smooth surface parameterized by r(u, v) = ⟨x(u,v), y(u,v), z(u,v)⟩ for (u, v) ∈ D ⊆ ℝ², and let f : ℝ³ → ℝ be a continuous scalar field defined on S. The surface integral of f over S is defined by pulling back the computation to the parameter domain D and introducing the area element.

SURFACE INTEGRAL OF A SCALAR FIELD

∬_S f(x, y, z) dS = ∬_D f(r(u, v)) |r_u × r_v| du dv

Here, r_u = ∂r/∂u and r_v = ∂r/∂v are tangent vectors, and |r_u × r_v| is the magnitude of their cross product, which equals the area scaling factor dS/(du dv).

Special Case: Graph of a Function z = g(x, y)

When the surface S is given explicitly as z = g(x, y) over a region D in the xy-plane, we can use the natural parameterization r(x, y) = ⟨x, y, g(x, y)⟩. In this case, the tangent vectors are r_x = ⟨1, 0, g_x⟩ and r_y = ⟨0, 1, g_y⟩, and their cross product yields ⟨−g_x, −g_y, 1⟩. Taking the magnitude gives the familiar formula below.

SURFACE INTEGRAL OVER A GRAPH z = g(x, y)

∬_S f dS = ∬_D f(x, y, g(x,y)) √(1 + (∂g/∂x)² + (∂g/∂y)²) dx dy

The factor √(1 + g_x² + g_y²) is recognizable as the surface area element for a graph. When f = 1, this reduces to the surface area formula from Calculus 3.

Special Case: Surfaces of Revolution

SURFACE OF REVOLUTION ABOUT THE x-AXIS

∬_S f dS = ∫_a^b ∫_0^{2π} f(x, h(x)cos θ, h(x)sin θ) · h(x) √(1 + (h'(x))²) dθ dx

Here h(x) ≥ 0 is the radius function generating the surface by rotation, and the parameterization is r(x, θ) = ⟨x, h(x)cos θ, h(x)sin θ⟩. The factor h(x)√(1 + (h'(x))²) is the area element for this geometry.

💡

Deriving dS via the Cross Product

To compute |r_u × r_v|, form the 3 × 3 determinant with i, j, k in the first row, the components of r_u in the second row, and the components of r_v in the third row. The resulting vector has components that are 2 × 2 minors of the Jacobian matrix [r_u | r_v], and its magnitude gives the area element.

SECTION 5

Common Surfaces & Their Area Elements

In practice, many surface integrals arise on standard geometric surfaces—spheres, cylinders, cones, and paraboloids. For each of these, there is a natural parameterization and a known formula for the area element dS. Having these at your fingertips dramatically speeds up computation, since you can skip the cross-product calculation and proceed directly to the double integral over the parameter domain.

Reference card for four common surfaces. Each panel shows the surface geometry, its standard parameterization, the area element dS, and the total surface area formula. Memorizing the area elements for the sphere (a² sin φ dφ dθ), cylinder (a dθ dz), and cone (r√2 dr dθ) will save significant time on exams.

Surface	Parameterization	dS	Parameter Domain
Sphere x² + y² + z² = a²	r(φ, θ) = ⟨a sin φ cos θ, a sin φ sin θ, a cos φ⟩	a² sin φ dφ dθ	φ ∈ [0, π], θ ∈ [0, 2π]
Cylinder x² + y² = a²	r(θ, z) = ⟨a cos θ, a sin θ, z⟩	a dθ dz	θ ∈ [0, 2π], z ∈ [0, h]
Cone z = √(x² + y²)	r(r, θ) = ⟨r cos θ, r sin θ, r⟩	r√2 dr dθ	r ∈ [0, h], θ ∈ [0, 2π]
Paraboloid z = x² + y²	r(r, θ) = ⟨r cos θ, r sin θ, r²⟩	r√(1 + 4r²) dr dθ	r ∈ [0, √c], θ ∈ [0, 2π]
Graph z = g(x, y)	r(x, y) = ⟨x, y, g(x, y)⟩	√(1 + $g_x$ ² + $g_y$ ²) dx dy	(x, y) ∈ D ⊆ ℝ²

Summary of standard parameterizations and area elements for common surfaces

SECTION 6

Worked Example

Let us compute the surface integral of f(x, y, z) = x² + y² + z² over the upper hemisphere S of radius 2 centered at the origin, defined by x² + y² + z² = 4 with z ≥ 0. Physically, this might represent integrating the square of the distance from the origin over a hemispherical shell—a calculation that arises naturally in moment-of-inertia problems.

∬_S (x² + y² + z²) dS over the upper hemisphere of radius 2

1
Step 1 — Parameterize the Surface

Use standard spherical coordinates with a = 2. Set r(φ, θ) = ⟨2 sin φ cos θ, 2 sin φ sin θ, 2 cos φ⟩ where φ ∈ [0, π/2] (upper hemisphere) and θ ∈ [0, 2π]. The area element for a sphere of radius a is dS = a² sin φ dφ dθ = 4 sin φ dφ dθ.

dS = 4 sin φ dφ dθ, with φ ∈ [0, π/2], θ ∈ [0, 2π]
2
Step 2 — Evaluate f on the Surface

On the sphere of radius 2, every point satisfies x² + y² + z² = 4. Therefore f(r(φ, θ)) = (2 sin φ cos θ)² + (2 sin φ sin θ)² + (2 cos φ)² = 4 sin²φ cos²θ + 4 sin²φ sin²θ + 4 cos²φ = 4(sin²φ + cos²φ) = 4. The function f is constant on the sphere, which greatly simplifies the integral.

f(r(φ, θ)) = 4 (constant on the sphere)
3
Step 3 — Set Up the Double Integral

Substituting into the surface integral formula: ∬_S f dS = ∫₀²π ∫₀^{π/2} 4 · 4 sin φ dφ dθ = 16 ∫₀²π ∫₀^{π/2} sin φ dφ dθ. The integrand factors into a product of functions of φ and θ alone, so the iterated integral separates.

∬_S f dS = 16 ∫₀²π dθ · ∫₀^{π/2} sin φ dφ
4
Step 4 — Evaluate the Inner Integral (φ)

∫₀^{π/2} sin φ dφ = [−cos φ]₀^{π/2} = −cos(π/2) + cos(0) = 0 + 1 = 1.

∫₀^{π/2} sin φ dφ = 1
5
Step 5 — Evaluate the Outer Integral (θ) and Combine

∫₀²π dθ = 2π. Multiplying all factors: ∬_S f dS = 16 × 2π × 1 = 32π. As a sanity check, note that the surface area of the upper hemisphere is 2π(2)² = 8π, and f = 4 is constant on S, so the integral should be 4 × 8π = 32π, which confirms our result.

∬_S (x² + y² + z²) dS = 32π

💡

Why the Shortcut Worked

The function f(x, y, z) = x² + y² + z² is constant on any sphere centered at the origin, since every point on such a sphere has the same distance from the center. Whenever the scalar field is constant on the surface—or can be recognized as such—the surface integral reduces to (constant value) × (surface area). Always check for this simplification before diving into computation.

SECTION 7

Comparison with Related Integral Types

Surface integrals of scalar fields sit within a larger ecosystem of integral types in multivariable calculus. Understanding where they fit—and how they differ from line integrals, double integrals, and surface integrals of vector fields—helps clarify both the purpose and the mechanics of each. The table below provides a side-by-side comparison of these integral types along several key dimensions.

Feature	Double Integral ∬_D f dA	Line Integral ∫_C f ds	Surface Integral ∬_S f dS
Domain	Flat region D ⊆ ℝ²	Curve C in ℝ² or ℝ³	Curved surface S in ℝ³
Parameters	2 (x, y)	1 (t)	2 (u, v)
Element	dA = dx dy	ds = \|r'(t)\| dt	dS = \| $r_u$ × $r_v$ \| du dv
Geometric role	Area in the plane	Arc length	Surface area
Typical output	Volume under a surface, total mass of flat plate	Total mass of wire, work (scalar version)	Total mass of shell, total charge on surface
Orientation needed?	No	No (for scalar); Yes (for vector)	No (scalar); Yes (vector/flux)

Comparison of double integrals, line integrals, and surface integrals of scalar fields

✦

KEY TAKEAWAY

The scalar surface integral ∬_S f dS is the two-dimensional analogue of the scalar line integral ∫_C f ds. In both cases, you are integrating a scalar function weighted by a geometric element (arc length ds or surface area dS) that corrects for the curvature and stretching of the domain. The key distinction is dimensionality: a line integral accumulates along a one-dimensional curve, while a surface integral accumulates over a two-dimensional surface, requiring two parameters and a cross product in place of a single derivative magnitude.

SECTION 8

Connection to Vector Surface Integrals & Theorems

The scalar surface integral studied in this lesson is the foundation upon which several deeper constructs in vector calculus are built. The most immediate extension is the surface integral of a vector field (the flux integral), which replaces the scalar f with the dot product F · n̂, where n̂ is the unit outward normal to the surface. This vector version requires an oriented surface and leads directly to the Divergence Theorem (Gauss's theorem) and Stokes' Theorem, the capstone results of multivariable calculus.

Aspect	Scalar Surface Integral ∬_S f dS	Vector Surface Integral ∬_S F · dS
Integrand	Scalar function f(x, y, z)	Vector field F · n̂ (normal component)
Area element	dS = \| $r_u$ × $r_v$ \| du dv (unsigned)	dS = ( $r_u$ × $r_v$ ) du dv (signed, vector)
Orientation required?	No — magnitude makes it positive	Yes — sign depends on chosen normal direction
Result type	Scalar (mass, charge, etc.)	Scalar (net flux through S)
Connection to major theorem	Appears in surface area computations and mass calculations	Central to Divergence Theorem and Stokes' Theorem

Scalar vs. vector surface integrals

Looking ahead, the Divergence Theorem states that ∬_S F · dS = ∭_E ∇ · F dV, converting a surface integral of a vector field over a closed surface to a triple integral over the enclosed volume. Stokes' Theorem relates ∬_S (∇ × F) · dS = ∮_C F · dr, linking a surface integral of the curl to a line integral around the boundary. Mastering the scalar surface integral is the essential prerequisite for both of these powerful theorems, since the flux integral builds directly on the same parameterization and cross-product machinery developed here.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain why the factor |r_u × r_v| is necessary in the surface integral formula. What would go wrong if we simply integrated f(r(u, v)) du dv without this factor?

Show answer

The factor |

r_u

r_v

| accounts for the distortion introduced by the parameterization. Different parameterizations stretch or compress the parameter domain differently relative to the actual surface area, so omitting this factor would produce a result that depends on the particular parameterization chosen rather than on the geometry of the surface itself. Geometrically, a small rectangle du × dv in the parameter plane maps to a curvilinear parallelogram on the surface whose area is approximately |

r_u

r_v

| du dv, not simply du dv. Without this correction, regions where the surface is steeply inclined would be underweighted, and the integral would not represent a genuine surface-area-weighted average.

Question 2: BASIC CALCULATION

Compute ∬_S z dS where S is the portion of the plane z = 2 − x − y lying in the first octant (x ≥ 0, y ≥ 0, z ≥ 0).

Show answer

The surface is a graph z = g(x, y) = 2 − x − y with

g_x

= −1 and

g_y

= −1, so dS = √(1 + 1 + 1) dx dy = √3 dx dy. The first-octant constraint z ≥ 0 gives x + y ≤ 2 with x ≥ 0, y ≥ 0. The integral becomes √3 ∫₀² ∫₀^{2−x} (2 − x − y) dy dx. Evaluating the inner integral: ∫₀^{2−x} (2 − x − y) dy = (2 − x)(2 − x) − (2 − x)²/2 = (2 − x)²/2. Then √3 ∫₀² (2 − x)²/2 dx = (√3/2)[−(2 − x)³/3]₀² = (√3/2)(8/3) = 4√3/3.

Question 3: INTERMEDIATE

Evaluate ∬_S (x² + y²) dS where S is the lateral surface of the cylinder x² + y² = 9, 0 ≤ z ≤ 4.

Show answer

Parameterize the cylinder with r(θ, z) = ⟨3 cos θ, 3 sin θ, z⟩ for θ ∈ [0, 2π] and z ∈ [0, 4]. The area element for a cylinder of radius a = 3 is dS = 3 dθ dz. On the surface, x² + y² = 9 cos²θ + 9 sin²θ = 9 (constant). Thus ∬_S (x² + y²) dS = ∫₀^{2π} ∫₀⁴ 9 · 3 dz dθ = 27 ∫₀^{2π} dθ · ∫₀⁴ dz = 27 × 2π × 4 = 216π. The integral simplifies because x² + y² is constant on the cylinder, a pattern analogous to the hemisphere example where x² + y² + z² was constant on the sphere.

Question 4: APPLIED

A thin hemispherical shell of radius R has a surface mass density σ(x, y, z) = σ₀(1 + z/R), where σ₀ is a positive constant and z ≥ 0. Find the total mass of the shell.

Show answer

The total mass is M = ∬_S σ dS. Parameterize the upper hemisphere: r(φ, θ) = ⟨R sin φ cos θ, R sin φ sin θ, R cos φ⟩ with φ ∈ [0, π/2], θ ∈ [0, 2π], and dS = R² sin φ dφ dθ. On this surface, z = R cos φ, so σ = σ₀(1 + cos φ). Thus M = ∫₀^{2π} ∫₀^{π/2} σ₀(1 + cos φ) R² sin φ dφ dθ = 2πσ₀R² ∫₀^{π/2} (sin φ + sin φ cos φ) dφ. Evaluating: ∫₀^{π/2} sin φ dφ = 1 and ∫₀^{π/2} sin φ cos φ dφ = 1/2. Therefore M = 2πσ₀R²(1 + 1/2) = 3πσ₀R². The mass exceeds the uniform-density result of 2πσ₀R² because the density increases toward the pole (z = R).

Question 5: CRITICAL THINKING

Let S be a smooth, closed surface enclosing a region E in ℝ³, and let f(x, y, z) = 1. Show that the surface integral ∬_S 1 dS gives the surface area of S. Then argue that this value is always strictly greater than the surface area of any sphere enclosing the same volume as E, unless S is itself a sphere. (Hint: use the isoperimetric inequality in three dimensions.)

Show answer

When f = 1, the surface integral ∬_S 1 dS = ∬_D |

r_u

r_v

| du dv simply sums up infinitesimal area elements over the entire surface, which by definition is the surface area A(S). For the second part, the three-dimensional isoperimetric inequality states that among all closed surfaces enclosing a fixed volume V, the sphere achieves the minimum surface area: A(S) ≥

(36πV²)^{1/3}

, with equality if and only if S is a sphere of radius r =

(3V/(4π))^{1/3}

. This follows from the classical Brunn-Minkowski theory. Therefore, if S is not a sphere, its surface area ∬_S dS must be strictly greater than the surface area of the sphere enclosing the same volume. This result connects the purely computational notion of a surface integral to a deep geometric optimization principle.

SUMMARY

Lesson Summary

Summary

The surface integral of a scalar field extends double integration from flat regions to curved surfaces in ℝ³. Given a surface S parameterized by r(u, v) and a continuous scalar function f, the integral is computed as ∬_D f(r(u,v)) | $r_u$ × $r_v$ | du dv, where the cross-product magnitude serves as the area magnification factor that converts parameter-plane area into true surface area. For surfaces given as graphs z = g(x, y), this simplifies to the factor √(1 + $g_x$ ² + $g_y$ ²).

The construction is independent of the choice of parameterization and does not require an oriented surface, distinguishing it from the vector surface integral (flux integral) that you will study next. Common applications include computing the total mass of a thin shell, the total charge on a surface, and the surface area itself (when f = 1). Mastery of this integral is the essential prerequisite for understanding the Divergence Theorem and Stokes' Theorem.