# Scaling up biology

Original paper

Physics is a discipline that attempts to develop a unifying, mathematical framework for understanding diverse phenomena. It connects things as different as planets orbiting the sun and a ball thrown through the air by showing that both these motions come from a single equation . Living things do not seem to obey such simplicity, but hidden beneath all the diversity and complexity of life are remarkably universal patterns called scaling laws. In a landmark 1997 paper by Geoffrey West, James Brown, and Brian Enquist, a simple explanation is given for how all organisms, from fleas to whales to trees, can be thought of as non-linearly scaled versions of each other.

A scaling law tells you how a property of an object, say the rate at which energy is consumed by an organism (its metabolic rate), changes with the object’s size. Just by looking at the data, many quantities scale as a power law of the mass, $A \propto M^{\alpha}$    (Eq. 1)

where ? is some number that, from the data, always seems to be a multiple of 1/4 . West, Brown, and Enquist build a theory showing how biology could have come up with this 1/4 power law, but in this article, I’m just going to focus on one specific example. I’m going to walk through the author’s arguments for how the metabolic rate, the rate at which an organism consumes energy, scales with an exponent of 3/4. They show that it all comes up from some basic assumptions about the networks that distribute nutrients to your body — your circulatory system .

These networks are assumed to have two characteristics . First, they are space-filling fractals. Fractals are shapes made of smaller, repeating versions of themselves no matter how far you zoom into it. However, our fractal blood vessels can’t get arbitrarily small, they have a “terminal unit”— the capillary. The second assumption about these networks is that all terminal units are the same size, regardless of organism size. With these two assumptions, the authors are able to derive the 3/4 power law for metabolic rate. Figure 1: Cartoon of a mammalian circulatory system on the left, which can be represented as a branching network model on the right. Adapted from Figure 1 of the original paper.

First, let’s build up a picture of what these networks look like. Figure 1 shows how the circulatory system can be thought of as a network structured into N levels, where each level k has $N_k$ tubes. At each level, a tube breaks into a number ( $m_k$) of smaller tubes. Each one of these tubes is idealized as a perfect cylinder with length $l_k$ and radius $r_k$, as shown in Figure 2. Figure 2: Illustration of the different parameters that each tube on the kth level of the network has. Adapted from Figure 1 of the original paper

How does blood move through this network? Well, the blood flow rate at each level of the network must be equal to the blood flow rate at every other level. Otherwise, you would have the equivalent of traffic jams in your arteries. You don’t want those. If the blood flow speed through one tube in the kth level is $u_k$, the blood flow rate through the entire kth level is $\dot{Q}_k = N_k \pi r_k^2 u_k = N_{cap} \pi r_{cap}^2 u_{cap} = \dot{Q}_{cap}$    (Eq. 2)

Your metabolic rate, B, depends on the flow rate through your capillaries, $\dot{Q}_{cap}$, so the authors assume that the two are proportional to each other: $B \propto \dot{Q}_{cap}$. Because all terminal units are the same size, the only variable left in Eq. 2 to relate to an animal’s mass is $N_{cap}$. Assuming that B scales like $B \propto M^{\alpha}$, and the authors predict $N_{cap} \propto M^{\alpha}$    (Eq. 3) Figure 3: Schematic of a branching point along the network, illustrating the definitions of the ratios and . In this case, .

To figure out the value of the exponent $\alpha$, the key is to get $N_{cap}$, which depends on the size of the organism, in terms of the capillary dimensions $r_c$ and $l_c$, which do not. To do this, the authors use relations derived using the self-similar geometry of the fractal network. When a tube breaks into smaller tubes, it does so with a ratio between the successive radii, $\beta_k = r_{k+1} / r_k$, and another ratio between the successive lengths, $\gamma_k = l_{k+1}/l_k$. This is illustrated in Figure 3. Because the network is fractal, the number of tubes each branch breaks into, $m_k$, the ratio of radii, $\beta_k$, and the ratio of lengths, $\gamma_k$, are all assumed to be constant for every k, $\beta_k = \beta, \; \gamma_k = \gamma, \; m_k = m \;\; \forall k$

Since, at every level, each branch breaks into m smaller branches, the total number of capillaries (i.e. the number of branches at level N) is $m^N$. Plugging this into Eq. 3, $\alpha = \frac{N \ln(m)}{\ln(M/M_0)}$    (Eq. 4)

Where $M_0^{\alpha}$ is the proportionality constant between $N_{cap}$ and $M^{\alpha}$. Remember, we’re trying to show that $\alpha = 3/4$.

Now that $N_{cap}$ has been rewritten in terms of network properties, the authors next turn their attention to  another quantity that scales with the organism size — its mass, M. To do this, the authors use the empirical fact that the total volume of blood, $V_b$, is proportional to the total mass of the organism, $V_b \propto M$. The total volume of blood is given by: $V_b = \sum_{k=0}^N V_k N_k = \sum_{k=0}^N \pi r_k^2 l_k m^k \propto \left( \gamma \beta^2 \right)^{-N} \propto M$    (Eq. 5)

In the above equation, the first proportionality sign (summing the series) requires a calculation that’s given here. The main idea of this calculation is that, because the ratios $r_{k+1} / r_k$ and $l_{k+1}/l_k$ are each constant, the sum in Eq. 5 can be turned into a geometric series which can be summed analytically. Plugging the final proportionality from Eq. 5 into Eq. 4, $\alpha = - \frac{\ln(m)}{\ln(\gamma \beta^2)}$    (Eq. 6)

To make further progress, we have to know something about $\gamma$ and $\beta$. Every tube of the network gives nutrients to a group of cells. As every good physicist does, the authors will assume that this group of cells has the volume of a sphere with a diameter equal to the length of the tube. The volumes serviced by each successive level are approximately equal to each other, $4/3 \pi (l_{k+1} / 2)^3 N_{k+1} \approx 4/3 \pi (l_k / 2)^3 N_k$. From this, they get an expression for $\gamma$: $\gamma_k^3 \equiv \left(\frac{l_{k+1}}{l_k}\right)^3 \approx \frac{N_k}{N_{k+1}} = \frac{1}{m}$    (Eq. 7)

which means $\gamma \approx m^{-1/3}$

Now the authors move on to $\beta$. Earlier, I argued that the flow rate has to be the same from one level of the network to the next to avoid “traffic jams” of blood. Since the tubes are assumed to be perfect cylinders, this boils down to the idea that the cross-sectional area of a parent tube being equal to the total cross-sectional area of its daughter tubes, $\pi r_k^2 = \pi r_{k+1}^2 m$. From this, the authors find an expression for $\beta$: $\beta_k^2 \equiv \left( \frac{r_{k+1}}{r_k} \right)^2 = \frac{1}{m}$     (Eq. 8)

Similar to the expression for $\gamma$, this means $\beta \approx m^{-1/2}$

Plugging in the expressions for $\gamma$ and $\beta$ in terms of m, we finally arrive at our desired result: $\alpha = - \frac{\ln(m)}{\ln(\gamma \beta^2)} = - \frac{\ln(m)}{\ln(m^{-1/3}(m^{-1/2})^2)} = 3/4$    (Eq. 9)

What West and his colleagues have done is use the fact that all organisms have to deliver nutrients to their individual parts to derive a general, universal scaling law. The authors go on to show that when you add a pump to the system, such as our heart, the analysis may get more complicated, but the ultimate result remains unchanged. All living things, regardless of size, seem to have arrived at the same solution for nutrient supply, building systems that are space-filling, fractal, and have the same size “terminal units”. Turns out we’re not so different after all. $F = Gm_1 m_2 / r^2$. ^

 For example:

• $\alpha = 3/4$ for cross section area of aortas of mammals, tree trunk sizes
• $\alpha = -1/4$ for cellular metabolic rate, heartbeat rate, population growth
• $\alpha = 1/4$ for time of blood circulation, life span, embryonic growth rate ^

 All the arguments hold for other distribution systems, such as our pulmonary system, plant vascular systems, and insect respiratory systems. ^

 There’s an additional assumption that the network is designed to minimize energy, but that won’t come into play in the part of the author’s arguments that I will be presenting here. ^

## 3 Replies to “Scaling up biology”

1. Tot says:

“With these two assumptions, the authors are able to derive the 3/4 power law for metabolic rate.”
At first, I thought you were about to show that there is a power law for metabolic rate, and that the exponent is 3/4, whereas (if I understood correctly) you show the latter from the former. Maybe emphasizing this would make the post clearer.

1. Danny Seara says:

Agreed! The power law is assumed to hold (as is substantiated by mounds and mounds of data), and finding the specific power of 3/4 is the point of this exercise. I’ll make that clearer.
2. Graham Scarr says: