Environmental Chemistry as a Systems Science

Compartments, fluxes, residence times and box models

Lesson 4001 of 4,500 · Environmental Chemistry

Learning objectives

Introduction

A pollutant released from a chimney does not stay put. It mixes into the air, reacts, dissolves in rain, settles onto soil, is washed into rivers and may end up in fish or sediment. Environmental chemistry is therefore not the study of isolated reactions in a flask, but of chemistry happening inside a connected system. To make sense of it we divide the environment into compartments, measure the flows of matter between them and ask how long substances stay in each one. This page introduces that systems toolkit, which underpins everything else in this unit.

Core explanation

Compartments. The Earth system is conventionally divided into the atmosphere , hydrosphere (oceans, lakes, rivers, groundwater), lithosphere (rocks, soils, sediments) and biosphere (living organisms). Each compartment may be subdivided: the atmosphere into troposphere and stratosphere, the ocean into a well-mixed surface layer and the deep ocean. A compartment is useful when its contents can be treated as reasonably uniform.

Burden and fluxes. The amount of a substance in a compartment is its burden (or inventory), M, often given in teragrams (1 Tg = 10¹² g) or petagrams (1 Pg = 10¹⁵ g). Matter enters through sources (emissions, chemical production, transport in) and leaves through sinks (chemical loss, deposition, transport out). Each is a flux , F, with units of mass per time.

Mass balance. The fundamental equation of any box model is conservation of mass:

dM/dt = F in − F out

If inflow exceeds outflow the burden grows; if outflow exceeds inflow it falls.

Steady state. When F in = F out, dM/dt = 0 and the burden is constant. Many natural cycles are close to steady state on human timescales, but human activity often pushes them away from it — the rising atmospheric CO₂ burden is the classic example.

Residence time. The average time a molecule spends in a compartment is the residence time (or lifetime), τ:

τ = M / F out

At steady state this also equals M / F in. A substance with a short residence time (days) is patchy and responds quickly to changes in emissions; one with a long residence time (decades) becomes well mixed around the globe and responds slowly.

First-order removal. Very often the loss rate is proportional to the amount present: F out = kM, where k is a first-order rate constant. Then τ = 1/k, exactly as for a first-order reaction in kinetics. With a constant source S, the box obeys dM/dt = S − kM, whose solution approaches the steady-state burden M ss = S/k with an e-folding time of τ. After one residence time the burden has moved 63% of the way to its new steady state; after three residence times about 95%.

Coupled boxes. Real systems are modelled with several boxes linked by exchange fluxes, each often first-order in the donor reservoir. Such models reveal that the slowest step — for example exchange with the deep ocean — can control how the whole system responds.

Formulae

dM/dt = F in − F out. Residence time τ = M / F out. First-order loss: F out = kM, so τ = 1/k. Steady-state burden with constant source S: M ss = S/k = Sτ. Approach to steady state from zero: M(t) = M ss(1 − e^(−t/τ)).

Step-by-step reasoning

To build a one-box model for a substance:

1. Define the compartment and assume it is well mixed. 2. List all sources and add them to get F in. 3. List all sinks; express each as a rate constant times the burden if first-order. 4. Write dM/dt = F in − kM, where k is the sum of the individual sink rate constants. 5. Solve for steady state (M ss = F in/k) and the residence time (τ = 1/k). 6. Check the well-mixed assumption by comparing τ with the mixing time of the compartment.

Visual explanation

Draw a bath tub. The tap is the source flux, the plughole the sink flux and the water level the burden. The level stops changing when the tap and plughole flows match. If the plughole drains faster when the tub is fuller, the level always settles at a definite height set by the tap flow.

Real-world analogy

A school is a box whose pupils are the burden. If 200 pupils join each year and each stays for 5 years, the school settles at about 1000 pupils. Residence time is simply the population divided by the number leaving each year.

Real-world example

Atmospheric methane has a burden of about 5000 Tg and total sinks of roughly 550 Tg per year, giving a residence time of about 9 years. Because this is long compared with the roughly one-year time for air to mix between the hemispheres, methane is fairly evenly distributed around the globe, with only a small north–south gradient.

Why?

Why does the residence time control how evenly a substance is spread? A molecule can only be carried far from its source if it survives long enough. Species lasting hours stay near their sources, while species lasting years are mixed by winds throughout the troposphere before they are removed.

Common misconception

"If emissions stop, the pollutant disappears immediately." The burden decays with the residence time. A gas with a 50-year residence time still has about 37% of its burden remaining 50 years after emissions cease.

Worked example

Question: A lake contains 2.0 × 10⁹ m³ of water. Rivers bring in 4.0 × 10⁸ m³ per year and the same volume flows out. What is the water residence time? A soluble pollutant is added at a constant rate and removed only by outflow; how long until it reaches about 95% of its steady-state concentration?

Reasoning: τ = M / F = 2.0 × 10⁹ / 4.0 × 10⁸ = 5.0 years. The approach to steady state follows 1 − e^(−t/τ); 95% is reached when t ≈ 3τ.

Answer: The residence time is 5.0 years, and about 95% of steady state is reached after roughly 15 years.

Quick check

1. A reservoir holds 600 Tg of a gas and loses 60 Tg per year. What is the residence time of the gas? Answer: τ = 600 ÷ 60 = 10 years, assuming the reservoir is at steady state.

Exam focus

Always state the steady-state assumption when using τ = M/F. Be comfortable converting between Tg, Pg and moles, and remember that for first-order loss τ = 1/k. Examiners often ask how long it takes to approach a new steady state: use multiples of τ.

Advanced insight

In coupled multi-box systems a perturbation does not decay with a single time constant. Atmospheric CO₂ is a striking case: its individual molecules exchange with the ocean and biosphere within a few years, yet a pulse of extra CO₂ takes centuries to millennia to be fully removed. The molecular residence time and the adjustment time of a perturbation are different quantities.

Summary

Environmental chemistry treats the Earth as linked compartments exchanging matter through fluxes. Mass balance, dM/dt = F in − F out, governs every box. At steady state inflow equals outflow and the residence time τ = M/F. For first-order loss τ = 1/k, and a box approaches its new steady state over a few residence times. Residence time controls how widely a substance spreads.

Practice questions

1. Write the mass-balance equation for a single well-mixed box. Answer: dM/dt = F in − F out, where M is the burden and F the fluxes in and out. 2. A gas has a first-order loss rate constant of 0.25 per year. What is its residence time? Answer: τ = 1/k = 1/0.25 = 4 years. 3. A constant source of 30 Tg per year feeds a box with a residence time of 2 years. What is the steady-state burden? Answer: M ss = Sτ = 30 × 2 = 60 Tg. 4. Explain why a gas with a residence time of one day shows large spatial variations in concentration. Answer: It is removed before winds can mix it far from its sources, so concentrations are high near sources and low elsewhere. 5. After emissions of a first-order pollutant stop, what fraction of the original burden remains after two residence times? Answer: e^(−2) ≈ 0.14, so about 14% remains.