Module 4: Allometric Scaling

Across 22 orders of magnitude in body mass โ€” from a mycoplasma at\(10^{-16}\) kg to a blue whale at \(10^{5}\) kg โ€” metabolic rate scales with mass to an almost universal exponent of 3/4. Heart rates scale with exponent\(-1/4\); lifespans and gestation times scale with \(+1/4\). These quarter-power laws, empirical since Kleiber (1932) and derived from first principles by West, Brown & Enquist (1997), are consequences of the fractal branching transport networks introduced in Module 3. This module derives them and shows how they set hard limits on maximum body size and tree height.

1. Kleiber's Law

In 1932, Max Kleiber measured oxygen consumption in mammals from mice to steers and found:

\[ B = B_0 \, M^{3/4} \]

\(B\) = basal metabolic rate (W); \(M\) = body mass (kg).

Rubner (1883) had proposed \(B \propto M^{2/3}\) based on the idea that heat loss is proportional to surface area. But Kleiber's data clearly favoured 3/4. Hemmingsen (1960) extended the relation to unicellular organisms and plants, spanning 22 orders of magnitude of mass. The exponent 3/4 is the most robust single exponent in all of biology.

1.1 Why 3/4 and Not 2/3?

The surface-law argument fails because heat dissipation is not the binding constraint; resource delivery through a fractal vascular network is. The WBE derivation (Module 3 Section 4):

  1. Space-filling network: \(L_k \propto n^{-k/3}\).
  2. Minimum dissipation (area-preserving for pulsatile flow): \(r_k \propto n^{-k/3}\).
  3. Invariant terminal units: \(N_{\text{cap}} = n^N\) capillaries per body.

Body mass \(M \propto V_{\text{blood}} \propto N_{\text{cap}}^{4/3}\); basal rate\(B \propto N_{\text{cap}}\). Eliminating\(N_{\text{cap}}\) gives \(B \propto M^{3/4}\).

2. Quarter-Power Scaling Laws

Once metabolic rate scales as \(M^{3/4}\), many other biological quantities derive automatically. Here are the main ones.

2.1 Heart Rate

Cardiac output \(\dot Q\) (litres/min) is proportional to metabolic rate:\(\dot Q \propto M^{3/4}\). Stroke volume \(V_s \propto M\) (ventricle volume scales with body size). Therefore heart rate:

\[ f_{\text{heart}} = \frac{\dot Q}{V_s} \propto \frac{M^{3/4}}{M^{1}} = M^{-1/4} \]

Mouse: ~600 bpm; human: ~70 bpm; elephant: ~30 bpm; blue whale: ~8 bpm.

2.2 Lifespan

Biological time (development, gestation, lifespan) slows with body size as\(t_{\text{bio}} \propto M^{1/4}\). Lifespan, gestation time, time to maturity, and the period of the circadian-but-longer-scale dampings all follow this pattern.

2.3 The One-Billion-Heartbeat Rule

Total heartbeats in a lifetime:

\[ N_{\text{beats}} = f_{\text{heart}} \cdot t_{\text{life}} \propto M^{-1/4} \cdot M^{+1/4} = M^0 \]

Independent of mass โ€” approximately \(10^9\) beats for most mammals.

Remarkably, a mouse and an elephant have similar total heartbeat counts. Humans, thanks to medicine, exceed this baseline (\(\approx 2.5\text{--}3\times 10^9\)) but the pattern persists across species. The same holds for breaths.

2.4 Other Exponents

Aorta radius

M^(3/8)

Average vessel length

M^(1/4)

Blood volume

M^1 (isometric)

Number of capillaries

M^(3/4)

Respiratory rate

M^(-1/4)

Muscle power per mass

M^(-1/4)

Home range area

M^(7/8) - M^1

Brain mass (mammals)

M^(3/4)

3. Derivation from WBE

Let the network have \(N\) branching levels; each level has a branching ratio\(n\) (usually 2). Denote level-\(k\) radius and length by\(r_k, L_k\).

3.1 Space-Filling

Each terminal unit (capillary) services a volume of tissue \(v_N\). Total volume is\(V = N_{\text{cap}} v_N\) and \(v_k \propto L_k^3\) because each branch fills its own cube:

\[ \frac{L_{k+1}}{L_k} = \gamma = n^{-1/3} \]

3.2 Area-Preserving Pulsatile Flow

For large arteries where blood is pulsatile, minimising wave reflection (Womersley impedance matching) requires \(\sum r_{k+1}^2 = r_k^2\), i.e.\(n r_{k+1}^2 = r_k^2\), so:

\[ \frac{r_{k+1}}{r_k} = \beta = n^{-1/2} \]

(In the small viscous regime like capillaries, Murray's cubic law gives \(\beta = n^{-1/3}\).)

3.3 Blood Volume Scales as M

Total blood volume:

\[ V_b = \sum_{k=0}^{N} n^k \pi r_k^2 L_k = \pi r_0^2 L_0 \sum_{k=0}^{N} (n \beta^2 \gamma)^k \]

With \(\beta^2 = n^{-1}\) and \(\gamma = n^{-1/3}\) we have\(n\beta^2\gamma = n^{-1/3}\). Summing the geometric series:

\[ V_b \approx \pi r_0^2 L_0 \cdot \frac{n^{-N/3}-1}{n^{-1/3}-1} \cdot n^{N/3} \cdot n^{-N/3} \sim n^{-N/3} \text{ terms dominate.} \]

More carefully, dominated by last level: \(V_b \propto n^N \pi r_N^2 L_N\). Given\(r_N, L_N\) fixed by invariant terminals, \(V_b \propto N_{\text{cap}}\). But also \(r_0 = r_N \beta^{-N}\), \(L_0 = L_N \gamma^{-N}\), so\(r_0^2 L_0 = r_N^2 L_N n^{N(1+1/3)} = r_N^2 L_N n^{4N/3}\). Combined with\(V_b \propto r_0^2 L_0\) (volume dominated by aorta segment):

\[ V_b \propto n^{4N/3} = N_{\text{cap}}^{4/3} \]

Since mass \(M \propto V_b\) (isometry) and \(B \propto N_{\text{cap}}\)(each capillary supplies a fixed metabolic unit):

\[ \boxed{\; B \propto M^{3/4} \;} \]

4. Maximum Body Size: Galileo's Problem

In Discorsi (1638), Galileo reasoned that an animal scaled up by a factor\(k\) in linear size has volume (and weight) \(\propto k^3\)but cross-sectional bone area \(\propto k^2\); bone stress therefore grows\(\propto k\). A simple upper bound:

\[ \sigma \approx \frac{\rho g L^3}{L^2} = \rho g L \]

For cortical bone yield stress \(\sigma_{\max} \approx 200\text{ MPa}\) and\(\rho \approx 1000\,\text{kg/m}^3\), the limiting length is:

\[ L_{\max} = \frac{\sigma_{\max}}{\rho g} = \frac{2\times 10^8}{10^4} \approx 2\times 10^4\,\text{m} \]

An absurdly large upper bound because we assumed geometric similarity. In reality, land animals adopt elastic similarity (McMahon 1973): \(r \propto L^{3/2}\), giving \(\sigma \propto L^{1/2}\) โ€” a much slower rise. The largest known land animals (sauropods ~90 tonnes) are near the practical bone-locomotor limit, while aquatic animals like blue whales (150 t) escape bone constraints thanks to buoyancy.

4.1 Circulatory Limits

WBE gives metabolic rate \(B \propto M^{3/4}\), but mass-specific rate\(B/M \propto M^{-1/4}\). Oxygen uptake rate drops with size, yet so does the cell-specific demand. The real limit arises from capillary density: \(\rho_{\text{cap}} \propto M^{-1/4}\). Below a critical density, diffusion cannot meet demand, setting a theoretical maximum near\(10^5\) kg โ€” close to the blue whale.

5. Tree Height: The Xylem Cavitation Limit

Plants have no heart; they draw water upward via the cohesion-tension mechanism. The water column in xylem is under tension, and the tensile strength of liquid water (before cavitation and embolism) is the hard limit on tree height.

5.1 Hydrostatic Component

Lifting water of density \(\rho_w = 1000\,\text{kg/m}^3\) through height \(H\)requires a tension:

\[ P_{\text{hydro}} = -\rho_w g H \]

A 100 m tree requires ~0.98 MPa of tension just against gravity.

5.2 Frictional Losses

The Hagen-Poiseuille equation gives the frictional gradient:

\[ \left|\frac{dP}{dz}\right|_{\text{fric}} = \frac{8\mu Q}{\pi r^4} \approx 0.01\text{--}0.02\,\text{MPa/m} \]

For a 100 m redwood, friction adds another ~1-2 MPa, yielding total tension\(|P| \approx 2\text{--}3\,\text{MPa}\).

5.3 The ~130 m Limit

Koch et al. (2004) measured tension at the tops of the tallest Sequoia sempervirens(Hyperion, 115.9 m). They found tension limits near 2 MPa โ€” the onset of cavitation for coniferous xylem. Extrapolating models gives a theoretical upper bound of\(H_{\max} \approx 122\text{--}130\) m, consistent with the tallest recorded tree, Eucalyptus regnans at ~133 m (historical records).

See the cross-course/tree-biophysicsand/ecological-biochemistryfor in-depth treatment of xylem vulnerability curves, stomatal regulation, and forest biogeography.

6. The Kleiber Plot

Kleiber's Law: B = B_0 M^{3/4} Across 22 Orders of Magnitude10^-1610^-1210^-810^-410^010^410^-1610^-1410^-1210^-1010^-810^-610^-410^-210^010^210^410^610^8body mass M (kg, log scale)metabolic rate B (W, log scale)B ~ M^(3/4)MycoplasmaE. coliYeastAmoebaParameciumMiteShrimpMouseRatCatDogHumanHorseElephantBlue whaleObserved (Hemmingsen compilation)WBE prediction: alpha = 3/4(surface law alpha = 2/3 would be too flat)

7. Allometric Limits Visualisation

Limits on Life: Size Ceilings Across Kingdoms0 m10 m30 m100 m200 mcharacteristic length (m)Mouse0.1 m, 20 gHuman1.7 m, 70 kgElephant4 m, 4 tBlue whale30 m, 150 tSauropod30 m, 90 tRedwood116 m (Hyperion)~130 m: theoretical tree height limit (xylem cavitation)~60 m: tallest living animal (aquatic), ~15 m: tallest land (extinct sauropod)typical large mammal scaleSize ceilings are set by distinct mechanisms:Land animals: bone stress + muscle power (~100 tonnes)Aquatic animals: prey density + circulatory limitsTrees: xylem tensile strength & cavitationInsects: tracheal diffusion & oxygen (< 1 cm)

8. Simulation 1 โ€” Fitting Kleiber's Law

Fit \(B = B_0 M^\alpha\) to Hemmingsen's compilation spanning mycoplasma to blue whale. The empirical slope consistently lies between 0.72 and 0.78, consistent with the WBE theoretical prediction of 0.75 and inconsistent with the surface law 0.67.

Python
script.py82 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

9. Simulation 2 โ€” Quarter-Power Traits

Eight biological traits plotted on log-log axes against body mass. Heart rate, breath rate, and muscle power per mass all have slope \(-1/4\); lifespan and gestation have slope\(+1/4\). The product heart-rate ร— lifespan is nearly constant, giving the famous one-billion-beats rule.

Python
script.py72 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

10. Simulation 3 โ€” Maximum Size Constraints

Three constraints plotted: Galileo's bone stress (geometric vs elastic similarity), xylem cavitation for trees, and circulatory power for animals. Their intersections with material limits set the observed maximum sizes.

Python
script.py100 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Module Summary

Kleiber's law

B = B_0 M^(3/4) across 22 orders of magnitude.

WBE derivation

Space-filling + minimum dissipation + invariant terminals -> 3/4.

Heart rate

f_heart ~ M^(-1/4); stroke volume ~ M.

Lifespan

t_life ~ M^(+1/4); biological time slows with size.

1 billion beats

f_heart * t_life = M^0 constant across mammals.

Galileo's limit

sigma = rho g L; elastic similarity r ~ L^(3/2) -> sauropod-scale cap.

Tree height limit

~130 m from xylem cavitation at -2 MPa.

Cross-links

tree-biophysics, ecological-biochemistry courses.

References

  1. Kleiber, M. (1932). Body size and metabolism. Hilgardia, 6, 315โ€“353.
  2. Hemmingsen, A. M. (1960). Energy metabolism as related to body size and respiratory surfaces. Reports of the Steno Memorial Hospital, 9(2), 1โ€“110.
  3. West, G. B., Brown, J. H., & Enquist, B. J. (1997). A general model for the origin of allometric scaling laws in biology. Science, 276, 122โ€“126.
  4. Schmidt-Nielsen, K. (1984). Scaling: Why is Animal Size So Important? Cambridge University Press.
  5. Calder, W. A. (1984). Size, Function, and Life History. Harvard.
  6. McMahon, T. A. (1973). Size and shape in biology. Science, 179, 1201โ€“1204.
  7. Koch, G. W. et al. (2004). The limits to tree height. Nature, 428, 851โ€“854.
  8. West, G. B. (2017). Scale: The Universal Laws of Growth, Innovation, Sustainability, and the Pace of Life. Penguin.
  9. Lindstedt, S. L. & Calder, W. A. (1981). Body size, physiological time, and longevity of homeothermic animals. Q. Rev. Biol., 56, 1โ€“16.
  10. Savage, V. M. et al. (2004). The predominance of quarter-power scaling in biology. Functional Ecology, 18, 257โ€“282.