Factors associated with excess all-cause mortality in the first wave of the COVID-19 pandemic in the UK: A time series analysis using the Clinical Practice Research Datalink

BackgroundExcess mortality captures the total effect of the Coronavirus Disease 2019 (COVID-19) pandemic on mortality and is not affected by misspecification of cause of death. We aimed to describe how health and demographic factors were associated with excess mortality during, compared to before, the pandemic.Methods and findingsWe analysed a time series dataset including 9,635,613 adults (≥40 years old) registered at United Kingdom general practices contributing to the Clinical Practice Research Datalink. We extracted weekly numbers of deaths and numbers at risk between March 2015 and July 2020, stratified by individual-level factors. Excess mortality during Wave 1 of the UK pandemic (5 March to 27 May 2020) compared to the prepandemic period was estimated using seasonally adjusted negative binomial regression models. Relative rates (RRs) of death for a range of factors were estimated before and during Wave 1 by including interaction terms. We found that all-cause mortality increased by 43% (95% CI 40% to 47%) during Wave 1 compared with prepandemic. Changes to the RR of death associated with most sociodemographic and clinical characteristics were small during Wave 1 compared with prepandemic. However, the mortality RR associated with dementia markedly increased (RR for dementia versus no dementia prepandemic: 3.5, 95% CI 3.4 to 3.5; RR during Wave 1: 5.1, 4.9 to 5.3); a similar pattern was seen for learning disabilities (RR prepandemic: 3.6, 3.4 to 3.5; during Wave 1: 4.8, 4.4 to 5.3), for black or South Asian ethnicity compared to white, and for London compared to other regions. Relative risks for morbidities were stable in multiple sensitivity analyses. However, a limitation of the study is that we cannot assume that the risks observed during Wave 1 would apply to other waves due to changes in population behaviour, virus transmission, and risk perception.ConclusionsThe first wave of the UK COVID-19 pandemic appeared to amplify baseline mortality risk to approximately the same relative degree for most population subgroups. However, disproportionate increases in mortality were seen for those with dementia, learning disabilities, non-white ethnicity, or living in London.

Helen Strongman and colleagues investigate the health and demographic factors associated with excess mortality during, as compared to before, the COVID-19 pandemic.     

The 1918-1919 influenza pandemic in England and Wales: spatial patterns in transmissibility and mortality impact

Spatial variations in disease patterns of the 1918-1919 influenza pandemic remain poorly studied. We explored the association between influenza death rates, transmissibility and several geographical and demographic indicators for the autumn and winter waves of the 1918-1919 pandemic in cities, towns and rural areas of England and Wales. Average measures of transmissibility, estimated by the reproduction number, ranged between 1.3 and 1.9, depending on model assumptions and pandemic wave and showed little spatial variation. Death rates varied markedly with urbanization, with 30-40% higher rates in cities and towns compared with rural areas. In addition, death rates varied with population size across rural settings, where low population areas fared worse. By contrast, we found no association between transmissibility, death rates and indicators of population density and residential crowding. Further studies of the geographical mortality patterns associated with the 1918-1919 influenza pandemic may be useful for pandemic planning.     

Mid‐Holocene decline in African buffalos inferred from Bayesian coalescent‐based analyses of microsatellites and mitochondrial DNA

Genetic studies concerned with the demographic history of wildlife species can help elucidate the role of climate change and other forces such as human activity in shaping patterns of divergence and distribution. The African buffalo (Syncerus caffer) declined dramatically during the rinderpest pandemic in the late 1800s, but little is known about the earlier demographic history of the species. We analysed genetic variation at 17 microsatellite loci and a 302‐bp fragment of the mitochondrial DNA control region to infer past demographic changes in buffalo populations from East Africa. Two Bayesian coalescent‐based methods as well as traditional bottleneck tests were applied to infer detailed dynamics in buffalo demographic history. No clear genetic signature of population declines related to the rinderpest pandemic could be detected. However, Bayesian coalescent modelling detected a strong signal of African buffalo population declines in the order of 75–98%, starting in the mid‐Holocene (approximately 3–7000 years ago). The signature of decline was remarkably consistent using two different coalescent‐based methods and two types of molecular markers. Exploratory analyses involving various prior assumptions did not seriously affect the magnitude or timing of the inferred population decline. Climate data show that tropical Africa experienced a pronounced transition to a drier climate approximately 4500 years ago, concurrent with the buffalo decline. We therefore propose that the mid‐Holocene aridification of East Africa caused a major decline in the effective population size of the buffalo, a species reliant on moist savannah habitat for its existence.     

Mortality Burden of the 2009 A/H1N1 Influenza Pandemic in France: Comparison to Seasonal Influenza and the A/H3N2 Pandemic

Background

The mortality burden of the 2009 A/H1N1 pandemic remains unclear in many countries due to delays in reporting of death statistics. We estimate the age- and cause-specific excess mortality impact of the pandemic in France, relative to that of other countries and past epidemic and pandemic seasons.

Methods

We applied Serfling and Poisson excess mortality approaches to model weekly age- and cause-specific mortality rates from June 1969 through May 2010 in France. Indicators of influenza activity, time trends, and seasonal terms were included in the models. We also reviewed the literature for country-specific estimates of 2009 pandemic excess mortality rates to characterize geographical differences in the burden of this pandemic.

Results

The 2009 A/H1N1 pandemic was associated with 1.0 (95% Confidence Intervals (CI) 0.2–1.9) excess respiratory deaths per 100,000 population in France, compared to rates per 100,000 of 44 (95% CI 43–45) for the A/H3N2 pandemic and 2.9 (95% CI 2.3–3.7) for average inter-pandemic seasons. The 2009 A/H1N1 pandemic had a 10.6-fold higher impact than inter-pandemic seasons in people aged 5–24 years and 3.8-fold lower impact among people over 65 years.

Conclusions

The 2009 pandemic in France had low mortality impact in most age groups, relative to past influenza seasons, except in school-age children and young adults. The historical A/H3N2 pandemic was associated with much larger mortality impact than the 2009 pandemic, across all age groups and outcomes. Our 2009 pandemic excess mortality estimates for France fall within the range of previous estimates for high-income regions. Based on the analysis of several mortality outcomes and comparison with laboratory-confirmed 2009/H1N1 deaths, we conclude that cardio-respiratory and all-cause mortality lack precision to accurately measure the impact of this pandemic in high-income settings and that use of more specific mortality outcomes is important to obtain reliable age-specific estimates.     

Movements and source–sink dynamics of a Masai giraffe metapopulation

Spatial variation in habitat quality and anthropogenic factors, as well as social structure, can lead to spatially structured populations of animals. Demographic approaches can be used to improve our understanding of the dynamics of spatially structured populations and help identify subpopulations critical for the long-term persistence of regional metapopulations. We provide a regional metapopulation analysis to inform conservation management for Masai giraffes (Giraffa camelopardalis tippelskirchi) in five subpopulations defined by land management designations. We used data from an individual-based mark–recapture study to estimate subpopulation sizes, subpopulation growth rates, and movement probabilities among subpopulations. We assessed the source–sink structure of the study population by calculating source–sink statistics, and we created a female-based matrix metapopulation model composed of all subpopulations to examine how variation in demographic components of survival, reproduction, and movement affected metapopulation growth rate. Movement data indicated no subpopulation was completely isolated, but movement probabilities varied among subpopulations. Source–sink statistics and net flow of individuals indicated three subpopulations were sources, while two subpopulations were sinks. We found areas with higher wildlife protection efforts and fewer anthropogenic impacts were sources, and less-protected areas were identified as sinks. Our results highlight the importance of identifying source–sink dynamics among subpopulations for effective conservation planning and emphasize how protected areas can play an important role in sustaining metapopulations.     

DNA‐based approaches for evaluating historical demography in terrestrial vertebrates

Contemporary DNA sequences can provide information about the historical demography of a species. However, different molecular markers are informative under different circumstances. In particular, mitochondrial (mt)DNA is uniparentally inherited and haploid in most vertebrates and thus has a smaller effective population size than diploid, biparentally inherited nuclear (n)DNA. Here, we review the characteristics of mtDNA and nDNA in the context of historical demography. In particular, we address how their contrasting rates of evolution and sex‐biased dispersal can lead to different demographic inferences. We do so in the context of an extensive review of the vertebrate literature that describes the use of mtDNA and nDNA sequence data in demographic reconstruction. We discuss the effects of coalescence, effective population size, substitution rates, and sex‐biased dispersal on informative timeframes and expected patterns of genetic differentiation. We argue that mtDNA variationin species with male‐biased dispersal can imply deviations from neutrality that do not reflect actual population expansion or selection. By contrast, mtDNA can be more informative when coalescence has occurred within the recent past, which appears to be the case with many vertebrates. We also compare the application and interpretation of demographic and neutrality test statistics in historical demography studies. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 112 , 367–386.     

Universal risk phenotype of US counties for flu-like transmission to improve county-specific COVID-19 incidence forecasts

The spread of a communicable disease is a complex spatio-temporal process shaped by the specific transmission mechanism, and diverse factors including the behavior, socio-economic and demographic properties of the host population. While the key factors shaping transmission of influenza and COVID-19 are beginning to be broadly understood, making precise forecasts on case count and mortality is still difficult. In this study we introduce the concept of a universal geospatial risk phenotype of individual US counties facilitating flu-like transmission mechanisms. We call this the Universal Influenza-like Transmission (UnIT) score, which is computed as an information-theoretic divergence of the local incidence time series from an high-risk process of epidemic initiation, inferred from almost a decade of flu season incidence data gleaned from the diagnostic history of nearly a third of the US population. Despite being computed from the past seasonal flu incidence records, the UnIT score emerges as the dominant factor explaining incidence trends for the COVID-19 pandemic over putative demographic and socio-economic factors. The predictive ability of the UnIT score is further demonstrated via county-specific weekly case count forecasts which consistently outperform the state of the art models throughout the time-line of the COVID-19 pandemic. This study demonstrates that knowledge of past epidemics may be used to chart the course of future ones, if transmission mechanisms are broadly similar, despite distinct disease processes and causative pathogens.     

Demographic balancing of migrant and resident elk in a partially migratory population through forage–predation tradeoffs

Partial migration is widespread in ungulates, yet few studies have assessed demographic mechanisms for how these alternative strategies are maintained in populations. Over the past two decades the number of resident individuals of the Ya Ha Tinda elk herd near Banff National Park has been increasing proportionally despite an overall population decline. We compared demographic rates of migrant and resident elk to test for demographic mechanisms partial migration. We determined adult female survival for 132 elk, pregnancy rates for 150 female elk, and elk calf survival for 79 calves. Population vital rates were combined in Leslie‐matrix models to estimate demographic fitness, which we defined as the migration strategy‐specific population growth rate. We also tested for differences in factors influencing risk of mortality between migratory strategies for adult females using Cox‐proportional hazards regression and time‐varying covariates of exposure to forage biomass, wolf predation risk, and group size. Despite higher pregnancy rates and winter calf weights associated with higher forage quality, survival of migrant adult females and calves were lower than resident elk. Resident elk traded high quality food to reduce predation risk by selecting areas close to human activity, and by living in group sizes 20% larger than migrants. Thus, residents experienced higher adult female survival and calf survival, but lower pregnancy and calf weights. Cause‐specific mortality of migrants was dominated by wolf and grizzly bear mortality, whereas resident mortality was dominated by human hunting. Demographic differences translated into slightly higher (2–3%), but non‐significant, resident population growth rate compared to migrant elk, suggesting demographic balancing between resident strategies during our study. Despite statistical equivalence, our results are also consistent with slow long‐term declines in migrants because of high predation because of higher wolf‐caused mortality in migrants. These results emphasize that migrants and residents will make different tradeoffs between forage and risk may affect the demographic balance of partially migratory populations, which may explain recent declines in migratory behavior in many ungulate populations around the world.     

Were Equatorial Regions Less Affected by the 2009 Influenza Pandemic? The Brazilian Experience

Although it is in the Tropics where nearly half of the world population lives and infectious disease burden is highest, little is known about the impact of influenza pandemics in this area. We investigated the mortality impact of the 2009 influenza pandemic relative to mortality rates from various outcomes in pre-pandemic years throughout a wide range of latitudes encompassing the entire tropical, and part of the subtropical, zone of the Southern Hemisphere (+5(°)N to -35(°)S) by focusing on a country with relatively uniform health care, disease surveillance, immunization and mitigation policies: Brazil. To this end, we analyzed laboratory-confirmed deaths and vital statistics mortality beyond pre-pandemic levels for each Brazilian state. Pneumonia, influenza and respiratory mortality were significantly higher during the pandemic, affecting predominantly adults aged 25 to 65 years. Overall, there were 2,273 and 2,787 additional P&I- and respiratory deaths during the pandemic, corresponding to a 5.2% and 2.7% increase, respectively, over average pre-pandemic annual mortality. However, there was a marked spatial structure in mortality that was independent of socio-demographic indicators and inversely related with income: mortality was progressively lower towards equatorial regions, where low or no difference from pre-pandemic mortality levels was identified. Additionally, the onset of pandemic-associated mortality was progressively delayed in equatorial states. Unexpectedly, there was no additional mortality from circulatory causes. Comparing disease burden reliably across regions is critical in those areas marked by competing health priorities and limited resources. Our results suggest, however, that tropical regions of the Southern Hemisphere may have been disproportionally less affected by the pandemic, and that climate may have played a key role in this regard. These findings have a direct bearing on global estimates of pandemic burden and the assessment of the role of immunological, socioeconomic and environmental drivers of the transmissibility and severity of this pandemic.     

Unexpected Demography in the Recovery of an Endangered Primate Population

Assessments of the status of endangered species have focused on population sizes, often without knowledge of demographic and behavioral processes underlying population recovery. We analyzed demographic data from a 28-year study of a critically endangered primate, the northern muriqui, to investigate possible changes in demographic rates as this population recovered from near extirpation. As the population increased from 60 to nearly 300 individuals, its growth rate declined due to increased mortality and male-biased birth sex ratios; the increased mortality was not uniform across ages and sexes, and there has been a recent increase in mortality of prime-aged males. If not for a concurrent increase in fertility rates, the population would have stabilized at 200 individuals instead of continuing to grow. The unexpected increase in fertility rates and in adult male mortality can be attributed to the muriquis’ expansion of their habitat by spending more time on the ground. The demographic consequences of this behavioral shift must be incorporated into management tactics for this population and emphasize the importance of understanding demographic rates in the recovery of endangered species.     

LOCAL ADAPTATION WHEN COMPETITION DEPENDS ON PHENOTYPIC SIMILARITY

Recent work incorporating demographic–genetic interactions indicates the importance of population size, gene flow, and selection in influencing local adaptation. This work typically assumes that density‐dependent survival affects individuals equally, but individuals in natural population rarely compete equally. Among‐individual differences in resource use generate stronger competition between more similar phenotypes (frequency‐dependent competition) but it remains unclear how this additional form of selection changes the interactions between population size, gene flow, and local stabilizing selection. Here, we integrate migration–selection dynamics with frequency‐dependent competition. We developed a coupled demographic‐quantitative genetic model consisting of two patches connected by dispersal and subject to local stabilizing selection and competition. Our model shows that frequency‐dependent competition slightly increases local adaptation, greatly increases genetic variance within patches, and reduces the amount that migration depresses population size, despite the increased genetic variance load. The effects of frequency‐dependence depend on the strength of divergent selection, trait heritability, and when mortality occurs in the life cycle in relation to migration and reproduction. Essentially, frequency‐dependent competition reduces the density‐dependent interactions between migrants and residents, the extent to which depends on how different and common immigrants are compared to residents. Our results add new dynamics that illustrate how competition can alter the effects of gene flow and divergent selection on local adaptation and population carrying capacities.     

Challenges to assessing connectivity between massive populations of the Australian plague locust

Linking demographic and genetic dispersal measures is of fundamental importance for movement ecology and evolution. However, such integration can be difficult, particularly for highly fecund species that are often the target of management decisions guided by an understanding of population movement. Here, we present an example of how the influence of large population sizes can preclude genetic approaches from assessing demographic population structuring, even at a continental scale. The Australian plague locust, Chortoicetes terminifera, is a significant pest, with populations on the eastern and western sides of Australia having been monitored and managed independently to date. We used microsatellites to assess genetic variation in 12 C. terminifera population samples separated by up to 3000 km. Traditional summary statistics indicated high levels of genetic diversity and a surprising lack of population structure across the entire range. An approximate Bayesian computation treatment indicated that levels of genetic diversity in C. terminifera corresponded to effective population sizes conservatively composed of tens of thousands to several million individuals. We used these estimates and computer simulations to estimate the minimum rate of dispersal, m, that could account for the observed range-wide genetic homogeneity. The rate of dispersal between both sides of the Australian continent could be several orders of magnitude lower than that typically considered as required for the demographic connectivity of populations.     

Potential population-level effects of increased haulout-related mortality of Pacific walrus calves

Availability of summer sea ice has been decreasing in the Chukchi Sea during recent decades, and increasing numbers of Pacific walruses have begun using coastal haulouts in late summer during years when sea ice retreats beyond the continental shelf. Calves and yearlings are particularly susceptible to being crushed during disturbance events that cause the herd to panic and stampede at these large haulouts, but the potential population-level effects of this mortality are unknown. We used recent harvest data, along with previous assumptions about demographic parameters for this population, to estimate female population size and structure in 2009 and project these numbers forward using a range of assumptions about future harvests and haulout-related mortality that might result from increased use of coastal haulouts during late summer. We found that if demographic parameters were held constant, the levels of harvest that occurred during 1990–2008 would have allowed the population to grow during that period. Our projections indicate, however, that an increase in haulout-related mortality affecting only calves has a greater effect on the population than an equivalent increase in harvest-related mortality distributed among all age classes. Therefore, disturbance-related mortality of calves at coastal haulouts may have relatively important population consequences.     

On the Use of Published Demographic Data for Population-Level Risk Assessment in Birds

As the practice of using population models for wildlife risk assessment has become more common, so has the practice of using surrogate data, typically taken from the published scientific literature, as inputs for demographic models. This practice clearly exposes the user to inferential errors. However, it is likely to continue because demographic data are expensive to gather. We review potential errors associated with the use of previously published demographic data and how those errors propagate into the endpoints of demographic projection models. We suggest methods for inferring bias in model endpoints when multiple and opposing biases are present in the demographic input data. We provide an example using Eastern Meadowlarks (Sturnella magna), a common songbird in Midwestern grasslands and agro-ecosystems. We conclude with a brief review of methods that could improve inference made using published demographic data, including methods from life-history theory, meta-analysis, and Bayesian statistics.     

Stage-structured ontogeny in resource populations generates non-additive stabilizing and de-stabilizing forces in populations and communities

Demographic heterogeneity influences how populations respond to density dependent intraspecific competition and trophic interactions. Distinct stages across an organism's development, or ontogeny, are an important example of demographic heterogeneity. In consumer populations, ontogenetic stage structure has been shown to produce categorical differences in population dynamics, community dynamics and even species coexistence compared to models lacking explicit ontogeny. The study of consumer–resource interactions must also consider the ontogenetic stage structure of the resource itself, particularly plants, given their fundamental role at the basis of terrestrial food webs. We incorporate distinct ontogenetic stages of plants into an adaptable multi-stage consumer–resource modeling framework that facilitates studying how stage specific consumers shape trophic dynamics at low trophic levels. We describe the role of density dependent demographic rates in mediating the dynamics of stage-structured plant populations. We then investigate how these demographic rates interact with consumer pressure to influence stability and coexistence in multiple stage-specific consumer–resource interactions. Results detail how density dependent effects across distinct ontogenetic stages in plant development produce non-additivity in the drivers of dynamic stability both in single populations and in consumer–resource settings, challenging the ubiquity of certain traditional ecological dynamic paradigms. We also find categorical differences in the population variability induced by herbivores consuming separate plant stages. Consumer–resource models, such as plant–herbivore interactions, often average out demographic heterogeneity in populations. Here, we show that explicitly including plant demographic heterogeneity through ontogeny yields distinct dynamic expectations for both plants and herbivores compared to traditional consumer–resource formulations. Our results indicate that efforts to understand the demographic effect of herbivores on plant populations may need to also consider the effects of plant demographics on herbivores and the reciprocal relationship between them.     

Soil-related variations in the population dynamics of six dipterocarp tree species with strong habitat preferences

Differences in the density of conspecific tree individuals in response to environmental gradients are well documented for many tree species, but how such density differences are generated and maintained is poorly understood. We examined the segregation of six dipterocarp species among three soil types in the Pasoh tropical forest, Malaysia. We examined how individual performance and population dynamics changed across the soil types using 10-year demographic data to compare tree performance across soil types, and constructed population matrix models to analyze the population dynamics. Species showed only minor changes in mortality and juvenile growth across soil types, although recruitment differed greatly. Clear, interspecific demographic trade-offs between growth and mortality were found in all soil types. The relative trade-offs by a species did not differ substantially among the soil types. Population sizes were projected to remain stable in all soil types for all species with one exception. Our life-table response experiment demonstrated that the population dynamics of a species differed only subtly among soil types. Therefore, species with strong density differences across soil types do not necessarily differ greatly in their population dynamics across the soil types. In contrast, interspecific differences in population dynamics were large. The trade-off between mortality and growth led to a negative correlation between the contributions of mortality and growth to variations in the population growth rate (λ) and thus reduced their net contributions. Recruitment had little impact on the variation in λ. The combination of these factors resulted in little variation in λ among species.     

Demographic effects of road mortality on mammalian populations: a systematic review

In light of rapidly expanding road networks worldwide, there is increasing global awareness of the growing amount of mammalian roadkill. However, the ways in which road mortality affects the population dynamics of different species remains largely unclear. We aimed to categorise the demographic parameters in mammalian populations around the world that are directly or indirectly affected by road mortality, as well as identify the most effective study designs for quantifying population-level consequences of road mortality. We conducted a comprehensive systematic review to synthesise literature published between 2000 and 2021 and out of 11,238 unique studies returned, 83 studies were retained comprising 69 mammalian species and 150 populations. A bias towards research-intensive countries and larger mammals was apparent. Although searches were conducted in five languages, all studies meeting the inclusion criteria were in English. Relatively few studies (13.3%) provided relevant demographic context to roadkill figures, hampering understanding of the impacts on population persistence. We categorised five direct demographic parameters affected by road mortality: sex- and age-biased mortality, the percentage of a population killed on roads per year (values up to 50% were reported), the contribution of roadkill to total mortality rates (up to 80%), and roadkill during inter-patch or long-distance movements. Female-biased mortality may be more prevalent than previously recognised and is likely to be critical to population dynamics. Roadkill was the greatest source of mortality for 28% of studied populations and both additive and compensatory mechanisms to roadkill were found to occur, bringing varied challenges to conservation around roads. In addition, intra-specific population differences in demographic effects of road mortality were common. This highlights that the relative importance of road mortality is likely to be context specific as the road configuration and habitat quality surrounding a population can vary. Road ecology studies that collect data on key life parameters, such as age/stage/sex-specific survival and dispersal success, and that use a combination of methods are critical in understanding long-term impacts. Quantifying the demographic impacts of road mortality is an important yet complex consideration for proactive road management.     

A linear mixed model to estimate COVID-19-induced excess mortality

The Corona Virus Disease (COVID-19) pandemic has increased mortality in countries worldwide. To evaluate the impact of the pandemic on mortality, the use of excess mortality rather than reported COVID-19 deaths has been suggested. Excess mortality, however, requires estimation of mortality under nonpandemic conditions. Although many methods exist to forecast mortality, they are either complex to apply, require many sources of information, ignore serial correlation, and/or are influenced by historical excess mortality. We propose a linear mixed model that is easy to apply, requires only historical mortality data, allows for serial correlation, and down-weighs the influence of historical excess mortality. Appropriateness of the linear mixed model is evaluated with fit statistics and forecasting accuracy measures for Belgium and the Netherlands. Unlike the commonly used 5-year weekly average, the linear mixed model is forecasting the year-specific mortality, and as a result improves the estimation of excess mortality for Belgium and the Netherlands.     

Demographic,mechanistic and density-dependent determinants of population growth rate: a case study in an avian predator

Identifying the determinants of population growth rate is a central topic in population ecology. Three approaches (demographic, mechanistic and density-dependent) used historically to describe the determinants of population growth rate are here compared and combined for an avian predator, the barn owl (Tyto alba). The owl population remained approximately stable (r approximately 0) throughout the period from 1979 to 1991. There was no evidence of density dependence as assessed by goodness of fit to logistic population growth. The finite (lambda) and instantaneous (r) population growth rates were significantly positively related to food (field vole) availability. The demographic rates, annual adult mortality, juvenile mortality and annual fecundity were reported to be correlated with vole abundance. The best fit (R(2) = 0.82) numerical response of the owl population described a positive effect of food (field voles) and a negative additive effect of owl abundance on r. The numerical response of the barn owl population to food availability was estimated from both census and demographic data, with very similar results. Our analysis shows how the demographic and mechanistic determinants of population growth rate are linked; food availability determines demographic rates, and demographic rates determine population growth rate. The effects of food availability on population growth rate are modified by predator abundance.     

Population and prehistory I: Food-dependent population growth in constant environments

We present a demographic model that describes the feedbacks between food supply, human mortality and fertility rates, and labor availability in expanding populations, where arable land area is not limiting. This model provides a quantitative framework to describe how environment, technology, and culture interact to influence the fates of preindustrial agricultural populations. We present equilibrium conditions and derive approximations for the equilibrium population growth rate, food availability, and other food-dependent measures of population well-being. We examine how the approximations respond to environmental changes and to human choices, and find that the impact of environmental quality depends upon whether it manifests through agricultural yield or maximum (food-independent) survival rates. Human choices can complement or offset environmental effects: greater labor investments increase both population growth and well-being, and therefore can counteract lower agricultural yield, while fertility control decreases the growth rate but can increase or decrease well-being. Finally we establish equilibrium stability criteria, and argue that the potential for loss of local stability at low population growth rates could have important consequences for populations that suffer significant environmental or demographic shocks.