Philippa Kaur

Spatial regime shifts over time and space
To assess how the dominance of large predatory fish (Eurasian perch Perca fluviatilis, northern pike Esox lucius) and three-spined stickleback (Gasterosteus aculeatus) changed over time and space, we collated juvenile fish abundance data from 13073 samplings conducted during 39 years (1979–2017) in 486 bays along a 1200 km stretch of the Swedish Baltic Sea coast. We used juvenile fish surveys because they (i) include stickleback, (ii) were conducted along the entire coast, including the outer archipelago, and (iii) capture patterns of recruitment failure; a proposed driver of local perch and pike stock decline34,35,37,45. The samplings were conducted by various monitoring programs and research projects to quantify fish recruitment. Much of the data was extracted from the Swedish national database for coastal fish (http://www.slu.se/kul). Other fish species also occurred in the data, but we—like others25,34,45—focused on the most common and strongly interacting species.
The timing and placement of most of the fish surveys was not chosen with this study in mind. We therefore included data from as many surveys as possible, ensuring that the dataset covered (i) gradients in distance to open sea and wave exposure, (ii) the entire Swedish east coast, and (iii) the longest time period possible (Supplementary Fig. 7). Nearly all sampling (94%) was conducted during July-September. To achieve the best possible spatial coverage we also included some bays only sampled in October and (for one bay) June. Initial exploration of data from bays sampled monthly from June to October suggested there were no large differences in the October and June data. Nine of the 486 bays occurred much further into the archipelago (49–67 km) than the rest (20 km from the open sea prior to 1995, we only used data from 1995–2017. Finally, we used a general linear multiple regression model to test how sampling year, wave exposure (two levels) and their interaction explained the maximum distance.
The bimodality in relative predator dominance seen over time and space (Figs. 1–3) could in theory be caused by variability in local abiotic conditions, and not by predator–prey interactions. We therefore tested whether any of four abiotic conditions estimated locally—water surface temperature, salinity, turbidity, water depth—could explain the variability in deviance residuals from the binomial glm (Supplementary Fig. 8), using a regular linear model (after assessing model assumptions and ruling out multicollinearity; see above). Salinity and turbidity (but not depth and temperature) had statistically significant but weak influences (R2 = 0.11), and the bimodality clearly remained (Supplementary Fig. 9).
Finally, we explored how time, distance to open sea, wave exposure, latitude and their two- and three-way interactions influenced log-transformed abundances of (i) predators (perch and pike pooled) and (ii) stickleback, using general linear models. We identified the most parsimonious models as outlined above.
Temporal regime shift at Forsmark: using the 34-year Forsmark time series, we first tested for temporal breakpoints in logit-transformed relative predator dominance data using change-point detection (strucchange) for linear models47. This method estimates the optimal number and (if identified) position of breakpoints using the Bayesian Information Criterion (BIC). Second, we explored what temporal change(s) in perch and stickleback abundances that preceded the shift. Because of highly non-linear patterns we modeled the temporal changes using generalized additive models (GAM) as implemented in the mgvc package68.
Importance of predator–prey reversal for fish recruitment: to assess the relative importance of predator–prey reversal for perch, pike and stickleback recruitment, we used statistical model selection based on path analyses; a form of structural equation modeling that can be used to tease apart direct vs. indirect (mediated) relationships between multiple ( >2) variables, and thereby assess the relative importance of direct vs. indirect relationships in systems48,69. Initial data exploration using multiple regression showed that one bay was a clear outlier due to 0 juvenile perch and pike, generating (i) too high leverage (influence on statistical relationships), (ii) heteroscedasticity and (iii) non-normally distributed errors. Since we suspected that juveniles had already migrated out of this bay, the bay was excluded (resulting in N = 31). Removing this statistical outlier resulted in that the model fulfilled test assumptions and the overall fit more than doubled (adjusted R2 increased from 0.17 to 0.37).
Based on ecological knowledge of the study system, we expressed 14 multivariate hypotheses of the direct and indirect drivers of perch and stickleback recruitment, as graphical network models of interacting paths36,49. Due to the relatively low sample size we restricted the number of paths to 7. The two simplest models assumed that perch+pike and stickleback juvenile abundance in summer (i.e. recruitment) was influenced by adult abundance in spring (stock recruitment) and cumulative cover of rooted vegetation, while adult abundance was explained by spring cumulative cover of all vegetation species36,46, and distance to open sea or wave exposure46. The more complex models included combinations of known predator–prey interactions: perch and pike controlling adult stickleback in spring through predation36, stickleback feeding on juvenile perch and pike45, and stickleback competing with juvenile perch for zooplankton prey37. We then analyzed each model using piecewise path analysis as implemented in the piecewiseSEM package48. First, we tested the goodness of fit of each model to the data using Shipley’s test of directional separation (D-sep)70. If missing paths were identified, they were included in a new model. For the models that fitted the data (p  > 0.05) we used the Akaike’s Information Criterion corrected for small samples (AICc) (calculated using Shipley’s general approach to calculate AIC for path analysis71) to compare relative model fit. A summary of all candidate models and details of the best-fitting model are presented in Supplementary Tables 1 and 2, respectively. The strength of paths in the best-fitting model are presented using standardized path coefficients, which (based on the best-fitting model, see Fig. 3) mean that 1 SD increase in pooled adult perch and pike abundance reduces adult stickleback abundance by 0.53 SD. We also calculated the amount of variation (R2) in adult stickleback abundance, pooled juvenile perch and pike, and juvenile stickleback abundance, that was explained by the paths. Finally, we tested whether the relative predator dominance of adult and juvenile fish in this smaller dataset (N = 31) was also uni- or bimodal, using Hartigan’s dip test (for details, see above).
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article. More

A convincing argument that weather influences COVID-19 can be formulated in three parts: (1) experimental data suggest SARS-CoV-2 persistence on surfaces or in the air is sensitive to temperature, humidity, and ultraviolet light; (2) other environmentally sensitive respiratory viruses are seasonal, and more common in winter; and therefore, (3) climatic effects could be protective over space (hot, dry places might have less transmission) and time (summer might see reduced transmission compared to winter). All three are plausible and are generally consistent, but in many places (including, and especially, on social media), the basic premise of each has been communicated to the public, and policymakers, in a way that obscures key nuance and creates false confidence.
Experimental evidence shows that SARS-CoV-2 is environmentally sensitive
Like other viruses with a lipid envelope, SARS-CoV-2 is probably sensitive to temperature, humidity, and solar radiation; this affects its ability to persist on surfaces and in air, and might have subtle impacts on transmission. But the finer details of microbiology are often lost, leading to false confidence in how lab studies could scale up to the real world. For example, studies showing that germicidal ultraviolet radiation in hospitals and laboratories (ultraviolet C (UV-C) wavelengths) kills the virus have been misconstrued as evidence that sunlight (a mix of UV-A and UV-B) would effectively neutralize the virus in outdoor public spaces, possibly at a scale detectable from case data1. Newer experimental evidence supports the hypothesis that sunlight might also have an effect on SARS-CoV-22, although a global study found only a 1% reduction in transmission linked to environmental UV radiation3. In the real world, these effects will be slight, and unlikely to set hard limits on transmission anywhere in the world4.
Other upper respiratory tract infections are seasonal, with declines in warmer months
Influenza, the common cold, and other respiratory infections show seasonal transmission that coincides with changes in temperature, humidity, and solar radiation. But seasonal epidemics are also a product of the transmissibility of a virus, the initial susceptibility of a population, and the degree and nature of immunity conferred by infections. In basic epidemiological models, stable “oscillations” like seasonal epidemic waves usually require some degree of immunity5,6; at the start of a pandemic, when transmissibility is high and immunity is low, even strong environmental drivers are unlikely to curb transmission. Previous influenza pandemics show the importance of this nuance7. “Seasonal” versus “pandemic” influenza refers not just to different epidemic phases, but entirely different viral strains, and population susceptibility to pandemic strains starts high enough for rapid spread regardless of the season. During the first wave of the 2009 A/H1N1 pandemic, epidemic growth was still possible in August, the most environmentally unfavorable point in the year, with immunity under 20%; months later, immunity—and consequently environmental sensitivity—may eventually have been high enough to lead to a winter-driven third wave8. Scientists anticipate a similar pattern for COVID-19: while the virus could develop seasonal oscillations if it becomes endemic9 (i.e., if pandemic control fails in the long term), current susceptibility to SARS-CoV-2 is high enough that summer weather is unlikely to be protective10.
Environmental drivers could plausibly create seasonal or geographic differences in COVID-19 outbreak intensity
But those drivers’ impacts are heavily confounded by immunity, interventions, human behavior, and other details that are usually left out of models, leading to potentially spurious conclusions. For example, most available contact tracing data indicate that the proportion of indoor transmission is high11,12, a pattern likely caused by a combination of social contact patterns (including both number, intensity, and duration of contacts), air circulation, and potential weather drivers like sunlight or humidity. However, when studies attempt to model links between temperature and transmission, they almost always use gridded climate data or local weather data that represents the outdoors, and is unrepresentative of indoor conditions a virus particle or aerosolized cloud would actually experience in most transmission events.
Perhaps the biggest confounder, social behavior is environmentally driven and seasonal, but is rarely weighed alongside environmental and immunity drivers as a hypothesis for why infectious diseases show seasonality13,14. For example, school terms are seasonal and have a marked influence on social mixing patterns relevant to influenza transmission, even in pandemics15,16,17. Without individual-level transmission data, it can be difficult to distinguish direct biological impacts of weather from behaviorally mediated seasonality, and in some cases, the two are blurred (e.g., vitamin D levels are driven by both weather and seasonal behavior). Confusing the two could easily lead to spurious predictions. If behavioral patterns become unpredictable—either because of externalities like social distancing restrictions or a feedback loop between science and public risk perceptions around seasonality—attempts at forecasting the pandemic based on environmental seasonality will only become more unreliable. More

NEWS Q&A
27 August 2020

The spill released a new type of low-sulfur fuel, and its ecological effects aren’t well studied, says environment advocate Jaqueline Sauzier.

Dyani Lewis

Search for this author in:

Around 1,000 tonnes of oil have leaked from the MV Wakashio off Mauritius since 6 August.Credit: Pierre Dalais/EPA-EFE/Shutterstock

When the cargo ship MV Wakashio ran aground on a coral reef on the southeast tip of Mauritius, in the Indian Ocean almost exactly a month ago, it unleashed a vast oil spill. The Japanese-owned vessel held 200 tonnes of diesel and 3,900 tonnes of fuel oil, an estimated 1,000 tonnes of which leaked into the sea when the ship’s hull cracked on 6 August. It is the first reported spill of a new type of low-sulfur fuel that has been introduced to reduce air pollution. The spill has left a 15-kilometre stretch of the coastline — an internationally recognized biodiversity hotspot — smeared with oil.
Jacqueline Sauzier, president of the non-profit Mauritius Marine Conservation Society in Phoenix, has been helping with volunteer efforts to contain the spill. She spoke to Nature about how the clean-up is progressing.
What has been the response to the spill?
Mauritius is not geared up to deal with a catastrophe of this size, so other countries have sent experts to help. A French team arrived first, from the nearby island Réunion, to erect ocean booms — floating structures that contain the spill. The United Nations sent a team including experts in oil spills and crisis management. They’ve been working with communities, the private sector and the government to coordinate clean-up efforts. Marine ecologists and others have arrived from Japan and the United Kingdom.
Mauritians were also very proactive. In one weekend, we made nearly 80 kilometres of make-shift ocean booms out of cane trash — the leftover leaves and waste from sugar-cane processing — to contain the oil. Empty bottles were put in the middle of the booms to make them float, and anchors were attached to keep them from drifting away with the current.
For ten days, people worked night and day to contain as much oil as possible so that it wouldn’t reach the shoreline, where it is more difficult to clean. We managed to contain and remove nearly 75% of the spilled oil. Only a small amount reached the shore. But there’s still the issue of water-soluble chemicals that come from the oil, but dissolve into the water and therefore aren’t scooped out with the oil that sits on the water’s surface.
What ecosystems have been affected?
When you look at images in the media, it can feel like the whole of Mauritius is under oil. But the oil reached only 15 kilometres of the 350-kilometre shoreline, so it could have been much worse.
Unfortunately, there are a lot of environmentally sensitive areas in the region affected. The ship ran aground off Pointe d’Esny and just to the north of Blue Bay Marine Park. These sites are listed under the Ramsar Convention on Wetlands of International Importance as biodiversity hotspots. Ocean currents carried the oil northwards, so fortunately there’s none in the Blue Bay Marine Park, but the mangroves on the shoreline north of Pointe d’Esny have been covered. This will definitely have an impact, because mangroves are the nursery of the marine environment.
The Île aux Aigrettes, a small island near the wreck, has also been affected. The island is home to vulnerable pink pigeons (Nesoenas mayeri) and other native birds, and Telfair’s skink (Leiolopisma telfairii). The Mauritian Wildlife Foundation in Port Louis was already working to restore the island’s unique plants and remove invasive species. The oil didn’t go onto the island itself, but chemicals might have seeped into the corals and fumes from the spill could also have an impact.
Two rivers open into the bay where the oil spill is. The brackish water at the mouths of the rivers is an important ecosystem, and the oil has managed to go up parts of the rivers. The oil slick also floated above a large and rare area of seagrass, which is home to seahorses. Although the oil didn’t touch the seagrass, we fear that chemicals in the water could reach them.

Jacqueline Sauzier is president of the Mauritius Marine Conservation Society in Phoenix.Credit: Jacqueline Sauzier

Are there particular species affected?
It is not one species that could be at risk. It is the whole ecosystem, because of the dispersal of water-soluble chemicals in the water. Filter feeders, such as corals and crustaceans and molluscs, are probably the first to be impacted. We haven’t seen lots of animals dying, but we will need to monitor for signs.
Bad weather over the past two weeks has also forced the ship against the coral reef. That pushed a lot of sand and broken coral over the reef into the lagoon, creating a sand bar just inside the reef. That could change the currents in the lagoon and will have an impact on coral growth.
The social impacts are also a big concern for us. Fishing communities living in the region cannot fish anymore, because the fish that have been caught contain high levels of arsenic.
Something that is also concerning is that we don’t know the possible long-term effects. The oil is a new low-sulfur fuel oil that is being introduced to reduce air pollution. This is the first time that type of oil has spilled, so there have been no long-term studies on the impacts.
What steps are being taken now?
As soon as the ship grounded, people started monitoring the quality of the water. So we have this baseline from before the spill and we know the target that we have to reach for remediation of the water.
Oil-spill experts are formulating a plan to clean the shoreline properly. The impact on the mangroves could be worse if the cleaning is done badly. It could also push chemicals into the sand, which could be released in warm weather a year or two from now.
The front part of the vessel has been tugged away to be sunk along the shipping route. This was the least bad option. The rear is still on the reef. It has been cleaned of fuel, but rust and paint could still cause damage. It’s also falling apart, which can break the coral
Can future spills be avoided?
We were lucky this time that the spill was small and the boat grounded when it did. First, if it had happened in April, we wouldn’t have been able to go out, because we were locked down because of COVID. Second, the cane harvest started at the end of June. So if the wreck had happened earlier, we wouldn’t have had the cane trash readily available to contain the oil.
The spill has opened people’s eyes. Mauritius lies near a shipping highway. Around 2,500 large vessels pass close to Mauritius every month. It is difficult to describe how very big the Wakashio is. From the shore, it’s as if something foreign is sitting on your veranda. But this is the third stranding that we have had in ten years. Each time, there is uproar from the community, saying large vessels are much too close to the island.
It is not the first time, but I really hope it is the last.

doi: 10.1038/d41586-020-02446-7

This interview has been edited for length and clarity.

Latest on:

Biodiversity

An essential round-up of science news, opinion and analysis, delivered to your inbox every weekday. More

Assembling and monitoring the synthetic bacterial communities
We selected three bacterial strains, a Chryseobacterium sp., a Staphylococcus sp., and a Bacillus sp., from a large collection of soil isolates, and we screened them with fluorescence-independent flow cytometry as we did previously [17]. We measured at an acquisition speed of 14 μl min−1 for 2 min per sample and setting a threshold of 10,000 regarding the height signal of the front scatter (FSC-H). Differently than in our previous work, here we acquired all scattering profiles based on growth assays at 30 °C that was the temperature at which we performed all the experiments. In addition, we screened the growing cultures with a temporal resolution of 20, rather than 30, min. We recorded significant interactions among the strains by comparing their single and mixed growth profiles at 30 °C (Dataset 1). These interactions were mainly positive, similarly to what we found previously at other temperatures [17]. We performed all the related growth assays in biological triplicates. All flow cytometry data are available in .fcs format online (http://flowrepository.org) under the “FR-FCM-Z25Q” identifier. Henceforth, when referring to the experiments we use the term “population” to describe the cells of a given strain within a flask at a given assay and the term “community” to describe the total bacterial cells within a flask at a given assay.
Quantification of “background noise”
Our flow cytometry method for screening the mixed bacterial cultures has an accuracy of 97% for sample densities above 105 cells ml−1 [17]. However, at lower densities sampling errors and instrument inconsistencies become increasingly important because the signal-to-noise ratio drops. This can result in substantially different counts among identical samples and can artificially inflate the observed variability. Thus, it was essential to quantify this “noise” before performing the main experiments and subtract it from the observed variability when quantifying drift.
To that end, we made a series of separate experiments to quantify “noise”. In these experiments, we mixed overnight cultures of the three strains in all the seven possible combinations and in final cell densities ranging from 1.6 × 104 to 6.3 × 107 strain−1 ml−1 (corresponding to the expected range of cell densities in the main experiment, see below), and we measured repeated aliquots from the same flask to determine the coefficient of variation (CV—Fig. 1a). We treated the samples in exactly the same way as in the main experiments to include the effect of sampling errors in our calculations. We acquired in total of 99 triplicate measurements of 1–3 populations for a total of 148 observations (Fig. 1a, Dataset 2). We hypothesized that the level of “noise” should be inversely related to the cell density of the sample, because the signal-to-noise ratio decreases at low cell densities in the flow cytometer. Accordingly, we fit different functions for the dependency of CV to cell density (Supplementary Table S1). Finally, we calculated the 99.5% confidence intervals of the best-fitting function (i.e., Michaelis-Menten) using the confint function of the MASS package [18] in R [19], and we defined the false discovery rate based on the number of observations that were above the upper 99.5% confidence interval (Supplementary Fig. S1). Finally, we verified that the levels of noise determined in this study are similar to the variability recorded from technical replicate samples taken during our previous experiment where we used the same bacterial system and instrument with identical settings [17] (Supplementary Fig. S2).
Main experiments
To quantify drift, we monitored the changes in population densities across identical starting communities incubated under the same environmental conditions (Fig. 1b). To that end, we mixed the three strains in all seven possible combinations, i.e., three monocultures, three mixed cultures of two strains and a mixed culture of all three strains, and in three different starting total cell densities (5 × 104 cells ml−1, 105 cells ml−1, and 106 cells ml−1). To perform each growth assay, we first inoculated the respective strains from overnight pure cultures in a single flask containing 300 ml of Luria–Bertani medium (Sigma). To reach the desired starting total cell density, we estimated the cell density of the overnight pure cultures with flow cytometry [17] and we inoculated the respective volume. We then mixed the culture thoroughly and we immediately split the volume equally into three flasks. We next sampled 500 μl from each flask and we compared the variability in the bacterial populations across the three flasks to the expected “background noise” for the same cell density. In specific, we examined the CVs of the bacterial populations and their z-scores compared to the “background noise”, i.e., how many standard deviations an observed CV differs from the expected “background noise” CV at a given cell density. If the observed z-scores were larger than 2 (95% CI), we aborted the given experiment because it indicated that we introduced variability when we mixed and split the cultures and thus the starting cultures could not be considered identical; this happened in ~50% of the cases. If the observed z-scores were lower than 2, indicating that the recorded variability was not statistically different or was less than the expected variability based on the “background noise”, we proceeded with the experiment, incubating the three flasks in the same chamber (New Brunswick Innova 42R, Eppendorf) at 30 °C and with shaking at 80 rpm.
We recorded the starting densities (Dataset 3, z-scores between −6.68 and 1.17) and the densities every 20 min until the end of the fourth hour of incubation starting from the 60th minute. To detect and quantify drift, our main assumption was that any larger-than-expected deviations in the population densities of identical starting communities incubated under the same environmental conditions could only be because of drift. Thus, we compared the observed CVs to the expected CVs based on the “background noise” by deploying the z-score. We quantified drift using two different thresholds:
1.
The “upper threshold” that focused on excluding false-positive observations. In this quantification, we used a cutoff significance level of z  > 3 (99.5% CI) and we ignored the lowest 17.57% of positive observations (i.e., 15 observations, corresponding to the FDR level of the “background noise”) to minimize the detection of false positives.

2.
The “mean threshold” that focused on excluding false negatives and increasing detectability. In this quantification, we used a cutoff significance level of z  > 0, meaning that we scored any observation greater than the mean noise function as positive.

The “upper threshold” quantification probably overestimates drift by taking into account only the highest among the recorded CV values while the “mean threshold” quantification underestimates drift by taking into account some low CV values that are very close to the noise levels. Thus, the “upper threshold” and “mean threshold” quantifications do not represent the true levels of drift (which are hard to define whatsoever in the presence of noise) but they rather represent the upper and lower boundaries within which the true levels of drift lie.
Estimating potential growth variability due to temperature differences within the incubation chamber
To ensure that the recorded variability in the population counts was not due to slight differences in temperature within the incubation chamber, we estimated the potential variability that could have resulted if each strain grew within the extremes of the recorded temperatures in the chamber. For that, we first measured the temperature within each flask at each experimental time point, five times per flask, using a digital immersion thermometer with an accuracy of 0.1 °C. The temperature varied by 0.15 °C ± 0.08 °C on average and by 0.28 °C at maximum. We then calculated the growth rates of each of the strains under the recorded temperature extremes at each time point by interpolating from previously recorded growth rates [17]. We interpolated both with respect to time and with respect to temperature because the previous data were recorded at intervals of 30 min and 0.5 °C (the latter at a range of 25–42 °C). Finally, for each strain, we calculated how the CV in the hypothetical population densities would increase if the strains were constantly growing within the recorded temperature extremes for the duration of the experiment and if the CV was calculated from three observations (like in the real experiments) of the resulting population density distributions. We note that with this analysis, we probably overestimated the hypothetical increase in population densities because we used growth rates from mixed cultures that were generally higher than those in monocultures because of the positive interactions among the strains (Dataset 1).
Simulations
To simulate drift in complex bacterial communities, we used in silico communities with diversity and abundance distributions similar to nature [20] where drift acts with magnitude according to our experimental data. A conceptual flowchart of the simulations can be found in Supplementary Information (Supplementary Fig. S3). Each simulation involved a metacommunity of 100 communities that were connected with dispersal and that initially contained 2000 species each.
We simulated dispersal occurring in a unidirectional way within a closed system; individuals from community n disperse to the community (n + 1) and individuals from community 100 disperse back to community 1. The strength of dispersal equaled to the percentage of individuals that disperse to the respective community and it varied between 2 and 20%. Our aim in modeling dispersal in this way was to create a setting where habitat fragmentation was high and therefore drift’s importance is expected to be more pronounced [21, 22], and where there was no gain or loss of individuals from outside the metacommunity.
We simulated selection as differences in the growth rates among species within a community. The growth rates were distributed normally with a mean of 1 (resembling systems at their carrying capacity) and with a standard deviation between 0.071 and 0.167. Growth rates changed at every generation by being re-drawn from the same distribution in an effort to represent fluctuating habitats where a given species is not always favored or disfavored. Therefore, in our simulations, the standard deviation in the growth rates represents the strength of selection, because the higher it is the bigger are the differences in the growth rates in a community and the changes in the growth of a species from generation to generation. The distribution of the abundances in a community at time zero was log-normal (mean = 4, sd = 1.1) and the distribution of the abundances of a given species across all communities was normal with a standard deviation equal to the strength of selection.
The metacommunity grew for 1000 generations under given dispersal and selection conditions with drift, where drift changed the assigned growth rates at every generation according to a distribution based on the defined threshold from the experimental data (“upper” or “mean” threshold). In parallel, an initially identical metacommunity grew under the same dispersal and selection conditions but without drift, meaning that the assigned growth rates at every generation did not change further. More details and an example on how we modeled changes in growth rates due to drift are presented in the Supplementary Text in Supplementary Information.
For a given generation, we calculated the effect of drift by comparing a given community in the metacommunity that grew under drift to the same community in the metacommunity that grew without drift. In specific, we examined the β-diversity by means of the Bray–Curtis (BC) community similarity and the differences in species richness and in Pielou’s evenness among drift-impacted and drift-free communities, calculating the metacommunity-wise mean and standard deviation on all these properties. Moreover, we kept track of the extinct species at the end of each simulation and we mapped their initial relative abundances, but here we report the metacommunity-wise median because the distribution of the relative abundances is skewed (Supplementary Fig. S4). We ran simulations under 50 different scenarios resulting from five levels of selection strength over ten levels of dispersal rate. To estimate the effect of drift on Bray–Curtis similarity in metacommunities with increasing number of species, we ran the same simulation at the highest selection and lowest dispersal levels, at intermediate selection and dispersal and at the lowest selection and highest dispersal, but we changed the number of species; we ran the simulation three times in metacommunities of 500, 1000, 2000, 4000, 6000, 8000, and 10,000 species. We performed all simulations in R. All code is available on GitHub (https://github.com/sfodel/Drift).
Reported β-diversity in stochastically assembled communities in nature
To compare our simulation results with the results from environmental surveys regarding the β-diversity in stochastically assembled communities, we searched for related studies using the following two criteria: (1) the study cites the works of Stegen and colleagues [12, 13], where the term “undominated” community assembly is presented formally for microbial ecology, (2) the study reports data on the range of the observed β-diversity in terms of Bray–Curtis dissimilarity (or similarity) in stochastically assembled communities, or this range can be inferred from the data presented in that study. More

Amos W (2009) Heterozygosity and mutation rate: evidence for an interaction and its implications. Bioessays 32:82–90
Google Scholar 

Amos W, Kosanović D, Eriksson A (2015) Inter-allelic interactions play a major role in microsatellite evolution. Proc R Soc B 282:20152125
PubMed  Google Scholar 

Amos W (2016) Heterozygosity increases microsatellite mutation rate. Biol Lett 12:20150929
PubMed  PubMed Central  Google Scholar 

Azevedo L, Serrano C, Amorim A, Cooper DN (2015) Trans-species polymorphism in humans and the great apes is generally maintained by balancing selection that modulates the host immune response. Hum Genomics 9:21
PubMed  PubMed Central  Google Scholar 

Beye M, Hasselmann M, Fondrk MK, Page RE, Omholt SW (2003) The gene csd is the primary signal for sexual development in the honeybee and encodes an SR-type protein. Cell 114:419–429
CAS  PubMed  Google Scholar 

Beye M, Seelmann C, Gempe T, Hasselmann M, Vekemans X, Fondrk MK et al. (2013) Gradual molecular evolution of a sex determination switch through incomplete penetrance of femaleness. Curr Biol 23:2559–2564
CAS  PubMed  Google Scholar 

Bezabih G, Adgaba N, Hepburn HR, Pirk CWW (2014) The territorial invasion of Apis florea in Africa. Afr Entomol 22:888–890
Google Scholar 

Biewer M, Lechner S, Hasselmann M (2016) Similar but not the same: insights into the evolutionary history of paralogous sex-determining genes of the dwarf honey bee Apis florea. Heredity 116:12–22
CAS  PubMed  Google Scholar 

Chapuis MP, Plantamp C, Streiff R, Blondin L, Piou C (2015) Microsatellite evolutionary rate and pattern in Schistocerca gregaria inferred from direct observation of germline mutations. Mol Ecol 24:6107–6119
CAS  PubMed  Google Scholar 

Chookajorn T, Kachroo A, Ripoll DR, Clark AG, Nasrallah JB (2004) Specificity determinants and diversification of the Brassica self-incompatibility pollen ligand. Proc Natl Acad Sci USA 101:911–917
CAS  PubMed  Google Scholar 

Cho S, Huang ZY, Green DR, Smith DR, Zhang J (2006) Evolution of the complementary sex-determination gene of honey bees: balancing selection and trans-species polymorphisms. Genome Res 16:1366–1375
CAS  PubMed  PubMed Central  Google Scholar 

Clarke BC (1979) The evolution of genetic diversity. Proc R Soc B 205:453–474
CAS  Google Scholar 

Ding G, Xu H, Oldroyd BP, Gloag RS (2017) Extreme polyandry aids the establishment of invasive populations of a social insect. Heredity 119:381–387
CAS  PubMed  PubMed Central  Google Scholar 

Ellegren H (2004) Microsatellites: simple sequences with complex evolution. Nat Rev Genet 5:435–445
CAS  PubMed  Google Scholar 

Fay JC, Wu CI (2000) Hitchhiking under positive Darwinian selection. Genetics 155:1405–1413
CAS  PubMed  PubMed Central  Google Scholar 

Fu YX (1997) Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147:915–925
CAS  PubMed  PubMed Central  Google Scholar 

Glémin S, Gaude T, Guillemin ML, Lourmas M, Olivieri I, Mignot A (2005) Balancing selection in the wild: testing population genetics theory of self-incompatibility in the rare species Brassica insularis. Genetics 171:279–289
PubMed  PubMed Central  Google Scholar 

Gervais CE, Castric V, Ressayre A, Billiard S (2011) Origin and diversification dynamics of self-incompatibility haplotypes. Genetics 188:625–636
CAS  PubMed  PubMed Central  Google Scholar 

Gibbs MJ, Armstrong JS, Gibbs AJ (2000) Sister-scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics 16:573–582
CAS  PubMed  Google Scholar 

Gloag R, Christie JR, Ding G, Stephens R, Buchmann G, Oldroyd BP (2019) Workers’ sons rescue genetic diversity at the sex locus in an invasive honey bee population. Mol Ecol 28:1585–1592
PubMed  Google Scholar 

Gloag R, Ding G, Christie JR, Buchmann G, Beekman M, Oldroyd BP (2017) An invasive social insect overcomes genetic load at the sex locus. Nat Ecol Evol 1:1–6
Google Scholar 

Hasselmann M, Beye M (2004) Signatures of selection among sex-determining alleles of the honey bee. Proc Natl Acad Sci USA 101:4888–4893
CAS  PubMed  Google Scholar 

Hasselmann M, Vekemans X, Pflugfelder J, Koeniger N, Koeniger G, Tingek S et al. (2008) Evidence for convergent nucleotide evolution and high allelic turnover rates at the complementary sex determiner gene of Western and Asian honeybees. Mol Biol Evol 25:696–708
CAS  PubMed  Google Scholar 

Hedrick PW (2005) A standardized genetic differentiation measure. Evolution 59:1633–1638
CAS  PubMed  Google Scholar 

Koch V, Nissen I, Schmitt BD, Beye M (2014) Independent evolutionary origin of fem paralogous genes and complementary sex determination in hymenopteran insects. PLoS ONE 9:e91883
PubMed  PubMed Central  Google Scholar 

Koetz AH (2013) Ecology, behaviour and control of Apis cerana with a focus on relevance to the Australian incursion. Insects 4:558–592
PubMed  PubMed Central  Google Scholar 

Kosakovsky Pond SL, Posada D, Gravenor MB, Woelk CH, Frost SD (2006) Automated phylogenetic detection of recombination using a genetic algorithm. Mol Biol Evol 23:1891–1901
PubMed  Google Scholar 

Kumar S, Stecher G, Li M, Knyaz C, Tamura K (2018) MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol 35:1547–1549
CAS  PubMed  PubMed Central  Google Scholar 

Kurbalija Novičić Z, Sayadi A, Jelić M, Arnqvist G (2020) Negative frequency dependent selection contributes to the maintenance of a global polymorphism in mitochondrial DNA. BMC Evol Biol 20:20. https://doi.org/10.1186/s12862-020-1581-2
CAS  Article  PubMed  PubMed Central  Google Scholar 

Lechner S, Ferretti L, Schöning C, Kinuthia W, Willemsen D, Hasselmann M (2014) Nucleotide variability at its limit? Insights into the number and evolutionary dynamics of the sex-determining specificities of the honey bee Apis mellifera. Mol Biol Evol 31:272–287
CAS  PubMed  Google Scholar 

Levin BR (1988) Frequency dependent selection in bacterial populations. Philos Trans R Soc Lond B 319:459–472
CAS  Google Scholar 

Li Y-C, Korol AB, Fahima T, Nevo E (2004) Microsatellites within genes: structure, function, and evolution. Mol Biol Evol 21:991–1007
CAS  PubMed  Google Scholar 

Lynch M (2015) Feedforward loop for diversity. Nature 523:414–416
CAS  PubMed  PubMed Central  Google Scholar 

Mackensen O (1951) Viability and sex determination in the honeybee (Apis mellifera L.). Genetics 36:500–509
CAS  PubMed  PubMed Central  Google Scholar 

Martin DP, Murrell B, Golden M, Khoosal A, Muhire B (2015) RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol 1:vev003
PubMed  PubMed Central  Google Scholar 

Martin DP, Posada D, Crandall KA, Williamson C (2005) A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. AIDS Res Hum Retroviruses 21:98–102
CAS  PubMed  Google Scholar 

May G, Matzke E (1995) Recombination and variation at the A mating-type of Coprinus cinereus. Mol Biol Evol 12:794–802
CAS  Google Scholar 

Moritz RFA, Haddad N, Bataieneh A, Shalmon B, Hefetz A (2010) Invasion of the dwarf honeybee Apis florea into the near East. Biol Invasions 12:1093–1099
Google Scholar 

Moritz RFA, Härtel S, Neumann P (2005) Global invasions of the western honeybee (Apis mellifera) and the consequences for biodiversity. Écoscience 12:289–301
Google Scholar 

Muirhead CA (2001) Consequences of population structure on genes under balancing selection. Evolution 55:1532–1541
CAS  PubMed  Google Scholar 

Muirhead CA, Glass NL, Slatkin M (2002) Multilocus self-recognition systems in fungi as a cause of trans-species polymorphism. Genetics 161:633–641
CAS  PubMed  PubMed Central  Google Scholar 

Nasrallah JB (1997) Evolution of the Brassica self-incompatibility locus: a look into S-locus gene polymorphisms. Proc Natl Acad Sci USA 94:9516–9519
CAS  PubMed  Google Scholar 

Nasrallah JB, Kao TH, Chen CH, Goldberg ML, Nasrallah ME (1987) Amino-acid sequence of glycoproteins encoded by three alleles of the S locus of Brassica oleracea. Nature 326:617–619
CAS  Google Scholar 

Oldroyd BP, Wongsiri S (2006) Asian honey bees. Biology, conservation and human interactions. Harvard University Press, Cambridge
Google Scholar 

Posada D, Crandall KA (2001) Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc Natl Acad Sci USA 98:13757–13762
CAS  PubMed  Google Scholar 

Radloff SE, Hepburn C, Hepburn HR, Fuchs S, Hadisoesilo S, Tan K et al. (2010) Population structure and classification of Apis cerana. Apidologie 41:589–601
Google Scholar 

Raper JR (1966) Sexuality in fungi. (Book reviews: genetics of sexuality in higher fungi). Science 154:758
Google Scholar 

Richman AD, Kohn JR (1999) Self-incompatibility alleles from Physalis: implications for historical inference from balanced genetic polymorphisms. Proc Natl Acad Sci USA 96:168–172
CAS  PubMed  Google Scholar 

Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE et al. (2017) DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol Biol Evol 34:3299–3302
CAS  PubMed  Google Scholar 

Schierup MH, Bechsgaard JS, Christiansen FB (2008) Selection at work in self-incompatible Arabidopsis lyrata. II. Spatial distribution of S haplotypes in Iceland. Genetics 180:1051–1059
PubMed  PubMed Central  Google Scholar 

Schierup MH, Vekemans X, Charlesworth D (2000) The effect of subdivision on variation at multi-allelic loci under balancing selection. Genet Res 76:51–62
CAS  PubMed  Google Scholar 

Smith DR (2011) Asian Honeybees and Mitochondrial DNA. In: Hepburn HR, Radloff SE (eds) Honeybees of Asia. Springer-Verlag Berlin, Heidelberg, p 69–93
Google Scholar 

Solignac M, Vautrin D, Loiseau A, Mougel F, Baudry E, Estoup A et al. (2003) Five hundred and fifty microsatellite markers for the study of the honeybee (Apis mellifera L.) genome. Mol Ecol Notes 3:307–311
CAS  Google Scholar 

Tajima F (1989) The effect of change in population size on DNA polymorphism. Genetics 123:597–601
CAS  PubMed  PubMed Central  Google Scholar 

Takahashi J, Shimizu S, Koyama S, Kimura K, Shimizu I, Yoshida T (2009) Variable microsatellite loci isolated from the Asian honeybee, Apis cerana (Hymenoptera; Apidae). Mol Ecol Resour 9:819–821
CAS  PubMed  Google Scholar 

Takahata N, Nei M (1990) Allelic genealogy under overdominant and frequency-dependent selection and polymorphism of major histocompatibility complex loci. Genetics 124:967–978
CAS  PubMed  PubMed Central  Google Scholar 

Wallberg A, Han F, Wellhagen G, Dahle B, Kawata M, Haddad N et al. (2014) A worldwide survey of genome sequence variation provides insight into the evolutionary history of the honeybee Apis mellifera. Nat Genet 46:1081–1088
CAS  PubMed  Google Scholar 

Walsh PS, Metzger DA, Higuchi R (1991) Chelex 100 as amedium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques 10:506–513
CAS  PubMed  Google Scholar 

Wang J (2015) Does GST underestimate genetic differentiation from marker data? Mol Ecol 24:3546–3558
CAS  PubMed  Google Scholar 

Weaver S, Shank SD, Spielman SJ, Li M, Muse SV, Kosakovsky Pond SL (2018) Datamonkey 2.0: a modern web application for characterizing selective and other evolutionary processes. Mol Biol Evol 35:773–777
CAS  PubMed  PubMed Central  Google Scholar 

Woyke J (1963) What happens to diploid drone larvae in honeybee colony? J Apic Res 2:73–75
Google Scholar 

Wright S (1939) The distribution of self-sterility alleles in populations. Genetics 24:538–552
CAS  PubMed  PubMed Central  Google Scholar 

Wu Q, Han TS, Chen X, Chen JF, Zou YP, Li ZW et al. (2017) Long-term balancing selection contributes to adaptation in Arabidopsis and its relatives. Genome Biol 18:217
PubMed  PubMed Central  Google Scholar 

Yang S, Wang L, Huang J, Zang X, Yuan Y, Chen JQ et al. (2015) Parent-progeny sequencing indicates higher mutation rates in heterozygotes. Nature 523:463–467
CAS  PubMed  Google Scholar 

Yokoyama S, Nei M (1979) Population dynamics of sex-determining alleles in honey bees and self-incompatibility alleles in plants. Genetics 91:609–626
CAS  PubMed  PubMed Central  Google Scholar 

Zareba J, Blazej P, Laszkiewicz A, Sniezewski L, Majkowski M, Janik S et al. (2017) Uneven distribution of complementary sex determiner (csd) alleles in Apis mellifera population. Sci Rep. 7:2317
PubMed  PubMed Central  Google Scholar  More

For predators, climate change-induced shifts in prey numbers, behaviours and spatial or temporal locations can be a major threat to food security. For predators that hoard prey to ensure survival through harsh winters, climate variation can have a doubled effect — influencing both food capture and store stability. Although northern latitude autumn and winter temperatures have increased dramatically in recent years, the effects of climate on foraging and storing throughout winter remain understudied.

Credit: Szymon Bartosz / Alamy Stock Photo

Giulia Masoero at the University of Turku, Finland, and colleagues analysed the impact of climate on Eurasian pygmy owl (Glaucidium passerinum) food-hoarding behaviour across 16 years. They found increased freeze–thaw frequency, lower winter precipitation and deeper snow cover were linked to greater hoard consumption. Higher autumn precipitation and an early hoarding start led to food rot, which reduced female owl recapture (indicating death or emigration).

Although owls delayed hoarding in autumns with fewer freeze–thaw events, suggesting some potential for climate change adaptation, the study indicates that altered climates can decrease owl overwinter survival, which may in turn have vast impacts on the boreal food web.

Author information

Affiliations

Nature Climate Change
Tegan Armarego-Marriott

Authors
Tegan Armarego-Marriott

Corresponding author
Correspondence to Tegan Armarego-Marriott.

Rights and permissions

About this article

Cite this article
Armarego-Marriott, T. Owls’ hoards rot. Nat. Clim. Chang. (2020). https://doi.org/10.1038/s41558-020-0903-0
Download citation More

1.
Hastorf, C. A. The Social Archaeology of Food: Thinking About Eating From Prehistory to the Present (Cambridge University Press, Cambridge, 2017).
Google Scholar 
2.
Lévi-Strauss, C. The culinary triangle. Partisan Review33, 586–595 (1966).
Google Scholar 

3.
Wrangham, R. W. Catching Fire: How Cooking Made Us Human (Basic Books Inc, New York, 2009).
Google Scholar 

4.
Roffet-Salque, M. et al. From the inside out: upscaling organic residue analyses of archaeological ceramics. J. Archaeol. Sci.16, 627–640 (2017).
Google Scholar 

5.
Appadurai, A. Gastro-Politics in Hindu South Asia. Am. Ethnol.8, 494–511 (1981).
Google Scholar 

6.
Douglas, M. Deciphering a meal. In Food and Culture: A Reader (eds Counihan, C. & Van Esterik, P.) 44–53 (Routledge, New York, 2008).
Google Scholar 

7.
Twiss, K. The Archaeology of Food and Identity (Center for Archaeological Investigations, Southern Illinois University at Carbondale, Carbondale, 2007).
Google Scholar 

8.
Rozin, P. Psychobiological perspectives on food preferences and avoidances. In Food and Evolution (eds Harris, M. & Ross, E. B.) 181–205 (Temple University Press, Philadelphia, 1987).
Google Scholar 

9.
Condamin, J., Formenti, F., Metais, M. O., Michel, M. & Blond, P. The application of gas chromatography to the tracing of oil in ancient amphorae. Archaeometry18, 195–201 (1976).
Google Scholar 

10.
Passi, S., Rothschild-Boros, M. C., Fasella, P., Nazzaro-Porro, M. & Whitehouse, D. An application of high performance liquid chromatography to analysis of lipids in archaeological samples. J. Lipid Res.22, 778–784 (1981).
PubMed  CAS  Google Scholar 

11.
Hastorf, C. A. & DeNiro, M. J. Reconstruction of prehistoric plant production and cooking practices by a new isotopic method. Nature315, 489–491 (1985).
ADS  CAS  Google Scholar 

12.
Patrick, M., de Koning, A. J. & Smith, A. B. Gas liquid chromatographic analysis of fatty acids in food residues from ceramics found in the Southwestern Cape, South Africa. Archaeometry27, 231–236 (1985).
CAS  Google Scholar 

13.
Evershed, R. P., Heron, C. & Goad, L. J. Epicuticular wax components preserved in potsherds as chemical indicators of leafy vegetables in ancient diets. Antiquity65, 540–544 (1991).
Google Scholar 

14.
Morton, J. D. & Schwarcz, H. P. Palaeodietary implications from stable isotopic analysis of residues on prehistoric Ontario ceramics. J. Archaeol. Sci.31, 503–517 (2004).
Google Scholar 

15.
Heron, C., Evershed, R. P. & Goad, L. J. Effects of migration of soil lipids on organic residues associated with buried potsherds. J. Archaeol. Sci.18, 641–659 (1991).
Google Scholar 

16.
Dudd, S. N., Evershed, R. P. & Gibson, A. M. Evidence for varying patterns of exploitation of animal products in different prehistoric pottery traditions based on lipids preserved in surface and absorbed residues. J. Archaeol. Sci.26, 1473–1482 (1999).
Google Scholar 

17.
Hart, J. P., Lovis, W. A., Schulenberg, J. K. & Urquhart, G. R. Paleodietary implications from stable carbon isotope analysis of experimental cooking residues. J. Archaeol. Sci.34, 804–813 (2007).
Google Scholar 

18.
Hansel, F. A., Copley, M. S., Madureira, L. A. S. & Evershed, R. P. Thermally produced ω-(o-alkylphenyl)alkanoic acids provide evidence for the processing of marine products in archaeological pottery vessels. Tetrahedron Lett.45, 2999–3002 (2004).
CAS  Google Scholar 

19.
Hansel, F. A. & Evershed, R. P. Formation of dihydroxy acids from Z-monounsaturated alkenoic acids and their use as biomarkers for the processing of marine commodities in archaeological pottery vessels. Tetrahedron Lett.50, 5562–5564 (2009).
CAS  Google Scholar 

20.
Evershed, R. P. Organic residue analysis in archaeology: the archaeological biomarker revolution. Archaeometry50, 895–924 (2008).
CAS  Google Scholar 

21.
Copley, M. S., Hansel, F. A., Sadr, K. & Evershed, R. P. Organic residue evidence for the processing of marine animal products in pottery vessels from the pre-colonial archaeological site of Kasteelberg D east, South Africa : research article. S. Afr. J. Sci.100, 279–283 (2004).
CAS  Google Scholar 

22.
DeNiro, M. J. & Epstein, S. Influence of diet on the distribution of carbon isotopes in animals. Geochim. Cosmochim. Acta42, 495–506 (1978).
ADS  CAS  Google Scholar 

23.
DeNiro, M. J. & Epstein, S. Influence of diet on the distribution of nitrogen isotopes in animals. Geochim. Cosmochim. Acta45, 341–351 (1981).
ADS  CAS  Google Scholar 

24.
Heron, C. & Craig, O. E. Aquatic resources in foodcrusts: identification and implication. Radiocarbon57, 707–719 (2015).
Google Scholar 

25.
Cramp, L. J. E. et al. Immediate replacement of fishing with dairying by the earliest farmers of the northeast Atlantic archipelagos. Proc. R. Soc. B281, 20132372 (2014).
PubMed  Google Scholar 

26.
Evershed, R. P. Biomolecular archaeology and lipids. World Archaeol.25, 74–93 (1993).
PubMed  CAS  Google Scholar 

27.
Evershed, R. P., Heron, C., Charters, S. & Goad, L. J. The survival of food residues: new methods of analysis, interpretation and application. Proc. Br. Acad.77, 187–208 (1992).
Google Scholar 

28.
Copley, M. S. et al. Direct chemical evidence for widespread dairying in prehistoric Britain. PNAS100, 1524–1529 (2003).
ADS  PubMed  CAS  Google Scholar 

29.
Evershed, R. P., Arnot, K. I., Collister, J., Eglinton, G. & Charters, S. Application of isotope ratio monitoring gas chromatography–mass spectrometry to the analysis of organic residues of archaeological origin. Analyst119, 909–914 (1994).
ADS  CAS  Google Scholar 

30.
Evershed, R. P. et al. Lipids as carriers of anthropogenic signals from prehistory. Philos. Trans. R. Soc. Lond. B354, 19–31 (1999).
CAS  Google Scholar 

31.
Dunne, J., Mercuri, A. M., Evershed, R. P., Bruni, S. & di Lernia, S. Earliest direct evidence of plant processing in prehistoric Saharan pottery. Nat. Plants3, 1–6 (2016).
Google Scholar 

32.
Roffet-Salque, M. et al. Widespread exploitation of the honeybee by early Neolithic farmers. Nature527, 226–230 (2015).
ADS  PubMed  CAS  Google Scholar 

33.
Hammann, S. & Cramp, L. J. E. Towards the detection of dietary cereal processing through absorbed lipid biomarkers in archaeological pottery. J. Archaeol. Sci.93, 74–81 (2018).
CAS  Google Scholar 

34.
Evershed, R. P., Copley, M. S., Dickson, L. & Hansel, F. A. Experimental evidence for the processing of marine animal products and other commodities containing polyunsaturated fatty acids in pottery vessels. Archaeometry50, 101–113 (2008).
CAS  Google Scholar 

35.
Craig, O. E. et al. Molecular and isotopic demonstration of the processing of aquatic products in northern European prehistoric pottery. Archaeometry49, 135–152 (2007).
CAS  Google Scholar 

36.
Craig, O. E. et al. Earliest evidence for the use of pottery. Nature496, 351–354 (2013).
ADS  PubMed  CAS  Google Scholar 

37.
Robson, H. K. et al. Diet, cuisine and consumption practices of the first farmers in the southeastern Baltic. Archaeol. Anthropol. Sci.11, 4011–4024 (2019).
PubMed  PubMed Central  Google Scholar 

38.
Oras, E. et al. The adoption of pottery by north-east European hunter-gatherers: evidence from lipid residue analysis. J. Archaeol. Sci.78, 112–119 (2017).
CAS  Google Scholar 

39.
Admiraal, M., Lucquin, A., von Tersch, M., Jordan, P. D. & Craig, O. E. Investigating the function of prehistoric stone bowls and griddle stones in the Aleutian Islands by lipid residue analysis. Quatern. Res.91, 1003–1015 (2019).
ADS  Google Scholar 

40.
Spangenberg, J. E., Jacomet, S. & Schibler, J. Chemical analyses of organic residues in archaeological pottery from Arbon Bleiche 3, Switzerland—evidence for dairying in the late Neolithic. J. Archaeol. Sci.33, 1–13 (2006).
Google Scholar 

41.
Guiry, E. et al. Differentiating salmonid migratory ecotypes through stable isotope analysis of collagen: archaeological and ecological applications. PLoS ONE15, e0232180 (2020).
PubMed  PubMed Central  CAS  Google Scholar 

42.
Mukherjee, A. J., Gibson, A. M. & Evershed, R. P. Trends in pig product processing at British Neolithic Grooved Ware sites traced through organic residues in potsherds. J. Archaeol. Sci.35, 2059–2073 (2008).
Google Scholar 

43.
Shoda, S., Lucquin, A., Ahn, J., Hwang, C. & Craig, O. E. Pottery use by early Holocene hunter-gatherers of the Korean peninsula closely linked with the exploitation of marine resources. Quatern. Sci. Rev.170, 164–173 (2017).
ADS  Google Scholar 

44.
Evershed, R. P. Experimental approaches to the interpretation of absorbed organic residues in archaeological ceramics. World Archaeol.40, 26–47 (2008).
Google Scholar 

45.
Grace, R. Review article use-wear analysis: the state of the art. Archaeometry38, 209–229 (1996).
Google Scholar 

46.
Oudemans, T. F. Applying organic residue analysis in ceramic studies in archaeology: a functional approach. Leiden J. Pottery Stud.23, 5–20 (2007).
Google Scholar 

47.
Skibo, J. M. Pottery use-alteration analysis. In Use-Wear and Residue Analysis in Archaeology (eds Marreiros, J. M. et al.) 189–198 (Springer International Publishing, New York, 2015).
Google Scholar 

48.
Marino, B. D. & Deniro, M. J. Isotopic analysis of archaeobotanicals to reconstruct past climates: effects of activities associated with food preparation on carbon, hydrogen and oxygen isotope ratios of plant cellulose. J. Archaeol. Sci.14, 537–548 (1987).
Google Scholar 

49.
Lovis, W. A., Urquhart, G. R., Raviele, M. E. & Hart, J. P. Hardwood ash nixtamalization may lead to false negatives for the presence of maize by depleting bulk δ13C in carbonized residues. J. Archaeol. Sci.38, 2726–2730 (2011).
Google Scholar 

50.
Hart, J. P., Urquhart, G. R., Feranec, R. S. & Lovis, W. A. Non-linear relationship between bulk δ13C and percent maize in carbonized cooking residues and the potential of false-negatives in detecting maize. J. Archaeol. Sci.36, 2206–2212 (2009).
Google Scholar 

51.
Warinner, C. & Tuross, N. Alkaline cooking and stable isotope tissue-diet spacing in swine: archaeological implications. J. Archaeol. Sci.36, 1690–1697 (2009).
Google Scholar 

52.
Warinner, C. Life and Death at Teposcolula Yucundaa: Mortuary, Archaeogenetic, and Isotopic Investigations of the Early Colonial Period in Mexico (Harvard University, New York, 2010).
Google Scholar 

53.
Renson, V. et al. Origin and diet of inhabitants of the Pacific Coast of Southern Mexico during the Classic Period—Sr, C and N isotopes. J. Archaeol. Sci. Rep.27, 101981 (2019).
Google Scholar 

54.
Szpak, P., Millaire, J.-F., White, C. D. & Longstaffe, F. J. Influence of seabird guano and camelid dung fertilization on the nitrogen isotopic composition of field-grown maize (Zea mays). J. Archaeol. Sci.39, 3721–3740 (2012).
CAS  Google Scholar 

55.
Wang, C. et al. Aridity threshold in controlling ecosystem nitrogen cycling in arid and semi-arid grasslands. Nat. Commun.5, 1–8 (2014).
ADS  Google Scholar 

56.
Yoneyama, T., Ito, O. & Engelaar, W. M. H. G. Uptake, metabolism and distribution of nitrogen in crop plants traced by enriched and natural 15N: Progress over the last 30 years. Phytochem. Rev.2, 121–132 (2003).
CAS  Google Scholar 

57.
Bocherens, H. & Drucker, D. Trophic level isotopic enrichment of carbon and nitrogen in bone collagen: case studies from recent and ancient terrestrial ecosystems. Int. J. Osteoarchaeol.13, 46–53 (2003).
Google Scholar 

58.
Minagawa, M. & Wada, E. Stepwise enrichment of 15N along food chains: further evidence and the relation between d15N and animal age. Geochim. Cosmochim. Acta48, 1135–1140 (1984).
ADS  CAS  Google Scholar 

59.
DeNiro, M. J. & Hastorf, C. A. Alteration of 15N14N and 13C12C ratios of plant matter during the initial stages of diagenesis: studies utilizing archaeological specimens from Peru. Geochim. Cosmochim. Acta49, 97–115 (1985).
ADS  CAS  Google Scholar 

60.
Fraser, R. A. et al. Assessing natural variation and the effects of charring, burial and pre-treatment on the stable carbon and nitrogen isotope values of archaeobotanical cereals and pulses. J. Archaeol. Sci.40, 4754–4766 (2013).
CAS  Google Scholar 

61.
Phillips, D. L., & Gregg, J. W. Source partitioning using stable isotopes: coping with too many sources. Oecologia136, 261–269 (2003).
ADS  PubMed  Google Scholar 

62.
Drieu, L. et al. Influence of porosity on lipid preservation in the wall of archaeological pottery. Archaeometry61, 1081–1096 (2019).
CAS  Google Scholar 

63.
Charters, S. Chemical Analysis of Absorbed Lipids and Laboratory Simulation Experiments to Interpret Archaeological Pottery Vessel Contents and Use (University of Bristol, Bristol, 1996).
Google Scholar 

64.
Bogaard, A., Heaton, T. H. E., Poulton, P. & Merbach, I. The impact of manuring on nitrogen isotope ratios in cereals: archaeological implications for reconstruction of diet and crop management practices. J. Archaeol. Sci.34, 335–343 (2007).
Google Scholar 

65.
Styring, A. K. et al. The effect of charring and burial on the biochemical composition of cereal grains: investigating the integrity of archaeological plant material. J. Archaeol. Sci.40, 4767–4779 (2013).
CAS  Google Scholar 

66.
Szpak, P. & Chiou, K. L. A comparison of nitrogen isotope compositions of charred and desiccated botanical remains from northern Peru. Veget Hist Archaeobot https://doi.org/10.1007/s00334-019-00761-2 (2019).
Article  Google Scholar 

67.
Casanova, E. et al. Accurate compound-specific 14 C dating of archaeological pottery vessels. Nature580, 506–510 (2020).
ADS  PubMed  CAS  Google Scholar 

68.
Boivin, N. L. et al. Ecological consequences of human niche construction: examining long-term anthropogenic shaping of global species distributions. Proc. Natl. Acad. Sci. USA113, 6388–6396 (2016).
PubMed  CAS  Google Scholar 

69.
Ceballos, G. et al. Accelerated modern human–induced species losses: Entering the sixth mass extinction. Sci. Adv.1, e1400253 (2015).
ADS  PubMed  PubMed Central  Google Scholar 

70.
Foley, S. F. et al. The Palaeoanthropocene: the beginnings of anthropogenic environmental change. Anthropocene3, 83–88 (2013).
Google Scholar 

71.
Tilman, D. & Lehman, C. Human-caused environmental change: Impacts on plant diversity and evolution. PNAS98, 5433–5440 (2001).
ADS  PubMed  CAS  Google Scholar 

72.
Wilk, R. Consumption, human needs, and global environmental change. Glob. Environ. Change12, 5–13 (2002).
Google Scholar 

73.
Correa-Ascencio, M. & Evershed, R. P. High throughput screening of organic residues in archaeological potsherds using direct acidified methanol extraction. Anal. Methods6, 1330–1340 (2014).
CAS  Google Scholar  More

1.
Turunen, J., Tomppo, E., Tolonen, K. & Reinikainen, A. Estimating carbon accumulation rates of undrained mires in Finland–application to boreal and subarctic regions. Holocene 12, 69–80 (2002).
ADS  Google Scholar 
2.
Gorham, E. Northern Peatlands: role in the carbon cycle and probable responses to climatic warming. Ecol. Appl. 1, 182–195 (1991).
PubMed  Google Scholar 

3.
Loisel, J. et al. A database and synthesis of northern peatland soil properties and Holocene carbon and nitrogen accumulation. Holocene 24, 1028–1042 (2014).
ADS  Google Scholar 

4.
Yu, Z., Loisel, J., Brosseau, D. P., Beilman, D. W. & Hunt, S. J. Global peatland dynamics since the Last Glacial Maximum. Geophys. Res. Lett. 37, L13402 (2010).
ADS  Google Scholar 

5.
Batjes, N. H. Total carbon and nitrogen in the soils of the world. Eur. J. Soil Sci. 65, 10–21 (2014).
CAS  Google Scholar 

6.
Dorrepaal, E. et al. Carbon respiration from subsurface peat accelerated by climate warming in the subarctic. Nature 460, 616–619 (2009).
ADS  CAS  Google Scholar 

7.
Nijp, J. J. et al. Rain events decrease boreal peatland net CO2 uptake through reduced light availability. Glob. Change Biol. 21, 2309–2320 (2015).
ADS  Google Scholar 

8.
Petrescu, A. M. R. et al. The uncertain climate footprint of wetlands under human pressure. Proc. Natl Acad. Sci. USA. 112, 4594–4599 (2015).
ADS  PubMed  Google Scholar 

9.
Wu, J. & Roulet, N. T. Climate change reduces the capacity of northern peatlands to absorb the atmospheric carbon dioxide: the different responses of bogs and fens. Glob. Biogeochem. Cycles 28, 1005–1024 (2014).
ADS  CAS  Google Scholar 

10.
Wang, X. et al. Soil respiration under climate warming: differential response of heterotrophic and autotrophic respiration. Glob. Change Biol. 20, 3229–3237 (2014).
ADS  Google Scholar 

11.
Mäkiranta, P. et al. Indirect regulation of heterotrophic peat soil respiration by water level via microbial community structure and temperature sensitivity. Soil Biol. Biochem. 41, 695–703 (2009).
Google Scholar 

12.
Bond-Lamberty, B., Wang, C. & Gower, S. T. A global relationship between the heterotrophic and autotrophic components of soil respiration? Glob. Change Biol. 10, 1756–1766 (2004).
ADS  Google Scholar 

13.
Järveoja, J., Nilsson, M. B., Gažovič, M., Crill, P. M. & Peichl, M. Partitioning of the net CO2 exchange using an automated chamber system reveals plant phenology as key control of production and respiration fluxes in a boreal peatland. Glob. Change Biol. 24, 3436–3451 (2018).
ADS  Google Scholar 

14.
Lloyd, J. & Taylor, J. A. On the temperature dependence of soil respiration. Funct. Ecol. 8, 315 (1994).
Google Scholar 

15.
van’t Hoff, J. H. Lectures on theoretical and physical chemistry: chemical dynamics. Part I (Edward Arnold, London, 1898).

16.
Arrhenius, S. Uber die Reaktionsgeschwindigkeit bei der Inversion von Rohrzucker durch Sauren. Z. f.ür. Phys. Chem. 4, 226–248 (1889).
Google Scholar 

17.
Reichstein, M. et al. On the separation of net ecosystem exchange into assimilation and ecosystem respiration: review and improved algorithm. Glob. Change Biol. 11, 1424–1439 (2005).
ADS  Google Scholar 

18.
Lasslop, G. et al. Separation of net ecosystem exchange into assimilation and respiration using a light response curve approach: critical issues and global evaluation. Glob. Change Biol. 16, 187–208 (2010).
ADS  Google Scholar 

19.
Lasslop, G. et al. On the choice of the driving temperature for eddy-covariance carbon dioxide flux partitioning. Biogeosciences 9, 5243–5259 (2012).
ADS  CAS  Google Scholar 

20.
Wohlfahrt, G. & Galvagno, M. Revisiting the choice of the driving temperature for eddy covariance CO2 flux partitioning. Agric. Meteorol. 237–238, 135–142 (2017).
Google Scholar 

21.
Phillips, S. C. et al. Interannual, seasonal, and diel variation in soil respiration relative to ecosystem respiration at a wetland to upland slope at Harvard Forest. J. Geophys. Res. Biogeosciences 115, G02019 (2010).

22.
Savage, K., Davidson, E. A. & Tang, J. Diel patterns of autotrophic and heterotrophic respiration among phenological stages. Glob. Change Biol. 19, 1151–1159 (2013).
ADS  CAS  Google Scholar 

23.
Thorne, R., Khomik, M., Hayman, E. & Arain, A. Response of soil CO2 efflux to shelterwood harvesting in a mature temperate pine forest. Forests 11, 304 (2020).
Google Scholar 

24.
Carbone, M. S., Winston, G. C. & Trumbore, S. E. Soil respiration in perennial grass and shrub ecosystems: linking environmental controls with plant and microbial sources on seasonal and diel timescales. J. Geophys. Res. Biogeosciences 113, G02022 (2008).

25.
Keane, J. Ben & Ineson, P. Technical note: differences in the diurnal pattern of soil respiration under adjacent Miscanthus × giganteus and barley crops reveal potential flaws in accepted sampling strategies. Biogeosciences 14, 1181–1187 (2017).
ADS  CAS  Google Scholar 

26.
Wutzler, T. et al. Basic and extensible post-processing of eddy covariance flux data with REddyProc. Biogeosciences 15, 5015–5030 (2018).
ADS  CAS  Google Scholar 

27.
Luo, Y. & Zhou, X. Soil respiration and the environment. (Academic Press, An Imprint of Elsevier Science, London, 2006).

28.
Hoffmann, M. et al. Automated modeling of ecosystem CO2 fluxes based on periodic closed chamber measurements: a standardized conceptual and practical approach. Agric. Meteorol. 200, 30–45 (2015).
Google Scholar 

29.
Rochette, P. & Hutchinson, G. Measurement of soil respiration in situ: chamber techniques. in Micrometeorology in agricultural systems, Agron. Monogr. 47 (ASA, CSSA and SSSA, Madison, WI 2005).

30.
Peichl, M. et al. A 12-year record reveals pre-growing season temperature and water table level threshold effects on the net carbon dioxide exchange in a boreal fen. Environ. Res. Lett. 9, 055006 (2014).
ADS  Google Scholar 

31.
Nilsson, M. et al. Contemporary carbon accumulation in a boreal oligotrophic minerogenic mire – a significant sink after accounting for all C-fluxes. Glob. Change Biol. 14, 2317–2332 (2008).
ADS  Google Scholar 

32.
Qiu, C. et al. ORCHIDEE-PEAT (revision 4596), a model for northern peatland CO2, water, and energy fluxes on daily to annual scales. Geosci. Model Dev. 11, 497–519 (2018).
ADS  CAS  Google Scholar 

33.
Abdalla, M. et al. Simulation of CO2 and attribution analysis at Six European Peatland sites using the ECOSSE model. Water, Air, Soil Pollut. 225, 2182 (2014).
ADS  Google Scholar 

34.
Metzger, C., Nilsson, M. B., Peichl, M. & Jansson, P.-E. Parameter interactions and sensitivity analysis for modelling carbon heat and water fluxes in a natural peatland, using CoupModel v5. Geosci. Model Dev. 9, 4313–4338 (2016).
ADS  CAS  Google Scholar 

35.
Jung, M. et al. Global patterns of land-atmosphere fluxes of carbon dioxide, latent heat, and sensible heat derived from eddy covariance, satellite, and meteorological observations. J. Geophys. Res. 116, G00J07 (2011).
Google Scholar 

36.
Ai, J. et al. MODIS-based estimates of global terrestrial ecosystem respiration. J. Geophys. Res. Biogeosciences 123, 326–352 (2018).
ADS  Google Scholar 

37.
Xiao, J. et al. Remote sensing of the terrestrial carbon cycle: a review of advances over 50 years. Remote Sens. Environ. 233, 111383 (2019).

38.
Badawy, B., Arora, V. K., Melton, J. R. & Nassar, R. Modeling the diurnal variability of respiratory fluxes in the Canadian Terrestrial Ecosystem Model (CTEM). J. Adv. Model. Earth Syst. 8, 614–633 (2016).
ADS  Google Scholar 

39.
Wu, Y., Verseghy, D. L. & Melton, J. R. Integrating peatlands into the coupled Canadian Land Surface Scheme (CLASS) v3.6 and the Canadian Terrestrial Ecosystem Model (CTEM) v2.0. Geosci. Model Dev. 9, 2639–2663 (2016).
ADS  CAS  Google Scholar 

40.
Waddington, J. M., Rotenberg, P. A. & Warren, F. J. Peat CO2 production in a natural and cutover peatland: implications for restoration. Biogeochemistry 54, 115–130 (2001).
CAS  Google Scholar 

41.
Glatzel, S., Basiliko, N. & Moore, T. Carbon dioxide and methane production potentials of peats from natural, harvested and restored sites, eastern Québec, Canada. Wetlands 24, 261–267 (2004).
Google Scholar 

42.
Hoyos-Santillan, J. et al. Quality not quantity: organic matter composition controls of CO2 and CH4 fluxes in neotropical peat profiles. Soil Biol. Biochem. 103, 86–96 (2016).
CAS  Google Scholar 

43.
Moore, T. R. & Dalva, M. Methane and carbon dioxide exchange potentials of peat soils in aerobic and anaerobic laboratory incubations. Soil Biol. Biochem. 29, 1157–1164 (1997).
CAS  Google Scholar 

44.
Nilsson, M. & Öquist, M. Partitioning litter mass loss into carbon dioxide and methane in peatland ecosystems. in Carbon Cycling in Northern Peatlands (American Geophysical Union, Washington, DC, 2009).

45.
Blodau, C., Basiliko, N. & Moore, T. R. Carbon turnover in peatland mesocosms exposed to different water table levels. Biogeochemistry 67, 331–351 (2004).
CAS  Google Scholar 

46.
D’Angelo, B. et al. Soil temperature synchronisation improves representation of diel variability of ecosystem respiration in Sphagnum peatlands. Agric. Meteorol. 223, 95–102 (2016).
Google Scholar 

47.
Phillips, C. L., Nickerson, N., Risk, D. & Bond, B. J. Interpreting diel hysteresis between soil respiration and temperature. Glob. Change Biol. 17, 515–527 (2011).
ADS  Google Scholar 

48.
Chapin, F. S. III, Matson, P. A. & Mooney, H. A. Carbon input to terrestrial ecosystems. in principles of terrestrial ecosystem. Ecology. (Springer, New York, NY, 2002).
Google Scholar 

49.
Vargas, R. & Allen, M. F. Environmental controls and the influence of vegetation type, fine roots and rhizomorphs on diel and seasonal variation in soil respiration. N. Phytol. 179, 460–471 (2008).
CAS  Google Scholar 

50.
Bahn, M., Schmitt, M., Siegwolf, R., Richter, A. & Brüggemann, N. Does photosynthesis affect grassland soil-respired CO2 and its carbon isotope composition on a diurnal timescale? N. Phytol. 182, 451–460 (2009).
CAS  Google Scholar 

51.
Laine, A. M. et al. Abundance and composition of plant biomass as potential controls for mire net ecosytem CO2 exchange. Botany 90, 63–74 (2012).
CAS  Google Scholar 

52.
Goulden, M. L. & Crill, P. M. Automated measurements of CO2 exchange at the moss surface of a black spruce forest. Tree Physiol. 17, 537–542 (1997).
CAS  PubMed  Google Scholar 

53.
Bubier, J., Crill, P., Mosedale, A., Frolking, S. & Linder, E. Peatland responses to varying interannual moisture conditions as measured by automatic CO2 chambers. Glob. Biogeochem. Cycles 17, 1066 (2003).
ADS  Google Scholar 

54.
Bond-Lamberty, B., Bronson, D., Bladyka, E. & Gower, S. T. A comparison of trenched plot techniques for partitioning soil respiration. Soil Biol. Biochem. 43, 2108–2114 (2011).
CAS  Google Scholar 

55.
Lai, D. Y. F., Roulet, N. T., Humphreys, E. R., Moore, T. R. & Dalva, M. The effect of atmospheric turbulence and chamber deployment period on autochamber CO2 and CH4 flux measurements in an ombrotrophic peatland. Biogeosciences 9, 3305–3322 (2012).
ADS  CAS  Google Scholar 

56.
Brændholt, A., Larsen, K. S., Ibrom, A. & Pilegaard, K. Overestimation of closed-chamber soil CO2 effluxes at low atmospheric turbulence. Biogeosciences 14, 1603–1616 (2017).
ADS  Google Scholar 

57.
Peichl, M., Sonnentag, O. & Nilsson, M. B. Bringing color into the picture: using digital repeat photography to investigate phenology controls of the carbon dioxide exchange in a boreal mire. Ecosystems 18, 115–131 (2015).
CAS  Google Scholar 

58.
Grinsted, A., Moore, J. C. & Jevrejeva, S. Application of the cross wavelet transform and wavelet coherence to geophysical time series. Nonlinear Process. Geophys. 11, 561–566 (2004).
ADS  Google Scholar 

59.
Welch, P. The use of fast Fourier transform for the estimation of power spectra: a method based on time averaging over short, modified periodograms. IEEE Trans. Audio Electroacoust. 15, 70–73 (1967).
ADS  Google Scholar  More