View source for Age-structured population

The mathematical modeling of age-structured populations has been initiated by the pioneering work of A.J. Lotka and A.C. McKendrick in the years 1920–1940 (see [[#References|[a6]]], [[#References|[a7]]], [[#References|[a8]]]). For a long time, the interest in age-structure was restricted mainly to demography, but in recent years a new impulse to the research came from other fields and nowadays (1990s) age-structure plays a fundamental role in ecology, epidemiology, cell growth, etc.. As a matter of fact, age is one of the most natural and important parameters structuring a population, because many internal variables, at the level of the single individual, are strictly depending on the age, so that differing ages imply different reproduction and survival capacities and, also, different behaviours.

The basic linear model, usually named the Lotka–McKendrick model, describes the growth of a single population under the assumption that age is the only source of difference among individuals. It is the equivalent of Malthusian growth for a homogeneous population. The state of the population is described by its age-distribution <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a1201701.png" />, where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a1201702.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a1201703.png" /> denote age and time, respectively; namely

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a1201704.png" /></td> </tr></table>

is the number of individuals of the population that at time <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a1201705.png" /> have ages in the interval <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a1201706.png" />. The dynamics is based on the following [[Cauchy problem|Cauchy problem]]

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a1201707.png" /></td> <td valign="top" style="width:5%;text-align:right;">(a1)</td></tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a1201708.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a1201709.png" /> are the age-specific fertility and mortality rate, respectively. Thus,

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017010.png" /></td> </tr></table>

is the (total) birth rate at time <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017011.png" /> and

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017012.png" /></td> </tr></table>

is the survival probability, i.e. the probability that an individual will survive until age <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017013.png" />.

In problem (a1), the first equation describes the aging process, while the boundary condition (the second equation) gives the renewal of the population through the input of newborns. This latter condition is a non-local one and makes the problem related to integral equations of Volterra type (cf. also [[Volterra equation|Volterra equation]]). In fact, the solution of (a1) has the form

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017014.png" /></td> <td valign="top" style="width:5%;text-align:right;">(a2)</td></tr></table>

where the birth rate <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017015.png" /> is the solution of the Volterra equation of the first kind

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017016.png" /></td> <td valign="top" style="width:5%;text-align:right;">(a3)</td></tr></table>

and

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017017.png" /></td> <td valign="top" style="width:5%;text-align:right;">(a4)</td></tr></table>

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017018.png" /></td> </tr></table>

A fundamental result concerning the asymptotic behaviour of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017019.png" /> states that

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017020.png" /></td> <td valign="top" style="width:5%;text-align:right;">(a5)</td></tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017021.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017022.png" />. Here, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017023.png" /> is the only real root of the so-called Lotka characteristic equation

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017024.png" /></td> <td valign="top" style="width:5%;text-align:right;">(a6)</td></tr></table>

any other root being non-real and with real part strictly less than <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017025.png" />. The particular solution <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017026.png" />, which is the asymptotic limit of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017027.png" /> in the sense stated by (a5), is called a persistent solution and corresponds to a stable age-profile <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017028.png" />. The basic parameter <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017029.png" /> is called the intrinsic Malthusian parameter and is strictly related to the so-called net reproduction rate

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017030.png" /></td> </tr></table>

which represents the number of newborns that a single individuals produces in his lifetime; thus, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017031.png" /> (respectively, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017032.png" /> or <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017033.png" />) when <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017034.png" /> (respectively, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017035.png" /> or <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017036.png" />). The asymptotic behaviour stated in (a5) is also referred to as the ergodic property, according to which the population forgets its initial age distribution, eventually attaining the stable profile <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017037.png" />.

Non-linear models, taking into account logistic effects and other mechanisms of the inter-specific interaction may be conceived in a large variety of ways describing how the age distribution influences the vital rates. A fairly general and often-used model is the following version of a model first proposed by M.E. Gurtin and R.C. MacCamy in [[#References|[a3]]]

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017038.png" /></td> <td valign="top" style="width:5%;text-align:right;">(a7)</td></tr></table>

Here, fertility and mortality depend upon the population through a variable <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017039.png" /> which is a weighted selection of the distribution. Stationary states for problem (a7) have the form

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017040.png" /></td> <td valign="top" style="width:5%;text-align:right;">(a8)</td></tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017041.png" /> is a solution of the equation

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017042.png" /></td> <td valign="top" style="width:5%;text-align:right;">(a9)</td></tr></table>

This equation means that at the equilibrium size <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017043.png" /> the net reproduction rate of the population must be equal to <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017044.png" />. The stability analysis of stationary states involves the investigation of a transcendental characteristic equation. Depending upon its own structure, the model may show different behaviour, including existence of stable periodic solutions.

Significant specific examples of (a7) are the birth-logistic model

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017045.png" /></td> </tr></table>

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017046.png" /></td> </tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017047.png" /> is a decreasing function, and the cannibalism model

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017048.png" /></td> </tr></table>

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017049.png" /></td> </tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/a/a120/a120170/a12017050.png" /> is increasing.

Beyond single species models, age structure is considered in many multi-population systems. In the domain of demography, two-sex models are of special interest and in the context of ecology the different cases of the inter-specific interaction are studied.

A well-investigated field is epidemics modeling, where age structure is taken into account for the description of diseases whose transmission strongly depends on particular age classes; typical examples are, of course, childhood diseases, as opposed to sexual transmitted ones. The analysis of the models shows that age structure may be responsible for periodicity in the outbreak of the disease.

In the same context of epidemics, a different kind of age structure is also considered; in fact, the so-called class age, i.e. the time elapsed since an individual has become infected, plays a role when variable infectiveness has to be considered. The description of the HIV/AIDS epidemic, with its long incubation period, requires that both chronological age and class age be taken into account.

The mathematical methods involved in the study of age-structured populations include: direct methods, which are widely used, and specific techniques from the theory of Volterra integral and integro-differential equations; however, in recent years, a functional-analytic approach based on the theory of semi-groups of operators (cf. also [[Semi-group of operators|Semi-group of operators]]) has also provided a natural and powerful framework for the research in this field (see [[#References|[a2]]]).

====References====
<table><TR><TD valign="top">[a1]</TD> <TD valign="top">  A.J. Coale,   "The growth and structure of human populations. A mathematical investigation" , Princeton Univ. Press  (1972)</TD></TR><TR><TD valign="top">[a2]</TD> <TD valign="top">  G. Webb,   "Theory of nonlinear age-dependent population dynamics" , M. Dekker  (1985)</TD></TR><TR><TD valign="top">[a3]</TD> <TD valign="top">  M.E. Gurtin,   R.C. MacCamy,   "Nonlinear age-dependent population dynamics"  ''Arch. Rat. Mech. Anal.'' , '''54'''  (1974)  pp. 281–300</TD></TR><TR><TD valign="top">[a4]</TD> <TD valign="top">  J. Impagliazzo,   "Deterministic aspects of mathematical demography" , ''Biomathematics'' , '''13''' , Springer  (1985)</TD></TR><TR><TD valign="top">[a5]</TD> <TD valign="top">  N. Keyfitz,   "Introduction to the mathematics of populations with revisions" , Addison-Wesley  (1977)</TD></TR><TR><TD valign="top">[a6]</TD> <TD valign="top">  A.J. Lotka,   "On a integral equation in population analysis"  ''Ann. Math. Stat.'' , '''10'''  (1939)  pp. 1–25</TD></TR><TR><TD valign="top">[a7]</TD> <TD valign="top">  A.J. Lotka,   "The stability of the normal age distribution"  ''Proc. Nat. Acad. Sci.'' , '''8'''  (1922)  pp. 339–345</TD></TR><TR><TD valign="top">[a8]</TD> <TD valign="top">  A.C. McKendrick,   "Applications of mathematics to medical problems"  ''Proc. Edinburgh Math. Soc.'' , '''44'''  (1926)  pp. 98–130</TD></TR><TR><TD valign="top">[a9]</TD> <TD valign="top">  J.A.J. Metz,   O. Diekmann,   "The Dynamics of physiologically structured populations" , ''Lecture Notes Biomath.'' , '''68''' , Springer  (1986)</TD></TR><TR><TD valign="top">[a10]</TD> <TD valign="top">  J. Song,   J. Yu,   "Population system control" , Springer  (1980)</TD></TR></table>