From Water Flow to Pollutant Load
A water-quality measurement such as \(3\) milligrams per litre describes how concentrated a pollutant is. It does not, by itself, tell us how much pollutant enters a reservoir, treatment plant, or other system. The amount entering also depends on how much water flows and for how long.
The calculation has three distinct quantities: concentration, mass rate, and accumulated load. Keeping them separate prevents a common modeling error: treating a concentration value as though it were already a mass or a mass rate.
The Mass Rate from Flow and Concentration
Let \(Q(t)\) be the water flow rate at time \(t\), measured as volume per unit time. Let \(C(t)\) be the pollutant concentration in that flowing water, measured as pollutant mass per unit volume. If the concentration represents the water passing through the measurement cross-section at time \(t\), then multiplying the two quantities gives the pollutant mass entering per unit time.
For example, if \(Q\) is in litres per minute and \(C\) is in milligrams per litre, then \(R\) is in milligrams per minute: the litre units cancel. If flow is in cubic metres per minute while concentration is in milligrams per litre, the units do not cancel until the flow is converted, since \(1\text{ m}^3=1000\text{ L}\).
Let \(M(t)\) denote the total pollutant mass that has entered between a chosen starting time \(t_0\) and time \(t\). When the mass rate is continuous, the entering load over \([t_0,t_1]\) is the integral of that rate:
The units provide a useful check: mass per time multiplied by time gives mass. The load is not concentration, and it is not the mass rate at the endpoint. It is the total mass accumulated over the stated interval.
Why the Product and Integral Give the Load
Consider a short time interval of length \(\Delta t\) near time \(t\). During that interval, approximately \(Q(t)\Delta t\) units of water pass the measurement point. At concentration \(C(t)\), the pollutant mass in that water is approximately \(C(t)Q(t)\Delta t\). Thus the mass entering during the short interval is approximately \(R(t)\Delta t\), where \(R(t)=Q(t)C(t)\).
Divide a longer interval into short subintervals. The total entering mass is approximated by adding the pollutant mass from each subinterval. As the subinterval lengths shrink, these sums approach the definite integral of the mass rate. This gives the accumulation formula when the rate is integrable; continuity of \(Q\) and \(C\) on the interval is a standard sufficient condition.
This reasoning also explains why the two time-varying quantities must be multiplied before integration. In general, integrating flow and concentration separately and then multiplying the results does not give the accumulated pollutant mass.
Worked Example: Flow and Concentration Both Vary
Worked Example: Flow and Concentration Both Vary
Suppose the flow into a system over the first \(10\) minutes is \(Q(t)=100+10t\) litres per minute, and the concentration is \(C(t)=2+0.2t\) milligrams per litre, where \(t\) is measured in minutes. Find the total pollutant load entering in those \(10\) minutes.
First calculate the mass rate. The litres cancel, so the result is in milligrams per minute:
Integrate this mass rate from \(t=0\) to \(t=10\):
The exact load is \(\frac{14000}{3}\) milligrams, or approximately \(4.67\) grams. The antiderivative check is \(\frac{d}{dt}(200t+20t^2+\frac{2}{3}t^3)=200+40t+2t^2=R(t)\), confirming the integrated mass-rate expression.
Unit Conversion Must Come Before Multiplication
The product rule for the physical quantities is simple, but incompatible units can hide a factor of \(1000\). Before multiplying, express flow and concentration using the same volume unit. The next example makes the conversion explicit.
Worked Example: Converting Cubic Metres to Litres
Suppose \(t\) is in minutes, \(Q(t)=0.4+0.02t\) cubic metres per minute, and \(C(t)=5+t\) milligrams per litre, for \(0\leq t\leq10\). Find the entering load.
Convert the flow to litres per minute using \(1\text{ m}^3=1000\text{ L}\):
Now multiply by the concentration. The result is a mass rate in milligrams per minute:
Integrate over the ten-minute interval:
The entering load is \(\frac{155000}{3}\) milligrams, approximately \(51.67\) grams. Differentiating the antiderivative gives \(2000+500t+20t^2\), exactly the mass rate found from the converted flow and concentration. Without the factor of \(1000\), the numerical result would be too small by that factor.
Operating Regimes and Piecewise Loads
Flow or concentration may change abruptly when a system changes operating regime. If the mass rate is defined by different formulas on different intervals, calculate the load on each interval and add the results. A single endpoint rate cannot replace this interval-by-interval accumulation.
Worked Example: Load Across Two Operating Regimes
A system receives water for five minutes. During the first two minutes, the flow is \(8\) litres per second and the concentration is \(1.5\) milligrams per litre. During the next three minutes, the flow is \(5\) litres per second and the concentration is \(4\) milligrams per litre. Find the total entering load. Use seconds for time.
In the first regime, the mass rate is \(8(1.5)=12\) milligrams per second. Two minutes is \(120\) seconds, so the load in this regime is \(12(120)=1440\) milligrams. In the second regime, the mass rate is \(5(4)=20\) milligrams per second. Three minutes is \(180\) seconds, so its load is \(20(180)=3600\) milligrams. Therefore:
One check is to write the cumulative load as a piecewise function, with \(t\) in seconds from the start of the interval:
At the changeover, both expressions give \(1440\) milligrams, so the cumulative load is continuous. At \(t=300\), the second expression gives \(1440+20(180)=5040\) milligrams, matching the sum of the two regime loads. On each interval, differentiating \(M(t)\) gives the corresponding mass rate: \(12\) or \(20\) milligrams per second.
What the Load Does—and Does Not—Measure
The integral \(\int Q(t)C(t)\,dt\) measures pollutant mass entering through the specified flow over the chosen interval. It is a gross entering load. It does not, by itself, determine how much pollutant remains in the system at the end: that may also depend on pollutant leaving, chemical or biological changes, and other sources or sinks. Those effects belong in a broader conservation balance.
The formula also depends on using a concentration representative of the water flow being measured. If the flow and concentration vary across a channel or pipe, a single concentration value must represent the flowing water appropriately for the product \(QC\) to represent the mass rate. In an application, the meaning and measurement basis of the reported concentration matter.
A Reliable Calculation Sequence
For a pollutant-load model, keep the physical interpretation visible at every step:
Specify the start and end times, and make sure the time unit used in the integral matches the rate units.
Convert flow or concentration so that their volume units cancel when multiplied.
Compute \(R(t)=Q(t)C(t)\), including any required conversion factor.
Integrate the mass rate over the interval, or integrate separately on each operating regime and add the loads.
The result should have units of pollutant mass and should describe entering load, not concentration or endpoint mass rate.
The central modeling idea is that water carries pollutant mass. Flow tells how much water passes per unit time; concentration tells how much pollutant is in each unit of water. Their product gives the pollutant mass rate, and its accumulation over time gives the total entering load.
Check Your Understanding
Use units and physical meaning, not just algebra, to answer these questions.
- If flow is measured in litres per minute and concentration in milligrams per litre, what are the units of their product?
- Why must the flow and concentration be multiplied before integrating when both vary with time?
- A flow is reported in cubic metres per minute and a concentration in milligrams per litre. What conversion is needed before forming the mass rate?
- If a system operates under two constant flow-concentration regimes, how can you calculate the total entering load?
- Does the entering load alone determine the pollutant mass remaining in a system at the end of the interval? Explain briefly.
