THE MORPHOLOGY OF GEOGRAPHICAL CHk~GE
W.R. Tobler
September 1975 WP-75-112
Working Papers are not intended for distribution outside of IIASA, and are solely for discussion and infor- mation purposes. The views expressed
are those of the author, and do not necessarily reflect those of IIASA.
The Morphology of Geographical Change W.R. Tobler
Geography is the study of the processes which result in changes of phenomena on the surface of the earth, and thus i t is useful to be able to describe change in an analytical manner. Here a single geographical variable is assumed to vary with time and this change is modelled in a number of simple ways. The objective is to see if this excursion is of use in the geographical study of growth poles. The topic is not one with which I am excessively familiar.
The most elementary case is to consider, in a rather mechanical fashion, a scalar variable and its changes through
time. This variable may appear in three guises, which how- ever may be used fairly interchangeably:
(I) as a function, G(x,y)
(II) as a finite set of real nu~oers, Gij (III) as a finite set of real numbers, Gk
In (I), x,y are considered as latitude and longitude coor- dinates (although no explicit recognition will be given to the roundness of the earth). G represents the value of the phenomenon of concern at the location x,y and exists at every location in the domain of interest. For example, G might represent the number of people residing within a circle of one kilometer radius centered at x,y; or i t might represent the average annual income of these people; or the percent of them which are unemployed, etc.
In (II), i and j are indices describing geographical locations when the data are assembled on a chessboard-like lattice. An interpretation is that G .. is a sample obtained
1J
from G(x,y) by multiplication with the brush function, the two-dimensional equivalent to Dirac's comb function. An
alternate interpretation is provided by grid coordinate data collection schemes as recently promulgated in Sweden and else- where.
The interpretation for (III) is similar, but the sub- script now refers to a bounded two-dimensional region of irregular shape~ a country or a census tract for example.
A conversion might be that
G(x,y)
=
Gk if and only if x,y is in region k.
Figure 1 shows such a function.
Figure 1. A piecewise continuous function z
=
f(x,y) ~ 1970 population density by states.The single symbol G, with or without subscripts, will refer to the state of the phenomenon at all locations at a single time period. The state of the phenomenon at a dif- ferent time will be distinguished by an asterisk~ G* for a later time period, G** for two periods later, G-* for an earlier time period, etc. Time will here be considered either continuous or in equal discrete steps.
- 3 -
Geographical Change Models
Regional policy is, in one elementary interpretation, concerned with the morphology of the function G(x,y); whether its variance is too large or too small; whether i t has too many peaks or too few; whether
J J
G(x,y) dxdy is too small or too large; and so on to include all of the static,descriptive, questions which one might ask of a social topography, as well as the normative desirata which one might envision. Of course, the concerns of regional policy are at least vector valued and not scalars, and should be more related to people than places. In any event, the dynamics can be viewed as an empirical question in which one asks how the morphology is changing. Let
be the temporal rate of change of the phenomenon of interest.
A number of specific models come to mind, e.g., dGdt(x,y)
=
R(x,y)where R(x,y) is a two-dimensional continuous random variable.
The change at each location is independent of that at all other locations, and is random. Presumably one needs to reject this model with one's data before going further, but I do not ever recall seeing such a test in the literature.
The next simplest model would appear to be that of constant regional growth:
dG(x,y) - constant
- c r t - -
Every place grows (or shrinks) by the same amount. If the measure G is in dollars this might reflect a static situation with inflation. Observationally, of course, the more usual case is that in which the rate of change depends on where one is located. We can now speculate on the form of this function.
One, still simple, model is that of proportional growth:
dGdt(x,y)
=
kG(x,y)where k is the same for all locations. Combining the three foregoing models would yield
dG(x,y)dt -- c + kG(x,y) + R(x,y)
This looks like a simple linear regression. A more compli- cated version, which however does not seem well specified, would be to allow spatial variations in both constants; that
is, c
=
C(x,y), k=
K(x,y).Continuing,
dGdt(x,y)
=
k(l - a) G(x,y)The alpha can be interpreted as interesting from a control
could be a constant or could describes logistic growth.
a target, and this might be theoretic viewpoint. Alpha vary spatially, a = A(x,y).
Alternate hypotheses might include
dGdt(x,y)
in which the rate of change is proportional to the spatial gradient, or
- 5 -
in which the temporal rate of change is proportional to the spatial rate of change of the gradient, or
d 2
G
=
k(a2G2 + a 2
G2) --2(x,y)
dt
ax ay
where the temporal acceleration of growth is to the spatial curvature of the phenomenon.
linear partial differential equations could special cases of
proportional These last two be considered
n dUG
L
£ -u(x,y) =u=l u dt
m v=l
L
with £u = Lu(x,y) and k v = Kv(x,y) to make for further generality.
The motivation for such models can be made clearer.
Consider a system of three adjacent spatial cells with populations G
l , G
2, and G 3:
1 2 3
Suppose further that the population of any cell is reduced by out-migration to neighbors and is increased by in-migration from neighbors. Let the proportion migrating from cell 1 to cell 2 be W21' with comparable interpretations for W23' W32' W12. If there are no other ways of changing the system, then
at some later time the population G
2
of cell 2 will bein which everything is accounted for. Now assruue (but only for expository convenience) spatial isotropy, namely,
w
2l= w
23= w
32=
W12=
W. As a consequenceThe change in population due to migration is then
G
2 -
G2=
W(Gl - 2G2 + G3), but this has occurred over some time period, say 6t. Consequently,G* - G
dG ~ 2 2 - W(G
at
~ 6t - 1 - 2G 2 + G3)Now consider the derivatives of a function f(x) at a finite set of points
Define x_I to be2 (f(x_l ) + f(x O))/2 and x+! to be2
(f(x
o)
+ f(xl ))/2, then the derivative at these points is approximated by- 7 -
where X
o -
x_I=
x+ I new function f' (x) =- X
o
= 6x. Thus we can consider thedf .
dx at the pOlnts x_~ and x+!:
d2 derivative of f' at x
o'
or ---2 .ax
f above values one findsis an approximation to the By substitution using the
Changing the labels to go back to the G's, we have
So that, approximately
dG dt
By analogous arguments in an orthogonal direction, a second spatial dimension can be added to obtain
which is one of the equations exhibited earlier, and i t no longer seems so mysterious and implausible. A comparable, but of course different, discussion could be provided for some of the other equations. The assumptions were, one recalls, that changes in the system occur only through migration, and spatial isotropy. The first situation is easily accounted for by adding terms for new entrants and departures, viz.
~~(x,y) ~ kG) + : ) ) + (B(X'Y) - 0(X'Y)) G(x,y)
The second assumption is· most easily removed by reverting to the equation
which has now to be completed by adding the births and deaths
Adding more spatial cells to this three cell scheme is clearly not difficult in an accounting sense. More interesting perhaps is to interpret k
=
K(x,y) in- 9 -
dG = k(a 2
G + a 2
G) + (8 + 8) G dt ax 2 ay2
This must surely mean that the earlier constant W varies from place to place. If one were to write out the entire system of equations, one accounting equation for each place,
G*2
G* - G
3 3
H12
then the isotropic assumption would yield a constant W matrix.
An intermediate form is obtained if the matrix is a band- matrix with entries some distance from the diagonal equal to
zero. This can be used to imply that there is no migration between distant cells. A convenient isotropic assumption is then that
G;. = 1)
p=+n p=-n
L
q=+n q=-n
L
W G. . + (8 .. - 8 .. ) G ..
pq l+p, )+q 1) 1) 1)
where the summation is over the n2 cells which make up the influence field around each cell.
A geographical growth pole is supposed to be a place from which growth effects spread. Re-examine the diffusion equation
dG
dt + (8 - 8) G
The first term on the right describes the spread effect,
migration to neighboring places, and the second term describes the "natural increase" and decrease. Suppose we remove the spread effects, subtracting from both sides
and now rewrite this in finite difference form as
G'!'.
1)
- L L
p q
W G. .
pq l+p,)+q = (S .. - 0 .. ) G
1) 1)
(wherein Wij
=
1 - W21 - W23 of the earlier equations). That is, if we separate out the spread effects, we can estimate the spatial growth multiplier in the "natural increase" term.
It, of course, can also be less than one, an inhibiting effect may exist. There seems to be hope that this can be done
empirically using data from two time periods. If one writes out the wave equation
in finite difference form then the second derivative term
on the left hand side means that we need data from an addition time period; i.e., the model is one with temporal lags, just as the presence of spatial derivatives calls our attention to the spatial lags.
An alternate type of model would be the purely spatial function
n n
G*(x,y) - G(x,y)
= L L
i=O j=O
i j a .. x y
1)
- 11 -
using an algebraic polynomial, or some other two-dimensional form (spherical harmonics, for example). Here the low-order terms (linear, quadratic) would indicate a spatial trend in the growth. A geographical trend, from East to West say, might then indicate the effects of an external influence.
Other than such effects the model is more descriptive than interesting.
Casetti and Semple have proposed an alternate model as follows
G* - G ~
k k
n
i~O
L
a. exp (-B.r.
k)
1 1 1
where r
ik
~
[(Xi - xk)2 - (Yi - Yk)2J! is the separation of the location i from the location k. Basically, the modelsays that there are n + 1 points around which growth declines.
The magnitude a. and the rate of decline
B.
are different for1 1
each of the several centers. The rate of decline can be made to vary in different directions around each of the centers by putting a separate Fourier series in the exponents, as has been done by Medvedkov. One of the difficulties of this model is its non-linearity, although an iterative computer program is available.
Suppose that a government decides to allocate funds according to F(x,y). Then one might postulate that
dG _dt - kF(x,y)
Suppose that the allocation is proportional to the population P(x,y), then
dG _dt - kP(x,y)
and if the population is the phenomenon of concern
dP kP(x,y) dt
one of the models encountered earlier. We must always
subtract the source of the funds too, the taxes T(x,y), and
f I
F(x,y) dxdy= J J
T(x,y) dxdyAssuming negligible costs for the bureaucracy. Thus
would be one model. What is now needed is a more process- oriented morphology.