Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">The Invisible Cliff: Abrupt Imposition of Malthusian Equilibrium in a Natural-Fertility, Agrarian Society</title><author><name sortKey="Puleston, Cedric" sort="Puleston, Cedric" uniqKey="Puleston C" first="Cedric" last="Puleston">Cedric Puleston</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Anthropology, University of California Davis, Davis, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Tuljapurkar, Shripad" sort="Tuljapurkar, Shripad" uniqKey="Tuljapurkar S" first="Shripad" last="Tuljapurkar">Shripad Tuljapurkar</name><affiliation><nlm:aff id="aff2"><addr-line>Department of Biology, Stanford University, Stanford, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Winterhalder, Bruce" sort="Winterhalder, Bruce" uniqKey="Winterhalder B" first="Bruce" last="Winterhalder">Bruce Winterhalder</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Anthropology, University of California Davis, Davis, California, United States of America</addr-line></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">24498131</idno><idno type="pmc">3909123</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3909123</idno><idno type="RBID">PMC:3909123</idno><idno type="doi">10.1371/journal.pone.0087541</idno><date when="2014">2014</date><idno type="wicri:Area/Pmc/Corpus">000376</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000376</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">The Invisible Cliff: Abrupt Imposition of Malthusian Equilibrium in a Natural-Fertility, Agrarian Society</title><author><name sortKey="Puleston, Cedric" sort="Puleston, Cedric" uniqKey="Puleston C" first="Cedric" last="Puleston">Cedric Puleston</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Anthropology, University of California Davis, Davis, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Tuljapurkar, Shripad" sort="Tuljapurkar, Shripad" uniqKey="Tuljapurkar S" first="Shripad" last="Tuljapurkar">Shripad Tuljapurkar</name><affiliation><nlm:aff id="aff2"><addr-line>Department of Biology, Stanford University, Stanford, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Winterhalder, Bruce" sort="Winterhalder, Bruce" uniqKey="Winterhalder B" first="Bruce" last="Winterhalder">Bruce Winterhalder</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Anthropology, University of California Davis, Davis, California, United States of America</addr-line></nlm:aff></affiliation></author></analytic><series><title level="j">PLoS ONE</title><idno type="eISSN">1932-6203</idno><imprint><date when="2014">2014</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>Analysis of a natural fertility agrarian society with a multi-variate model of population ecology isolates three distinct phases of population growth following settlement of a new habitat: (1) a sometimes lengthy <italic>copial</italic> phase of surplus food production and constant vital rates; (2) a brief <italic>transition</italic> phase in which food shortages rapidly cause increased mortality and lessened fertility; and (3) a <italic>Malthusian</italic> phase of indefinite length in which vital rates and quality of life are depressed, sometimes strikingly so. Copial phase duration declines with increases in the size of the founding group, maximum life expectancy and fertility; it increases with habitat area and yield per hectare; and, it is unaffected by the sensitivity of vital rates to hunger. Transition phase duration is unaffected by size of founding population and area of settlement; it declines with yield, life expectancy, fertility and the sensitivity of vital rates to hunger. We characterize the transition phase as the <italic>Malthusian transition interval (MTI)</italic>, in order to highlight how little time populations generally have to adjust. Under food-limited density dependence, the copial phase passes quickly to an equilibrium of grim Malthusian constraints, in the manner of a runner dashing over an invisible cliff. The three-phase pattern diverges from widely held intuitions based on standard Lotka-Verhulst approaches to population regulation, with implications for the analysis of socio-cultural evolution, agricultural intensification, bioarchaeological interpretation of food stress in prehistoric societies, and state-level collapse.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Wood, Jw" uniqKey="Wood J">JW Wood</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Lee, Rd" uniqKey="Lee R">RD Lee</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Holling, Cs" uniqKey="Holling C">CS Holling</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lee, Ct" uniqKey="Lee C">CT Lee</name></author><author><name sortKey="Tuljapurkar, S" uniqKey="Tuljapurkar S">S Tuljapurkar</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Puleston, Co" uniqKey="Puleston C">CO Puleston</name></author><author><name sortKey="Tuljapurkar, S" uniqKey="Tuljapurkar S">S Tuljapurkar</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lee, Ct" uniqKey="Lee C">CT Lee</name></author><author><name sortKey="Puleston, Co" uniqKey="Puleston C">CO Puleston</name></author><author><name sortKey="Tuljapurkar, S" uniqKey="Tuljapurkar S">S Tuljapurkar</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Read, Dw" uniqKey="Read D">DW Read</name></author><author><name sortKey="Leblanc, Sa" uniqKey="Leblanc S">SA LeBlanc</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Winterhalder, B" uniqKey="Winterhalder B">B Winterhalder</name></author><author><name sortKey="Baillargeon, W" uniqKey="Baillargeon W">W Baillargeon</name></author><author><name sortKey="Cappelletto, F" uniqKey="Cappelletto F">F Cappelletto</name></author><author><name sortKey="Daniel Jr, Ir" uniqKey="Daniel Jr I">IR Daniel Jr</name></author><author><name sortKey="Prescott, C" uniqKey="Prescott C">C Prescott</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Neubert, Mg" uniqKey="Neubert M">MG Neubert</name></author><author><name sortKey="Caswell, H" uniqKey="Caswell H">H Caswell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gaillard, J M" uniqKey="Gaillard J">J-M Gaillard</name></author><author><name sortKey="Festa Bianchet, M" uniqKey="Festa Bianchet M">M Festa-Bianchet</name></author><author><name sortKey="Yoccoz, Ng" uniqKey="Yoccoz N">NG Yoccoz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tsoularis, A" uniqKey="Tsoularis A">A Tsoularis</name></author><author><name sortKey="Wallace, J" uniqKey="Wallace J">J Wallace</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sibly, Rm" uniqKey="Sibly R">RM Sibly</name></author><author><name sortKey="Barker, D" uniqKey="Barker D">D Barker</name></author><author><name sortKey="Denham, Mc" uniqKey="Denham M">MC Denham</name></author><author><name sortKey="Hone, J" uniqKey="Hone J">J Hone</name></author><author><name sortKey="Pagel, M" uniqKey="Pagel M">M Pagel</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Tainter, Ja" uniqKey="Tainter J">JA Tainter</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Middleton, Gd" uniqKey="Middleton G">GD Middleton</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Turner Ii, Bl" uniqKey="Turner Ii B">BL Turner II</name></author><author><name sortKey="Sabloff, Ja" uniqKey="Sabloff J">JA Sabloff</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Webster, D" uniqKey="Webster D">D Webster</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hamblin, Rl" uniqKey="Hamblin R">RL Hamblin</name></author><author><name sortKey="Pitcher, Bl" uniqKey="Pitcher B">BL Pitcher</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Turchin, P" uniqKey="Turchin P">P Turchin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hunt, Bg" uniqKey="Hunt B">BG Hunt</name></author><author><name sortKey="Elliott, Ti" uniqKey="Elliott T">TI Elliott</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Drysdale, R" uniqKey="Drysdale R">R Drysdale</name></author><author><name sortKey="Zanchetta, G" uniqKey="Zanchetta G">G Zanchetta</name></author><author><name sortKey="Hellstrom, J" uniqKey="Hellstrom J">J Hellstrom</name></author><author><name sortKey="Maas, R" uniqKey="Maas R">R Maas</name></author><author><name sortKey="Fallick, A" uniqKey="Fallick A">A Fallick</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gill, Rb" uniqKey="Gill R">RB Gill</name></author><author><name sortKey="Mayewski, Pa" uniqKey="Mayewski P">PA Mayewski</name></author><author><name sortKey="Nyberg, J" uniqKey="Nyberg J">J Nyberg</name></author><author><name sortKey="Haug, Gh" uniqKey="Haug G">GH Haug</name></author><author><name sortKey="Peterson, Lc" uniqKey="Peterson L">LC Peterson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Peterson, Lc" uniqKey="Peterson L">LC Peterson</name></author><author><name sortKey="Haug, Gh" uniqKey="Haug G">GH Haug</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shaowu, W" uniqKey="Shaowu W">W Shaowu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cullen, Hm" uniqKey="Cullen H">HM Cullen</name></author><author><name sortKey="Demenocal, Pb" uniqKey="Demenocal P">PB deMenocal</name></author><author><name sortKey="Hemming, S" uniqKey="Hemming S">S Hemming</name></author><author><name sortKey="Hemming, G" uniqKey="Hemming G">G Hemming</name></author><author><name sortKey="Brown, Fh" uniqKey="Brown F">FH Brown</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Staubwasser, M" uniqKey="Staubwasser M">M Staubwasser</name></author><author><name sortKey="Sirocko, F" uniqKey="Sirocko F">F Sirocko</name></author><author><name sortKey="Grootes, Pm" uniqKey="Grootes P">PM Grootes</name></author><author><name sortKey="Segl, M" uniqKey="Segl M">M Segl</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kennett, Dj" uniqKey="Kennett D">DJ Kennett</name></author><author><name sortKey="Breitenbach, Sfm" uniqKey="Breitenbach S">SFM Breitenbach</name></author><author><name sortKey="Aquino, Vv" uniqKey="Aquino V">VV Aquino</name></author><author><name sortKey="Asmerom, Y" uniqKey="Asmerom Y">Y Asmerom</name></author><author><name sortKey="Awe, J" uniqKey="Awe J">J Awe</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wright, Le" uniqKey="Wright L">LE Wright</name></author><author><name sortKey="White, Cd" uniqKey="White C">CD White</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, J" uniqKey="Liu J">J Liu</name></author><author><name sortKey="Dietz, T" uniqKey="Dietz T">T Dietz</name></author><author><name sortKey="Carpenter, Sr" uniqKey="Carpenter S">SR Carpenter</name></author><author><name sortKey="Alberti, M" uniqKey="Alberti M">M Alberti</name></author><author><name sortKey="Folke, C" uniqKey="Folke C">C Folke</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="O Sullivan, P" uniqKey="O Sullivan P">P O'Sullivan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Turchin, P" uniqKey="Turchin P">P Turchin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Turchin, P" uniqKey="Turchin P">P Turchin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Diamond, Jm" uniqKey="Diamond J">JM Diamond</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kirch, Pv" uniqKey="Kirch P">PV Kirch</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Johnston, Kj" uniqKey="Johnston K">KJ Johnston</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS One</journal-id><journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">plosone</journal-id><journal-title-group><journal-title>PLoS ONE</journal-title></journal-title-group><issn pub-type="epub">1932-6203</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">24498131</article-id><article-id pub-id-type="pmc">3909123</article-id><article-id pub-id-type="publisher-id">PONE-D-13-31216</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0087541</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Agriculture</subject><subj-group><subject>Agricultural Production</subject></subj-group><subj-group><subject>Sustainable Agriculture</subject></subj-group></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biology</subject><subj-group><subject>Computational Biology</subject><subj-group><subject>Population Modeling</subject></subj-group></subj-group><subj-group><subject>Ecology</subject><subj-group><subject>Agroecology</subject><subj-group><subject>Agro-Population Ecology</subject><subject>Agronomic Ecology</subject><subject>Ecological Political Economy</subject></subj-group></subj-group><subj-group><subject>Ecological Metrics</subject><subj-group><subject>Carrying Capacity</subject><subject>Extinction Risk</subject><subject>Maximum Sustainable Yield</subject><subject>Population Growth</subject><subject>Population Size</subject></subj-group></subj-group><subj-group><subject>Behavioral Ecology</subject></subj-group></subj-group><subj-group><subject>Population Biology</subject><subj-group><subject>Population Dynamics</subject><subj-group><subject>Predator-Prey Dynamics</subject></subj-group></subj-group><subj-group><subject>Population Metrics</subject><subj-group><subject>Birth Rate</subject><subject>Carrying Capacity</subject><subject>Death Rate</subject><subject>Population Growth</subject><subject>Population Size</subject></subj-group></subj-group><subj-group><subject>Population Ecology</subject><subject>Population Modeling</subject></subj-group></subj-group><subj-group><subject>Theoretical Biology</subject></subj-group></subj-group></article-categories><title-group><article-title>The Invisible Cliff: Abrupt Imposition of Malthusian Equilibrium in a Natural-Fertility, Agrarian Society</article-title><alt-title alt-title-type="running-head">Abrupt Imposition of Malthusian Equilibrium</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Puleston</surname><given-names>Cedric</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author"><name><surname>Tuljapurkar</surname><given-names>Shripad</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Winterhalder</surname><given-names>Bruce</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Anthropology, University of California Davis, Davis, California, United States of America</addr-line></aff><aff id="aff2"><label>2</label><addr-line>Department of Biology, Stanford University, Stanford, California, United States of America</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>O'Connor</surname><given-names>Kathleen A.</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1"><addr-line>University of Washington, United States of America</addr-line></aff><author-notes><corresp id="cor1">* E-mail: <email>puleston@ucdavis.edu</email></corresp><fn fn-type="conflict"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="con"><p>Conceived and designed the experiments: CP ST BW. Performed the experiments: CP ST BW. Analyzed the data: CP ST BW. Wrote the paper: CP BW.</p></fn></author-notes><pub-date pub-type="collection"><year>2014</year></pub-date><pub-date pub-type="epub"><day>31</day><month>1</month><year>2014</year></pub-date><volume>9</volume><issue>1</issue><elocation-id>e87541</elocation-id><history><date date-type="received"><day>29</day><month>7</month><year>2013</year></date><date date-type="accepted"><day>30</day><month>12</month><year>2013</year></date></history><permissions><copyright-year>2014</copyright-year><copyright-holder>Puleston et al</copyright-holder><license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p></license></permissions><abstract><p>Analysis of a natural fertility agrarian society with a multi-variate model of population ecology isolates three distinct phases of population growth following settlement of a new habitat: (1) a sometimes lengthy <italic>copial</italic> phase of surplus food production and constant vital rates; (2) a brief <italic>transition</italic> phase in which food shortages rapidly cause increased mortality and lessened fertility; and (3) a <italic>Malthusian</italic> phase of indefinite length in which vital rates and quality of life are depressed, sometimes strikingly so. Copial phase duration declines with increases in the size of the founding group, maximum life expectancy and fertility; it increases with habitat area and yield per hectare; and, it is unaffected by the sensitivity of vital rates to hunger. Transition phase duration is unaffected by size of founding population and area of settlement; it declines with yield, life expectancy, fertility and the sensitivity of vital rates to hunger. We characterize the transition phase as the <italic>Malthusian transition interval (MTI)</italic>, in order to highlight how little time populations generally have to adjust. Under food-limited density dependence, the copial phase passes quickly to an equilibrium of grim Malthusian constraints, in the manner of a runner dashing over an invisible cliff. The three-phase pattern diverges from widely held intuitions based on standard Lotka-Verhulst approaches to population regulation, with implications for the analysis of socio-cultural evolution, agricultural intensification, bioarchaeological interpretation of food stress in prehistoric societies, and state-level collapse.</p></abstract><funding-group><funding-statement>This work is supported by the NSF (Division of Behavioral and Cognitive Sciences) Human Systems Dynamics, Collaborative Research proposal, Development and Resilience of Complex Socioeconomic Systems: A Theoretical Model and Case Study from the Maya Lowlands [Proposal# 0827275]. Website: <ext-link ext-link-type="uri" xlink:href="http://www.nsf.gov/div/index.jsp?div=BCS">http://www.nsf.gov/div/index.jsp?div=BCS</ext-link>. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement></funding-group><counts><page-count count="13"></page-count></counts></article-meta></front><body><sec id="s1"><title>Introduction</title><p>Malthus <xref rid="pone.0087541-Malthus1" ref-type="bibr">[1]</xref>, Boserup <xref rid="pone.0087541-Boserup1" ref-type="bibr">[2]</xref>, Wood <xref rid="pone.0087541-Wood1" ref-type="bibr">[3]</xref>, Lee <xref rid="pone.0087541-Lee1" ref-type="bibr">[4]</xref>,<xref rid="pone.0087541-Lee2" ref-type="bibr">[5]</xref> and others have given us riches of partial insights and unresolved puzzles on the topics of population, environment, economics and social evolution. Despite an extensive intellectual history and acknowledged importance, our methods for analyzing population and ecology are not well developed. We commonly fall back on the insights of a few pioneers and analytical models of too few parameters. Many of us have formed our basic intuitions from conceptual bases not much advanced over those developed in the 18<sup>th</sup> and early 19<sup>th</sup> century by Verhulst (1804–1849), Lotka (1880–1949), and Volterra (1860–1940) (see <xref rid="pone.0087541-Hutchinson1" ref-type="bibr">[6]</xref>). Our models of these evolutionary and historical processes often are built from too few variables to accurately represent their dynamics <xref rid="pone.0087541-Kohler1" ref-type="bibr">[7]</xref>,<xref rid="pone.0087541-Turchin1" ref-type="bibr">[8]</xref>,<xref rid="pone.0087541-Wright1" ref-type="bibr">[9]</xref>. For example, these approaches begin by defining carrying capacity as a parameter instead of treating it as a variable dependent on features of the population and its environment. This makes it difficult to analyze realistic population dynamics and regulation. Holling <xref rid="pone.0087541-Holling1" ref-type="bibr">[10]</xref> and others have made headway into developing simple yet powerful models of forager behavior and predator-prey dynamics in animal populations, but these are only distantly applicable to human agriculturalists, and have nothing to say about the links between food, labor, social and economic behavior and demography.</p><p>These complex interactions were considered in schematic form in a classic work by Wood <xref rid="pone.0087541-Wood1" ref-type="bibr">[3]</xref> and later incorporated into a modeling framework called “food-limited demography,” developed in Lee and Tuljapurkar <xref rid="pone.0087541-Lee3" ref-type="bibr">[11]</xref>, Puleston and Tuljapurkar <xref rid="pone.0087541-Puleston1" ref-type="bibr">[12]</xref> and Lee <italic>et al.</italic><xref rid="pone.0087541-Lee4" ref-type="bibr">[13]</xref>. We extend it here to examine the nature and timing of the trajectories of food-limited agricultural populations. We examine the interplay between the size and demographic characteristics of a population and the production and consumption of its food supply. The variables and dynamics we seek to understand rest at the heart of problems from the prehistoric origins of agriculture and the state, to contemporary issues of population and environmental degradation. The framework iteratively assesses the effects of (a) age-sex distribution and population size on a society's labor availability and food requirements; (b) environmental factors affecting yield of resources, given labor; (c) consumption requirements of the population, given its size and age composition; and, (d) the effect of food availability on vital rates of fertility and mortality. From these linked interactions flows a revised schedule of fertility and mortality, leading to a new age distribution and population size, completing the logical circuit of the model. How much the population grows or shrinks depends on how much food it can grow using the available labor pool, technology and skills, and the agrarian landscape it occupies. A well-fed population can grow rapidly, but a hungry one will decline from elevated risk of death and lowered fertility.</p><p>Here we apply the model to examine the developmental course of an agrarian society, from initial settlement to equilibrium, as a means of illuminating relationships between population history, environment and density dependent controls. This paper anticipates a second <xref rid="pone.0087541-Winterhalder1" ref-type="bibr">[14]</xref> in which we explore the interaction of agrarian production, peasant livelihood and the state-level extraction of revenue in the form of goods or labor.</p><p>We envision a population colonizing an empty habitat and show that the resulting population trajectory is characterized by three stages: (i) a relatively long period of food abundance (<italic>copial phase</italic>), followed by (ii) a generally very short <italic>transition phase</italic>, to a (iii) <italic>Malthusian phase</italic> representing a stable equilibrium of indefinite length. We estimate the relative length of these phases as a function of basic environmental and demographic properties. Of special interest is the transition phase (ii), the period initiated by the first impacts of food shortfalls on vital rates. We define the term <italic>Malthusian transition interval (MTI)</italic> to characterize the span between the first experience of hunger and when the population effectively reaches its Malthusian equilibrium. This period represents the interval in which the effects of population pressure have become obvious, but in which there still may be room to make behavioral adaptations to forestall a density-dependent equilibrium. We bound the end of this interval at the arbitrary point when the population's doubling time has slowed to 1000 years. The population is effectively stationary at this point, even though it will continue to approach equilibrium for an infinite span. We show that the <italic>MTI</italic> typically is quite short and, counter-intuitively, that it is shorter the better are the copial phase circumstances of the population. Higher environmental yield, longer lifespan, lower background mortality and higher fertility all act to impose a more abrupt, shorter and, we presume for the populations experiencing them, more painful passage to the full effects of Malthusian food limitations. Like a runner moving quickly over flat ground, the faster the pace, the less time to react when the road begins to fall away beneath her.</p><p>We find population ecology dynamics that differ in important ways from those generated by Lotka-Verhulst approaches. Citing examples from Maya prehistory, we show how the model can be used to evaluate hypotheses to the effect that labor shortages, peasant revolt, or cyclic structural instabilities contributed to state level collapse. We likewise show how bioarchaeologists who adopt Lotka-Verhulst expectations are likely to make incorrect empirical inferences about the impact of population on societal development.</p></sec><sec id="s2"><title>Model</title><sec id="s2a"><title>The Conceptual Model</title><p>We focus on food- and space-limited population dynamics in a constant environment, and on colonization of an empty habitat. Our results would pertain as well to populations recovering from severe declines in their numbers, or to those suddenly released from extant density dependent constraints by environmental changes or technological innovations that significantly elevate subsistence production.</p><p>The model is depicted in <xref ref-type="fig" rid="pone-0087541-g001">Figure 1</xref>; parameters are defined, assigned ranges and units and given brief descriptions in <xref ref-type="table" rid="pone-0087541-t001">Table 1</xref>. We work through the conceptual and mathematical relationships of the model, and then turn attention to the adjustable parameters that allow us to make analytical and interpretive use of it. <xref ref-type="sec" rid="s2">Model</xref> mechanics are represented by the components within the dashed rectangle of <xref ref-type="fig" rid="pone-0087541-g001">Figure 1</xref>.</p><fig id="pone-0087541-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0087541.g001</object-id><label>Figure 1</label><caption><title>Population Ecology in an Agrarian Society.</title><p>Components within the dashed rectangle make up the model itself; solid arrows depict the causal relationships that link them. Model parameters, those that must be set and can be manipulated for analytical purposes, are shown outside the dashed box. Dotted arrows show their point of operation in the logical structure of the model. Time is indexed implicitly in all model terms; <italic>x</italic> indexes age.</p></caption><graphic xlink:href="pone.0087541.g001"></graphic></fig><table-wrap id="pone-0087541-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0087541.t001</object-id><label>Table 1</label><caption><title>Model components and parameters.</title></caption><alternatives><graphic id="pone-0087541-t001-1" xlink:href="pone.0087541.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Parameter</td><td align="left" rowspan="1" colspan="1">Interpretation</td><td align="left" rowspan="1" colspan="1">Baseline value & unit</td><td align="left" rowspan="1" colspan="1">Socio-environmental features represented in this parameter</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1"><italic>Y</italic></td><td align="left" rowspan="1" colspan="1">yield/area</td><td align="left" rowspan="1" colspan="1">21,000 kcal/ha/day</td><td align="left" rowspan="1" colspan="1">Edaphic conditions, climate, cultigens, technology</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>H</italic></td><td align="left" rowspan="1" colspan="1">longest age-specific agricultural work day</td><td align="left" rowspan="1" colspan="1">5 hr/day/ind (10 hr/day/worker)</td><td align="left" rowspan="1" colspan="1">Climate, physiological stress, socio-economic factors, socio-cultural belief and practice concerning agrarian labor</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>k</italic></td><td align="left" rowspan="1" colspan="1">area worked/hr</td><td align="left" rowspan="1" colspan="1">0.0944 ha-day/worker/hr</td><td align="left" rowspan="1" colspan="1">Topography, technology, energy subsidies</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>φ<sub>x</sub></italic></td><td align="left" rowspan="1" colspan="1">proportion of hours, <italic>H</italic>, worked by age</td><td align="left" rowspan="1" colspan="1">0≤<italic>φ<sub>x</sub></italic>≤1; <italic>φ</italic><sub>0</sub> = 0.738</td><td align="left" rowspan="1" colspan="1">Age-specific physiological and skill-related constraints on work; socio-cultural belief and practice related to labor</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>m<sub>x</sub></italic><sub>(0)</sub></td><td align="left" rowspan="1" colspan="1">maximum, age-specific fertility, <italic>E</italic>≥1</td><td align="left" rowspan="1" colspan="1">daughters born to a woman from age <italic>x</italic> to <italic>x</italic>+1</td><td align="left" rowspan="1" colspan="1">Life history context; socio-cultural belief and practice</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>p<sub>x</sub></italic><sub>(0)</sub></td><td align="left" rowspan="1" colspan="1">maximum, age-specific survival, <italic>E</italic>≥1</td><td align="left" rowspan="1" colspan="1">probability of survival from age <italic>x</italic> to <italic>x</italic>+1</td><td align="left" rowspan="1" colspan="1">Physiological and epidemiological environment; socio-cultural belief and practice</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>J</italic></td><td align="left" rowspan="1" colspan="1">Maximum, age-specific kcal requirement</td><td align="left" rowspan="1" colspan="1">2785 kcal/day</td><td align="left" rowspan="1" colspan="1">Body mass, climatic conditions, metabolic costs of maintenance and reproduction</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>ρ<sub>x</sub></italic></td><td align="left" rowspan="1" colspan="1">Proportion of consumption, <italic>J</italic>, by age</td><td align="left" rowspan="1" colspan="1">0≤<italic>ρ<sub>x</sub></italic>≤; <italic>ρ</italic><sub>0</sub> = 0.832</td><td align="left" rowspan="1" colspan="1">Age-specific nutritional needs; socio-cultural belief and practice related to food distribution</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>A<sub>m</sub></italic></td><td align="left" rowspan="1" colspan="1">Arable land available</td><td align="left" rowspan="1" colspan="1">1000 ha</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>F</italic></td><td align="left" rowspan="1" colspan="1">Fraction of <italic>A<sub>m</sub></italic> in cultivation</td><td align="left" rowspan="1" colspan="1">0≤<italic>F</italic>≤1</td><td align="left" rowspan="1" colspan="1">A function of population density and interference competition, decreasing returns to marginal land, etc.</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>α<sub>x</sub></italic></td><td align="left" rowspan="1" colspan="1">Elasticity at <italic>E</italic> = 1, at age <italic>x</italic>, for survival rate, <italic>p<sub>x</sub></italic></td><td align="left" rowspan="1" colspan="1">0.00279≤<italic>α<sub>x</sub></italic>≤0.156; <italic>α<sub>25</sub></italic> = 0.00464 at <italic>E</italic> = 1</td><td align="left" rowspan="1" colspan="1">Climate, other environmental factors affecting age-specific survival under conditions of food shortfalls from optimum; physical environment, disease, cultural practices, technology</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>γ</italic></td><td align="left" rowspan="1" colspan="1">Elasticity at <italic>E</italic> = 1, at any age for fertility, <italic>m<sub>x</sub></italic></td><td align="left" rowspan="1" colspan="1"><italic>γ</italic> = 0.135 at <italic>E</italic> = 1</td><td align="left" rowspan="1" colspan="1">Any socio-environmental factor affecting age-specific fertility under conditions of food shortfalls from optimum physical environment, disease, cultural practices, technology</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>E</italic></td><td align="left" rowspan="1" colspan="1">Food ratio</td><td align="left" rowspan="1" colspan="1"><italic>E<sub>t</sub></italic><sub> = 1</sub> = 2.99; <inline-formula><inline-graphic xlink:href="pone.0087541.e001.jpg"></inline-graphic></inline-formula> = 0.6683</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>r</italic>, <italic>b</italic>, <italic>d</italic></td><td align="left" rowspan="1" colspan="1">per capita reproductive (<italic>r</italic>), birth (<italic>b</italic>) and death (<italic>d</italic>) rate</td><td align="left" rowspan="1" colspan="1"><italic>r</italic><sub>0</sub> = 0.0176; <italic>b</italic><sub>0</sub> = 0.0369; <italic>d</italic><sub>0</sub> = 0.0192</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>TFR<sub>d</sub></italic></td><td align="left" rowspan="1" colspan="1">Total fertility rate, in daughters</td><td align="left" rowspan="1" colspan="1"><italic>TFR<sub>d</sub></italic> = 2.44</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>e</italic><sub>0</sub></td><td align="left" rowspan="1" colspan="1">Life expectancy, from birth (age 0)</td><td align="left" rowspan="1" colspan="1"><italic>e</italic><sub>0</sub> = 45 yr</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>N<sub>t</sub></italic></td><td align="left" rowspan="1" colspan="1">Population size at time <italic>t</italic></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td></tr></tbody></table></alternatives></table-wrap><p>The model tracks the age structure and size of the population (“Population age structure; total size, <italic>N</italic>”). We conceptualize all individuals as female and make adjustments to their fertility schedule in such a ways as to include an equal number of phantom males. The number of individuals in each annual age class determines “Age specific labor availability” and “Age specific food requirements.” Available labor determines total production of food according the specific environmental and socio-economic parameters that characterize the production system. Production is then summed as “Food kcal available.” Age specific food requirements likewise are summed, generating the “Baseline kcal requirement,” baseline in the sense that it is the minimum kilocalories required for all individuals in the population to achieve their optimal vital rates.</p><p>The “Food ratio (<italic>E</italic>)” characterizes our approach as one of food-limited demography. If the age-specific labor supply is large enough and environmental conditions accommodate sufficient production, then more food will be available from subsistence production than is required for the population to meet its baseline requirements and <italic>E</italic>≥1. If not, <italic>E</italic><1, and the population exists in a state of suboptimal diet measured as kcal deprivation relative to the baseline. More simply, <italic>E</italic><1 represents a state of hunger.</p><p>The operational effect of this condition is shown in the component “Age-specific fertility (<italic>m<sub>x</sub></italic>) & survivorship (<italic>l<sub>x</sub></italic>).” If <italic>E</italic>≥1, then fertility and survivorship are at their baseline maxima. As <italic>E</italic> drops below 1, food deprivation depresses fertility and elevates mortality. A gamma probability density function represents the response of fertility (<italic>γ</italic>) and age-specific survival (<italic>α<sub>x</sub></italic>) to shortfalls. These curves are parameterized from empirical data (discussed in <xref rid="pone.0087541-Lee3" ref-type="bibr">[11]</xref>) showing that sensitivity to kcal shortfalls and hunger is initially small but becomes quite large as food availability falls toward 50% of baseline level.</p><p>Continuing the cycle, age-specific fertility and mortality then have their actuarial consequences in an updated “Population age structure; total size, <italic>N</italic>,” completing an iteration of the model. The time index advances one step, a new year begins.</p><p>For ease of comprehension we have described the model as a loop of causal interactions. In fact, it is represented and iterated or solved as a family of projection matrices (see <xref rid="pone.0087541-Caswell1" ref-type="bibr">[15]</xref>), with linear algebra and occasionally the software environment of MatLab® doing the heavy work.</p><p>The components outside of the dashed rectangle are the system parameters. Where these require a specific value we have chosen figures representative of prehistoric subsistence populations engaged in agriculture (<xref ref-type="table" rid="pone-0087541-t001">Table 1</xref> and see <xref rid="pone.0087541-Lee3" ref-type="bibr">[11]</xref>). The exact value chosen is less important than the opportunity to examine the effect of that parameter varied over a range of potentially realistic values. Thus, rather than asking “What is the equilibrium population density of humans if yield (<italic>γ</italic>) = 21,000 kcal/ha/day?,” we might ask, “What is the effect of doubling yield on population density or life expectancy?” We use the model to make <italic>structural</italic> predictions <xref rid="pone.0087541-Winterhalder2" ref-type="bibr">[16]</xref> of the form, “If <italic>x</italic> gets larger, what is the direction and relative magnitude of the effect on <italic>y</italic>?”</p><p><italic>H</italic> (<xref ref-type="fig" rid="pone-0087541-g001">Figure 1</xref>) is the number of hours of productive subsistence work per day by the age class or classes with the greatest daily labor investment in producing food. The parameter <italic>φ<sub>x</sub></italic> determines the relative work investment by any age class <italic>x</italic>, scaled to 0≤<italic>φ<sub>x</sub></italic>≤1. By definition, <italic>φ<sub>x</sub></italic> = 1 for individuals working <italic>H</italic> hours; <italic>φ<sub>x</sub></italic> = 0 for complete dependents such as infants and the very elderly. Together, <italic>H</italic> and <italic>φ<sub>x</sub></italic> give us control over the age specific subsistence work effort of the population (“Age-specific labor availability”). We define <italic>φ</italic><sub>0</sub> as the sum over all age classes of each <italic>φ<sub>x</sub></italic> times the fraction of the copial population in that class. This tells us how many “full-time” equivalent workers there are per person and, as such, is a measure of the dependency ratio.</p><p>The “Age-specific food requirements” (kcals) are handled similarly. <italic>J</italic> is the requirement of the age class or classes with the greatest food needs and <italic>ρ<sub>x</sub></italic> scales the requirement for all other age classes, 0≤<italic>ρ<sub>x</sub></italic>≤1. We determine <italic>ρ<sub>x</sub></italic> according to estimates from representative modern populations; it is independent of the level of agricultural labor assigned an age class, <italic>φ<sub>x</sub></italic>. <italic>ρ</italic><sub>0</sub> measures the equivalent number of most energy costly individuals per capita in the copial phase. Food is allocated across age classes in equal proportion to need.</p><p>The next group of parameters sets the environmental conditions of production. <italic>A<sub>m</sub></italic> gives the total area (ha) of arable land available. We assume it to be of uniform quality with an average annual food yield (<italic>γ</italic>) measured in kcal/ha/day. <italic>k</italic> is a constant that converts labor time into area cultivated.</p><p><italic>F</italic> captures the difficulty of putting into production increasingly marginal plots as the population grows. At its simplest, <italic>F</italic> would be a straight line of positive slope, each new farmer bringing into cultivation the same increment of land until the habitat was fully allocated and more agriculturalists do not lead to an increase in food supply. We assume that interference competition results in a steadily decreasing efficiency of agricultural expansion by employing a continuous exponential function that is nearly linear when the population is small. Concave downward, its slope approaches an asymptote of <italic>F</italic> = 1 as <italic>N</italic> becomes large.</p><p>Finally, we use a gamma probability density function (<italic>pdf</italic>) (“Vital rate elasticities”) to determine how vital rates of age-specific survival and age-specific fertility respond to food shortfalls. The sigmoid form suggests that modest kcal shortfalls of <italic>E</italic><1 have a minimal effect, but that more significant food deprivation (0.5<<italic>E</italic><0.8) has a rapidly accelerating impact on survival and fertility. The slope of the gamma pdf determines the local <italic>elasticity</italic> or responsiveness of survival and fertility to food shortages (see <xref rid="pone.0087541-Caswell1" ref-type="bibr">[15]</xref>). The shape of the gamma pdfs has been calibrated to empirical observations <xref rid="pone.0087541-Lee3" ref-type="bibr">[11]</xref>,<xref rid="pone.0087541-Bengtsson1" ref-type="bibr">[17]</xref>. These elasticities are bio-cultural in that they are determined by socio-economic factors and choice as well as by hormonal and metabolic responses to food deprivation.</p><p>Linkages between nutrition, fertility and mortality are complex and the relationships we employ may not hold in all cases, particularly as they include behavioral components expected to vary across societies. We find that the broad conclusions of the analysis remain robust even if fertility rates are entirely unresponsive to food availability. In this case equilibrium is achieved by elevated mortality alone. Of course because fertility is a parameter in our analytical expressions, it is possible to consider the results of varying the sensitivity of fertility across a wider range than we have done here.</p><p>The food-limited framework has advantages over the Lotka-Verhulst approach. Most significant, “carrying capacity” need not be set <italic>a priori</italic>. Instead, the number of individuals at equilibrium is a function of the interaction between labor, environment, food availability and vital rates. By representing population dynamics in this way, we can examine their dependence on factors not present in the logistic framework. Similar advantages arise from deriving growth rate from underlying processes of fertility and mortality. A central variable of the model, the food ratio, drives vital demographic rates but also serves as a measure of human welfare, for which there is no analog in Lotka-Verhulst. Finally, the food-limited approach allows the consideration of any number of novel socio-economic features, such as the effects of the method and severity of taxation <xref rid="pone.0087541-Winterhalder1" ref-type="bibr">[14]</xref>, innovation that increases either yield or worker efficiency <xref rid="pone.0087541-Puleston1" ref-type="bibr">[12]</xref>, and storage <xref rid="pone.0087541-Winterhalder3" ref-type="bibr">[18]</xref>.</p><p>The most important cost of this approach: its complexity means we often must forsake exact analytical solutions for approximations and simulation. Another is that the large number of parameters require significant effort to tailor the approach to new scenarios. It could also be argued that that the model is not complex enough. It relies on numerous simplifying assumptions, and factors that do not appear may in fact be quite important. Our goal is to better understand the interactions of the features that are included, not to discount any that might be missing. We prize tractability for the analytical insights it affords us, even at the cost of some realism.</p></sec><sec id="s2b"><title>The Model as Mathematics</title><p>We solve for copial phase dynamics analytically. To examine dynamics in the remaining phases we iterate the population forward one year at a time, beginning with a small founder group and baseline rates of survival and fertility, indicating that enough food is available to avoid discounting vital rates due to hunger. The food-dependent vital rates are converted to a Leslie matrix to calculate the next year's population from the current population vector. This is equivalent to multiplying the age-specific probabilities of survival by the appropriate numbers of people in each age group going forward. Newborns are created by doing the same with the age-specific fertility rates.</p><p>In some cases we calculate an equilibrium state from initial parameters, regardless of the pathway of transitory states, using the method of comparative statics. This entails solving for the value of <italic>E</italic> at which the population equilibrates. This must be determined numerically. The point at which the population just replaces itself each year is the value of <italic>E</italic> for which the dominant eigenvalue of the food-dependent Leslie matrix is 1. The associated eigenvector is proportional to the equilibrium population structure. Sections 2 and 3 of <xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref> give more details.</p></sec></sec><sec sec-type="results" id="s3"><title>Results</title><sec id="s3a"><title>Growth to Equilibrium and Its Consequences</title><sec id="s3a1"><title>A Standard Temporal Pathway</title><p>In <xref ref-type="fig" rid="pone-0087541-g002">Figure 2</xref> a population of 20 settles in a rich but unoccupied landscape (<italic>A<sub>m</sub></italic> = 1000 ha; <italic>Y</italic> = 21,000 kcal/ha/day). Their total fertility rate (daughters only) is <italic>TFR<sub>d</sub></italic> = 2.44, and their lifespan is <italic>e</italic><sub>0</sub> = 45 years. At the time of arrival the food ratio, <italic>E</italic>, is just below 3. The population can produce three times the food required to optimize its vital rates. Consequently, the growth rate is at its baseline maximum, <italic>r</italic> = 0.0176. From this founding event, the population grows and the food ratio declines. For 346 years the birth rate remains high and undiminishing, the death rate low and not increasing. At year 347 the food ratio declines below <italic>E</italic> = 1 and the growth rate plummets as vital rates react to food shortfalls; total population growth decelerates. Forty years later, suboptimal food availability drives the net growth rate <italic>r</italic> approximately to zero by substantially degrading survival and fertility. All three variables of <xref ref-type="fig" rid="pone-0087541-g002">Figure 2</xref> settle asymptotically on the value that they will maintain indefinitely. There are <inline-formula><inline-graphic xlink:href="pone.0087541.e002.jpg"></inline-graphic></inline-formula> = 13,509 individuals, consuming <inline-formula><inline-graphic xlink:href="pone.0087541.e003.jpg"></inline-graphic></inline-formula> = 0.668 of their optimal diet and achieving an average expected lifespan <inline-formula><inline-graphic xlink:href="pone.0087541.e004.jpg"></inline-graphic></inline-formula> = 30.0 yr.</p><fig id="pone-0087541-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0087541.g002</object-id><label>Figure 2</label><caption><title>Growth in a Space-Limited Agrarian Population.</title><p>From top to bottom the graphs show concordant changes in population size (<italic>N</italic>), food availability (<italic>E</italic>) and population growth rate (<italic>r</italic>). At time zero a group of 20 people (<italic>N</italic><sub>0</sub> = 20) colonize a productive but space- and resource-limited agricultural landscape. High food availability (E>1) and an open frontier for settlement allow the population to grow at its maximum unconstrained rate for 346 years, when it becomes sufficiently large that the environment is filled and can no longer provide the amount of food energy optimal for baseline vital rates. As food availability falls below <italic>E</italic> = 1, convergence between per capita fertility and mortality rates causes a decline in growth rate (<italic>r</italic>) and the population asymptotically approaches an equilibrium. We describe this time course in terms of a <italic>copial</italic> phase, in which the population is demographically unconstrained and growing at this maximum rate; a short <italic>transitional</italic> phase, in which fertility declines and mortality increases; and, a <italic>Malthusian</italic> phase in which fertility matches mortality and growth ceases. Model parameters initially set to: <italic>TFR<sub>d</sub></italic> (measured in daughters) = 2.44; average lifespan <italic>e</italic><sub>0</sub> = 45 years; agricultural area <italic>A<sub>m</sub></italic> = 1000 ha; yield <italic>Y</italic> = 21,000 Kcal/ha/day; youths start work at age 10 and adults end working life at 65.</p></caption><graphic xlink:href="pone.0087541.g002"></graphic></fig><p>We divide this temporal path into three phases: (1) a <italic>copial phase</italic>, the period of potential surplus production between initial settlement and the point at which suboptimal food intake (hunger) is first experienced (i.e., <italic>E</italic><1); (2) a <italic>transition phase</italic>, which represents a period of increasing food-limited constraint on reproduction and survival; and (3) an <italic>equilibrium or Malthusian phase</italic>, in which food-limitations have caused rates of birth and death to converge, stabilizing population size. Separation between phases (1) and (2) is precise at <italic>E</italic> = 1. For the purposes of calculating the <italic>MTI</italic> we set the dividing line between phases (2) and (3) at a doubling time of 1000 years, the point at which growth virtually has stopped.</p></sec><sec id="s3a2"><title>Malthusian Transition Interval (MTI)</title><p>Much of our analysis concerns the effect of environmental and social factors on the relative duration of the first two of these temporal phases, distinguished as a long and bucolic period of relative plenty and a much shorter period of rapid adjustment to a final period of severe Malthusian constraints based in hunger. To focus this effort, we introduce the concept of <italic>Malthusian transition interval (MTI)</italic> to characterize the transition phase.</p><p>Consider this scenario. A small population arrives as first colonists in a region with well-defined geographical boundaries (e.g., an island or valley). The climate, the state of its soil, and the technology, skills and cultivars available to the population are propitious. The founding group begins with a potential surplus of food (<inline-formula><inline-graphic xlink:href="pone.0087541.e005.jpg"></inline-graphic></inline-formula>>1). The population grows exponentially at a constant rate (<italic>r</italic><sub>0</sub>) determined by its baseline fertility and mortality. Its frontier of unused agricultural land diminishes. Interference competition begins to constrain per capita access to land, captured by the flattening slope of <italic>F</italic> (the fraction of arable land that is in cultivation). Surplus declines to and then drops below <italic>E</italic> = 1. As per capita rates of death and birth converge due to food shortfalls, population growth stops.</p><p>The start of the transition phase has a sudden and palpable impact on human quality of life (<xref ref-type="fig" rid="pone-0087541-g003">Figure 3</xref>). Life expectancy is diminished and infant-youth mortality is elevated. Total fertility rate falls to replacement. These changes likely would be unexpected and they happen quickly. Total fertility rate, average life span and probability of survival to age five are high and constant throughout the copial phase, but once <italic>E</italic> drops below 1 they plummet toward their equilibrium values. In model analyses with realistic parameters the copial phase typically lasts between 300 and 400 years (12–16, 25-year generations). The transitional phase and <italic>MTI</italic> are much shorter, often less than 50–60 years (two-three generations).</p><fig id="pone-0087541-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0087541.g003</object-id><label>Figure 3</label><caption><title>Malthusian transition interval (<italic>MTI</italic>) and Quality of Life Demographic Variables.</title><p>(A) Food availability declines from <italic>E</italic> = 2.99, to <italic>E</italic> = 1 at year 347. It stabilizes at <italic>E</italic> = 0.668 approximately 51 years later. Coincident with the beginning of the transition phase: (B) total fertility rate (<italic>TFR<sub>d</sub></italic>) declines sharply from 2.44 to replacement, 2.10; (C) average life span falls from <italic>e</italic><sub>0</sub> = 45 to <inline-formula><inline-graphic xlink:href="pone.0087541.e006.jpg"></inline-graphic></inline-formula> = 30.0 years; and (D) the probability of a youth surviving to age five falls from 0.773 to 0.652. Although the production surplus declines continuously and gradually, Malthusian constraints impose themselves abruptly as variables representing demographic quality of life plummet.</p></caption><graphic xlink:href="pone.0087541.g003"></graphic></fig><p><italic>MTI</italic> highlights the implications of these two time frames. Density dependence strikes with little warning. It sharply and quickly curtails human quality of life by elevating hunger and infant mortality, and by lowering life expectancy and fertility. Between the initial signal of a problem and its full impact, the opportunity for adaptation to avert the full force of the Malthusian outcome is short. We believe this has implications for the effects of population ecology on socio-economic history.</p></sec></sec><sec id="s3b"><title>Factors Affecting Phase Duration</title><p>The effects of environment, work, consumption, and demography on the duration of the copial and transitional phases of population growth are not easy to predict. Generally this transition occurs rapidly, but in a narrow range of specific circumstances it can be prolonged. We examine the consequences of altering initial colony size (<italic>N<sub>t</sub></italic><sub> = 0</sub>), the area of the habitat being settled <italic>A<sub>m</sub></italic>, environmental yield (<italic>Y</italic>), baseline life expectancy (<italic>e</italic><sub>0</sub>), total fertility rate (<italic>TFR<sub>d</sub></italic>), and the elasticity of fertility and age-specific survival at <italic>E</italic> = 1 (<italic>γ</italic>, <italic>α<sub>x</sub></italic>) (see also supplemental materials in <xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref>).</p><sec id="s3b1"><title>Founding Population Size</title><p>A larger initial settlement group shortens the copial phase; however, the duration of the transitional phase or <italic>MTI</italic> is unaffected. We demonstrate this and subsequent results in two formats. <xref ref-type="fig" rid="pone-0087541-g004">Figure 4</xref> shows it in terms of the vital rates that drive population growth; Table S1a (<xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref>; Section 1) shows the duration of the copial phase declines with increasing size of the founding group, whereas the duration of the transitional phase is constant over <italic>N</italic><sub>0</sub>. During the copial phase a constant per capita birth rate exceeds a constant per capita death rate by a wide margin (<xref ref-type="fig" rid="pone-0087541-g004">Figure 4</xref>); <italic>r</italic> is relatively large. As <italic>E</italic> drops below 1, initiating the transitional phase, these rates begin rapidly to converge. They overshoot slightly and then settle on the equalized value that determines replacement without growth, the Malthusian equilibrium. Per capita birth rate declines from 0.037 to 0.033 (a drop of 11%), whereas per capita death rate climbs from 0.019 to 0.033 (an increase of 71%). Our parameterization of vital rate elasticities often results in mortality being more responsive to food deprivation than is fertility. <inline-formula><inline-graphic xlink:href="pone.0087541.e007.jpg"></inline-graphic></inline-formula> = 13,509, life expectancy <inline-formula><inline-graphic xlink:href="pone.0087541.e008.jpg"></inline-graphic></inline-formula> has fallen from 45.0 to 30.0, and <inline-formula><inline-graphic xlink:href="pone.0087541.e009.jpg"></inline-graphic></inline-formula> = 0.668, indicating that at equilibrium the population is consuming 67% of its demographically optimal diet.</p><fig id="pone-0087541-g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0087541.g004</object-id><label>Figure 4</label><caption><title>Effect of Settlement Population Size <italic>N</italic><sub>0</sub> on the Time Course of Vital Rates.</title><p>Increasing the size of the founding group by factors of 10, from 2 to 20 to 200, decreases the length of the copial phase by a constant 131 years each step, a log-linear relationship. During the copial phase per capita births and death rates, hence the reproductive rate, <italic>r</italic>, are unchanging. The transition phase in which birth and death rates converge (at <italic>b</italic> = <italic>d</italic> = 0.033) arrives sooner with larger settlement, but it otherwise is of the same form and duration for each scenario. At Malthusian equilibrium, the population numbers <inline-formula><inline-graphic xlink:href="pone.0087541.e010.jpg"></inline-graphic></inline-formula> = 13,509 individuals, who have a life expectancy of <inline-formula><inline-graphic xlink:href="pone.0087541.e011.jpg"></inline-graphic></inline-formula> = 30.0 yrs; the food ratio is <inline-formula><inline-graphic xlink:href="pone.0087541.e012.jpg"></inline-graphic></inline-formula> = 0.668.</p></caption><graphic xlink:href="pone.0087541.g004"></graphic></fig></sec><sec id="s3b2"><title>Area of Habitat Available for Settlement</title><p>As the available habitat grows larger, the copial phase increases in duration. The duration of the relatively short transitional phase is unaffected by area (<xref ref-type="fig" rid="pone-0087541-g005">Figure 5</xref>; <xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref> Table S1b). Increasing the agricultural land available (<italic>A<sub>m</sub></italic>) reduces the initial the density of farmers, lowering interference competition among them. Populations able to expand over more arable land will take longer to feel the density dependent effects of hunger and diminishing returns on labor. As a result, high marginal productivities are prolonged and the decline of surplus (<italic>E</italic>≥1) is slowed. <italic>MTI</italic> however is independent of the amount of arable land (<xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref> Table S1b). The point at which <italic>E</italic>>1 is associated with a particular population <italic>density</italic> and the effect of density on the timing of transition phase approach to equilibrium is invariant across different sized habitats.</p><fig id="pone-0087541-g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0087541.g005</object-id><label>Figure 5</label><caption><title>Effect of Agricultural Area (<italic>A</italic><sub>m</sub>) on the Time Course of Vital Rates.</title><p>Increasing the maximum area available for agricultural in-fill (<italic>A</italic><sub>m</sub>) by factors of 10, from 100 to 1000 to 10,000 ha, increases the length of the copial phase by a constant 131 years each increment, a log-linear relationship. Habitat area has no effect on transition phase length; the passage time from E = 1 to the neighborhood of the Malthusian equilibrium is identical in all cases. At equilibrium, food availability (<inline-formula><inline-graphic xlink:href="pone.0087541.e013.jpg"></inline-graphic></inline-formula> = 0.668) and lifespan (<inline-formula><inline-graphic xlink:href="pone.0087541.e014.jpg"></inline-graphic></inline-formula> = 30) are unaffected by area. By contrast, total population size increases proportionally: <inline-formula><inline-graphic xlink:href="pone.0087541.e015.jpg"></inline-graphic></inline-formula> = 1351; <inline-formula><inline-graphic xlink:href="pone.0087541.e016.jpg"></inline-graphic></inline-formula> = 13,509; <inline-formula><inline-graphic xlink:href="pone.0087541.e017.jpg"></inline-graphic></inline-formula> = 135,090.</p></caption><graphic xlink:href="pone.0087541.g005"></graphic></fig><p>Population grows at a constant exponential rate in the copial phase, ensuring that the effect of agricultural area on in-fill is log-linear. For our parameter set, increasing arable area by a factor of ten (from 100 ha to 1000 ha) adds 131 years to the copial phase. A second 10-fold increase adds another 131 years. The same log-linear effect is evident for founding group size (<xref ref-type="fig" rid="pone-0087541-g004">Figure 4</xref>).</p></sec><sec id="s3b3"><title>Habitat Yield</title><p>Increasing the yield or productivity of the environment independent of its area lengthens the copial phase and it shortens the transitional phase (<xref ref-type="fig" rid="pone-0087541-g006">Figure 6</xref>; <xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref> Table S1c). The extension of the copial phase is a result of the greater capacity of the system to produce a food surplus (i.e., a greater <italic>E<sub>t</sub></italic><sub> = 0</sub>). Maximum growth rates can be sustained for a long period and the size, and hence density, of the population will be larger at the point that <italic>E</italic> first falls below 1. Thus even though the equilibrium population size increases approximately linearly with yield <xref rid="pone.0087541-Puleston1" ref-type="bibr">[12]</xref>, the <italic>MTI</italic> is an expression of the exponential growth equation (see <xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref> Equation S3.8). A richer yield per hectare increases the maximum value of <italic>E</italic>, and as a consequence, populations in high-yield environments pass through the transitional phase more quickly.</p><fig id="pone-0087541-g006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0087541.g006</object-id><label>Figure 6</label><caption><title>Effect of Yield (<italic>Y</italic>) on the Time Course of Vital Rates.</title><p>The solid, low-yield, curves represent 7000 kcal/ha/day; the dashed, moderate yield curves a 3-fold increase to 21,000 kcal/ha/day; and, the dotted, high-yield curves another 3-fold increase to 63,000 kcal/ha/day. Because the low yield situation is below adequate kcal intake at settlement, there is no copial phase as fertility falls and mortality rises smoothly with further declines in food availability. The approach to equilibrium is slow, the transition phase extended. The moderate- and high-yield cases allow the population to grow at its maximum, constant rate until food availability drops below <italic>E</italic> = 1. They exhibit a lengthy copial phase. The moderate-yield population experiences hunger first, but in each case the approach to equilibrium is rapid once <italic>E</italic> has fallen below 1. At equilibrium, <inline-formula><inline-graphic xlink:href="pone.0087541.e018.jpg"></inline-graphic></inline-formula> = 0.668 and <inline-formula><inline-graphic xlink:href="pone.0087541.e019.jpg"></inline-graphic></inline-formula> = 30.0 for each level of yield. Equilibrium population density, however, increases with yield; <inline-formula><inline-graphic xlink:href="pone.0087541.e020.jpg"></inline-graphic></inline-formula> = 2799; <inline-formula><inline-graphic xlink:href="pone.0087541.e021.jpg"></inline-graphic></inline-formula> = 13,509; <inline-formula><inline-graphic xlink:href="pone.0087541.e022.jpg"></inline-graphic></inline-formula> = 27,287.</p></caption><graphic xlink:href="pone.0087541.g006"></graphic></fig><p>In a very low yield environment, a population experiences hunger early and at low density, with the result that it grows more slowly toward its equilibrium size. Settings of poor agricultural yields are unique in that <italic>MTI</italic> may be relatively long and may exceed the copial phase in duration (<xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref> Table S1c).</p><p>These relationships are evident in vital rates (<xref ref-type="fig" rid="pone-0087541-g006">Figure 6</xref>). Low yielding environments may have no copial phase or a quite short one matched to a relatively long and gradual transitional phase witnessing the slow convergence of birth and death rates. Moderate to high yielding environments revert to the more common pattern in which there is a copial phase significantly longer than the transitional one and, although delayed, the transition happens with more typical abruptness. Intuition might suggest that adding arable land to a society's territory would have an effect similar to that of increasing the yield of a smaller area. This is not the case. While larger area (<italic>A<sub>m</sub></italic>) and yield (<italic>Y</italic>) both increase copial phase duration, only the latter has an effect on the transitional phase. The analytical approximation of <italic>MTI</italic> (<xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref> Equation S3.8) tracks the trajectory of population density, which is a function of yield, but not area.</p><p>Comparison of <xref ref-type="fig" rid="pone-0087541-g004">Figures 4</xref>, <xref ref-type="fig" rid="pone-0087541-g005">5</xref> and <xref ref-type="fig" rid="pone-0087541-g006">6</xref> shows that births and deaths equilibrate at the same value (0.033) for all treatments involving these three parameters (<italic>N<sub>t</sub></italic><sub> = 0</sub>, <italic>A<sub>m</sub></italic> & <italic>Y</italic>). Life expectancy at equilibrium likewise is the same across all cases, <inline-formula><inline-graphic xlink:href="pone.0087541.e023.jpg"></inline-graphic></inline-formula> = 30.0, as is the hunger level, <inline-formula><inline-graphic xlink:href="pone.0087541.e024.jpg"></inline-graphic></inline-formula> = 0.668. The size of the population at equilibrium does not change as a result of different settlement numbers (<italic>N<sub>t</sub></italic><sub> = 0</sub>); however, it increases as a function of area (<italic>A<sub>m</sub></italic>) and yield (<italic>Y</italic>). The Malthusian phase life history experience of individuals is unaffected by the size of the founding group, the area settled or the yield of the environment, parameters which affect the manner in which the equilibrium is approached but not its demographic or nutritional character. Among the other inferences possible, high yields do not improve quality of life at Malthusian equilibrium.</p></sec><sec id="s3b4"><title>Background Mortality</title><p>From settlement and environmental factors we turn to vital rate parameters, the first of which is density independent, background mortality, implemented as an increase in the baseline expected lifespan. Increasing background mortality causes a pronounced increase in the length of the copial phase (<xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref> Table S1d). In an environment in which even a well-fed individual (<italic>E</italic>≥1) lives only 30 years on average, the effective growth rate is so small that it takes 1204 years to fill the habitat to <italic>E</italic> = 1 (<xref ref-type="fig" rid="pone-0087541-g007">Figure 7</xref>). If that individual expects to live to 45 or even 60 years on average, the population grows rapidly enough that it reaches <italic>E</italic> = 1 within 347 or 237 years, respectively. Similarly, as background mortality declines, so does <italic>MTI</italic>. Populations in environments we would consider propitious will, once they experience hunger, approach their Malthusian destiny more rapidly, with less time to react. <italic>MTI</italic> = 65 years for a baseline expected lifespan of <italic>e</italic><sub>0</sub> = 30 and drops to 48 years for baseline <italic>e</italic><sub>0</sub> = 60.</p><fig id="pone-0087541-g007" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0087541.g007</object-id><label>Figure 7</label><caption><title>Effect of Background Mortality, Measured as Life Expectancy (<italic>e</italic><sub>0</sub>), on the Time Course of Vital Rates.</title><p>The solid lines represent high background mortality corresponding to a short life expectancy of <italic>e</italic><sub>0</sub> = 30 years. The dashed curves depict a medium mortality, baseline life expectancy of <italic>e</italic><sub>0</sub> = 45 years; the dotted lines low background mortality, baseline expected lifespan of <italic>e</italic><sub>0</sub> = 60 years. In the high-mortality environment, the maximum growth rate is small and the copial phase lasts a lengthy 1204 years. Mortality starts high and it changes very little at equilibrium (<inline-formula><inline-graphic xlink:href="pone.0087541.e025.jpg"></inline-graphic></inline-formula> = 26.8) because equilibrium food availability is fairly high (<inline-formula><inline-graphic xlink:href="pone.0087541.e026.jpg"></inline-graphic></inline-formula> = 0.851). By contrast, the fast-growing, low-background-mortality population experiences a short copial phase of 237 years. Once it drops into hunger (<italic>E</italic> = 1) it approaches its Malthusian equilibrium very quickly, experiencing a large rise in mortality. At equilibrium, expected lifespan is longest in the most favorable environment (<inline-formula><inline-graphic xlink:href="pone.0087541.e027.jpg"></inline-graphic></inline-formula> = 34.0, having however dropped from <italic>e</italic><sub>0</sub> = 60). Hunger is acute in this case. In numbers: High mortality, Harsh Env.: <inline-formula><inline-graphic xlink:href="pone.0087541.e028.jpg"></inline-graphic></inline-formula> = 0.851; <inline-formula><inline-graphic xlink:href="pone.0087541.e029.jpg"></inline-graphic></inline-formula> = 26.8; <inline-formula><inline-graphic xlink:href="pone.0087541.e030.jpg"></inline-graphic></inline-formula> = 10,441. Moderate Mortality Env.: <inline-formula><inline-graphic xlink:href="pone.0087541.e031.jpg"></inline-graphic></inline-formula> = 0.668; <inline-formula><inline-graphic xlink:href="pone.0087541.e032.jpg"></inline-graphic></inline-formula> = 30.0; <inline-formula><inline-graphic xlink:href="pone.0087541.e033.jpg"></inline-graphic></inline-formula> = 13,509. Low Mortality, Propitious Env.: <inline-formula><inline-graphic xlink:href="pone.0087541.e034.jpg"></inline-graphic></inline-formula> = 0.562; <inline-formula><inline-graphic xlink:href="pone.0087541.e035.jpg"></inline-graphic></inline-formula> = 34.0; <inline-formula><inline-graphic xlink:href="pone.0087541.e036.jpg"></inline-graphic></inline-formula> = 16,075.</p></caption><graphic xlink:href="pone.0087541.g007"></graphic></fig><p>Expected lifespan is the first parameter we have examined for which the equilibrium values for birth and death rate and, as a consequence, other demographic quality of life variables differ among treatments. When background mortality is high, the equilibrium population reaches a density of <inline-formula><inline-graphic xlink:href="pone.0087541.e037.jpg"></inline-graphic></inline-formula> = 10,441, with an equilibrium life expectancy of <inline-formula><inline-graphic xlink:href="pone.0087541.e038.jpg"></inline-graphic></inline-formula> = 26.8 years, and a relatively low level of hunger, consuming 85.1% of its baseline kilocalorie diet. By contrast, when background mortality is low, the equilibrium population is much larger, 16,075. While it has a longer equilibrium expected lifespan (<inline-formula><inline-graphic xlink:href="pone.0087541.e039.jpg"></inline-graphic></inline-formula> = 34.0), it suffers hunger at the alarming level of 56.2% of its baseline diet. High background mortality results in a smaller population with shortened expected lifespans but less hunger at equilibrium for those individuals who survive their otherwise harsh conditions.</p><p>These results demonstrate that equilibrium conditions are set by the combination of density independent and density dependent mortality, offsetting fertility. As density independent mortality increases it takes less food-induced impedance to sustain the Malthusian equilibrium. Lessened background mortality generates a population abundant in numbers but poorly fed. Population control by strictly Malthusian means is a punishing situation either way.</p></sec><sec id="s3b5"><title>Fertility</title><p>The effect of different per capita fertility rates mirrors that of death rates, an expected result in that <italic>r</italic> = <italic>d</italic> – <italic>b</italic> (<xref ref-type="fig" rid="pone-0087541-g008">Figure 8</xref>; <xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref> Table S1e). Higher total fertility (<italic>TFR<sub>d</sub></italic>) results in a sharp decline in the duration of the copial phase and a less evident decline in the much shorter, transitional phase. High fertility populations approach the point of hunger quickly; they likewise move from <italic>E</italic> = 1 toward equilibrium at a faster rate than a population with lower baseline fertility.</p><fig id="pone-0087541-g008" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0087541.g008</object-id><label>Figure 8</label><caption><title>Effect of Total Fertility Rate (<italic>TFR<sub>d</sub></italic>) on the Time Course of Vital Rates.</title><p>Assessing <italic>TFR<sub>d</sub></italic> at <italic>E</italic> = 1, the low-fertility population's <italic>TFR<sub>d</sub></italic> is 2 daughters per mother (solid line), the moderate <italic>TFR<sub>d</sub></italic> is 3 (dashed line) and the high fertility population has a <italic>TFR<sub>d</sub></italic> of 4 (dotted line). Mortality tends to be more responsive than fertility in our parameterization of the model, so when fertility is high the gap between the vital rates is closed at equilibrium by a marked rise in morality. The result is higher fertility <italic>and</italic> mortality at equilibrium than in the other scenarios. Generally, lower maximum fertility also means lower per capita mortality at equilibrium, as shown by the per capita rate at which morality and fertility converge. The high-fertility population also grows quickly; it thus has a shorter copial phase, experiences hunger sooner, and has a steeper, more abrupt approach to equilibrium. It is the least well-fed at equilibrium and has a very short life expectancy. In numbers: High fertility population: <inline-formula><inline-graphic xlink:href="pone.0087541.e040.jpg"></inline-graphic></inline-formula> = 0.579; <inline-formula><inline-graphic xlink:href="pone.0087541.e041.jpg"></inline-graphic></inline-formula> = 20.4; <inline-formula><inline-graphic xlink:href="pone.0087541.e042.jpg"></inline-graphic></inline-formula> = 15,956. Moderate fertility population: <inline-formula><inline-graphic xlink:href="pone.0087541.e043.jpg"></inline-graphic></inline-formula> = 0.624; <inline-formula><inline-graphic xlink:href="pone.0087541.e044.jpg"></inline-graphic></inline-formula> = 25.4; <inline-formula><inline-graphic xlink:href="pone.0087541.e045.jpg"></inline-graphic></inline-formula> = 14,601. Low fertility population: <inline-formula><inline-graphic xlink:href="pone.0087541.e046.jpg"></inline-graphic></inline-formula> = 0.731; <inline-formula><inline-graphic xlink:href="pone.0087541.e047.jpg"></inline-graphic></inline-formula> = 35.3; <inline-formula><inline-graphic xlink:href="pone.0087541.e048.jpg"></inline-graphic></inline-formula> = 12,247.</p></caption><graphic xlink:href="pone.0087541.g008"></graphic></fig><p>A low initial fertility rate, set at a <italic>TFR<sub>d</sub></italic> of 2 daughters per woman, corresponds to greater food availability at equilibrium (<inline-formula><inline-graphic xlink:href="pone.0087541.e049.jpg"></inline-graphic></inline-formula> = 0.731), a longer average life span (<inline-formula><inline-graphic xlink:href="pone.0087541.e050.jpg"></inline-graphic></inline-formula> = 35.3), and a smaller population of <inline-formula><inline-graphic xlink:href="pone.0087541.e051.jpg"></inline-graphic></inline-formula> = 12,247 (<xref ref-type="fig" rid="pone-0087541-g008">Figure 8</xref>). Malthusian equilibrium is achieved by a relatively small increase in density dependent, per capita death rates (0.0199 to 0.0279) matched to a lesser drop in birth rates (0.0303 to 0.0279). By contrast, for the high fertility population (<italic>TFR<sub>d</sub></italic> = 4), the Malthusian state is one of greater hunger (<inline-formula><inline-graphic xlink:href="pone.0087541.e052.jpg"></inline-graphic></inline-formula> = 0.579), shorter expected life span (<inline-formula><inline-graphic xlink:href="pone.0087541.e053.jpg"></inline-graphic></inline-formula> = 20.4) and a larger population (<inline-formula><inline-graphic xlink:href="pone.0087541.e054.jpg"></inline-graphic></inline-formula> = 15,956). Per capita death rate at equilibrium has increased by a substantial 2.5 times, from 0.0194 to 0.0479, again reflecting the brutal mechanics of Malthusian controls.</p></sec><sec id="s3b6"><title>Elasticities of Survival and Fertility</title><p>Finally, we consider the responsiveness of age-specific vital rates to shortfalls in food availability (<xref ref-type="fig" rid="pone-0087541-g009">Figure 9</xref>; <xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref> Table S1f). Because density dependence operates only if <italic>E</italic><1, the elasticity of vital rates has no effect on the duration of the copial phase. However, increasing the elasticity of vital rates to food shortfalls shortens the transition phase (<xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref> Table S1f). We portray this result by multiplying elasticities in the vicinity of <italic>E</italic> = 1 by values from 0.5 to 2, a procedure that captures increasing sensitivity to a unit decline in food availability. Vital rate elasticity can be conceptualized as determining which of the per capita vital rates, death or birth, converges most quickly on the other and, consequently, which makes the greatest contribution to population equilibrium. Keeping in mind that these rates are falling, greater elasticity of survival rates enhances the impact of mortality (per capita death rate) and greater elasticity of fertility rates enhances the impact of fertility (<xref ref-type="fig" rid="pone-0087541-g009">Figure 9A, 9B</xref>) on the achieved equilibrium. For both parameters, this entails a more rapid transition to Malthusian equilibrium. Note, however, that the 4-fold change in vital rate elasticities depicted in <xref ref-type="fig" rid="pone-0087541-g009">Figure 9</xref> has a fairly modest impact on equilibrium values for the food ratio (<inline-formula><inline-graphic xlink:href="pone.0087541.e055.jpg"></inline-graphic></inline-formula>), life expectancy (<inline-formula><inline-graphic xlink:href="pone.0087541.e056.jpg"></inline-graphic></inline-formula>), and population size (<inline-formula><inline-graphic xlink:href="pone.0087541.e057.jpg"></inline-graphic></inline-formula>).</p><fig id="pone-0087541-g009" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0087541.g009</object-id><label>Figure 9</label><caption><title>Effect of Elasticity on the Time Course of Vital Rates.</title><p>(A) Elasticity of survival (mortality); (B) Elasticity of fertility. Elasticities play no role in the duration of the copial phase (much of which is elided in the graphs); they increase sensitivity of vital rates to hunger, shortening the duration of the transitional phase. Increasing the elasticity of mortality at <italic>E</italic> = 1 from half its empirical value (labeled “low” in the graphs) to its empirical value (normal) to double its empirical value (high) causes mortality to climb more quickly once <italic>E</italic> drops below 1, resulting in relatively high equilibration of per capita birth and death rates. Higher elasticity of fertility at <italic>E</italic> = 1 causes birth rates to drop more quickly with births and deaths. High elasticities are associated at equilibrium with more adequate Kcal intake, longer lifespan and lower population size. In numbers: Low survival elasticity: <inline-formula><inline-graphic xlink:href="pone.0087541.e058.jpg"></inline-graphic></inline-formula> = 0.622; <inline-formula><inline-graphic xlink:href="pone.0087541.e059.jpg"></inline-graphic></inline-formula> = 31.8; <inline-formula><inline-graphic xlink:href="pone.0087541.e060.jpg"></inline-graphic></inline-formula> = 14,530. Normal survival elasticity: <inline-formula><inline-graphic xlink:href="pone.0087541.e061.jpg"></inline-graphic></inline-formula> = 0.668; <inline-formula><inline-graphic xlink:href="pone.0087541.e062.jpg"></inline-graphic></inline-formula> = 30.0; <inline-formula><inline-graphic xlink:href="pone.0087541.e063.jpg"></inline-graphic></inline-formula> = 13,509. High survival elasticity: <inline-formula><inline-graphic xlink:href="pone.0087541.e064.jpg"></inline-graphic></inline-formula> = 0.721; <inline-formula><inline-graphic xlink:href="pone.0087541.e065.jpg"></inline-graphic></inline-formula> = 28.5; <inline-formula><inline-graphic xlink:href="pone.0087541.e066.jpg"></inline-graphic></inline-formula> = 12,478. Low fertility elasticity: <inline-formula><inline-graphic xlink:href="pone.0087541.e067.jpg"></inline-graphic></inline-formula> = 0.652; <inline-formula><inline-graphic xlink:href="pone.0087541.e068.jpg"></inline-graphic></inline-formula> = 28.3; <inline-formula><inline-graphic xlink:href="pone.0087541.e069.jpg"></inline-graphic></inline-formula> = 13,894. Normal fertility elasticity: <inline-formula><inline-graphic xlink:href="pone.0087541.e070.jpg"></inline-graphic></inline-formula> = 0.668; <inline-formula><inline-graphic xlink:href="pone.0087541.e071.jpg"></inline-graphic></inline-formula> = 30.0; <inline-formula><inline-graphic xlink:href="pone.0087541.e072.jpg"></inline-graphic></inline-formula> = 13,509. High fertility elasticity: <inline-formula><inline-graphic xlink:href="pone.0087541.e073.jpg"></inline-graphic></inline-formula> = 0.694; <inline-formula><inline-graphic xlink:href="pone.0087541.e074.jpg"></inline-graphic></inline-formula> = 32.3; <inline-formula><inline-graphic xlink:href="pone.0087541.e075.jpg"></inline-graphic></inline-formula> = 12,959.</p></caption><graphic xlink:href="pone.0087541.g009"></graphic></fig></sec></sec></sec><sec sec-type="discussion" id="s4"><title>Discussion</title><sec id="s4a"><title>Resource-limited Population Phases and Their Determinants</title><p>Naturalistic population growth in a space-limited agrarian population is characterized by copial, transitional and Malthusian phases. We introduced the term <italic>Malthusian transition interval (MTI)</italic> to characterize the brevity of the transition phase and its singular consequence: a rapid shift from full diet and positive demographic indicators to hunger and depressed demographic indicators (<xref ref-type="fig" rid="pone-0087541-g003">Figure 3</xref>). Malthusian constraints arrive abruptly, with little warning; once initiated, they settle onto a population rapidly. Their impact on quality of life variables can be quite severe; <italic>ceteris paribus</italic>, they last indefinitely. Ironically, the better the copial phase conditions that precede transition to Malthusian equilibrium – e.g., high fertility, low background mortality (longer lifespan) – the quicker and grimmer will be the plunge to the Malthusian condition when it arrives. When the pressure on basic demographic adjustment is greatest, the time available to respond adaptively is shortest.</p><p>Passage times through the copial and transitional phases depend on environmental and demographic parameters (<xref ref-type="supplementary-material" rid="pone.0087541.s001">File S1</xref> Table S1). Increased size of the founding group <italic>N<sub>t</sub></italic><sub> = 0</sub> shortens the copial phase; increased area available for settlement <italic>A<sub>m</sub></italic> lengthens it. In both cases the relationship is log-linear. Neither variable, however, affects the duration of the transition phase. Yield (<italic>Y</italic>), has opposing effects. Very low yield may eliminate the copial phase altogether; increasing yield lengthens the copial phase but it shortens the transitional phase. Longer baseline expected lifespan (<italic>e</italic><sub>0</sub>) shortens both the copial and transitional phases, the former to a much greater degree than the latter. Greater baseline fertility (<italic>TFR<sub>d</sub></italic>) has the same effect. Finally, increasing the elasticity of fertility and mortality in response to food shortfalls has no effect on the copial phase, but shortens the transition phase.</p><p>In all circumstances except that of very low yield (<italic>Y</italic>), the duration of the copial phase exceeds that of the transition phase by a wide margin. Realistic gamma function parameters for the elasticity of survival and fertility produce another result: for all of our variables and scenarios, increasing mortality contributes more to the convergence of vital rates and thus to equilibrium than does declining fertility (<xref ref-type="fig" rid="pone-0087541-g004">Figures 4</xref> through <xref ref-type="fig" rid="pone-0087541-g009">9</xref>). This is a direct reflection of the fact that elasticities of survival are larger than those of fertility near <italic>E</italic> = 1 <xref rid="pone.0087541-Lee3" ref-type="bibr">[11]</xref>.</p></sec><sec id="s4b"><title>Resource-Limited Population Ecology</title><p>Most previous attempts to define population growth rates in terms of resources either lack demographic complexity or are not developed into full dynamic models (e.g., <xref rid="pone.0087541-Bengtsson1" ref-type="bibr">[17]</xref>,<xref rid="pone.0087541-Caughley1" ref-type="bibr">[19]</xref>,<xref rid="pone.0087541-Read1" ref-type="bibr">[20]</xref>,<xref rid="pone.0087541-Scott1" ref-type="bibr">[21]</xref>,<xref rid="pone.0087541-Winterhalder4" ref-type="bibr">[22]</xref>). We seek more complex dynamics from the bottom up, by defining explicit links among environment, labor capacity, caloric requirements, competition, food availability and vital rates (<xref ref-type="fig" rid="pone-0087541-g001">Figure 1</xref>). Other models assume that there is a single “density dependent” vital rate (e.g., <xref rid="pone.0087541-Neubert1" ref-type="bibr">[23]</xref>), that vital rates respond in sequence from most to least sensitive <xref rid="pone.0087541-Gaillard1" ref-type="bibr">[24]</xref>, or that the population is made up of a single age class and the overall growth rate is the food-dependent variable (e.g., <xref rid="pone.0087541-Winterhalder4" ref-type="bibr">[22]</xref>). Our framework attains greater realism through the use of elasticities to quantify a proportional response of age-specific vital rates to food deficit.</p><p>In general, the literature on the population ecology of pre-modern societies takes as its inspiration the logistic model of Lotka-Verhulst <xref rid="pone.0087541-Hutchinson1" ref-type="bibr">[6]</xref>. The logistic calculates impedance to an underlying exponential growth rate, assuming that the growth rate slows from the very moment that population increase begins. The factor damping growth is either linear, as in Verhulst's original formulation, or curved, as in the theta logistic <xref rid="pone.0087541-Tsoularis1" ref-type="bibr">[25]</xref>. Equilibrium occurs when the population arrives at a carrying capacity assigned <italic>a priori</italic>. This pattern of impedance and point of carrying capacity are either deduced from first principles or derived statistically from fitting the model to empirical observations <xref rid="pone.0087541-Sibly1" ref-type="bibr">[26]</xref>.</p><p>We address two shortcomings of the Lotka-Verhulst approach. First, because the impact of the Lotka-Verhulst feedback term is gradual and is felt from the very beginning of growth, it necessarily entails a smooth approach to equilibrium over the whole course of a population's trajectory. There is nothing like a copial phase of unchanging maximal growth, no possibility of a prolonged surplus, and no sharp discontinuities in growth rates (cf. <xref ref-type="fig" rid="pone-0087541-g002">Figures 2</xref>, <xref ref-type="fig" rid="pone-0087541-g003">3</xref>), in short, no abruptness to the dynamics. The Lotka-Verhulst approach produces dynamic gradualism; our model allows for a more punctuated form of population history.</p><p>The second shortcoming is related. In the Lotka-Verhulst model, both <italic>r</italic> and carrying capacity (typically, <italic>K</italic>) are set as initial parameters. In our approach, these two variables are outcomes of lower-level variables, some of them easier to observe and assess empirically. This potentially increases the empirical reliability and the analytical flexibility of the model. Instead of a single value for the intrinsic rate of increase, <italic>r</italic>, we are able to derive growth from birth and death rates over age structure, represented as elastic responses to food intake. Instead of a single measure of carrying capacity, we build up to density dependent resource limitations from estimates of agricultural area, yield, age-specific work effort, and interference competition.</p></sec><sec id="s4c"><title>Social Evolution and the Question of State Collapse</title><p>Hypotheses that population has had a causal role in social evolution are prominent in both the popular <xref rid="pone.0087541-Diamond1" ref-type="bibr">[27]</xref> and the scientific literature <xref rid="pone.0087541-Tainter1" ref-type="bibr">[28]</xref>,<xref rid="pone.0087541-Middleton1" ref-type="bibr">[29]</xref>,<xref rid="pone.0087541-TurnerII1" ref-type="bibr">[30]</xref>. Population pressure is seen to compel new adaptive forms like agricultural production or state-level political organization, or it is seen to exceed the adaptive capacity of societies and lead to their decline or collapse due to environmental degradation or peasant revolt. Population competes with other explanatory candidates, among them environmental change, economic innovation, warfare, ideological development and the agency of individuals seeking prestige. Unifying frameworks are rare but needed, given the complexity of human society and the many factors that reasonably have been postulated to cause early civilizations to rise, persist and fail. By including in a single systematic framework, environment, labor productivity, nutritional requirements, age structure and vital rate demography our approach seeks to provide a more synthetic framework.</p><p>Interactions among population, labor and food availability are often implicated in collapse narratives, but the models on which these arguments are based may not be sufficient to the job. We focus here on the Maya case, because it is one of the best known cases in which advanced agrarian political systems failed. For example, Culbert <xref rid="pone.0087541-Culbert1" ref-type="bibr">[31]</xref> proposes that Classic Maya political society collapsed because of a labor shortage. Although approaching high levels of population, Maya society was vulnerable to diversion of workers from agriculture to warfare and the construction of monumental architecture. By contrast, a population ecology approach suggests that when the population is large and has limited ability to expand into agricultural frontiers there is no better time to divert workers from the fields, as their marginal productivity in agriculture is approaching zero. If relatively few workers were needed to construct the monuments and other public works of Maya polities <xref rid="pone.0087541-Webster1" ref-type="bibr">[32]</xref>, there would remain a pool available for diversion into warfare without much impact on agricultural production. Redirecting farmers to become soldiers who suffer elevated (density independent) mortality might actually <italic>stabilize</italic> subsistence production, as the non-warrior, producing population would suffer less hunger and enjoy a longer expected lifespan. Empirically, Culbert <xref rid="pone.0087541-Culbert1" ref-type="bibr">[31]</xref> might be shown to be correct, but counter-intuitive elements of resource-limited population ecology reduce the <italic>a priori</italic> plausibility of the proposed hypothesis.</p><p>Some analysts of Maya prehistory <xref rid="pone.0087541-Hamblin1" ref-type="bibr">[33]</xref> propose that the agrarian workforce caused the collapse of Maya polities by revolting against intolerable levels of exploitation. Population ecology insights can elevate or diminish the likelihood of such a hypothesis being correct. For instance, its plausibility is low if producers are in the copial phase of population growth, but greater if they already are suffering from the effects of a Malthusian equilibrium. High density independent background mortality would lower the plausibility of this hypothesis by elevating levels of food availability at equilibrium; low background mortality would make it more likely.</p><p>A central question in the collapse literature is whether we need look exclusively to exogenous factors for the causes of collapse. Societies have been challenged from the outside: Huns attack, droughts and pestilence ruin crops, foreigners introduce new diseases. However, endogenous practices also can contribute to sudden instability (<xref ref-type="fig" rid="pone-0087541-g003">Figure 3</xref>). Social practices affecting work, fertility and consumption, can result in severe dislocations, depending on where in its growth trajectory a population happens to be. Innocuous practices such as high fertility established early in a growth cycle, in times of plenty, may become maladaptive later in that cycle.</p><p>Turchin <xref rid="pone.0087541-Turchin2" ref-type="bibr">[34]</xref> finds evidence that agrarian societies exhibit low amplitude population cycles on a period of 200–300 years, a wave-length roughly consistent with our observation of growth to equilibrium over 6–12 generations. Turchin proposes a “demographic-structural” theory that envisions rapid population growth provoking the breakdown of central state authority, which in turn generates conditions (e.g., lawlessness, conflict, degradation of public infrastructure protecting health and food production) that eventually cause a population downturn. Lowered population levels are then followed by restoration of public order and investment, and a renewed cycle of growth. Temporal lags between the linked cycles of population and public order drive the out-of-phase historical observations of oscillations. Our emphasis on the punctuated character of Malthusian constraints is consistent with the demographic-structural approach, but would place greater weight on the structural impact of environmental and demographic variables. Even if central authority remains intact and the more violent expressions of public breakdown muted, it may prove difficult for a society to avoid a prolonged demographic downturn when faced with abrupt imposition of food-limited, Malthusian constraints.</p><p>Our reading of the state collapse literature suggests that there is an overwhelming preference for hypotheses based on external, often environmental causation <xref rid="pone.0087541-Middleton1" ref-type="bibr">[29]</xref>,<xref rid="pone.0087541-Hunt1" ref-type="bibr">[35]</xref>. Extreme drought <xref rid="pone.0087541-Drysdale1" ref-type="bibr">[36]</xref>,<xref rid="pone.0087541-Gill1" ref-type="bibr">[37]</xref>,<xref rid="pone.0087541-Peterson1" ref-type="bibr">[38]</xref>,<xref rid="pone.0087541-Shaowu1" ref-type="bibr">[39]</xref>,<xref rid="pone.0087541-Cullen1" ref-type="bibr">[40]</xref>,<xref rid="pone.0087541-Staubwasser1" ref-type="bibr">[41]</xref> is a favorite candidate. Population, if considered, is in the background, and is seen as a gradual process raising vulnerability to a sudden but exogenous insult. We see in this view the lingering influence of Lotka-Verhulst and the association of density dependent population growth with gradualism. But population ecology processes that are gradual in some respects (e.g., declining surplus, increasing density), can entail abrupt and dramatic changes in others: average life span, fertility and survival among them (<xref ref-type="fig" rid="pone-0087541-g003">Figure 3</xref>). As quality of life indicators these demographic variables measure potential for social disruption. An unexpected plunge over two generations from abundance to the grim reality of Malthusian limits acts with nearly the same suddenness and degree of impact as a short sequence of recurring droughts. A demographic crisis may be just as unexpected, impose levels of hunger (<italic>E</italic>) and hardship just as severe. <italic>MTI</italic>, the opportunity to forestall by adjustment, will be just as short.</p><p>Fortunately, we can predict how particular environmental or demographic circumstances will affect the timing and magnitude of punctuated demographic change. A large habitat area extends the period of copial abundance experienced by a founding population, but it does not prolong the descent to equilibrium once it starts. High yield likewise delays the eventual Malthusian reckoning, but it then shortens <italic>MTI</italic>. A low-yield environment is the only circumstance among the scenarios that we examine in which demographic variables signal an impending equilibrium gradually over the course of habitat in-fill (<xref ref-type="fig" rid="pone-0087541-g006">Figure 6</xref>). Demographically hospitable environments, measured by low background mortality and high fertility, shorten both the copial and transitional phases. As measured by hunger (<italic>E</italic><1), they impose especially severe food shortages once the Malthusian phase occurs. Across all environments, high population density is achieved at a high cost to quality of life and, across most all environments, low yield again being the exception, <italic>MTI</italic> is quite short relative to the length of the food-abundant phase. Each of these predictions tells us something about the likelihood that rapid demographic change lies behind population upheaval and socio-political collapse in particular cases.</p><p>An agrarian producing population may endure at a high density set by food-limited constraints even as the rapid decline in human welfare associated with the Malthusian transition phase causes, or abets the causes of, political collapse. This is consistent with current thinking in the Maya case, which entails the interaction of high population density, increasing warfare, and a gradual and long-standing environmental downturn <xref rid="pone.0087541-Kennett1" ref-type="bibr">[42]</xref>. Malthusian disruption to the social fabric and infrastructure of everyday life might well be enough to drive a highly interconnected civilization to fail, whether or not suffering other stresses.</p></sec><sec id="s4d"><title>Bioarchaeology and Population History</title><p>Anthropologists commonly use bioarchaeological assessments of human health and nutritional status from mortuary populations to assess whether or not prehistoric societies are suffering from food shortages. This methodology and interpretations resulting from it are then used to assess the likelihood that over-population or habitat degradation explain socio-political developments. Absent evidence of osteological insult during the lead-up to a phenomenon like social collapse, population-based explanations are considered to be refuted. For example, in their survey of bioarchaeological evidence from Maya prehistory, Wright and White <xref rid="pone.0087541-Wright2" ref-type="bibr">[43]</xref>:147 state that the ecological model for lowland collapse and abandonment “requires … consistent change in diet from Preclassic and Early Classic Periods to the Terminal Classic.” This imposes a requirement of gradual increase in measures of nutritional deficit over 600–800 years of expanding population density. Wright and White find instead that the isotopic and bioarchaeological evidence shows local diversity and no trends through this lengthy time span. They interpret this result as inconsistent with ecological or population imbalance models.</p><p>However, underlying this interpretation is the assumption that the stress of over-population must build gradually over a significant period of time prior to a dramatic event like socio-political collapse. This assumption would be consistent with a Lotka-Verhulst approach; it is not consistent with our model results showing a long period of relative abundance that ends abruptly with a rapid transition to Malthusian straits. Given the standard error ranges of archaeological dating, and the aggregate nature of mortuary populations, in all except the very low-yield cases evidence of nutritional shortfalls will appear in the archaeological record more-or-less coincident with their social consequences (e.g., collapse). They may be dismissed as an effect rather than a cause. Surplus will decline slowly in the long-term lead-up to disruptive demographic stress (<xref ref-type="fig" rid="pone-0087541-g003">Figure 3</xref>) – the granaries slowly will go empty and the rats become skinny and sparse – but in a growth trajectory we should not expect to find bioarchaeological indicators of stress in humans until the Malthusian crisis of density dependence is only a few generations from its full and crushing impact.</p></sec></sec><sec sec-type="conclusions" id="s5"><title>Conclusion</title><p>Societies are complex systems linking natural and social features <xref rid="pone.0087541-Liu1" ref-type="bibr">[44]</xref>; their study requires models that do the same, while drawing on comparative, historical and prehistorical records <xref rid="pone.0087541-OSullivan1" ref-type="bibr">[45]</xref>. The modeling effort we describe instantiates these goals in a population ecology approach that embodies the fundamental principles discussed by Turchin <xref rid="pone.0087541-Turchin3" ref-type="bibr">[46]</xref>: exponential growth, self-limitation and consumer-resource oscillation. With Turchin <xref rid="pone.0087541-Turchin4" ref-type="bibr">[47]</xref>, we take population ecology models to be tools for seeking understanding of regularities in history.</p><p>Colonizing populations or those recovering from an epidemic or other catastrophe so that they are faced with comeback in a relatively empty habitat will enjoy a lengthy period of comfortable growth before they hit a quick transition to deprivation of indefinite length. Except for a very brief transition, populations in which fertility and mortality are density dependent will exist in one of two states: comfortable or deprived. This runs counter to intuitions formed by familiarity with variants of the logistic equation.</p><p>We propose that growing populations in the copial phase have few means of detecting their fate more than several generations in advance of it being fully upon them. It lurks like an invisible cliff. Their <italic>Malthusian transition interval (MTI)</italic> is short. From initial settlement, the growth rate of colonizers remains constantly high and population size increases at an accelerating rate. Justifications for complacency abound. Surplus food availability declines but until food availability drops below the critical value of <italic>E</italic> = 1, quality of life variables like life span, infant survival and fertility are unaffected by habitat in-fill. Absent nutritional or demographic strain felt at the household level, swelling aggregate population size is unlikely to be perceived or monitored as a signal of pending trouble. Similarly, we do not expect individuals to anticipate crisis from a declining surplus; food availability in excess of need weakens the signal. Growth rate (<italic>r</italic>) by contrast reflects palpable and personal linkages among mortality, fertility and hunger, and it or its constituent vital rates should constitute a strong household signal of environmental difficulty. However, these are constant until a few generations before the Malthusian equilibrium. The most informative and directly reliable signals, those entailing hunger and acute changes in mortality and fertility, offer virtually no lead-time.</p><p>The brevity of the <italic>MTI</italic> limits options for response if adaptations require a long lead time or heavy investments of planning and implementation. It is difficult to envision that a society could successfully undertake the introduction of facilities like extensive irrigation canals, terracing or raised fields on the timespan of a couple of generations after perceiving their need, during a crisis. In a sanguine conclusion to his review of the “Archaeology of Overshoot and Collapse,” Tainter <xref rid="pone.0087541-Tainter1" ref-type="bibr">[28]</xref>:72 says:</p><disp-quote><p>“overshoot is teleological, as if humans could set a target for population or consumption…[it] denies the human capacity for flexible adjustments…It usually is possible to coax more resource production by applying capital and technology, increasing labor, applying energy subsidies, and making production more knowledge-intensive.”</p></disp-quote><p>Given realistic appraisal of how unexpectedly and quickly they may impose themselves, population crises may offer little opportunity for the kind of adjustments Tainter envisions.</p><p>These observations also suggest, however, that we should be cautious in extrapolating from the fate of places like Easter Island to our own situation <xref rid="pone.0087541-Diamond2" ref-type="bibr">[48]</xref>. The collapse of Pacific Island societies sometimes is portrayed as a failure to recognize the need for adaptive response. Plentiful signs were not read or were ignored, reasonable adjustments were not made, and social upheaval ensued when population finally exceeded resource capacities (cf. <xref rid="pone.0087541-Kirch1" ref-type="bibr">[49]</xref>). But this reading of the evidence entails a Lotka-Verhulst version of the world in which human failures are measured by long-standing neglect of distress signals. Our results suggest that the very salubrious Pacific Island setting is precisely the one in which the signal would anticipate problems by the slimmest of margins, perhaps coming too late for an effective response no matter how conscientious the population.</p><p>Many models that place population at the heart of social upheaval envision linked resource degradation as the ultimate cause. Overexploitation provokes deteriorating yields or, in the Maya case “a population collapse precipitated by agricultural failure” <xref rid="pone.0087541-Johnston1" ref-type="bibr">[50]</xref>:141. This combination of causes may well have been the case historically. But it is important to acknowledge that population can act abruptly and with severity even in the absence of an environmental crisis. In our analyses the environmental capacity to yield resources is not declining; our plunge to Malthusian constraints does not invoke factors like reduced yield, forest clearance or soil erosion. Rather, the fate of the population resides primarily in the responses of fertility and mortality in the face of calorie shortages. Population stress and collapse might arise from punctuated demographic failure quite apart from environmental manifestation of resource overexploitation. Abrupt imposition of Malthusian constraints and insufficient <italic>MTI</italic> may well be a stand-alone culprit in social evolution.</p></sec><sec sec-type="supplementary-material" id="s6"><title>Supporting Information</title><supplementary-material content-type="local-data" id="pone.0087541.s001"><label>File S1</label><caption><p><bold>Three sections of supporting material, including Table S1 and two sections of text.</bold></p><p>(DOC)</p></caption><media xlink:href="pone.0087541.s001.doc"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><ack><p>The authors would like to thank Sheryl Gerety, Michael Graves, Patrick Kirch and Thegn Ladefoged for helpful comments. Members of the UC Davis Human Behavioral Ecology (HBE) lab group have been generous in their insights and criticisms. We also thank Michael Gurven and an anonymous reviewer for comments and questions that led to substantial improvements in this manuscript.</p></ack><ref-list><title>References</title><ref id="pone.0087541-Malthus1"><label>1</label><mixed-citation publication-type="book">Malthus TR (2003 [1798]) An Essay on the Principle of Population: Text, Sources and Background, Criticism; Appleman P, editor. New York: W. W. Norton.</mixed-citation></ref><ref id="pone.0087541-Boserup1"><label>2</label><mixed-citation publication-type="book">Boserup E (1965) The Conditions of Agricultural Growth: The Economics of Agrarian Change under Population Pressure. New York: Aldine.</mixed-citation></ref><ref id="pone.0087541-Wood1"><label>3</label><mixed-citation publication-type="journal"><name><surname>Wood</surname><given-names>JW</given-names></name> (<year>1998</year>) <article-title>A theory of preindustrial population dynamics: Demography, economy, and well-being in Malthusian systems</article-title>. <source>Current Anthropology</source><volume>39</volume>: <fpage>99</fpage>–<lpage>135</lpage></mixed-citation></ref><ref id="pone.0087541-Lee1"><label>4</label><mixed-citation publication-type="book">Lee RD (1986) Malthus and Boserup: A dynamic synthesis. In: Coleman D, Schofield R, editors. The State of Population Theory. New York: Basil Blackwell. pp. 96–130.</mixed-citation></ref><ref id="pone.0087541-Lee2"><label>5</label><mixed-citation publication-type="journal"><name><surname>Lee</surname><given-names>RD</given-names></name> (<year>2003</year>) <article-title>Rethinking the evolutionary theory of aging: Transfers, not births, shape senescence in social species</article-title>. <source>Proceedings of the National Academy of Sciences USA</source><volume>100</volume>: <fpage>9637</fpage>–<lpage>9642</lpage></mixed-citation></ref><ref id="pone.0087541-Hutchinson1"><label>6</label><mixed-citation publication-type="book">Hutchinson GE (1978) An Introduction to Population Ecology. New Haven, CN: Yale University Press.</mixed-citation></ref><ref id="pone.0087541-Kohler1"><label>7</label><mixed-citation publication-type="book">Kohler TA, van der Leeuw SE (2007) Introduction: Historical socionatural systems and models. In: Kohler TA, van der Leeuw SE, editors. The Model-Based Archaeology of Socionatural Systems. Santa Fe: School for Advanced Research Press. pp. 1–12.</mixed-citation></ref><ref id="pone.0087541-Turchin1"><label>8</label><mixed-citation publication-type="book">Turchin P (2003) Historical Dynamics: Why States Rise and Fall. Princeton, NJ: Princeton University Press.</mixed-citation></ref><ref id="pone.0087541-Wright1"><label>9</label><mixed-citation publication-type="book">Wright HT (2007) Progress in cultural modeling. In: Kohler TA, van der Leeuw SE, editors. The Model-Based Archaeology of Socionatural Systems. Santa Fe: School for Advanced Research Press. pp. 229–232.</mixed-citation></ref><ref id="pone.0087541-Holling1"><label>10</label><mixed-citation publication-type="journal"><name><surname>Holling</surname><given-names>CS</given-names></name> (<year>1959</year>) <article-title>Some characteristics of simple types of predation and parasitism</article-title>. <source>Canadian Entomologist</source><volume>XCI</volume>: <fpage>385</fpage>–<lpage>398</lpage></mixed-citation></ref><ref id="pone.0087541-Lee3"><label>11</label><mixed-citation publication-type="journal"><name><surname>Lee</surname><given-names>CT</given-names></name>, <name><surname>Tuljapurkar</surname><given-names>S</given-names></name> (<year>2008</year>) <article-title>Population and prehistory I: Food-dependent population growth in constant environments</article-title>. <source>Theoretical Population Biology</source><volume>73</volume>: <fpage>473</fpage>–<lpage>482</lpage><pub-id pub-id-type="pmid">18439637</pub-id></mixed-citation></ref><ref id="pone.0087541-Puleston1"><label>12</label><mixed-citation publication-type="journal"><name><surname>Puleston</surname><given-names>CO</given-names></name>, <name><surname>Tuljapurkar</surname><given-names>S</given-names></name> (<year>2008</year>) <article-title>Population and prehistory II: Space-limited human populations in constant environments</article-title>. <source>Theoretical Population Biology</source><volume>74</volume>: <fpage>147</fpage>–<lpage>160</lpage><pub-id pub-id-type="pmid">18598711</pub-id></mixed-citation></ref><ref id="pone.0087541-Lee4"><label>13</label><mixed-citation publication-type="journal"><name><surname>Lee</surname><given-names>CT</given-names></name>, <name><surname>Puleston</surname><given-names>CO</given-names></name>, <name><surname>Tuljapurkar</surname><given-names>S</given-names></name> (<year>2009</year>) <article-title>Population and prehistory III: Food-dependent demography in variable environments</article-title>. <source>Theoretical Population Biology</source><volume>76</volume>: <fpage>179</fpage>–<lpage>188</lpage><pub-id pub-id-type="pmid">19540865</pub-id></mixed-citation></ref><ref id="pone.0087541-Winterhalder1"><label>14</label><mixed-citation publication-type="book">Winterhalder B, Puleston CO (n.d.) The exchequer's guide to revenue in the agrarian state (ms.). Davis, CA.</mixed-citation></ref><ref id="pone.0087541-Caswell1"><label>15</label><mixed-citation publication-type="book">Caswell H (2004) Matrix Population Models, 2nd edition. Sunderland, MA Sinauer.</mixed-citation></ref><ref id="pone.0087541-Winterhalder2"><label>16</label><mixed-citation publication-type="book">Winterhalder B (2002) Models. In: Hart JP, Terrell JE, editors. Darwin and Archaeology: A Handbook of Key Concepts. Westport, Connecticut: Bergin & Garvey. pp. 201–223.</mixed-citation></ref><ref id="pone.0087541-Bengtsson1"><label>17</label><mixed-citation publication-type="book">Bengtsson T, Campbell C, Lee JZ (2004) Life under Pressure: Mortality and Living Standards in Europe and Asia, 1700–1900. Cambridge, MA: MIT Press.</mixed-citation></ref><ref id="pone.0087541-Winterhalder3"><label>18</label><mixed-citation publication-type="book">Winterhalder B, Puleston C, Ross C (n.d.) Yield fluctuations, household-level food storage and population-level consequences in prehistoric, seasonal agrarian societies (ms.). Davis, CA.</mixed-citation></ref><ref id="pone.0087541-Caughley1"><label>19</label><mixed-citation publication-type="book">Caughley G (1976) Wildlife management and the dynamics of ungulate populations. In: Coaker TH, editor. Applied Ecology, Volume 1. London: Academic Press. pp. 183–246.</mixed-citation></ref><ref id="pone.0087541-Read1"><label>20</label><mixed-citation publication-type="journal"><name><surname>Read</surname><given-names>DW</given-names></name>, <name><surname>LeBlanc</surname><given-names>SA</given-names></name> (<year>2003</year>) <article-title>Population growth, carrying capacity, and conflict</article-title>. <source>Current Anthropology</source><volume>44</volume>: <fpage>59</fpage>–<lpage>85</lpage></mixed-citation></ref><ref id="pone.0087541-Scott1"><label>21</label><mixed-citation publication-type="book">Scott S, Duncan CJ (2003) Demography and Nutrition: Evidence from Historical and Contemporary Populations. New York: John Wiley and Sons.</mixed-citation></ref><ref id="pone.0087541-Winterhalder4"><label>22</label><mixed-citation publication-type="journal"><name><surname>Winterhalder</surname><given-names>B</given-names></name>, <name><surname>Baillargeon</surname><given-names>W</given-names></name>, <name><surname>Cappelletto</surname><given-names>F</given-names></name>, <name><surname>Daniel Jr</surname><given-names>IR</given-names></name>, <name><surname>Prescott</surname><given-names>C</given-names></name> (<year>1988</year>) <article-title>The population ecology of hunter-gatherers and their prey</article-title>. <source>Journal of Anthropological Archaeology</source><volume>7</volume>: <fpage>289</fpage>–<lpage>328</lpage></mixed-citation></ref><ref id="pone.0087541-Neubert1"><label>23</label><mixed-citation publication-type="journal"><name><surname>Neubert</surname><given-names>MG</given-names></name>, <name><surname>Caswell</surname><given-names>H</given-names></name> (<year>2000</year>) <article-title>Density-dependent vital rates and their population dynamic consequences</article-title>. <source>Journal of Mathematical Biology</source><volume>41</volume>: <fpage>103</fpage>–<lpage>121</lpage><pub-id pub-id-type="pmid">11039693</pub-id></mixed-citation></ref><ref id="pone.0087541-Gaillard1"><label>24</label><mixed-citation publication-type="journal"><name><surname>Gaillard</surname><given-names>J-M</given-names></name>, <name><surname>Festa-Bianchet</surname><given-names>M</given-names></name>, <name><surname>Yoccoz</surname><given-names>NG</given-names></name> (<year>1988</year>) <article-title>Population dynamics of large herbivores: Variable recruitment with constant adult survival</article-title>. <source>TREE</source><volume>13</volume>: <fpage>58</fpage>–<lpage>63</lpage><pub-id pub-id-type="pmid">21238201</pub-id></mixed-citation></ref><ref id="pone.0087541-Tsoularis1"><label>25</label><mixed-citation publication-type="journal"><name><surname>Tsoularis</surname><given-names>A</given-names></name>, <name><surname>Wallace</surname><given-names>J</given-names></name> (<year>2002</year>) <article-title>Analysis of logistic growth models</article-title>. <source>Mathematical Biosciences</source><volume>179</volume>: <fpage>21</fpage>–<lpage>55</lpage><pub-id pub-id-type="pmid">12047920</pub-id></mixed-citation></ref><ref id="pone.0087541-Sibly1"><label>26</label><mixed-citation publication-type="journal"><name><surname>Sibly</surname><given-names>RM</given-names></name>, <name><surname>Barker</surname><given-names>D</given-names></name>, <name><surname>Denham</surname><given-names>MC</given-names></name>, <name><surname>Hone</surname><given-names>J</given-names></name>, <name><surname>Pagel</surname><given-names>M</given-names></name> (<year>2005</year>) <article-title>On the regulation of populations of mammals, birds, fish, and insects</article-title>. <source>Science</source><volume>309</volume>: <fpage>607</fpage>–<lpage>610</lpage><pub-id pub-id-type="pmid">16040705</pub-id></mixed-citation></ref><ref id="pone.0087541-Diamond1"><label>27</label><mixed-citation publication-type="book">Diamond JM (2004) Collapse: How Societies Choose to Fail or Succeed. New York: Penguin.</mixed-citation></ref><ref id="pone.0087541-Tainter1"><label>28</label><mixed-citation publication-type="journal"><name><surname>Tainter</surname><given-names>JA</given-names></name> (<year>2006</year>) <article-title>Archaeology of overshoot and collapse</article-title>. <source>Annual Review of Anthropology</source><volume>35</volume>: <fpage>59</fpage>–<lpage>74</lpage></mixed-citation></ref><ref id="pone.0087541-Middleton1"><label>29</label><mixed-citation publication-type="journal"><name><surname>Middleton</surname><given-names>GD</given-names></name> (<year>2012</year>) <article-title>Nothing lasts forever: Environmental discourses on the collapse of past societies</article-title>. <source>Journal of Archaeological Research</source><volume>20</volume>: <fpage>257</fpage>–<lpage>307</lpage></mixed-citation></ref><ref id="pone.0087541-TurnerII1"><label>30</label><mixed-citation publication-type="journal"><name><surname>Turner II</surname><given-names>BL</given-names></name>, <name><surname>Sabloff</surname><given-names>JA</given-names></name> (<year>2012</year>) <article-title>Classic Period collapse of the Central Maya lowlands: Insights about human-environment relationships for sustainability</article-title>. <source>Proceedings of the National Academy of Sciences USA</source><volume>109</volume>: <fpage>13908</fpage>–<lpage>13914</lpage></mixed-citation></ref><ref id="pone.0087541-Culbert1"><label>31</label><mixed-citation publication-type="book">Culbert TP (1988) Collapse of Classic Maya civilization. In: Yoffee N, Cowgill GL, editors. The Collapse of Ancient States and Civilizations. Tucson, AZ: University of Arizona Press. pp. 69–101.</mixed-citation></ref><ref id="pone.0087541-Webster1"><label>32</label><mixed-citation publication-type="journal"><name><surname>Webster</surname><given-names>D</given-names></name> (<year>1985</year>) <article-title>Surplus, labor, and stress in Late Classic Maya society</article-title>. <source>Journal of Anthropological Research</source><volume>41</volume>: <fpage>375</fpage>–<lpage>399</lpage></mixed-citation></ref><ref id="pone.0087541-Hamblin1"><label>33</label><mixed-citation publication-type="journal"><name><surname>Hamblin</surname><given-names>RL</given-names></name>, <name><surname>Pitcher</surname><given-names>BL</given-names></name> (<year>1980</year>) <article-title>The Classic Maya collapse: Testing class conflict hypotheses</article-title>. <source>American Antiquity</source><volume>45</volume>: <fpage>246</fpage>–<lpage>267</lpage></mixed-citation></ref><ref id="pone.0087541-Turchin2"><label>34</label><mixed-citation publication-type="journal"><name><surname>Turchin</surname><given-names>P</given-names></name> (<year>2009</year>) <article-title>Long-term population cycles in human societies</article-title>. <source>The Year in Ecology and Conservation Biology, 2009: Annals of the New York Academy of Science</source><volume>1162</volume>: <fpage>1</fpage>–<lpage>17</lpage></mixed-citation></ref><ref id="pone.0087541-Hunt1"><label>35</label><mixed-citation publication-type="journal"><name><surname>Hunt</surname><given-names>BG</given-names></name>, <name><surname>Elliott</surname><given-names>TI</given-names></name> (<year>2005</year>) <article-title>A simulation of the climatic conditions associated with the collapse of the Maya civilization</article-title>. <source>Climate Change</source><volume>69</volume>: <fpage>393</fpage>–<lpage>407</lpage></mixed-citation></ref><ref id="pone.0087541-Drysdale1"><label>36</label><mixed-citation publication-type="journal"><name><surname>Drysdale</surname><given-names>R</given-names></name>, <name><surname>Zanchetta</surname><given-names>G</given-names></name>, <name><surname>Hellstrom</surname><given-names>J</given-names></name>, <name><surname>Maas</surname><given-names>R</given-names></name>, <name><surname>Fallick</surname><given-names>A</given-names></name>, <etal>et al</etal> (<year>2006</year>) <article-title>Late Holocene drought responsible for the collapse of Old World civilizations is recorded in an Italian cave flowstone</article-title>. <source>Geology</source><volume>34</volume>: <fpage>101</fpage>–<lpage>104</lpage></mixed-citation></ref><ref id="pone.0087541-Gill1"><label>37</label><mixed-citation publication-type="journal"><name><surname>Gill</surname><given-names>RB</given-names></name>, <name><surname>Mayewski</surname><given-names>PA</given-names></name>, <name><surname>Nyberg</surname><given-names>J</given-names></name>, <name><surname>Haug</surname><given-names>GH</given-names></name>, <name><surname>Peterson</surname><given-names>LC</given-names></name> (<year>2007</year>) <article-title>Drought and the Maya collapse</article-title>. <source>Ancient Mesoameria</source><volume>18</volume>: <fpage>283</fpage>–<lpage>302</lpage></mixed-citation></ref><ref id="pone.0087541-Peterson1"><label>38</label><mixed-citation publication-type="journal"><name><surname>Peterson</surname><given-names>LC</given-names></name>, <name><surname>Haug</surname><given-names>GH</given-names></name> (<year>2005</year>) <article-title>Climate and the collapse of Maya civilization</article-title>. <source>American Scientist</source><volume>93</volume>: <fpage>322</fpage>–<lpage>329</lpage></mixed-citation></ref><ref id="pone.0087541-Shaowu1"><label>39</label><mixed-citation publication-type="journal"><name><surname>Shaowu</surname><given-names>W</given-names></name> (<year>2005</year>) <article-title>Abrupt climate change and collapse of ancient civilizations at 2200BC - 2000BC</article-title>. <source>Progress in Natural Science</source><volume>15</volume>: <fpage>908</fpage>–<lpage>914</lpage></mixed-citation></ref><ref id="pone.0087541-Cullen1"><label>40</label><mixed-citation publication-type="journal"><name><surname>Cullen</surname><given-names>HM</given-names></name>, <name><surname>deMenocal</surname><given-names>PB</given-names></name>, <name><surname>Hemming</surname><given-names>S</given-names></name>, <name><surname>Hemming</surname><given-names>G</given-names></name>, <name><surname>Brown</surname><given-names>FH</given-names></name>, <etal>et al</etal> (<year>2010</year>) <article-title>Climate change and the collapse of the Akkadian empire: Evidence from the deep sea</article-title>. <source>Geology</source><volume>28</volume>: <fpage>379</fpage>–<lpage>382</lpage></mixed-citation></ref><ref id="pone.0087541-Staubwasser1"><label>41</label><mixed-citation publication-type="journal"><name><surname>Staubwasser</surname><given-names>M</given-names></name>, <name><surname>Sirocko</surname><given-names>F</given-names></name>, <name><surname>Grootes</surname><given-names>PM</given-names></name>, <name><surname>Segl</surname><given-names>M</given-names></name> (<year>2003</year>) <article-title>Climate change at the 4.2 ka BP termination of the Indus valley civilization and Holocene south Asian monsoon variability</article-title>. <source>Geophysical Research Letters</source><volume>30</volume>: <fpage>1425</fpage>–<lpage>1428</lpage></mixed-citation></ref><ref id="pone.0087541-Kennett1"><label>42</label><mixed-citation publication-type="journal"><name><surname>Kennett</surname><given-names>DJ</given-names></name>, <name><surname>Breitenbach</surname><given-names>SFM</given-names></name>, <name><surname>Aquino</surname><given-names>VV</given-names></name>, <name><surname>Asmerom</surname><given-names>Y</given-names></name>, <name><surname>Awe</surname><given-names>J</given-names></name>, <etal>et al</etal> (<year>2012</year>) <article-title>Development and disintegration of Maya political systems in response to climate change</article-title>. <source>Science</source><volume>388</volume>: <fpage>788</fpage>–<lpage>791</lpage><pub-id pub-id-type="pmid">23139330</pub-id></mixed-citation></ref><ref id="pone.0087541-Wright2"><label>43</label><mixed-citation publication-type="journal"><name><surname>Wright</surname><given-names>LE</given-names></name>, <name><surname>White</surname><given-names>CD</given-names></name> (<year>1996</year>) <article-title>Human biology in the Classic Maya collapse: Evidence from paleopathology and paleodiet</article-title>. <source>Journal of World Prehistory</source><volume>10</volume>: <fpage>147</fpage>–<lpage>198</lpage></mixed-citation></ref><ref id="pone.0087541-Liu1"><label>44</label><mixed-citation publication-type="journal"><name><surname>Liu</surname><given-names>J</given-names></name>, <name><surname>Dietz</surname><given-names>T</given-names></name>, <name><surname>Carpenter</surname><given-names>SR</given-names></name>, <name><surname>Alberti</surname><given-names>M</given-names></name>, <name><surname>Folke</surname><given-names>C</given-names></name>, <etal>et al</etal> (<year>2007</year>) <article-title>Complexity of coupled human and natural systems</article-title>. <source>Science</source><volume>317</volume>: <fpage>1513</fpage>–<lpage>1516</lpage><pub-id pub-id-type="pmid">17872436</pub-id></mixed-citation></ref><ref id="pone.0087541-OSullivan1"><label>45</label><mixed-citation publication-type="journal"><name><surname>O'Sullivan</surname><given-names>P</given-names></name> (<year>2008</year>) <article-title>The ‘collapse’ of civilizations: What palaeoenvironmental reconstruction cannot tell us, but anthropology can</article-title>. <source>The Holocene</source><volume>18</volume>: <fpage>45</fpage>–<lpage>55</lpage></mixed-citation></ref><ref id="pone.0087541-Turchin3"><label>46</label><mixed-citation publication-type="journal"><name><surname>Turchin</surname><given-names>P</given-names></name> (<year>2001</year>) <article-title>Does population ecology have general laws?</article-title><source>Oikos</source><volume>94</volume>: <fpage>17</fpage>–<lpage>26</lpage></mixed-citation></ref><ref id="pone.0087541-Turchin4"><label>47</label><mixed-citation publication-type="journal"><name><surname>Turchin</surname><given-names>P</given-names></name> (<year>2008</year>) <article-title>Arise ‘Cliodynamics’</article-title>. <source>Nature</source><volume>454</volume>: <fpage>34</fpage>–<lpage>35</lpage><pub-id pub-id-type="pmid">18596791</pub-id></mixed-citation></ref><ref id="pone.0087541-Diamond2"><label>48</label><mixed-citation publication-type="journal"><name><surname>Diamond</surname><given-names>JM</given-names></name> (<year>2004</year>) <article-title>Twilight at Easter</article-title>. <source>The New York Review of Books</source><volume>LI</volume>: <fpage>6</fpage>–<lpage>10</lpage></mixed-citation></ref><ref id="pone.0087541-Kirch1"><label>49</label><mixed-citation publication-type="journal"><name><surname>Kirch</surname><given-names>PV</given-names></name> (<year>1997</year>) <article-title>Microcosmic histories: Island perspectives on “Global” change</article-title>. <source>American Anthropologist</source><volume>99</volume>: <fpage>30</fpage>–<lpage>42</lpage></mixed-citation></ref><ref id="pone.0087541-Johnston1"><label>50</label><mixed-citation publication-type="journal"><name><surname>Johnston</surname><given-names>KJ</given-names></name> (<year>2003</year>) <article-title>The intensification of pre-industrial cereal agriculture in the tropics: Boserup, cultivation lengthening, and the Classic Maya</article-title>. <source>Journal of Anthropological Archaeology</source><volume>22</volume>: <fpage>126</fpage>–<lpage>161</lpage></mixed-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Archeologie/explor/PaleopathV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 0003760 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 0003760 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Archeologie
|area= PaleopathV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé=
|texte=
}}
| This area was generated with Dilib version V0.6.27. |  |