
==== Front
Nature
Nature
Nature
0028-0836
1476-4687
Nature Publishing Group UK London

38693262
7354
10.1038/s41586-024-07354-8
Article
Frequent disturbances enhanced the resilience of past human populations
http://orcid.org/0000-0003-4244-7495
Riris Philip priris@bournemouth.ac.uk

1
http://orcid.org/0000-0003-1368-6331
Silva Fabio 1
http://orcid.org/0000-0001-6727-5138
Crema Enrico 2
http://orcid.org/0000-0003-0758-5032
Palmisano Alessio 3
http://orcid.org/0000-0002-0789-3724
Robinson Erick 456
http://orcid.org/0000-0003-1260-7646
Siegel Peter E. 7
French Jennifer C. 8
Jørgensen Erlend Kirkeng 9
http://orcid.org/0000-0002-4333-1972
Maezumi Shira Yoshi 10
http://orcid.org/0000-0001-8293-8147
Solheim Steinar 11
Bates Jennifer 12
Davies Benjamin 13
Oh Yongje 12
http://orcid.org/0000-0002-3247-437X
Ren Xiaolin 14
1 https://ror.org/05wwcw481 grid.17236.31 0000 0001 0728 4630 Department of Archaeology and Anthropology, Bournemouth University, Poole, UK
2 https://ror.org/013meh722 grid.5335.0 0000 0001 2188 5934 Department of Archaeology, University of Cambridge, Cambridge, UK
3 https://ror.org/048tbm396 grid.7605.4 0000 0001 2336 6580 Department of Historical Studies, University of Turin, Torino, Italy
4 Native Environment Solutions, Boise, ID USA
5 https://ror.org/02vg22c33 grid.474431.1 0000 0004 0525 4843 Division of Atmospheric Sciences, Desert Research Institute, Reno, NV USA
6 https://ror.org/03efmqc40 grid.215654.1 0000 0001 2151 2636 School of Human Evolution and Social Change, Arizona State University, Tempe, AZ USA
7 https://ror.org/01nxc2t48 grid.260201.7 0000 0001 0745 9736 Department of Anthropology, Montclair State University, Montclair, NJ USA
8 https://ror.org/04xs57h96 grid.10025.36 0000 0004 1936 8470 Department of Archaeology, Classics, and Egyptology, University of Liverpool, Liverpool, UK
9 https://ror.org/02xhrye98 grid.436614.2 0000 0001 0730 2472 NIKU High North Department, Norwegian Institute for Cultural Heritage Research, Tromsø, Norway
10 https://ror.org/00js75b59 Department of Archaeology, Max Planck Institute of Geoanthropology, Jena, Germany
11 https://ror.org/01xtthb56 grid.5510.1 0000 0004 1936 8921 The Museum of Cultural History, University of Oslo, Oslo, Norway
12 https://ror.org/04h9pn542 grid.31501.36 0000 0004 0470 5905 Department of Archaeology and Art History, Seoul National University, Seoul, Republic of Korea
13 https://ror.org/05wvpxv85 grid.429997.8 0000 0004 1936 7531 Environmental Studies, Tufts University, Boston, MA USA
14 grid.9227.e 0000000119573309 Institute for the History of Natural Sciences, Chinese Academy of Sciences, Beijing, People’s Republic of China
1 5 2024
1 5 2024
2024
629 8013 837842
6 12 2023
26 3 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
The record of past human adaptations provides crucial lessons for guiding responses to crises in the future1–3. To date, there have been no systematic global comparisons of humans’ ability to absorb and recover from disturbances through time4,5. Here we synthesized resilience across a broad sample of prehistoric population time–frequency data, spanning 30,000 years of human history. Cross-sectional and longitudinal analyses of population decline show that frequent disturbances enhance a population’s capacity to resist and recover from later downturns. Land-use patterns are important mediators of the strength of this positive association: farming and herding societies are more vulnerable but also more resilient overall. The results show that important trade-offs exist when adopting new or alternative land-use strategies.

Analysis of population decline shows that frequent disturbances enhance a population’s capacity to resist and recover from downturns and that trade-offs exist when adopting new or alternative land-use strategies.

Subject terms

Environmental social sciences
Climate-change impacts
Social evolution
Climate-change adaptation
Archaeology
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcMain

Understanding the range of past responses of human societies to disturbances is a global priority across the social and natural sciences and will support the development of solutions to future crises1–3. Numerous case studies have addressed past cultural collapse, transformation and persistence, although how to best characterize these processes is a subject for debate6. A major unresolved issue is the lack of comparability between cases of population resilience in the historical sciences4,5. Few studies explicitly model impacts, recovery and adaptation, or formally account for long-term history, which contains important and previously overlooked variation within and between cultural or environmental settings. Furthermore, a tendency to narrowly focus on responses to extreme events in both natural and social systems7,8 may overemphasize local or short-term adaptive success at the expense of understanding large-scale or long-term vulnerabilities6,9. A well-known example is the shift to a narrow marine diet among the Greenland Norse that initially offset the short-term risk of crop failure yet heightened societal vulnerability to longer-term North Atlantic cooling10. Here, we establish a global comparative approach to long-term resilience to identify the factors that structure the response of prehistoric populations to disturbances through time. The approach measures population capacity to withstand changes, as well as the rate of recovery following a disturbance through the common proxy of radiocarbon time–frequency data11,12. Disturbances are the inferred drivers of relative reductions in population or archaeological activity in prehistory, which are described variously as recessions, downturns, busts, negative deviations or similar12–15 and form the focus of this study, using summed probability distributions (SPDs) of calibrated radiocarbon dates. SPDs function as an index of relative levels of human activity, or population change, over time16,17. Population downturns are defined as periods when SPDs are significantly below expected growth trajectories in response to disturbances. Our efforts focus on two key questions: (1) how quickly do past populations recover after downturns; and (2) what factors mediate past resistance and resilience to downturns?.

To quantify patterns in population resistance–resilience, we performed a meta-analysis of 16 published study regions that used archaeological radiocarbon data to reconstruct regional palaeodemographic trends (Supplementary Table 1). Our approach trades specificity for a large-scale, comparative perspective that still accounts for variation between cases to focus on the emergent properties of the statistical analysis. Studies were reviewed based on three criteria: evidence for significant downturns, their scope and the inclusion of radiocarbon datasets. A lack of any single criterion resulted in the exclusion of a study. Cases with no reported downturns were not included, nor were those whose scope was restricted to specific activities within a regional radiocarbon dataset, such as flint mining18. Where published data have been superseded by later compilations, dates were added from the People3k database19, a systematic compilation of cleaned radiocarbon dates, based on the geographical area of the original study. Our global sample of regions ranges from the Arctic to the tropics and spans 30,000 years of history (Fig. 1a). Population downturns were reproduced using our protocol (Methods) and resistance–resilience metrics (Fig. 1b) were collected on 154 periods of population downturn, with a median of 8.5 periods in each region (Fig. 1c and Supplementary Tables 1 and 2). The numerical metrics collectively capture the severity, chronology and frequency of periods of statistically significant population downturn. Disturbances were classified into both general categories and specific drivers according to the original studies and expert interpretation of regional social, cultural and environmental history (Extended Data Table 1). The broad category of land use and evidence for adaptive change during, or in the wake of, downturns were also recorded. The focus of this meta-analysis is to identify the factors governing the relative depth of downturns (resistance) and the rate of recovery at their conclusion (resilience). A suite of hierarchical linear mixed models was developed to test for significant associations between parameters while controlling for regional variability (Methods).Fig. 1 Map and summaries of all world regions and disturbance types included in this study.

a, World map of study regions made with Natural Earth. Time ranges in calendar years before present (cal bp) and 14C dates used in this meta-analysis are indicated for each region. b, Conceptual diagram of measuring resistance–resilience on palaeodemographic downturns against expected growth trajectories. Downturn A is longer and faster, with low resistance and low relative rate of recovery (resilience). Downturn B is shorter and slower, displaying higher resistance and higher resilience. Other combinations of high or low resistance and/or resilience are possible (Extended Data Fig. 4). Parameters b, x and e for the equations in Table 1 are indicated in white dots. c, Summary of disturbance types within three broad categories reported in published palaeodemographic studies. Listed proxies are examples drawn from the surveyed literature. ‘Unclear’ downturns (downturns with a lack of clear evidence for a driver) are not shown (Extended Data Table 1).

This approach to past resistance–resilience provides a glimpse into population-level responses to disturbances throughout human history. Results demonstrate that a single factor—the frequency of downturns—increases both the ability to withstand disturbances and to recover from them. Additionally, the frequency of downturns experienced by a given population is influenced by land use: agricultural and agropastoral populations experience significantly more downturns over time than other land-use patterns recorded during downturns (Extended Data Table 1). These findings collectively suggest that the global shift to food-producing economies during the Holocene (starting 11,700 calendar years before present) may have increased population vulnerability to disturbances, yet in the process enhanced their adaptive capacity through repeated exposure. Parallels to our long-term perspective on human population change in macroecology suggest that comparative approaches in the historical sciences have the potential to generate profound insights into past human–environmental relationships on a broad scale.

Cross-sectionally high resilience

The most common high-level driver of downturns, after those with a lack of evidence for a specific cause (unclear, n = 65), is environmental (n = 48, 31%), followed by mixed (cultural–environmental) (n = 33, 21%). The most common disturbance type is aridity in the environmental category (n = 31), followed by mobility (n = 20) in the cultural category. Only a third of recorded downturns have resistance values that drop more than 50% from pre-downturn activity levels (n = 53, 34.4%; Fig. 2a). Most downturns (n = 133, 86%) end before baseline relative population levels are attained; in other words, observed SPD values are lower at the end of most downturns than at the pre-disturbance reference value. Resilience is relatively high across all cases (median = 0.64, n = 154), with 40% (n = 63) still attaining 90% of pre-downturn conditions by their end (Fig. 2b). Full returns to SPD baseline conditions are frequently interrupted by subsequent downturns. Downturns associated with cultural drivers return the highest median resilience (0.74), whereas the median mixed (cultural–environmental) resistance is highest at 0.79. Resistance (0.65) and resilience are lowest among climate-driven downturns (0.57). However, we do not find support for significant differences in either metric across disturbance categories (analysis of variance, resistance: F = 0.541, P = 0.65, d.f. = 3, n = 154; resilience: F = 0.04, P = 0.98, d.f. = 3, n = 154).Fig. 2 Summaries of downturn numerical metrics showing the relationship between resistance, resilience and the duration and relative pace of downturns.

a,b, Distributions of resistance (a) and resilience (b) across disturbance categories (Extended Data Table 1). The lower and upper hinges correspond to the 25th and 75th percentiles The upper and lower whiskers extend from the hinges to 1.5 × the interquartile range (IQR). The thick black line represents the group median. Data beyond the whiskers are individually plotted outlying points. The dashed line represents the combined data median. c, Distribution of the duration of downturns and the time to SPD minima across all recorded downturns. d, Relative pace (overall duration normalized by time to minimum) approximates a normal distribution and enables comparison of the speed of downturns. Modelled downturn durations skew towards multidecadal and centennial timescales.

Global variation in recovery rate

Initial modelling indicates that the geographical location of downturns does not affect resistance, with the exceptions of significantly higher values in the Caribbean archipelago and the Italian peninsula. Conversely, there is support for significantly lower values for resilience in three regions: the Central China Plains, the Caribbean archipelago and the Korean Peninsula, when compared with all other regions (Extended Data Table 2). Further examination of these cases reveals that a large minority of downturns in these three regions return negative values of resilience (n = 11, 39%), which are produced when the population proxy exceeds the SPD value at the start of the downturn by its end (Extended Data Fig. 4, 12). Although the SPD population proxy remains below modelled expectations in all of these cases and therefore they are, in the strictest sense, downturns relative to the null models, the results nevertheless imply that populations in these regions were, on average, able to recover faster than the norm. Owing to the observed range of variation and its potential impact on the metrics, study region was retained as a random effect variable in the mixed-effect models.

Long-term downturns are the norm

The durations of population downturns tend towards centennial (100–500 years, n = 48, 31%) and decadal periods (less than or equal to 50 years, n = 47, 30.5%), with a median of 98 years across the sample. Downturns lasting longer than 500 years are a minority (n = 29). The time taken to reach SPD minima skews further towards decadal timescales (Fig. 2c). Only a single downturn that commences with the 8.2-thousand-years-ago event in the Near East20 has a time to minimum longer than a millennium (2,070 years). Both variables have a strong positive skew (duration = 2.84, time to minimum = 4.23). To control for the distribution and broad range of variation in these variables, the time to minimum was normalized by downturn duration to produce an index of its relative speed, which we term ‘pace’ (Fig. 2d). This transformation enables comparison of the variation between downturns, with higher numbers reflecting slower downturns (σ = 0.55, s.d. = 0.23) and lacking support for non-normally distributed values as in the time to minimum and duration variables (Shapiro–Wilk W = 0.98396, P = 0.07108). Consequently, relative pace was employed as a fixed-effect candidate in the mixed-effect modelling.

Land use mediates resilience

The frequency of downturns over time by region was estimated by normalizing the cumulative number of downturns in a region by their duration (Table 1). We transformed this to the logarithm of events per millennium to account for its strongly skewed distribution. The variable allows us to compare the rate at which downturns occur. It consistently displays the strongest relationship to both resistance and resilience (P < 0.001 in both cases) (Extended Data Table 3) and is the only fixed variable retained by the information criterion-based selection procedure.Table 1 Metrics and model parameters extracted from global case studies of past population downturns

Name	Description	Notes	
Parameters	
Overall duration	Length of downturn, in calendar years before present		
ndownturn	Cumulative number of downturns in a region		
Time to minimum	Time to SPD minimum value, in calendar years before present		
Tstart	Calendar year before present at start of downturn		
Tend	Calendar year before present at end of downturn		
b	Baseline SPD value at the start of a downturn	Fig. 1 and Extended Data Fig. 4	
x	Minimum SPD value during a downturn	
e	SPD value at the end of a downturn	
Independent	
Resistance	The depth of a downturn relative to baseline conditions. Range: 0–1	1−2×∣b−x∣∣b∣+∣b−x∣	
Resilience	Rate of recovery to baseline conditions, controlling for maximum impact of downturn. Range: −1–1	2×∣b−x∣∣b−x∣+∣b−e∣−1	
Dependent	
Geographical information	Latitude, longitude, world region.		
Relative pace	The time to SPD minimum normalized by downturn duration. Range: 0–1	timetominimumoverallduration	
Frequency of downturns	Cumulative number of downturns per region standardized by duration in calendar years, per millennium.	ndownturnduration×1,000	
Category	Driver of disturbance.	Extended Data Table 3	
Type	Type of disturbance.	
Land use	Dominant land-use pattern.	
Change	Cultural changes over the interval of the downturn.	Boolean variable, with specific responses noted separately	
Detailed descriptions of resistance–resilience metrics are in Methods.

Collectively, these results indicate that populations experience an enhanced ability to withstand disturbances as frequency increases, as well as to recover in the aftermath (Fig. 3a,b). Further examination of frequency of downturns shows that, among the recorded disturbance types, changes to mobility regimes (median frequency of downturns = 2.26, n = 20) and high environmental variability (median frequency of downturns = 2.19, n = 17) occur at the highest rate per millennium, whereas cooling occurs at the lowest rate (median frequency of downturns = 0.761, n = 4) (Fig. 3c). In terms of regional variation, the highest frequency of downturn is recorded in the South African Greater Cape Floristic Region (median frequency of downturns = 2.41, n = 17 over 9,950 years) and the lowest in the Korean Peninsula (median frequency of downturns = 0.58, n = 3 over 4,000 years).Fig. 3 Effect sizes of model terms and dominant land use during downturns over time.

a,b, Resilience (a) and resistance (b) are strongly influenced by the rate of downturn per millennium, with a stronger effect for resistance. Values are extracted from fitted models I and II. c, Standardized regression coefficients for significant (P < 0.05) mixed-model terms for n = 154 independent samples based on a two-sided test. Hunter-gatherer was set as the reference group for land use. d, Proportions of dominant land-use types during downturns, in 1,000-year time slices. Pleistocene downturns (before 11,000 calendar years before present, n = 24) have been combined. The 95% confidence intervals are indicated in a–c by shaded bands and error bars. *P < 0.05, **P < 0.01 and ***P < 0.001.

Treating the frequency of downturns as a response variable (Methods) revealed that agricultural and agropastoralist land-use patterns are associated with significantly higher rates of downturn compared with low-level food production, marine foraging or mixed subsistence. Results from this additional modelling exercise suggest that the frequency of downturns is likely to have an important effect on resistance and recovery, whereas the frequency of downturns itself covaries with the dominant pattern of land use and disturbance type in a given time and place. A larger sample size of downturns would increase the explanatory power of our approach and enhance our characterization of these relationships.

Discussion

This meta-analysis has examined the potential factors influencing resistance and resilience across a broad archaeological sample, and was controlled for regional variation. The frequency of downturns is the main determinant of the observed outcomes among the examined factors. The relationship between resistance and resilience displays variable rates, although these events all tend to unfold at multidecadal to centennial timescales. Well-known historical examples support this finding: the catastrophic depopulation of indigenous groups across the Americas took place over centuries21, and the collapse of imperial power in Western Rome was preceded by a long period of rural depopulation22. Systematic data on independent population downturns in prehistory are less common. However, what data are available23,24 corroborate this result: palaeodemographic downturns resolved in radiocarbon data tend to last at least one human generation but frequently much longer.

The frequency of downturns is associated with both the ability of past populations to withstand downturns and the rate of recovery following them across a broad sample of human populations. The results suggest the existence of a common mechanism among human populations that confers resilience to disturbances. The size of this interaction is greater for resistance (eta-squared (η2) = 0.46, P < 0.001; Fig. 3b) than resilience (η2 = 0.29, P < 0.001; Fig. 3a). In practical terms, this suggests that the ability to withstand downturns is distinct from the ability to recover in their wake. We note that this result does not imply a monocausal or deterministic relationship; there are likely to be several adaptive pathways that increase population resistance and resilience by means of the mechanism of increasing downturn frequency.

Further testing indicates that land-use practices associated with food production, such as farming and herding, are significantly and positively correlated with the frequency of downturns (Fig. 3c and Extended Data Table 3). From the early Holocene, the proportion of land use associated with food production in our sample of downturns also increased, as the aggregate global population became gradually more reliant on domesticated species for meeting subsistence needs (Fig. 3d). Collectively, these trends show that although populations generally increased resistance and resilience over time, the heightened rate of downturn over time is itself likely to be linked to the historical tendency towards food-producing subsistence systems. Current archaeological evidence does not indicate that this process was unidirectional or inevitable; foraging and food production are neither mutually exclusive nor in opposition. Our synthetic findings agree with specific cases showing that the behavioural and social changes that food production entailed had trade-offs in other arenas23.

Traditional agricultural or agropastoral practices include a diverse range of socio-ecological systems that are often highlighted as potential sources of inspiration for solutions to current sustainability, biodiversity and conservation challenges5,10,25,26. The results suggest that population downturns or collapse are an inherent property of these systems and a potential trade-off of promoting their use. Systematic reviews in disturbance ecology indicate that frequent natural disturbances enhance the long-term resilience of key ecosystem services and that localized subsystem declines are an important mechanism through which this occurs27. Our study provides insight into the possible existence of an analogous relationship within our sample of human populations; higher frequencies are strongly correlated with both smaller downturns and closer matches to pre-downturn values in the SPD proxies. We suggest that humanity’s overall constant long-term population growth28 may have been sustained due, in part, to the emergent positive feedback between vulnerability, resistance and recovery documented here.

Population decline has been termed an ‘inevitable’ feature of our species’ demographic dynamics29. In their systematic review of historical collapse and resilience dynamics, Cumming and Peterson1 list depopulation as a key metric or factor in every ancient socio-ecological system studied. We anticipate that this singular status will continue undiminished. Our contribution indicates that downturns play an important role in human population history by enhancing the resilience of survivor populations. We speculate that the creation of biased cultural transmission may be responsible; downturns provide critical opportunities for landscape learning and the strengthening of local-to-regional knowledge networks to propagate through a cultural system30,31. Population downturns have been identified as potential triggers of labour investment in infrastructure, social cohesion and technological advancement15. Together, these mechanisms have the potential to enhance the preferential transmission of knowledge and practices that promote future resistance or resilience10. Raising population thresholds by intensifying land use may also heighten the risk of more serious collapse in return for increasingly marginal benefits1,24,32. Other non-trivial and historically contingent factors that are likely to affect outcomes are the diversity and ecology of domesticated species assemblages, degree and type of political complexity, and settlement patterning in relation to the environment. Indicators such as the cessation of monument construction, loss of literacy or economic turmoil can provide additional insight into the consequences and/or potential drivers of population downturns. These potential causal links must be rigorously tested, however, as they are not easily disentangled. A realistic model for the generative mechanism underlying the resilience of human populations will therefore have to be multiscalar and sensitive to cascading effects to account for how various exogenous impacts unfold and endogenous strategies are developed to solve them. The approach used for our comparative analysis of palaeodemographic resistance–resilience, however, does not distinguish between these elements of the studied populations. Future research may translate between our work and the microscale, the patterns of which are only truly understandable within the kind of generalist, synthetic frame of reference that we provide.

This study finds parallels in macroecology, where analogous resistance–resilience outcomes have been suggested to only fully resolve at centennial timescales or above33,34. Archaeology is uniquely suited to examining past population history, and the dynamics that underlie these trends, at such long-term timescales35–37. Understanding past societies’ responses to crises is often explicitly motivated by the goal of applying learnings from the historical sciences to present-day policy and activism, contributing to the ultimate objective of fostering resilient adaptations for the future6. Most archaeological work on past resilience is historically particularistic4,9 and emphasizes the contingencies, decisions and practices that underwrote successful adaptations in specific times and places38,39. This specificity can be illuminating but, if the historical sciences are to play any role in fostering future resilience, improving our understanding of the processes and drivers that influence long-term, centennial-scale resilience is a necessary prerequisite5. Our approach has highlighted the global relationship between population change and frequency of disturbance over millennial timescales and applies across a broad geographical and chronological sample of past populations, including prehistoric examples that have been overlooked in systematic reviews of societal resilience more broadly. Improved clarity on the drivers of exposure frequency and type in the past will help to reveal the mechanism(s) behind the dynamics we describe and their potential limits, which is particularly important as environmental variability is predicted to increase in the future40,41. Archaeological and historical case studies have focused on the frequency of volcanism, warfare42 and hydroclimatic oscillations23,43, as well as the rate of abrupt or extreme events in general8. Comparable evidence on different categories of downturns is necessary. Synthesis of these or similar data in a comparative framework can provide important insights into the causal links between population resilience, risk of exposure and, ultimately, the ability to recover.

Methods

Calibration and aggregation

Archaeological radiocarbon dates were collated from published studies that previously adopted null hypothesis significance testing (NHST) approaches towards prehistoric demography. Our literature search identified 16 studies that collectively span six continents and approximately 30,000 radiocarbon years (Supplementary Tables 1 and 3–18). We applied a consistent protocol to the calibration of radiocarbon assays. The calibrate function within the R package rcarbon17 was used to convert 14C radiocarbon years to calendar years before present. The IntCal20 (ref. 44) and SHCal20 (ref. 45) curves were applied to dates in the northern and southern hemispheres, respectively. Calibrated ages are reported as the 95.4% (two-sigma, 2σ) age range. Laboratory codes and site identification numbers were appended to each calibrated date range and postcalibration distributions were not normalized. To account for between-site variation in sampling intensity, several dates from a single site that are within 50 radiocarbon years of each other were pooled (‘binned’) before aggregation into regional SPDs.

Bayesian model fitting

Our protocol aimed to replicate the results of the 16 identified case studies to the greatest extent possible. To reproduce statistically significant negative population events (‘downturns’) produced by rcarbon’s ‘modelTest()’ function in the original studies, while simultaneously addressing the well-known limitations of using summed probability distributions in NHST, we adopted an alternative, Bayesian modelling approach. Markov Chain Monte Carlo (MCMC) methods implemented in the nimbleCarbon package46 for radiocarbon data can obtain robust parameter estimates, accounting for radiocarbon measurement errors and sampling error simultaneously. Previously, this has been a major drawback of NHST approaches to aggregated radiocarbon data, with several alternatives proposed in the literature47–49. Using the MCMC-derived parameter estimates in posterior predictive checks enabled us to detect periods where expected growth trajectories were lower than the fitted model parameters and which were more robust and accurate than least-squares regression approaches applied directly on SPDs. The protocol produced outputs that are analogous to those in previous NHST studies (Extended Data Fig. 1).

We analysed regional SPDs separately by fitting identical bounded exponential growth models to each dataset. This common model is defined by three parameters: growth rate (r) and boundaries (a and b). A weakly informative exponentially distributed prior (λ = 1/000.4) was used for r to capture a broad range of potential growth rates among the cases. Parameters a and b were adopted from the original studies. Markov chain traceplots (Extended Data Figs. 2 and 3) evaluate chain mixing alongside model convergence (Gelman-Rubin Ȓ) and effective sample size diagnostics (Extended Data Table 5). Three chains of 50,000 iterations were run per region, with a burn-in of 5,000 iterations and a thinning interval of 2. To ensure comparability of results with published studies that used a logistic growth model as a null hypothesis, regional datasets were subset based on expert judgement at documented palaeodemographic transitions and treated as two separate exponential growth models. Subsetting was only performed on the Near East and Italy, Sicily and Sardinia datasets. Downturns adjacent to transition points were removed from the sample to avoid including data points introduced by the subsetting. Posterior predictive checks were executed using samples of parameter estimates obtained by the MCMC approach to simulate and back-calibrate a number of radiocarbon dates equal to the sample size. The procedure was repeated 1,000 times to derive critical envelopes.

Resilience metrics

The resilience metrics target periods when empirical SPD curves are below the 90% confidence envelopes of the fitted models, according to the posterior predictive checks (periods termed ‘downturns’). Extraction was performed using a bespoke R function modified from Riris and De Souza12, which is available at ref. 50. The principal response metrics in our analysis are resistance and resilience (Extended Data Fig. 4). Respectively, these metrics quantify the normalized response to downturns and the rate of recovery relative to baseline conditions. Resistance is measured on SPDs using two parameters: SPD values at the start of a downturn (b) and at downturn minima (x), whereas resilience is measured across the entire period of decline until its end (e) relative to the minimum and baseline (Table 1). Resistance ranges between 0 and 1, indicating a 100% change from baseline to no change. Resilience ranges between −1 and 1, with 1 indicating full recovery by the end of the downturn. Negative values of resilience indicate that the baseline value has been exceeded, although remaining outside the expectations of the null model. Finally, zero indicates no recovery11,51. Variations in the shape of the SPD can result under comparable scenarios (that is, with similar demographic dynamics, archaeological sample sizes and disturbances) and hence can produce different results in the metrics because of calibration effects despite their similarity. Although this issue may contribute to noise and error in measurement, it is nevertheless highly unlikely to be systematic or to correlate with other variables of interest. Finally, we conservatively only consider events greater than ten years in duration for our statistical modelling. These events are at an elevated risk of being artefacts of the null model rather than true downturns with an unclear driver.

We also collected information on the start and ends of downturns, their duration, elapsed time until SPD minima were reached, the cumulative number of a downturn and the frequency of downturns (Supplementary Table 2). The frequency of downturns is calculated on a per downturn basis within study regions, using the cumulative number of the downturn, normalized by its duration. We report frequency of downturns for each downturn as a logarithm per millennium.

Statistical modelling

The target of our comparative analysis is the resistance and resilience of human populations to disturbance as defined in each individual study. Our approach assumed that high values reflect resilient populations that successfully reestablish growth regimes after periods of decline related to disturbance events. We also assumed that downturns are randomly distributed in time and geographical space. To account for the influence of interregional and interevent cultural variation on outcomes, we drew on expert judgement and close readings of the published literature to record the broad category and specific type of disturbance during downturns, as well as the dominant land-use type and the nature of resulting socio-cultural changes, if any (Extended Data Table 1). These variables provided a control on whether a given population within a cultural system retained its identity and function over time, or whether system transformation and adaptive change is archaeologically evident.

Linear mixed-effect models were executed to evaluate the presence and strength of relationships between resistance, resilience, the recorded variables and case study locations. This analysis was performed using the cAIC4 and lme4 R packages52,53; scripts are available at ref. 50. Initial models were defined with resistance and resilience as response variables, with only case identifiers (region) as a random effect. As observed downturns were sequential within each case, the random effect controlled for potential pseudoreplication and avoided the need to weight the data by group size. A stepwise search using Akaike’s information criterion was implemented for investigating the information gain of including fixed effects in each model in turn. These candidate models were sequentially fitted using restricted maximum likelihood. Most fixed covariates (Extended Data Table 3) were left out of the final models. Region was retained as a random effect in all cases, to produce two models:resistance≈1∣region+frequencyofdownturn

resilience≈1∣region+frequencyofdownturn

Model output is summarized in Extended Data Table 3 and diagnostics are shown in Extended Data Figs. 5 and 6. We present standardized residuals, by region and in full, as well as leverage and Cook’s distance.

To further explore the relationship between rates of downturns and resistance and resilience, we performed an additional modelling exercise with the same random effect and full suite of fixed effects, with the frequency of downturn as the independent variable (Extended Data Table 4 and Extended Data Fig. 7).frequencyofdownturn≈1∣region+landuse+change+disturbancetype+pace

The effect sizes (standardized coefficients) of the significant model terms are plotted graphically in Fig. 3c and reported in full in Extended Data Table 4. We report effect sizes in the text as η2, that is, the total variance explained by differences between means.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Online content

Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at 10.1038/s41586-024-07354-8.

Supplementary information

Supplementary Information

Supplementary Tables 1–18.

Peer Review File

Supplementary Table 1 Summary of attributes of regions used in this study.

Supplementary Table 2 Resistance and resilience metrics extracted from periods of population downturn detected in posterior predictive checks.

Supplementary Table 3 Radiocarbon data for Arid Zone Australia.

Supplementary Table 4 Radiocarbon data for Caribbean Archipelago.

Supplementary Table 5 Radiocarbon data for China Central Plains.

Supplementary Table 6 Radiocarbon data for Circumpolar Norway.

Supplementary Table 7 Radiocarbon data for Greece.

Supplementary Table 8 Radiocarbon data for Highland and Coastal Peru.

Supplementary Table 9 Radiocarbon data for Italy, Sardinia and Sicily.

Supplementary Table 10 Radiocarbon data for Korean Peninsula.

Supplementary Table 11 Radiocarbon data for Near East.

Supplementary Table 12 Radiocarbon data for Southeastern Norway.

Supplementary Table 13 Radiocarbon data for South Africa (GCFR).

Supplementary Table 14 Radiocarbon data for South Africa (SRZ).

Supplementary Table 15 Radiocarbon data for Tropical Lowlands.

Supplementary Table 16 Radiocarbon data for Utah.

Supplementary Table 17 Radiocarbon data for Wyoming.

Supplementary Table 18 Radiocarbon data for Yukon.

Extended data figures and tables

Extended Data Fig. 1 Posterior predictive checks for 16 study regions.

Regions shaded in blue indicate period below modelled growth trajectories: downturns.

Extended Data Fig. 2 Traceplots of chain mixing for the MCMC of each study region and subset.

There is adequate mixing and convergence across chains.

Extended Data Fig. 3 Traceplots of chain mixing for the MCMC of each study region and subset.

There is adequate mixing and convergence across chains.

Extended Data Fig. 4 Resistance and Resilience as functions of variation in minima (x) and end-points (e) when baselines (b) are held constant (0.02).

(A) Resistance equals zero for two different values of x, when x = 0 and when x is two times b. Negative resistance is theoretically possible but does not occur on SPDs, (B) Resilience as a function of x, (C) Resilience with a varying end point (e), and constant start (b) and minimum (x = 0.01). Resilience equals 1 only when e = b, (D) Interpretative scatterplot of indicative resistance-resilience outcomes.

Extended Data Fig. 5 Model diagnostics for Model I (Resistance).

(A) Standardised residuals against fitted values for each study region. (B) All residuals versus fitted values. (C) Standardised residuals by study region for n = 154 independent samples across 16 regions. The lower and upper hinges correspond to the 25th and 75th percentiles The upper and lower whiskers extend from the hinges to 1.5 * IQR (where IQR is the inter-quartile range, distance between 25th and 75th percentiles). Data beyond the whiskers are individually plotted outlying points. The black points represent the group median. (D) Leverage and Cook’s Distance for 154 observations.

Extended Data Fig. 6 Model diagnostics for Model II (Resilience).

(A) Standardised residuals against fitted values for each study region. (B) All residuals versus fitted values. (C) Standardised residuals by study region for n = 154 independent samples across 16 regions. The lower and upper hinges correspond to the 25th and 75th percentiles The upper and lower whiskers extend from the hinges to 1.5 * IQR (where IQR is the inter-quartile range, distance between 25th and 75th percentiles). The black points represent the group median. Data beyond the whiskers are individually plotted outlying points. (D) Leverage and Cook’s Distance for 154 observations.

Extended Data Fig. 7 Model diagnostics for Model III.

(A) Standardised residuals against fitted values for each study region. (B) All residuals versus fitted values. (C) Standardised residuals by study region for n = 154 independent samples across 16 regions. The lower and upper hinges correspond to the 25th and 75th percentiles The upper and lower whiskers extend from the hinges to 1.5 * IQR (where IQR is the inter-quartile range, distance between 25th and 75th percentiles). The black points represent the group median. Data beyond the whiskers are individually plotted outlying points. (D) Leverage and Cook’s Distance for 154 observations.

Extended Data Table 1 Summary of disturbance categories and sub-types collected from published palaeodemographic studies

Extended Data Table 2 Initial mixed-effect model diagnostics for Resistance and Resilience with Study Region included as a random effect

Significant terms (p < 0.05) based on a two-sided test in bold.

Extended Data Table 3 Final mixed-effect model diagnostics for Resistance and Resilience with Study Region included as a random effect and Frequency of downturn as a fixed effect

Significant terms (p < 0.05) based on a two-sided test in bold.

Extended Data Table 4 Mixed-effect diagnostics for Rate of downturns with Study Region included as a random effect and Land Use, Change, Disturbance Type, and Pace as a fixed effects

Significant terms (p < 0.05) based on a two-sided test in bold.

Extended Data Table 5 Convergence (Gelman-Rubin Ȓ) and Effective Sample Size diagnostics for MCMC fits

The Gelman-Rubin convergence diagnostic provides a convergence summary based on multiple chains, in our case, three per MCMC fit.

Extended data

is available for this paper at 10.1038/s41586-024-07354-8.

Supplementary information

The online version contains supplementary material available at 10.1038/s41586-024-07354-8.

Acknowledgements

We thank P. Lockwood for comments on an earlier version of this work. We acknowledge support from the Arts and Humanities Research Council grant AH/X002217/1 (P.R.), Samsung Electronics grant A0342-20220007 (J.B.), Leverhulme Trust grant no. PLP-2019–304 (E.C.) and the Youth Innovation Promotion Association of the Chinese Academy of Sciences grant YIPA-CAS, 2022149 (X.R.).

Author contributions

P.R. and F.S. conceptualized the study. P.R., E.C., A.P. and F.S. developed the methodology. P.R., B.D., E.K.J., Y.O., A.P., X.R., E.R., P.E.S. and S.S. carried out the investigation. P.R., E.C. and F.S. did the analysis. P.R., J.C.F. and P.E.S. wrote the article. P.R., J.B., E.C., J.C.F., E.K.J., S.Y.M., A.P., E.R., P.E.S., F.S. and S.S. edited the article. P.R., S.Y.M. and E.R. carried out the visualization.

Peer review

Peer review information

Nature thanks Marko Porčić and John Haldon for their contribution to the peer review of this work. Peer reviewer reports are available.

Data availability

The data supporting the findings of this study are available via Zenodo at 10.5281/zenodo.10061467 (ref. 53).

Code availability

The code supporting the findings of this study is available via Zenodo at 10.5281/zenodo.10061467 (ref. 53).

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Cumming GS Peterson GD Unifying research on social–ecological resilience and collapse Trends Ecol. Evol. 2017 32 695 713 10.1016/j.tree.2017.06.014 28734593
Cumming, G. S. & Peterson, G. D. Unifying research on social–ecological resilience and collapse. Trends Ecol. Evol. 32, 695–713 (2017).28734593 10.1016/j.tree.2017.06.014
2. Haldon J History meets palaeoscience: consilience and collaboration in studying past societal responses to environmental change Proc. Natl Acad. Sci. USA 2018 115 3210 3218 10.1073/pnas.1716912115 29531084
Haldon, J. et al. History meets palaeoscience: consilience and collaboration in studying past societal responses to environmental change. Proc. Natl Acad. Sci. USA 115, 3210–3218 (2018).29531084 10.1073/pnas.1716912115
3. IPBES Global assessment report on biodiversity and ecosystem services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services. Zenodo10.5281/zenodo.3831673 (2019).
4. Bradtmöller M Grimm S Riel-Salvatore J Resilience theory in archaeological practice–an annotated review Quat. Int. 2017 446 3 16 10.1016/j.quaint.2016.10.002
Bradtmöller, M., Grimm, S. & Riel-Salvatore, J. Resilience theory in archaeological practice–an annotated review. Quat. Int. 446, 3–16 (2017).10.1016/j.quaint.2016.10.002
5. Silva F Developing transdisciplinary approaches to sustainability challenges: the need to model socio-environmental systems in the longue durée Sustainability 2022 14 10234 10.3390/su141610234
Silva, F. et al. Developing transdisciplinary approaches to sustainability challenges: the need to model socio-environmental systems in the longue durée. Sustainability 14, 10234 (2022).10.3390/su141610234
6. Degroot D Towards a rigorous understanding of societal responses to climate change Nature 2021 591 539 550 10.1038/s41586-021-03190-2 33762769
Degroot, D. et al. Towards a rigorous understanding of societal responses to climate change. Nature 591, 539–550 (2021).33762769 10.1038/s41586-021-03190-2
7. Broska LH Poganietz WR Vögele S Extreme events defined—a conceptual discussion applying a complex systems approach Futures 2020 1 102490 10.1016/j.futures.2019.102490
Broska, L. H., Poganietz, W. R. & Vögele, S. Extreme events defined—a conceptual discussion applying a complex systems approach. Futures 1, 102490 (2020).10.1016/j.futures.2019.102490
8. Pausas JG Leverkus AB Disturbance ecology in human societies People Nat. 2023 5 1082 1093 10.1002/pan3.10471
Pausas, J. G. & Leverkus, A. B. Disturbance ecology in human societies. People Nat. 5, 1082–1093 (2023).10.1002/pan3.10471
9. Middleton GD The show must go on: collapse, resilience, and transformation in 21st-century archaeology Rev. Anthropol. 2017 46 78 105 10.1080/00938157.2017.1343025
Middleton, G. D. The show must go on: collapse, resilience, and transformation in 21st-century archaeology. Rev. Anthropol. 46, 78–105 (2017).10.1080/00938157.2017.1343025
10. Jackson RC Dugmore AJ Riede F Rediscovering lessons of adaptation from the past Glob. Environ. Change 2018 52 58 65 10.1016/j.gloenvcha.2018.05.006
Jackson, R. C., Dugmore, A. J. & Riede, F. Rediscovering lessons of adaptation from the past. Glob. Environ. Change 52, 58–65 (2018).10.1016/j.gloenvcha.2018.05.006
11. Van Meerbeek K Jucker T Svenning JC Unifying the concepts of stability and resilience in ecology J. Ecol. 2021 109 3114 3132 10.1111/1365-2745.13651
Van Meerbeek, K., Jucker, T. & Svenning, J. C. Unifying the concepts of stability and resilience in ecology. J. Ecol. 109, 3114–3132 (2021).10.1111/1365-2745.13651
12. Riris P De Souza JG Formal tests for resistance-resilience in archaeological time series Front. Ecol. Evol. 2021 9 740629 10.3389/fevo.2021.740629
Riris, P. & De Souza, J. G. Formal tests for resistance-resilience in archaeological time series. Front. Ecol. Evol. 9, 740629 (2021).10.3389/fevo.2021.740629
13. Shennan S Regional population collapse followed initial agriculture booms in mid-Holocene Europe Nat. Commun. 2013 4 2486 10.1038/ncomms3486 24084891
Shennan, S. et al. Regional population collapse followed initial agriculture booms in mid-Holocene Europe. Nat. Commun. 4, 2486 (2013).24084891 10.1038/ncomms3486
14. Bevan A Holocene fluctuations in human population demonstrate repeated links to food production and climate Proc. Natl Acad. Sci. USA 2017 114 E10524 E10531 10.1073/pnas.1709190114 29158411
Bevan, A. et al. Holocene fluctuations in human population demonstrate repeated links to food production and climate. Proc. Natl Acad. Sci. USA 114, E10524–E10531 (2017).29158411 10.1073/pnas.1709190114
15. Freeman J Mauldin RP Whisenhunt M Hard RJ Anderies JM Repeated long-term population growth overshoots and recessions among hunter-gatherers The Holocene 2023 7 09596836231183072
Freeman, J., Mauldin, R. P., Whisenhunt, M., Hard, R. J. & Anderies, J. M. Repeated long-term population growth overshoots and recessions among hunter-gatherers. The Holocene 7, 09596836231183072 (2023).
16. Freeman J Byers DA Robinson E Kelly RL Culture process and the interpretation of radiocarbon data Radiocarbon 2018 60 453 467 10.1017/RDC.2017.124
Freeman, J., Byers, D. A., Robinson, E. & Kelly, R. L. Culture process and the interpretation of radiocarbon data. Radiocarbon 60, 453–467 (2018).10.1017/RDC.2017.124
17. Crema ER Bevan A Inference from large sets of radiocarbon dates: software and methods Radiocarbon 2021 63 23 39 10.1017/RDC.2020.95
Crema, E. R. & Bevan, A. Inference from large sets of radiocarbon dates: software and methods. Radiocarbon 63, 23–39 (2021).10.1017/RDC.2020.95
18. Schauer P Supply and demand in prehistory? Economics of Neolithic mining in northwest Europe J. Anthropol. Archaeol. 2019 54 149 160 10.1016/j.jaa.2019.03.001
Schauer, P. et al. Supply and demand in prehistory? Economics of Neolithic mining in northwest Europe. J. Anthropol. Archaeol. 54, 149–160 (2019).10.1016/j.jaa.2019.03.001
19. Bird D p3k14c, a synthetic global database of archaeological radiocarbon dates Sci. Data. 2022 9 27 27 10.1038/s41597-022-01118-7 35087092
Bird, D. et al. p3k14c, a synthetic global database of archaeological radiocarbon dates. Sci. Data. 27, 27 (2022).35087092 10.1038/s41597-022-01118-7
20. Palmisano A Lawrence D de Gruchy MW Bevan A Shennan S Holocene regional population dynamics and climatic trends in the Near East: a first comparison using archaeo-demographic proxies Quat. Sci. Rev. 2021 252 106739 10.1016/j.quascirev.2020.106739
Palmisano, A., Lawrence, D., de Gruchy, M. W., Bevan, A. & Shennan, S. Holocene regional population dynamics and climatic trends in the Near East: a first comparison using archaeo-demographic proxies. Quat. Sci. Rev. 252, 106739 (2021).10.1016/j.quascirev.2020.106739
21. Koch A Brierley C Maslin MM Lewis SL Earth system impacts of the European arrival and Great Dying in the Americas after 1492 Quat. Sci. Rev. 2019 207 13 36 10.1016/j.quascirev.2018.12.004
Koch, A., Brierley, C., Maslin, M. M. & Lewis, S. L. Earth system impacts of the European arrival and Great Dying in the Americas after 1492. Quat. Sci. Rev. 207, 13–36 (2019).10.1016/j.quascirev.2018.12.004
22. Storey, R. & Storey, G. R. Rome and the Classic Maya: Comparing the Slow Collapse of Civilizations (Routledge, 2017).
23. Finley JB Robinson E DeRose RJ Hora E Multidecadal climate variability and the florescence of Fremont societies in Eastern Utah Am. Antiq. 2020 85 93 112 10.1017/aaq.2019.79
Finley, J. B., Robinson, E., DeRose, R. J. & Hora, E. Multidecadal climate variability and the florescence of Fremont societies in Eastern Utah. Am. Antiq. 85, 93–112 (2020).10.1017/aaq.2019.79
24. Freeman J Landscape engineering impacts the long-term stability of agricultural populations Hum. Ecol. 2021 49 369 382 10.1007/s10745-021-00242-z
Freeman, J. et al. Landscape engineering impacts the long-term stability of agricultural populations. Hum. Ecol. 49, 369–382 (2021).10.1007/s10745-021-00242-z
25. Boivin N Crowther A Mobilizing the past to shape a better Anthropocene Nat. Ecol. Evol. 2021 5 273 284 10.1038/s41559-020-01361-4 33462488
Boivin, N. & Crowther, A. Mobilizing the past to shape a better Anthropocene. Nat. Ecol. Evol. 5, 273–284 (2021).33462488 10.1038/s41559-020-01361-4
26. Burke A The archaeology of climate change: the case for cultural diversity Proc. Natl Acad. Sci. USA 2021 118 e2108537118 10.1073/pnas.2108537118 34301875
Burke, A. et al. The archaeology of climate change: the case for cultural diversity. Proc. Natl Acad. Sci. USA 118, e2108537118 (2021).34301875 10.1073/pnas.2108537118
27. Seidl R Rammer W Spies TA Disturbance legacies increase the resilience of forest ecosystem structure, composition, and functioning Ecol. Appl. 2014 24 2063 2077 10.1890/14-0255.1 27053913
Seidl, R., Rammer, W. & Spies, T. A. Disturbance legacies increase the resilience of forest ecosystem structure, composition, and functioning. Ecol. Appl. 24, 2063–2077 (2014).27053913 10.1890/14-0255.1
28. Zahid HJ Robinson E Kelly RL Agriculture, population growth, and statistical analysis of the radiocarbon record Proc. Natl Acad. Sci. USA 2016 113 931 935 10.1073/pnas.1517650112 26699457
Zahid, H. J., Robinson, E. & Kelly, R. L. Agriculture, population growth, and statistical analysis of the radiocarbon record. Proc. Natl Acad. Sci. USA 113, 931–935 (2016).26699457 10.1073/pnas.1517650112
29. Shennan S Sear R Archaeology, demography and life history theory together can help us explain past and present population patterns Philos. Trans. R. Soc. Lond 2021 376 20190711 10.1098/rstb.2019.0711 33250034
Shennan, S. & Sear, R. Archaeology, demography and life history theory together can help us explain past and present population patterns. Philos. Trans. R. Soc. Lond. B 376, 20190711 (2021).33250034 10.1098/rstb.2019.0711
30. Rockman, M. in Macroevolution in Human Prehistory: Evolutionary Theory and Processual Archaeology (eds Prentiss, A., Kuijit, I. & Chatters, J. C.) 51–71 (Springer, 2009).
31. Galan J Landscape adaptation to climate change: local networks, social learning and co-creation processes for adaptive planning Glob. Environ. Change 2023 78 102627 10.1016/j.gloenvcha.2022.102627
Galan, J. et al. Landscape adaptation to climate change: local networks, social learning and co-creation processes for adaptive planning. Glob. Environ. Change 78, 102627 (2023).10.1016/j.gloenvcha.2022.102627
32. De Souza JG Climate change and cultural resilience in late pre-Columbian Amazonia Nat. Ecol. Evol. 2019 3 1007 1017 10.1038/s41559-019-0924-0 31209292
De Souza, J. G. et al. Climate change and cultural resilience in late pre-Columbian Amazonia. Nat. Ecol. Evol. 3, 1007–1017 (2019).31209292 10.1038/s41559-019-0924-0
33. Cole LE Bhagwat SA Willis KJ Recovery and resilience of tropical forests after disturbance Nat. Commun. 2014 20 3906 10.1038/ncomms4906
Cole, L. E., Bhagwat, S. A. & Willis, K. J. Recovery and resilience of tropical forests after disturbance. Nat. Commun. 20, 3906 (2014).10.1038/ncomms4906
34. Cant J Capdevila P Beger M Salguero‐Gómez R Recent exposure to environmental stochasticity does not determine the demographic resilience of natural populations Ecol. Lett. 2023 26 1186 1199 10.1111/ele.14234 37158011
Cant, J., Capdevila, P., Beger, M. & Salguero‐Gómez, R. Recent exposure to environmental stochasticity does not determine the demographic resilience of natural populations. Ecol. Lett. 26, 1186–1199 (2023).37158011 10.1111/ele.14234
35. Redman CL Resilience theory in archaeology Am. Anthropol. 2005 107 70 77 10.1525/aa.2005.107.1.070
Redman, C. L. Resilience theory in archaeology. Am. Anthropol. 107, 70–77 (2005).10.1525/aa.2005.107.1.070
36. French JC Riris P Fernandez-Lopez de Pablo J Lozano S Silva F A manifesto for palaeodemography in the twenty-first century Philos. Trans. R. Soc. Lond. B 2021 376 20190707 10.1098/rstb.2019.0707 33250019
French, J. C., Riris, P., Fernandez-Lopez de Pablo, J., Lozano, S. & Silva, F. A manifesto for palaeodemography in the twenty-first century. Philos. Trans. R. Soc. Lond. B 376, 20190707 (2021).33250019 10.1098/rstb.2019.0707
37. Freeman, J. et al. The long-term expansion and recession of human populations. Proc. Natl Acad. Sci. USA 121, e2312207121 (2024).
38. Allen KJ Coupled insights from the palaeoenvironmental, historical and archaeological archives to support social-ecological resilience and the sustainable development goals Environ. Res. Lett. 2022 17 055011 10.1088/1748-9326/ac6967
Allen, K. J. et al. Coupled insights from the palaeoenvironmental, historical and archaeological archives to support social-ecological resilience and the sustainable development goals. Environ. Res. Lett. 17, 055011 (2022).10.1088/1748-9326/ac6967
39. Schug GR Climate change, human health, and resilience in the Holocene Proc. Natl Acad. Sci. USA 2023 120 e2209472120 10.1073/pnas.2209472120 36649426
Schug, G. R. et al. Climate change, human health, and resilience in the Holocene. Proc. Natl Acad. Sci. USA 120, e2209472120 (2023).36649426 10.1073/pnas.2209472120
40. Wisner, B. G., Blaikie, P., Cannon, T. & Davis, I. At Risk: Natural Hazards, People’s Vulnerability and Disasters (Routledge, 2014).
41. Thornton PK Ericksen PJ Herrero M Challinor AJ Climate variability and vulnerability to climate change: a review Glob. Change Biol. 2014 20 3313 3328 10.1111/gcb.12581
Thornton, P. K., Ericksen, P. J., Herrero, M. & Challinor, A. J. Climate variability and vulnerability to climate change: a review. Glob. Change Biol. 20, 3313–3328 (2014).10.1111/gcb.12581
42. Gao C Volcanic climate impacts can act as ultimate and proximate causes of Chinese dynastic collapse Commun. Earth Environ. 2021 2 234 10.1038/s43247-021-00284-7
Gao, C. et al. Volcanic climate impacts can act as ultimate and proximate causes of Chinese dynastic collapse. Commun. Earth Environ. 2, 234 (2021).10.1038/s43247-021-00284-7
43. Douglas PM Demarest AA Brenner M Canuto MA Impacts of climate change on the collapse of lowland Maya civilization Annu. Rev. Earth Planet Sci. 2016 44 613 645 10.1146/annurev-earth-060115-012512
Douglas, P. M., Demarest, A. A., Brenner, M. & Canuto, M. A. Impacts of climate change on the collapse of lowland Maya civilization. Annu. Rev. Earth Planet Sci. 44, 613–645 (2016).10.1146/annurev-earth-060115-012512
44. Reimer PJ The IntCal20 northern hemisphere radiocarbon age calibration curve (0–55 cal kBP) Radiocarbon 2020 62 725 757 10.1017/RDC.2020.41
Reimer, P. J. et al. The IntCal20 northern hemisphere radiocarbon age calibration curve (0–55 cal kBP). Radiocarbon 62, 725–757 (2020).10.1017/RDC.2020.41
45. Hogg AG SHCal20 southern hemisphere calibration, 0–55,000 years cal BP Radiocarbon 2020 62 759 778 10.1017/RDC.2020.59
Hogg, A. G. et al. SHCal20 southern hemisphere calibration, 0–55,000 years cal BP. Radiocarbon 62, 759–778 (2020).10.1017/RDC.2020.59
46. Crema, E. R. nimbleCarbon (v.0.2.1): models and utility functions for Bayesian analyses of radiocarbon dates with NIMBLE. GitHub https://github.com/ercrema/nimbleCarbon (2022).
47. Carleton WC Evaluating Bayesian radiocarbon‐dated event count (REC) models for the study of long‐term human and environmental processes J. Quat. Sci. 2021 36 110 123 10.1002/jqs.3256
Carleton, W. C. Evaluating Bayesian radiocarbon‐dated event count (REC) models for the study of long‐term human and environmental processes. J. Quat. Sci. 36, 110–123 (2021).10.1002/jqs.3256
48. Timpson A Barberena R Thomas MG Méndez C Manning K Directly modelling population dynamics in the South American Arid Diagonal using 14C dates Philos. Trans. R. Soc. Lond. B 2021 376 20190723 10.1098/rstb.2019.0723 33250032
Timpson, A., Barberena, R., Thomas, M. G., Méndez, C. & Manning, K. Directly modelling population dynamics in the South American Arid Diagonal using 14C dates. Philos. Trans. R. Soc. Lond. B 376, 20190723 (2021).33250032 10.1098/rstb.2019.0723
49. Crema ER Statistical inference of prehistoric demography from frequency distributions of radiocarbon dates: a review and a guide for the perplexed J. Archaeol. Method Theory 2022 29 1387 1418 10.1007/s10816-022-09559-5
Crema, E. R. Statistical inference of prehistoric demography from frequency distributions of radiocarbon dates: a review and a guide for the perplexed. J. Archaeol. Method Theory 29, 1387–1418 (2022).10.1007/s10816-022-09559-5
50. Riris, P. Data and scripts for the paper ‘Frequent disturbances enhance the resistance and recovery of past human populations’. Zenodo10.5281/zenodo.10061467 (2023).
51. Orwin KH Wardle DA New indices for quantifying the resistance and resilience of soil biota to exogenous disturbances Soil Biol. Biochem. 2004 36 1907 1912 10.1016/j.soilbio.2004.04.036
Orwin, K. H. & Wardle, D. A. New indices for quantifying the resistance and resilience of soil biota to exogenous disturbances. Soil Biol. Biochem. 36, 1907–1912 (2004).10.1016/j.soilbio.2004.04.036
52. Bates D Mächler M Bolker B Walker S Fitting linear mixed-effects models using lme4 J. Stat. Softw. 2015 67 1 48 10.18637/jss.v067.i01
Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).10.18637/jss.v067.i01
53. Säfken B Rügamer D Kneib T Greven S Conditional model selection in mixed-effects models with cAIC4 J. Stat. Softw 2021 99 1 30 10.18637/jss.v099.i08
Säfken, B., Rügamer, D., Kneib, T. & Greven, S. Conditional model selection in mixed-effects models with cAIC4. J. Stat. Softw. 99, 1–30 (2021).10.18637/jss.v099.i08
