applet-magic.comThayer WatkinsSilicon Valley & Tornado Alley USA |
---|

in Biological Distribution Systems |

Geoffrey B. West along with his collaborators James H. Brown and Brian J. Enquist in
the Biology Department of the University of New Mexico at Albuquerque developed a model of the distribution systems of plants and animals that
explains allometric scaling laws. A distribution system such as the circulatory system that
carries blood to the body cells in animals consists of a branching network of tubes starting
with the aorta at the heart which feeds into successively smaller tubes until the smallest
capillaries are reached. The model was published in *Science* Magazine with
the title, "A General Model for the Origin of Allometric Scaling Laws in Biology," Although it
was published in 1997 it is still not as well known as it should be. This material is an
explanation of the model. The original article is a bit hard to follow in crucial spots. However
even if there are some flaws of expostion in the article the model is sound and clearly the
right overall approach to modeling such biological systems. To make reference easier the
model will be referred to as the West Model.

The model treats a system such as the cariovascular circulatory system of animals as a branching network of smaller and smaller tubes that has some of the characteristics of a fractal but is not actually a fractal. The acknowledges three explicit assumptions but there are a number of implicit assumptions utilized in the development. The explicit assumptions are:

- The branching network is fractal-like and volume-filling.
- The minimal elements in the network are the same size for all creatures. This means that not only the radii the same but also length and other characteristics.
- The structure of the network minimizes fluid resistance.

In the course of the analysis a number of other assumptions were made. Some of those assumption made were rather artifical and could be replaced with equivalent more plausible assumptions. The implicit assumptions will be pointed out in the explanation of the model and a summary presented after the model has been covered.

Distribution systems the model applies to include the bronchial tubes in the lungs of animals and plant vascular systems as well as the blood circulatory systems of animals. They have a tree-like structure that terminates in some minimal size element. This is called fractal-like because a true fractal would have no terminal elements. Blood circulation system of animals branches out until the capillaries are reached then the collection system goes from capillaries to veins. The circulatory system is thus two fractal-like systems back-to-back. The exposition of the West model is concerned just with the distribution system but the model applies equally well to the collection system.

Let the branching network have K levels starting with level zero. In the West article the number
of levels was denoted as N but N was also used in reference to the number of elements in a
level. It is less confusing to let the number of levels be denoted by K (uppercase). Lower case
k is used as the index of the levels.
Let the length and radius
of the tubes in the k-th level be denoted as *l*_{k} and r_{k}, respectively.
The pressure drop over the length of the level k tubes will be denoted as Δp_{k},
and the cross-section average flow rate as u_{k}. At each level there is a branching
number for the number of smaller tubes each tube branches into. This is denoted as n_{k}.

The total number of tubes in the k-th level, N_{k}, is then given by the product of
the branching numbers of the previous levels; i.e.,

The total flow of fluid through the k-th level, Q_{k}, is
N_{k}πr<_{k}²u_{k}. Since
the fluid is conserved

The volume of fluid in a k-th level element, V_{k}, is given by

But West implicitly assumes that the ratio of radius to length is the same in all levels. This is not explicitly stated in the 1997 paper but is equivalent to what is stated. Let d be this common ratio. Then

At this point in the analysis West limits consideration to networks in which
the branching number is the same for all levels; i.e., n_{k}=n for all
k. West makes two additional assumptions. One, that the volume of fluid in one
level is the same as in the previous level. Two, that the cross-sectional area
of the branches of a tube is the same as that of the tube.

Thus, it is assumed that

for all k, which implies that

N

and hence

(

and therefore

γ =

From the assumption that the cross-section area of the n branches of a tube is equal to the area of that tube it follows that

and therefore

r

This common ratio of successive radii West calls ß.

These two assumptions are reasonable if one assumes the average fluid
velocity is constant from one level to the next. West makes these assumptions
and then *derives* the constancy of fluid velocity. It is more sensible to
make the constancy of the fluid velocity the assumption and derive the
area-preserving and volume-preserving assumptions from it.

The total volume of fluid in the distribution system is the sum of the volumes in each level; i.e.,

= πr

= πr

= πr

Note that πr_{0}²*l*_{0}=V_{0} so the above can
be expressed as

which for nγß²<1 and K large reduces to

V

West presumes that nγß²<1 but without supporting arguments. The parameters γ and ß are less than unity but n is greater than unity. In fact, West's assumption of area preservation between levels implies that nß²=1 and hence nγß² is really just γ and hence is necessarily less than unity.

Note that the volume of the terminal elements is given by

so

V

and, taking into account that nß²=1

V

West presumes that the total fluid volume needs to be proportional to body
mass. This is plausible in that each unit of body mass needs a definite quantity
of nutrients and oxygen. However changes in fluid velocity could compensate for
reduced fluid volume. Under West's assumption and the above equation it follows
that (γß²)^{-K} is proportional to body mass.

The number of terminal network elements N_{K} is given by

The metabolism rate corresponds to the total fluid flow rate, which under the previous assumptions is constant for all levels in the distribution network. The geometric characteristics of the terminal elements of the distribution network are assumed to be the same at all creature scales. This means that the total fluid flow in the terminal elements would be proportional to the number of terminal elements. Thus, if B represents the metabolic rate

M = c

where c

These relationship provide a link between B and M through the variable K; i.e.,

ln(M) = ln(c

It was previously found that γ is proportional to n^{-1/3}

and ß
to n^{-1/2} so (γß²) is proportional to
n^{-1/3-2(1/2)}=n^{-4/3}. This means that

Solving this last for K*ln(n) in terms of M and substituting the result into the equation for ln(B) gives B proportional to M to the (3/4) power.

The size characteristics of the terminal elements are scale invariant so the dependencies of the sizes of the other elements on scale are easily determined; i.e.,

Since ß = n

r

r

The total number of terminal elements N_{K} was previously found to
be proportional to the metabolic rate B, which in turn was found proportional to
M^{3/4}. Therefore r_{0} is proportional to
(M^{3/4})^{1/2}=M^{3/8}.

Likewise

and γ=n

and thus

B

(To be continued.)

**The energy minimization component of the model is dealt with
elsewhere. That part of the model is not a success but its failure does not detract from the
quality of the model.
**

applet-magic HOME PAGE OF Thayer
Watkins |