Population

Crude death rate (CDR)

From population.bas:

"There are three numbers that capture fertility:  TFR, CBR, the fertility distribution. There are three numbers that caputure mortality:  LIFEXP, CDR, mortality distribution. They will not be completely compatible.  The distributions must be used to handle cohort dynamics.  But they may not be available for all countries (e.g. Taiwan, Palestine), whereas the others are.  For many years we used CBR and CDR to adjust the distributions and then calculated TFR and LIFEXP from them.  Over time, TFR and LIFEXP have become the key variables that dynamically drive the population model, and good data are now available. We have therefore switched over to them as the bases for initial conditions of the model, calculating CBR and CDR from them, giving us also POPR." <header><hgroup>

Demography/Population

Overview

The most recent and complete demographic model documentation is available on Pardee's website. Although the text in this interactive system is, for some IFs models, often significantly out of date, you may still find the basic description useful to you.

The population sub-model of IFs uses the cohort component analysis approach of many population models, including the studies done by the United Nations (United Nations, 1956 and 1977). The structure of the IFs population model drew initially on the World Integrated Model (WIM) or the second generation Mesarovic-Pestel Model (Hughes, 1980), but has changed much over time.  In particular, José Solórzano and Randall Kuhn have made many contributions to its development.

The approach relies upon age, fertility, and mortality distributions for each country/region with 22 cohorts: one for infants, 20 of five-year size, and one for all individuals of age 100 or older. A major advantage of five-year cohorts is that data sources generally present demographic data in that form. Ideally, however, the cohort size should correspond to the model time step so as to avoid "numerical diffusion," the propagation of change from a five-year cohort to an adjoining cohort in a single year. To prevent such numerical diffusion, IFs actually runs an age distribution with 100 single-year cohorts and advances that over time, collapsing to 22 cohorts only for the calculations of births and deaths.

Because extensions of life expectancy are occurring steadily and there is at least the possibility of substantial breakthroughs, the IFs project has also created the option of extending the number of cohorts from 22 up to as many as 42 (allowing separate representation of age categories up to 200+).  The capability is normally turned off, but instructions for turning on extended aging can be found here. <header><hgroup>

Structure and Agent System: Demographic

System/Subsystem	Demographic
Organizing Structure	Cohort-component
Stocks	Population by age-sex
Flows	Birth, death, migration
Key Aggregate  Relationships  (illustrative, not comprehensive)	Life expectancy (from health model)
Key Agent-Class Behavioral  Relationships (illustrative, not comprehensive)	Household fertility and migration

Humans as individuals within households interact in larger demographic systems or structures. The computer model should represent the behavior of such households, such as decisions to have children or to emigrate. And it should represent the larger demographic structures that incorporate the decisions of millions of such households. A typical approach to representing such demographic systems is through age-sex cohort distributions (see the figure below showing an example from the model). IFs also uses fertility and mortality distributions by age and sex and tracks migration across countries.

Demographers have widely accepted the representation of demographic systems and the development of demographic models with cohort-component structures. In fact, the United Nations, the U.S. Census Bureau, and the International Institute for Applied Systems Analysis (IIASA), pre-eminent demographic forecasting institutions, all use cohort-component modeling (O’Neill and Balk 2001)

Dominant Relations: Population

 The dominant population (POP) equation is a simple addition of births (BIRTHS) at the bottom of the cohort distribution, subtraction of deaths (DEATHS) from each population cohort, and advance of people to the next cohort over time.

The following key dynamics are directly linked to the Dominant Relations:

Births are primarily a function of the total fertility rate (TFR), which in the longer term responds especially to education level of the adult population. The model user has direct control over TFR with a multiplier (tfrm ), but also much control for low fertility countries with a parameter specifying long-term stabilization level and lower boundary for fertility (tfrmin ). There is also a secular trend reduction in fertility (controlled by ttfrr ).

Deaths are primarily a function of life expectancy (LIFEXP), itself computed within the IFs health model where, like fertility, it responds in the long run to adult education and also to GDP per capita and technology change. The model user has direct control over all deaths with a mortality multiplier (mortm ) and over those specific to a cause of health with an alternative multiplier (hlmortm ). There is also a secular trend reduction in mortality (controlled by tmortr ).

The larger demographic model in combination with the health model provides representation of and control over migration; the fertility impact of infant mortality and contraception use rates; and the mortality impact of many factors including undernutrition, smoking rates, and indoor air pollution from open burning of solid fuels.

Demographic Flow Charts

Overview

The demographic model represents the population of each geographic unit in terms of 22 cohorts (infants, five-year intervals up to age 99, and those aged 100 and older), separately for females and males. An age distribution records the population in each cohort and sex category. The sum across all cohorts in the age distribution and both sexes is the total population. A fertility distribution determines births, which are added to the bottom of the age distribution, while a mortality distribution determines deaths, which are subtracted from the appropriate cohort of the age distribution.

Those who might like to turn on the extension of age-cohort representation, to as many as 42, can do so by making changes in the IFsInit table of the IFsInit.mdb file.  Specifically, the  NCohorts field can be changed to as many cohorts as 42 and the NAges field can be changed up to 200.  Registering these changes requires a rebuild of the Base Case (see documentation of Extended Features).

The population model is central to many broader dynamics of IFs. Two key feedback loops drive its own dynamics. The first is a positive feedback loop around fertility, linking population and births (causing population to drive exponentially upward if nothing else changes), while the second is a negative loop around mortality, linking population and deaths (causing population to decline). This second loop actually runs through the health model of IFs where deaths are computed (switching the control parameter hlmodelsw from 1 to 0 would, however, cause the model to revert to an earlier formulation in which life expectancy was computed as function of GDP per capita and controlled the death rate and deaths; it would turn off the health model's impact).  A Malthusian variation of the negative feedback loop involving deaths may be of interest to those who believe that food supplies do or will play an important role in population dynamics (as they clearly do in countries with very low nutritional levels) by raising mortality rates, especially of children. See the topic on nutrition.  Whether population rises or falls depends on the relative strength of those two loops.

The easiest and most often used scenario handles for the population model are a multiplier on the total fertility rate (the number of children borne by an average woman in a lifetime), namely tfrm , a multiplier on the total mortality rate, mortm , and a multiplier on mortality by cause, hlmortm .

A large number of indicators are calculated in IFs from the age distribution:

• the median age of the country's population (POPMEDAGE)
• population aged 15 to 65 (POP15TO65)
• population above age 65 (POPGT65)
• population below age 15 (POPLE15)
• population pre working age (POPPREWORK), controlled by the parameter specifying the work starting age (workageentry )
• population post working age (POPRETIRED), controlled by the parameter specifying the retirement age (workageretire )
• population within the working years (POPWORKING)
• the potential support ratio, or the population from 15 to 64 over that above 65 (POTSUPRAT)
• an indicator of the youth bulge or the population from 15 to 29 as a portion of that 15 and above (YTHBULGE)

In addition there are a number of indicators calculated from the size of country populations:

• the growth rate of population (POPR)
• the world population (WPOP)
• growth in the world population (WPOPR)

More description is available on the dynamics around fertility and mortality as well as several specialized topics on topics such as nutrition levels and migration.

The central indicator of fertility is the total fertility rate (TFR), the number of children that the average woman will bear throughout her life. Fertility generally decreases in the long run as deep or distal drivers such as GDP per capita (from the economic model) or formal years of education of adults (EDYRSAG15 from the education model) increase; our own analysis suggested the use of education years, a result that Angeles (2010) reinforced.

In addition there are more proximate drivers, some of which can change more rapidly than GDP per capita or education levels and thereby affect fertility rate.  We represent two, namely infant mortality rates (INFMOR) and contraception usage (CONTRUSE).  The health model of IFs determines infant mortality rates.  Sudden change in those do not, however, immediately affect fertility rates and we smooth changes in rates so as to introduce an approximately 10-year lag, consistent again with the findings of Angeles (2010). Based on cross-sectional analysis the population model of IFs forecasts the percentage of population using modern contraception as a function of GDP per capita at purchasing power parity (GDPPCP).  We found, however, that there was additional growth over time and the parameter for time-related usage growth (tconr ) controls that. The user can also change contraception use via an exogenous multiplier (contrusm ).

Although those three distal and proximate drivers substantially determine the forecasts of fertility rate, there are several additional elements that influence it.  First, we calculate the historical growth rate of TFR (TFRgr) and use that internal variable in the first few years so as to maintain an inertial pattern of change in TFR consistent with history; we phase out that inertial element in favor of endogenously computed factors over a 10-year period.  Second, we have used another time dependent parameter (ttfrr ) to allow  introduction of  somewhat faster or slower growth rates in TFR.  Mostly we have used that as a tuning parameter to adjust our long-term global population forecasts to be more consistent with those of others such as the UN Population Division or the International Institute of Applied Systems Analysis.  Normally that parameter is very small or zero.  Third, we provide the user with a direct multiplier on the fertility rate (tfrm ).  And finally, not knowing what the long-term minimum fertility rate might be in a world where for many countries rates have fallen very substantially below replacement rates, we provide such a minimum (tfrmin ).  Since some countries are below most expected minimums and therefore below common values of that parameter, we phase that minimum in over time with a convergence parameter (tfrconv ), which serves double duty by also marking the number of years of convergence of TFR itself to the values that the function with distal and proximate drivers produces. 

Mortality Detail

The current default representation of mortality and life expectancy in IFs relies entirely on the health model so please see its documentation.  That model computes deaths by age and sex and uses those to compute total deaths (DEATHS) as well deaths by category of cause (DEATHCAT).  It also computes life expectancy (LIFEXP) and infant mortality (INFMOR), variables of importance to the population model.  Two parameters in the health model allow multiplicative intervention with respect to total deaths (mortm ) and those by cause (hlmortm ).

There is, however, a legacy representation of mortality in IFs (available if the health model is switched off with hlmodelsw =0) that reverses that logic and uses the model's calculation of an initial estimate of life expectancy to drive mortality by age and sex (not cause). Life expectancy normally increases and mortality normally decreases as GDP per capita rises (see the economic model) or as the income share of the poorest 20% of the population increases. In the legacy representation, the initial calculation of life expectancy is imposed on an initial mortality distribution that provides a country-specific age and sex profile of mortality.

A number factors then further affect and alter the mortality distribution in the legacy mortality structure. These include deaths related to warfare (CIVDM), to AIDS, and possibly to starvation (via infant mortality, because it is primarily the very young who are at risk). In addition, the user of the model may introduce greater or lesser mortality via a mortality multiplier.

At the same time, however, increases in life expectancy shift the mortality distribution from its initial condition towards an ultimate life (survivor) table as life expectancy approaches that built into the ultimate life table (approximately 85 in the 1998 revision of the UN population tables).

Once the mortality distribution adjustments are made and deaths can be computed from it, the life expectancy is recomputed. <header><hgroup>

Nutrition

</hgroup></header> As noted in the overview of demographics, there is an element of mortality calculations that has historically captured considerable attention from many of those interested in long-term population forecasting, going back to at least Malthus's elaboration of the negative feedback loop linking undernutrition or starvation to higher death rates.  Donella Meadows, et al. (1972) popularized this in their discussion of The Limits to Growth.

In the current default representation of the model, with the health model turned on, that health model computes the rate of undernutrition among children (MALNCHP), the resultant total numbers of children undernourished (MALNCHIL), and the deaths associated with undernutrition.  Undernutrition in the health model is a function of calories, but also of access to improved sanitation and clean water.  Health interventions, including those to reduce diarrhea, can supplement greater access to calories to reduce the undernutrition and associated mortality.

In the legacy version of the population model there is a cruder and more overtly Malthusian representation.  A comparison of calories available with those needed can generate starvation deaths. Historical and contemporary data do not exist, however, to support calculation of starvation (the deaths of those who are severely malnourished are normally attributed to various diseases that prey on them, such as diarrhea among children).  For this reason and because of the considerably greater sophistication of the health model, we recommend leaving the health model engaged.