Chapter 20 Limits Calculus Accumulation and Differential Equations

Year 1 · Mathematics & Statistics · G117-L2

Learn limits, derivatives, integrals and elementary differential equations through a worked sediment-accumulation model with explicit assumptions.

A repeatable workflow from a testable question to an evidence-based interpretation. Original Geology Hub teaching diagram.

Learning objectives: Interpret a derivative as rate, use an integral for accumulation, recognize elementary differential equations, and distinguish a model from measured history.

A limit makes local change precise

A function connects an input with an output, such as time with sediment thickness. A limit describes the value approached as the input approaches a specified point. Continuity means the function’s value agrees with that limiting behaviour where the condition applies. Real geological records can contain abrupt boundaries or missing intervals, so a continuous mathematical curve is often an approximation over a chosen scale.

The derivative is the limiting rate of change of a function. On a graph it is the local tangent slope. If thickness h is measured in metres and time t in years, dh/dt has units of metres per year. An average rate over an interval is not necessarily the rate at every instant. The distinction matters when a short flood deposit is averaged across a long hiatus.

Differentiation rules organize useful models

The power rule gives d(tⁿ)/dt = ntⁿ⁻¹ in the usual domain. Constant factors can be carried through differentiation. Products, quotients and compositions require their own rules; the chain rule accounts for an inner variable changing with another variable. Exponential and logarithmic functions appear naturally in growth, decay and scaling models. Always retain the dimensions of coefficients so a symbolic expression remains physically interpretable.

Integration reverses differentiation in the sense of finding antiderivatives and also calculates accumulated quantities through definite integrals. An indefinite integral includes an arbitrary constant; a definite integral over stated limits gives a value. Substitution simplifies an integral by changing variables, while integration by parts reorganizes products using the product rule. These techniques are tools, not evidence that a selected geological model is correct.

Limits distinguish local behaviour from a single value

A limit asks what output a function approaches as its input approaches a point. A function can approach a value even if no measurement exists exactly at that point. Continuity additionally requires the defined value to agree with the limit. In a geological model, a smooth accumulation curve may be appropriate within one interval but inappropriate across an erosion event. Mathematical continuity is therefore a modelling condition to test against the geological record rather than an assumption every sedimentary succession must satisfy.

Differentiation follows the structure of a function

The power rule handles algebraic powers; exponentials and logarithms have their own derivatives. For a dimensionless variable x, d(e^x)/dx = e^x and d(ln x)/dx = 1/x where x is positive. A product needs the product rule, and a nested function needs the chain rule. If a model contains e^−kt, its time derivative includes the factor −k. Omitting that factor loses both the correct rate and the inverse-time units carried by k.

Integration recovers accumulation and introduces conditions

Substitution simplifies a composite expression by choosing a new variable and transforming its differential consistently. Integration by parts rearranges a product according to ∫u dv = uv − ∫v du. A definite integral includes limits and returns an accumulated value; an indefinite integral returns a family with an arbitrary constant. Integrating a constant discharge of 3 cubic metres per second for 20 seconds gives 60 cubic metres. The result is volume, not a new discharge, because integration multiplies the rate’s units by time.

Order and degree describe an equation differently

The order is the highest derivative present. When an equation is polynomial in its derivatives, degree describes the power of the highest-order derivative after the equation has the appropriate polynomial form. The model dN/dt = −kN is first order and first degree. Separation gives dN/N = −k dt and then an exponential solution for positive N. An initial amount selects one member of that solution family. The equation does not establish whether an actual geological system follows the assumed proportional-loss mechanism.

Second-order models need two independent conditions

A constant-coefficient model x″ + ω²x = 0 has sine and cosine solutions for a fixed positive angular frequency ω. Initial position and initial velocity determine their coefficients. Adding a damping term changes the behaviour, and adding an external forcing term can produce another response. Such models connect calculus with oscillating instruments. They illustrate why knowing an equation’s general solution is insufficient: the physical interpretation also requires units, initial conditions and a justification for the forces or rates represented.

Worked accumulation model

Let a fictional depositional thickness follow h(t) = 2t + 0.1t², with t in thousands of years and h in metres. The coefficients therefore carry units consistent with those definitions. Differentiation gives dh/dt = 2 + 0.2t metres per thousand years. At t = 5, the model rate is 3 metres per thousand years. Thickness at that time is 12.5 metres.

The average rate from zero to five thousand years is 12.5/5 = 2.5 metres per thousand years, different from the final instantaneous rate. Integrating the rate from zero to five reproduces 12.5 metres of accumulation. This agreement checks the algebra. It does not account for compaction, erosion or gaps in deposition, which would need additional assumptions or measurements.

Differential equations link state and rate

A differential equation contains a function and its derivatives. The model dN/dt = −kN states that loss rate is proportional to the amount remaining, with positive k. Its solution N = N₀e⁻ᵏᵗ requires an initial amount N₀. Separation of variables solves this first-order example. The condition k has inverse-time units ensures that the exponent is dimensionless.

Second-order equations can describe an acceleration or restoring process. For example, x″ + ω²x = 0 models ideal undamped oscillation, with solutions determined by initial position and velocity. A real instrument can require damping and forcing terms. Stating an equation’s order and coefficients is useful, but interpreting the boundary or initial conditions is equally important.

For practice, draw thickness, rate and accumulated-area sketches for the worked model. Mark the point at five thousand years and label each vertical axis with its different units. Your final result should explain in words why a steeper thickness curve implies a greater rate and why integration adds a history of rates into a net quantity.

Self-check

What are the units of dh/dt if h is metres and t is years?

Metres per year.

Why does an indefinite integral include a constant?

Differentiation removes a constant, so antiderivatives form a family until an additional condition is given.

Does algebraic consistency validate a geological model?

No. Independent observations must test whether its assumptions fit the geological system.

Sources and further study

Topic framework: the supplied University of Peshawar BS Geology curriculum, session 2004–05 onward. This independently authored learning resource is not affiliated with the university and does not confer a degree or accreditation.

Chapter 20 multiple choice questions

Choose one best answer for each question. Open each answer panel after making your choice.

A model gives thickness h(t) = 2t + 0.1t², with h in metres and t in thousands of years. What rate dh/dt is predicted at t = 5?
1. 2 m per thousand years
2. 3 m per thousand years
3. 12.5 m per thousand years
4. 25 m per thousand years
Answer and explanation
B. 3 m per thousand years The derivative is 2 + 0.2t, which equals 3 at t = 5.
For that same model, what is the average accumulation rate from t = 0 to t = 5?
1. 3 m per thousand years
2. 12.5 m per thousand years
3. 0.5 m per thousand years
4. 2.5 m per thousand years
Answer and explanation
D. 2.5 m per thousand years Thickness increases by 12.5 m over five thousand-year units, giving 12.5/5 = 2.5.
An accumulation rate is measured in metres per year. What does its definite integral over a time interval represent?
1. An instantaneous rate in metres per year squared
2. A dimensionless probability
3. Net accumulated thickness in metres
4. A temperature change in degrees
Answer and explanation
C. Net accumulated thickness in metres Integrating a rate over time gives the accumulated quantity, with rate units multiplied by time.
A loss model is dN/dt = −kN with k positive. Which additional information selects one solution from the family N = N₀e⁻ᵏᵗ?
1. An initial amount N₀
2. Only the sign of k, without an initial condition
3. Only the units of N, without its starting value
4. Only the knowledge that N decreases, without its starting value
Answer and explanation
A. An initial amount N₀ An initial condition fixes the otherwise undetermined solution constant.
A sediment model differentiates and integrates consistently but ignores erosion and compaction. What is the strongest conclusion?
1. The algebra proves a complete depositional history.
2. The model automatically measures all missing time.
3. Correct integration makes compaction physically impossible.
4. Its algebra is consistent, but observations must still test its geological assumptions.
Answer and explanation
D. Its algebra is consistent, but observations must still test its geological assumptions. Mathematical consistency is necessary for calculation but does not validate the physical model.

All 84 chapters

Search This Blog

Geology Hub