Philippa Kaur

Jiang, Z. H. Bamboo and Rattan in the World (China Forest Publishing House, 2007).
Google Scholar 
Zhou, B. Z., Fu, M. Y., Xie, J. Z., Yang, X. S. & Li, Z. C. Ecological functions of bamboo forest: Research and application. J. For. Res. 16, 143–147 (2005).Article 

Google Scholar 
Su, W., Fan, S., Zhao, J. & Cai, C. Effects of various fertilization placements on the fate of urea-15N in moso bamboo forests. For. Ecol. Manag. 453, 117632 (2019).Article 

Google Scholar 
Zhao, J. et al. Ammonia volatilization and nitrogen runoff losses from moso bamboo forests under different fertilization practices. Can. J. For. Res. 49(3), 213–220 (2019).CAS 
Article 

Google Scholar 
Lima, R. A. F., Rother, D. C., Muler, A. E., Lepsch, I. F. & Rodrigues, R. R. Bamboo overabundance alters forest structure and dynamics in the Atlantic Forest hotspot. Biol. Conserv. 147(1), 32–39 (2012).Article 

Google Scholar 
Kobayashi, K., Kitayama, K. & Onoda, Y. A. A simple method to estimate the rate of the bamboo expansion based on one-time measurement of spatial distribution of culms. Ecol. Res. 33(6), 1137–1143 (2018).CAS 
Article 

Google Scholar 
Xu, Q. F. et al. Rapid bamboo invasion (expansion) and its effects on biodiversity and soil processes. Glob. Ecol. Conserv. 21, e00787 (2020).Article 

Google Scholar 
Isagi, Y. & Torii, A. Range expansion and its mechanisms in a naturalized bamboo species, Phyllostachys pubescens, Japna. J. Sustain. Forest. 6(1–2), 127–141 (1997).Article 

Google Scholar 
Dong, C. L. et al. Effect of new rhizome growth on the fringe of the forest of Phyllostachys heterocycla cv. pubescens by different measure. J. Anhui Agric. Univ. 27(2), 150–153 (2000) (In Chinese with English abstract).
Google Scholar 
Bai, S. B. et al. Plant species diversity and dynamics in forests invaded by Moso bamboo (Phyllostachys edulis) in Tianmu Mountain Nature Reserve. Biodivers. Sci. 21(3), 288–295 (2013) (In Chinese with English abstract).Article 

Google Scholar 
Lin, Q. Q., Wang, B., Ma, Y. D., Wu, C. Y. & Zhao, M. S. Effects of Phyllostachys pubescens forest expansion on biodiversity in Tianmu Mountain Nature Reserve. J. Northeast Forest Univ. 42(9), 43–47 (2014) (In Chinese with English abstract).
Google Scholar 
Okutomi, K., Shinoda, S. & Fukuda, H. Causal analysis of the invasion of broad-leaved forest by bamboo in Japan. J. Veg. Sci. 7(5), 723–728 (1996).Article 

Google Scholar 
Ouyang, M. et al. Effects of the expansion of Phyllostachys edulis on species composition, structure and diversity of the secondary evergreen broad-leaved forests. Biodivers. Sci. 24(6), 649–657 (2016) (In Chinese with English abstract).Article 

Google Scholar 
Larpkern, P., Moe, S. R. & Totland, Ø. The effects of environmental variables and human disturbance on woody species richness and diversity in a bamboo-deciduous forest in northeastern Thailand. Ecol. Res. 24(1), 147–156 (2009).Article 

Google Scholar 
Larpkern, P., Moe, S. R. & Totland, Ø. Bamboo dominance reduces tree regeneration in a disturbed tropical forest. Oecologia 165(1), 161–168 (2011).ADS 
PubMed 
Article 

Google Scholar 
Griscom, B. W. & Ashton, M. S. A self-perpetuating bamboo disturbance cycle in a neotropical forest. J. Trop. Ecol. 22(5), 587–597 (2006).Article 

Google Scholar 
Yin, J. et al. Abandonment lead to structural degradation and changes in carbon allocation patterns in Moso bamboo forests. For. Ecol. Manag. 449, 117449 (2019).Article 

Google Scholar 
Suzuki, S. & Nakagoshi, N. Expansion of bamboo forests caused by reduced bamboo-shoot harvest under different natural and artificial conditions. Ecol. Res. 23(4), 641–647 (2008).Article 

Google Scholar 
Cai, L., Zhang, R. L., Li, C. F. & Ding, Y. A method to inhabit the expansion of Phyllostachys pubescens stands based on the analysis of underground rhizome. J. Northeast Forest Univ. 31(5), 68–70 (2003) (In Chinese with English abstract).
Google Scholar 
Rice, E. L. Allelopathy (Academic Press, 1984).
Google Scholar 
Huang, W. et al. Allelopathic effects of Cinnamomum septentrionale leaf litter on Eucalyptus grandis saplings. Glob. Ecol. Conserv. 21, e00872 (2020).Article 

Google Scholar 
Turk, M. A. & Tawaha, A. M. Allelopathic effect of black mustard (Brassica nigra L.) on germination and growth of wild oat (Avena fatua L.). Crop Prot. 22(4), 673–677 (2003).Article 

Google Scholar 
Kong, C. H., Li, H. B., Hu, F., Xu, X. H. & Wang, P. Allelochemicals released by rice roots and residues in soil. Plant Soil 288(1–2), 47–56 (2006).CAS 
Article 

Google Scholar 
Duke, S. O. Weeding with allelochemicals and allelopathy-a commentary. Pest Manag. Sci. 63(4), 307–307 (2007).CAS 
PubMed 
Article 

Google Scholar 
Braine, J. W., Curcio, G. R., Wachowicz, C. M. & Hansel, F. A. Allelopathic effects of Araucaria angustifolia needle extracts in the growth of Lactuca sativa seeds. J. For. Res. 17(5), 440–445 (2012).CAS 
Article 

Google Scholar 
Soltys, D., Krasuska, U., Bogatek, R. & Gniazdowska, A. Allelochemicals as bioherbicides-present and perspectives. In Herbicides-current research and case studies in use (eds Price, A. J. & Kelton, J. A.) (IntechOpen, 2013).
Google Scholar 
Qin, J. H. et al. Allelopathic effects of the different allelochemical pathways of sesame extracts. J. Foshan Univ. 31(4), 1–5 (2013) (In Chinese with English abstract).
Google Scholar 
Khasabulli, B. D., Musyimi, D. M., George, O. & Gichuhi, M. N. Allelopathic effect of Bidens Pilosa on seed germination and growth of Amaranthus Dubius. J. Asian Sci. Res. 8(3), 103–112 (2018).
Google Scholar 
Boter, M. et al. An integrative approach to analyze seed germination in Brassica napus. Front. Plant Sci. 10, 1342 (2019).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Saha, D., Marble, S. C. & Pearson, B. J. Allelopathic effects of common landscape and nursery mulch materials on weed control. Front. Plant Sci. 9, 733 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Bachheti, A., Sharma, A., Bachheti, R. K., Husen, A. & Pandey, D. P. Plant allelochemicals and their various applications. In Co-Evolution of Secondary Metabolites Reference Series in Phytochemistry (eds Mérillon, J. M. & Ramawat, K.) (Springer, 2020).
Google Scholar 
Duary, B. Effect of leaf extract of sesame (Sesamum indicum L.) on germination and seedling growth of blackgram (Vigna mungo L.) and rice (Oryza sativa L.). Allelopathy J. 10(2), 153–156 (2002).
Google Scholar 
Soleymani, A. & Shahrajabian, M. H. Study of allelopathic effects of sesame (Sesamum indicum) on canola (Brassica napus) growth and germination. Intl. J. Agri. Crop Sci. 4(4), 183–186 (2012).
Google Scholar 
Gorai, M., Aloui, W. E., Yang, X. & Neffati, M. Toward understanding the ecological role of mucilage in seed germination of a desert shrub Henophyton desert: Interactive effects of temperature, salinity and osmotic stress. Plant Soil 374(1–2), 727–738 (2014).CAS 
Article 

Google Scholar 
Wang, C., Wu, B. & Jiang, K. Allelopathic effects of Canada goldenrod leaf extracts on the seed germination and seedling growth of lettuce reinforced under salt stress. Ecotoxicology 28, 103–116 (2019).CAS 
PubMed 
Article 

Google Scholar 
Wang, X. L. et al. Allelopathic effects of exotic mangrove species Laguncularia racemosa on Bruguiera gymnorhiza. J. Xiamen Univ. 56(3), 339–345 (2017) (In Chinese with English abstract).
Google Scholar 
Shah, A. N. et al. Allelopathic influence of sesame and green gram intercrops on cotton in a replacement series. Clean: Soil, Air, Water 45(1), 1–10 (2017).
Google Scholar 
Amare, T. Allelopathic effect of aqueous extracts of parthenium (Parthenium hysterophorus L.) parts on seed germination and seedling growth of maize (Zea mays L.). J. Agric. Crop 4(12), 157–163 (2018).
Google Scholar 
Yan, X. F., Du, Q., Fang, S. & Zhou, L. B. Allelopathic effects of water extraction of Rhus typhina on Zea mays seeds germination. Seed 29(3), 15–18 (2010) (In Chinese with English abstract).
Google Scholar 
Yan, X. F., Zhou, Y. F. & Du, Q. Allelopathic effects of water extraction from root and leaf litter of Rhus typhina on the germination of wheat seeds. Seed 30(5), 17–20 (2011) (In Chinese with English abstract).
Google Scholar 
Wang, X. et al. Allelopathic effects of aqueous leaf extracts from four shrub specious on seed germination and initial growth of Amygdalus pedunculata Pall. Forests 9, 711 (2018).Article 

Google Scholar 
Alencar, N. L. M. et al. Ultrastructural and biochemical changes induced by salt stress in Jatropha curcas seeds during germination and seedling development. Funct. Plant Biol. 42(9), 865–874 (2015).CAS 
PubMed 
Article 

Google Scholar 
Lozano-Isla, F., Campos, M. L. O., Endres, L., Bezerra-Neto, E. & Pompelli, M. F. Effects of seed storage time and salt stress on the germination of Jatropha curcas L. Ind. Crop Prod. 118, 214–224 (2018).CAS 
Article 

Google Scholar 
Wu, J. R., Chen, Z. Q. & Peng, S. L. Allelopathic potential of invasive weeds: Alternanthera philoxeroide, Ipomoea cairica and Spartina alterniflora. Allelopathy J. 17(2), 279–285 (2006).
Google Scholar 
Sahu, A. & Devkota, A. Allelopathic effects of aqueous extract of leaves of Mikania micrantha H.B.K. on seed germination and seedling growth of Oryza sativa L. and Raphanus sativus L. Sci. World 11(11), 70–77 (2013).Article 

Google Scholar 
Gatti, A. B., Ferreira, A. G., Arduin, M. & Perez, S. C. G. D. A. Allelopathic effects of aqueous extracts of Artistolochia esperanzae O.Kuntze on development of Sesamum indicum L. seedlings. Acta Bot. Bras. 24(2), 454–461 (2010).Article 

Google Scholar 
Hou, Y. P. et al. Effects of litter from dominant tree species on seed germination and seedling growth of exotic plant Rhus typhina in hilly areas in Shandong peninsula. Sci. Silvae Sin. 52(6), 28–34 (2016) (In Chinese with English abstract).
Google Scholar 
Jiang, Z. et al. Effects of root exudates from Picea asperata seedlings on the seed germination and seedling growth of two herb species. Sci. Silvae Sin. 55(6), 160–166 (2019) (In Chinese with English abstract).
Google Scholar 
Hagan, D. L., Jose, S. & Lin, C. Allelopathic exudates of cogongrass (Imperata cylindrical): Implications for the performance of native pine savanna plant species in the Southeastern US. J. Chem. Ecol. 39, 312–322 (2013).CAS 
PubMed 
Article 

Google Scholar 
Cheng, F. & Cheng, Z. Research progress on the use of plant allelopathy in agriculture and the physiological and ecological mechanisms of allelopathy. Front. Plant Sci. 6, 1020 (2015).PubMed 
PubMed Central 

Google Scholar 
Bogatek, R., Gniazdowska, A., Zakrzewska, W., Oracz, K. & Gawronski, S. W. Allelopathic effects of sunflower extracts on mustard seed germination and seedling growth. Biol. Plantarum 50(1), 156–158 (2006).Article 

Google Scholar 
Politycka, B. Peroxidase activity and lipid peroxidation in roots of cucumber seedlings influenced by derivatives of cinnamic and benzoic acids. Acta Physiol. Plant. 18(4), 365–370 (1996).CAS 

Google Scholar 
Williamson, G. B. & Richardson, D. Bioassays for allelopathy: Measuring treatment responses with independent controls. J. Chem. Ecol. 14(1), 181–187 (1988).Article 

Google Scholar  More

In the version of this article initially published, there were mistakes in affiliations 1, 2 and 6. The corrected affiliations should read as follows: 1. Department of Earth System Science, Ministry of Education Key Laboratory for Earth System Modeling, Institute for Global Change Studies, Tsinghua University, Beijing, China; 2. Ministry of Education Ecological Field Station for East Asian Migratory Birds, Department of Earth System Science, Tsinghua University, Beijing, China; 6. Department of Geography, Department of Earth Sciences, and Institute for Climate and Carbon Neutrality, The University of Hong Kong, Hong Kong, China. The affiliations have been corrected in the HTML and PDF versions of the article. More

The quality of global maps can be assessed in different ways. One way is global assessment where a single statistic is chosen to summarize the quality of the entire map: the map accuracy. For a categorical variable, this can be the probability that for a randomly chosen location on the map, the map value corresponds to the true value. For a continuous variable, it can be the RMSE, describing for a randomly chosen location on the map the expected difference between the mapped value and the true value. When a probability sample, such as a completely spatially random sample, is available for the area for which a global assessment is needed, then map accuracy can be estimated model-free (also called design-based, e.g., by using the unweighted sample mean in case of a completely spatially random sample). This circumvents modeling of spatial correlation because observations are independent by design6,9. This approach is called model-free because no model needs to be assumed about the distribution or correlation of the data: the only source of randomness is the random selection of sample units from a target population. If a probability sample is not available this approach cannot be used, and automatically the accuracy assessment approach becomes model-based10, which involves modeling a spatial process by assuming distributions and taking spatial correlations into account, and choosing estimation methods accordingly.Using naive random n-fold or leave-one-out cross-validation methods (or a simple random train-test split) to assess global model quality (usually equated with map accuracy) makes sense when the data are independent and identically distributed. When this is not the case, dependencies between nearby samples, e.g., in a spatial cluster, are ignored and result in biased, overly optimistic model assessment, as shown in, e.g., Ploton et al.5. Alternative cross-validation approaches such as spatial cross-validation5,11 that control for such dependencies are the only way to overcome this bias. Different spatial cross-validation strategies have been developed in the past few years, all aiming at creating independence between cross-validation folds5,11,12,13. Cross-validation creates prediction situations artificially by leaving out data points and predicting their value from the remaining points. If the aim is to assess the accuracy of a global map, the prediction situations created need to resemble those encountered while predicting the global map from the reference data (see Fig. 1 and discussions in Milà et al.14). This occurs naturally when reference data were obtained by (completely spatially random) probability sampling, but in other cases, this has to be forced for instance by controlling spatial distances (spatial cross-validation). Such forcing, however, is only possible when the distances in space that need to be resembled are available in the reference data. In the extreme case where all reference data come from a single cluster, this is impossible. When all reference data come from a small number of clusters, larger distances are available between clusters but do not provide substantial independent information about variation associated with these distances. Lack of information about larger distances means that we cannot assess the quality of predictions associated with such distances and cannot properly estimate global quality measures. Alternative approaches such as experiments with synthetic data15 or a validation using independent data at a higher level of integration16 would then be options to support confidence in the predictions.Another way of accuracy assessment is local assessment: for every location, a quality measure is reported, again as probability or prediction error. Such a local assessment predicts how close the map value is to newly observed values at particular locations. If the measurement error is quantified explicitly, a smoother, measurement-error-free value may be predicted10. If the model accounts for change of support10,17, predictions errors may refer to average values over larger areas such as 1 × 1, 5 × 5, or 10 × 10 km grid cells. Examples of local assessment in the context of global ecological mapping are modeled prediction errors using Quantile Regression Forests18 or mapped variance of predictions made by ensembles1,2. Neither of these examples quantifies spatial correlation or measurement error, or addresses change of support, as it is known from other modeling frameworks19. By omitting to model the spatial process, the local accuracy estimates as presented in the global studies that motivated this comment are disputable.The difference between global and local assessment is striking, in particular for global maps. A global, single number averages out all variability in prediction errors, and obscures any differences, e.g., between continents or climate zones. It is of little value for interpreting the quality of the map for particular regions. More

The data used here provides evidence that particular relationships may determine sex-specific patterns of witchcraft accusation. Cases where women were targeted frequently came from affinal kin, while those directed at men were often from unrelated individuals and blood relatives. Most previous research on factors that determine the sex of accused ‘witches’ has largely consisted of qualitative studies of a single society or a few societies48, or historical studies that have not tested for correlations49. Our findings, in support of the overarching hypothesis that accusations may be driven by various forms of competition, can be tentatively aligned with evolutionary literature on patterns of intrasexual and kin competition, intersexual conflict and polygamous mating30,31,50.Men were more often accused than women in our sample, although we did not have a prediction in relation to this. But the finding suggests how overall patterns of competition within relationships may contribute to societal ‘phenotypes’ of witches as male or female. The ethnography of the Ndembele perhaps indicates why women were less frequently targeted in Bantu societies: ‘in a case of witchcraft, the complainant is actuated by caprice, jealousy or pique; and the defendant is a person of wealth or popularity, and is always a man, for the women have neither wealth nor honor worth coveting’51.Our predictions about how the sex of accused ‘witches’ might be associated with particular relationship categories were supported. The majority of accusations targeting men came from unrelated individuals, which is unsurprising, as inclusive fitness52 would not mitigate the effects of competition between them. Blood relatives were the next most common relationship category directing accusations at men. This aligns with more recent studies indicating that witchcraft fears between family members are significant in parts of Africa, to the extent that they can be construed as ‘the dark side of kinship’53. In evolutionary terms, kin may compete with one another in environments where resources are limited30,31,50 and in societies with patrilineal inheritance related males, and particularly brothers, compete for resources in order to marry31. This aligns with an ethnographic observation that among the Banyoro witchcraft accusations often occurred between brothers over inheritance, but not between brothers and sisters, whose interests did not conflict21. The situations relating to accusations of men were also often connected to the acquisition of wealth and status, such as rivalry over village headmanships32, power struggles between a chief’s counsellors54 or disputes over inheritance55. These connections can be found in more recent contexts such as twentieth century Ghana, where notions of obtaining political power and wealth through occult means involving human sacrifice were pervasive56.Accusations of women were more likely to come from affines. Husbands were the largest category of affinal kin to accuse women (Supplementary Fig. 2). The higher rate of accusations from husbands to wives than wives to husbands aligns with evolutionary perspectives suggesting male coercion of females is a strategy to maximize male reproductive success39,41. Accusations of wives who were suspected of being unfaithful can be interpreted as a strategy for reducing investment in unrelated offspring35,41. In a case from the Shona a woman gave birth to a stillborn child. This was attributed to an affair before marriage, and was followed by divorce and the repayment of bridewealth to her husband, who commented she was ‘a witch, a woman who had killed her own child’48. Other ethnographic accounts suggest accusations of wives by husbands were an attempt to gain control within the marital relationship55.A significant number of accusations of women by affinal kin were from co-wives in polygynous marriages, and these were often notably associated with jealousy connected to a husband’s attention and investment32. Evolutionary models predict competition for reproductive resources would occur among co-resident breeding women57, as has been found to occur among the Mosuo of southwest China58. In the patrilocal social systems that are predominant in our sample, women disperse at marriage and are isolated from kin, so conflict may be more extreme30. This is consistent with ethnographic observations reporting that the relationship between co-wives in polygynous marriages was often (although not always) marked by conflict, and liable to produce witchcraft accusations38,59.There were accusations of women from other categories of their affinal kin (Supplementary Fig. 2). These again may result from competition for a husband’s time and resources between his kin and wife. New wives may be vulnerable in environments where they enter their husband’s families as unrelated strangers, and are potentially expendable, at least before the arrival of offspring. Some accounts of accusations indicate that accusations of wives by in-laws in patrilocal households are common29.Accusations directed at elderly individuals targeted women more often than men. This may form part of a broader pattern of geronticide: societies close to subsistence-level are documented as sometimes accepting the abandonment or killing of elderly people19,60. In modern Tanzania, ‘witches’ are mostly post-reproductive women, who are more likely to be murdered in periods of income shock19. This is also the case in contemporary Ghana, where accusations are frequently directed at middle-aged or elderly women, whose families may subsequently cease to provide them with financial or material assistance61. In our sample, elderly women may have been targeted more frequently as a result of longer female lifespans: in a polygynous society, men may marry younger women, so wives would be widowed at an earlier age than husbands. Among the Bantu, older men were accused, but some were possibly protected by their status.Accusers’ payoffs from accusations are not always explicit but they can be inferred. The most common outcome of accusations in our sample was that accused ‘witches’ were exiled from their communities or forced to move from where they were living. This would mean resources and cooperative assistance they would have used became available to their accusers or others nearby. Where the accused acquires a negative reputation, which was the second most common outcome, there may be a subtle removal of benefits, which may be preferred to direct ‘punishment’ as it is less costly62. Accusers’ gains need not be direct, as harming behaviours may reduce the overall pressure of competition in an environment28. 8% of accusations in the sample resulted in the acquisition of either resources or political positions from the accused, or in preventing the accused from acquiring them. Where the accused were penalised in other ways, such performing ceremonies to reconcile with accusers, this is perhaps akin to classic cooperation models involving the punishment of defectors (although the accused may not actually be uncooperative)11, providing accusers with subordinate partners who offer fitness benefits to avoid more serious allegations63. Where an accusation does not ‘stick’, ethnographic accounts sometimes indicate it was reversed through divination or ordeal54. In other cases, for various reasons accusations are short-lived and forgotten about4. Finally, although not tested in this dataset, accusers may gain informal prestige and dominance, an outcome analogous to competitive punishment63.Not all of the cases in our dataset support the hypothesis that witchcraft accusations are a mechanism for competition. There is a significant proportion where the accusation of a particular individual appears to be incidental, or dependent a on circumstantial association between the ‘witch’ and a negative event. Such accusations are unlikely to provide accusers with a competitive advantage. There are several possible explanations for such cases. They are in line with the hypothesis that witchcraft belief arises from attempts to identify the cause of an impactful misfortune3,4. Cultural evolutionary explanations of witchcraft beliefs suggest that they are a maladaptive attempt to explain misfortune. Although it is inaccurate, belief in witches is maintained through bias and selective inattention to evidence that would otherwise counter it64. Alternatively this could be viewed under the contention that superstitious beliefs (or errors in attributing cause and effect) are broadly adaptive if they occasionally lead individuals to acts which provide them with fitness benefits65.Although witchcraft accusations may be a mechanism for mitigating the damage to accusers’ reputations in harmful competitive acts, as with any behavioural strategy it is not without risks. Accusers may suffer costs in the form subsequent reputational damage or counter-accusations, as with punishment63, depending on factors such the level of support for an accusation by other members of the community.One limitation of our dataset is that it contains realized allegations of witchcraft, that cannot be tested against baseline population measures. We could not examine the risk that a particular individual, such as an elderly woman, would be accused. Instead, the analysis shows the odds, given an accusation occurred, that the ‘witch’ was male or female, given certain predictors. For example, if the accused was elderly, there are increased odds they were female rather than male.A dataset using historic witchcraft cases is almost certainly affected by selection bias. Cases with sensational outcomes are more likely to be reported, and cases that are dismissed or where the accused removes themselves from their accusers are liable to be overlooked19. Most incidents in our sample were reported anecdotally. Obtaining a random sample of witchcraft accusations within a population is challenging, if not impossible1,66. Attempts to systematically collect cases within a given location and timeframe cannot guarantee that all are brought to the attention of researchers19. Comparative studies of this kind usually use all the data that is available and control for confounding effects. Our sensitivity analyses suggest the large number of accusations of men in the dataset probably reflects patterns of accusations in these societies, rather than male-focused bias from ethnographers. There are many accounts of cultures where witches are predominantly male33,34,49. But the accuracy of historic ethnographic accounts cannot be verified, especially in relation to one-off events such as witchcraft accusations, just as it is unclear how much uncertainty there is in the ethnographic record overall67. Ethnographers may not always have noted the characteristics of the individuals involved, or there may be times where they were mistaken in reporting the circumstances surrounding an accusation. There are several explanations for cases where the identities of accusers or purported victims of witchcraft were not reported. Not all cases had identifiable ‘victims’, for example when the accused was thought to have used witchcraft to promote their own success, or ethnographers could not denote the relationship between the accused and their accusers when suspicions of witchcraft were communicated through general gossip. In a small number of cases, ethnographer perspectives on accusations (and possible inability to access further information) are salient, as they may ascribe more importance to one relationship over another in reporting a case, such as a witch’s envy of their victim, or a witch’s argument with an accuser.However, it is likely that ethnographers were for the most part accurate in documenting variables of interest such as the sex of an accused individual and their relationships with accusers. There is less certainty in relation to the situation connected to an accusation, especially taking into recent research that indicates the prevalence of phenomena such as the misperception of causation68,69. Our attempts to account for such possibilities with sensitivity analyses and meta-data on the production of ethnographies cannot conclusively provide reassurance that bias has not affected results, and so this section of the analysis should be treated with caution and regarded as exploratory. The situations documented in our study do however align with accounts of accusations from more contemporary observers and studies from different geographic locations, suggesting that similar causes of accusations arise convergently in different societies. For example in modern contexts accusations have led to accusers gaining land or property in India6 and cessation of the obligation to provide material and financial assistance to elderly relatives in Ghana61. One advantage of our cross-cultural data being drawn from numerous ethnographies is that it is not reliant on the perspective of one individual, meaning that random perceptual error or individual (as opposed to cultural) bias is more likely to be mitigated in the results than would be the case in the study of a single culture by one ethnographer.As a further limitation, we were reliant on accessible ethnographic records from the best-documented societies. Although selection bias in favour of better described societies is present in our sample, this should not impact the main aim of this research, which is to understand the determinants of witchcraft accusations being directed at male or female targets.Overall our findings may indicate allegations of witchcraft stem from diverse forms of competition between individuals. This aligns with evolutionary approaches to competition and conflict. Accusations may provide fitness benefits by allowing individuals to target competitors, but the exact form and direction of competition is determined by aspects of socio-ecology. This in turn influences which sex is most likely to be accused and the overall portrayal of witches in a society. Accusations may be more likely to occur in some relationships rather than others, when there is a gain for the accuser, as in disputes over inheritance and property, or where another individual may pose a threat, or by simply reducing numbers of competitors. The success of witchcraft accusations in removing competitors and their flexibility as an adaptive strategy may explain their widespread distribution. More

Data collection of flowering and fruiting phenologyMonthly reproductive phenology data recorded over 35 years (from April 1976 to September 2010) were collected from the Bulletin Fenologi Biji Benih dan Anak Benih (Bulletin of Seed and Seedling Phenology), which was deposited at the FRIM library. The bulletin reported seed and seedling availabilities and the flowering and fruiting phenology of trees at several research stations in Malaysia. The present study collected flowering and fruiting records of trees grown in FRIM arboretums located approximately 12 km northwest of Kuala Lumpur, Malaysia (latitude 3°24 ‘N, longitude 101°63 ‘E, elevation 80 m). There are both dipterocarp and non-dipterocarp arboretums in FRIM, both of which were founded in 1929. These arboretums preserve and maintain living trees for research and other purposes. Each month, three research staff members of FRIM with sufficient phenology monitoring training made observations with binoculars to record the presence of flowers and fruits on trees of each species on the forest floor from April 1976 to September 2010. The phenological status of the trees was recorded as flowering during the developmental stages from flower budding to blooming and as fruiting during the developmental stages from the occurrence of immature fruit to fruit ripening. Because only one or two individuals per species are grown at the FRIM arboretums, the flowering and fruiting phenology were monitored using these individuals. The resultant flowering and fruiting phenology data included a time series of binary data (1 for presence and 0 for absence) with a length of 417 months.The original data included 112 dipterocarp and 240 non-dipterocarp species. We excluded 17 dipterocarps and 125 non-dipterocarp species based on the following five criteria for data accuracy.

1.

Percentage of missing values is ≤50%: If the monthly flowering or fruiting phenology data of a species included a substantially large number of missing values ( >50%), the species was excluded.

2.

Stable flowering period: We considered an observation to be unreliable if the flowering period was significantly different among flowering events (if the coefficient of variation in the flowering period was larger or equal to 1.0).

3.

Flowering period is shorter than or equal to 12 months: we considered an observation to be unreliable if the flowering period was longer than 12 months because it was unlikely that the same tree would flower continuously for longer than 1 year.

4.

The flowering and fruiting frequencies were not significantly different between the first and second half of the census period: when the flowering frequency was zero for the first half of the observation period but was larger than 0.1 for the second half of the observation period, or when the flowering frequency was zero for the second half of the observation period but was larger than 0.1 for the first half of the observation period, we removed these species because data are not reliable (e.g., physiological conditions may have changed significantly). We adopted the same criteria for the fruiting phenology data.

5.

We removed overlapping species, herb species, and specimens with unknown species names.

After removing unreliable species based on the five criteria explained above, we obtained 95 dipterocarp and 115 non-dipterocarp species (Supplementary Data 1). We used these species for further analyses. It is unlikely that our final data includes trees that were replaced by young trees during the census period because newly planted seedlings do not flower over 20–30 years until they are fully grown to the reproductive stage ( >20–30 cm DBH)45.Detection of seasonality in reproductive phenologyTo compare the flowering and fruiting phenology seasonality among different families, nine families that included at least five species were used. The number of flowering or fruiting events was counted for each month from January to December during a census, and then the frequency distribution was drawn as a histogram. Similarly, we also generated a histogram for the seed dispersal month, which was calculated as the month when fruiting ended (i.e., when the binary fruiting phenology data changed from one to zero).Classification of phenological patternsTo classify the phenological patterns, we performed time-series clustering using the R package TSclust46 with the hierarchical clustering method based on the Dynamic Time Warping distance of the flowering phenology data of each species. For this analysis, time points at which there were missing values for at least one species were excluded. Because of the large number of missing values in non-Dipterocarpaceae species, we performed time-series clustering only for the Dipterocarpaceae species based on 394 time points in total. The number of phenological clusters was estimated based on AIC, as explained below.Climate dataDaily minimum, mean, and maximum temperatures and precipitation data monitored at the FRIM KEPONG (3° 14’ N, 101° 42’ E, elevation 97 m) weather station were provided by the Malaysian Meteorological Department. We used the daily minimum temperature for our analysis because there were fewer missing values compared to the numbers of missing daily mean and daily maximum temperature values. The periods in which climate data were available were from 1 March 1973 to 31 March 1996, and from 23 July 1997 to 20 April 2005. We removed periods in which there were missing values spanning longer than 5 days. When the range of missing values spanned a period shorter than 3 days, we approximated these missing values using the mean minimum temperatures recorded on the adjacent three days. Although solar radiation data were not available for our study, the use of precipitation is sufficient for model fitting because there is a significant negative correlation between solar radiation and precipitation in Southeast Asia47.Climate data generated by GCMsAs the future climate inputs, we used bias-corrected climate input data from 1 January 2050 to 31 December 2099, with a daily temporal resolution and a 0.5° spatial resolution, provided by the ISI-MIP project48; these data are based on the Coupled Model Intercomparison Project Phase 5 outputs from three GCMs: GFDL–ESM2M, IPSL–CM5A-LR, and MIROC5. To compare the flowering phenology between 1976–1996 and 2050–2099, bias-corrected GCM data from 1 May 1976 to 31 March 1996, were also used. This period (1 May 1976–31 March 1996) is consistent with the period used for model fitting. We selected daily minimum temperature and precipitation time series from the 0.5° grid cells corresponding to the study site for phenology monitoring at FRIM. To compare flowering phenology among regions, we also used the same set of data from three other regions in Southeast Asia: Trang Province in Thailand (7° 4’ N, 99° 47’ E), Lambir Hills National Park in Malaysia (4° 2’ N, 113° 50’ E), and central Kalimantan in Indonesia (0° 06’ S, 114° 0’ E). Because the study site in FRIM was not in the center of a 0.5° grid cell, we interpolated the data using four grid cells in the vicinity of the observation site. We used the weighted average according to the distance between each observation site and the center of each corresponding grid cell.Although the climate input data provided by ISI-MIP were already bias-corrected, we conducted additional bias correction at FRIM using a historical scenario for each GCM data set and the observed weather data from 1 January 1976 to 31 December 2004 based on previously presented protocol49. We did not implement any bias correction for the frequency of dry days or precipitation intensity of wet days49 because we only focused on the average precipitation.The variances in the annual fluctuation of the monthly mean precipitation were not the same between the observation data and historical GCM runs at FRIM. For all three GCMs (GFDL–ESM2M, IPSL–CM5A-LR, and MIROC5), the variances in the yearly fluctuation output by the GCMs tended to be larger than that of the observed data at the FRIM KEPONG weather station during winter and spring. On the other hand, during summer and fall, the variances output by the GCMs tended to be smaller than that of the observed data. These biases could not be corrected using the previous method49. Therefore, we conducted the following bias correction for these data:$${p}_{i,m,y}^{{{{{{rm{GCM}}}}}}* }={r}_{i,m,y}^{{{{{{rm{GCM}}}}}}}cdot left[{F}_{Gamma }^{-1}left({F}_{Gamma }left({delta }_{m,y}^{{{{{{rm{GCM}}}}}}}|{k}_{m,y},{theta }_{m,y}right)|{k}_{m,y}^{* },{theta }_{m,y}^{* }right)cdot {rho }_{m,y}^{{{{{{rm{GCM}}}}}}}right],$$
(1)
where ({p}_{i,m,y}^{{{{{{rm{GCM}}}}}}* }) is the bias-corrected precipitation value of the target GCM at year y, month m, and date i. In the equation, ({r}_{i,m,y}^{{{{{{rm{GCM}}}}}}}) is the ratio of the precipitation value of the GCM relative to the monthly mean value. Then, the following equation is used:$${r}_{i,m,y}^{{{{{{rm{GCM}}}}}}}=frac{{p}_{i,m,y}^{{{{{{rm{GCM}}}}}}}}{{bar{p}}_{m,y}^{{{{{{rm{GCM}}}}}}}},$$
(2)
where ({p}_{i,m,y}^{{{{{{rm{GCM}}}}}}}) is the precipitation value (not bias-corrected) of the GCM at year (y), month (m), and date i and ({bar{p}}_{m,y}^{{{{{{rm{GCM}}}}}}}) is the monthly mean precipitation value of the GCM at year (y) and month (m). In Eq. 1, ({F}_{Gamma }) represents the cumulative distribution function of a gamma distribution, ({F}_{Gamma }^{-1}) represents the inverse function of the cumulative distribution function of the gamma distribution, and ({k}_{m,y}) and ({theta }_{m,y}) are the shape parameters. In Eq. 1, ({delta }_{m,y}^{{{{{{rm{GCM}}}}}}}) indicates the deviation of the monthly mean from the normal climate value of the corresponding period, and this value is calculated as follows:$${delta }_{m,y}^{{{{{{rm{GCM}}}}}}}=frac{{bar{p}}_{m,y}^{{{{{{rm{GCM}}}}}}}}{{rho }_{m,y}^{{{{{{rm{GCM}}}}}}}},$$
(3)
where ({rho }_{m,y}^{{{{{{rm{GCM}}}}}}}) is the normal climate value during the target period. In this method, we defined the normal climate value as the mean of the monthly mean precipitation values over 31 years.$${rho }_{m,y}^{{{{{{rm{GCM}}}}}}}=frac{1}{31}mathop{sum }limits_{j=y-15}^{y+15}{bar{p}}_{m,j}^{{{{{{rm{GCM}}}}}}}.$$
(4)
When the mean of a gamma distribution is fixed at one, the shape parameters are represented as follows:$${k}_{m,y}=frac{1}{Vleft({delta }_{m,y}^{{{{{{rm{GCM}}}}}}}right)},$$
(5)
$${theta }_{m,y}=frac{1}{{k}_{m,y}},$$
(6)
where (Vleft({delta }_{m,y}^{{{{{{rm{GCM}}}}}}}right)) indicates the variance in ({delta }_{m,y}^{{{{{{rm{GCM}}}}}}}) at month (m) over 31 years.In this method, we assumed that the ({delta }_{m,y}^{{{{{{rm{GCM}}}}}}}) value follows a gamma distribution and that the ratio of the variance of ({delta }_{m,y}^{{{{{{rm{GCM}}}}}}}) to the variance of ({delta }_{m,y}^{{{{{{rm{obs}}}}}}}) is maintained even in the future scenario. Here, ({delta }_{m,y}^{{{{{{rm{obs}}}}}}}) represents the deviation of the monthly mean in the observation data from the normal climate value.$${delta }_{m,y}^{{{{{{rm{obs}}}}}}}=frac{{bar{p}}_{m,y}^{{{{{{rm{obs}}}}}}}}{{rho }_{m}^{{{{{{rm{obs}}}}}}}},$$
(7)
$${rho }_{m}^{{{{{{rm{obs}}}}}}}=frac{1}{28}mathop{sum }limits_{j=1976}^{2004}{bar{p}}_{m,y}^{{{{{{rm{obs}}}}}}}.$$
(8)
In the above equations, ({bar{p}}_{m,y}^{{{{{{rm{obs}}}}}}}) indicates the monthly mean precipitation value in the observed data. As mentioned above, because we assume that the ratio of the variance in ({delta }_{m,y}^{{{{{{rm{GCM}}}}}}}) to the variance in ({delta }_{m,y}^{{{{{{rm{obs}}}}}}}) is maintained, ({k}_{m,y}^{* }) and ({theta }_{m,y}^{* }) are calculated as follows:$${k}_{m,y}^{* }=frac{{k}_{m,y}}{alpha },$$
(9)
$${theta }_{m,y}^{* }=frac{1}{{k}_{m,y}^{* }},$$
(10)
where$$alpha =frac{Vleft({delta }_{m,y}^{{{{{{{rm{GCM}}}}}}}^{{{{{{rm{h}}}}}}}}right)}{Vleft({delta }_{m,y}^{{{{{{rm{obs}}}}}}}right)}.$$
(11)
In Eq. 11, ({delta }_{m,y}^{{{{{{{rm{GCM}}}}}}}^{{{{{{rm{h}}}}}}}}) is the deviation of the monthly mean of the historical GCM precipitation data from the normal climate value. Here, we defined the normal climate value as the average monthly mean during 1976–2004.The method proposed here is an original bias correction method, but the above equations are easily derived if we assume that the ({delta }_{m,y}^{{{{{{rm{GCM}}}}}}}) value follows a gamma distribution and that the ratio of the variance in ({delta }_{m,y}^{{{{{{rm{GCM}}}}}}}) to the variance in ({delta }_{m,y}^{{{{{{rm{obs}}}}}}}) is maintained even in the future scenario. Notably, because we combined this method with the bias correction method described previously49, Eq. 2 should be expressed as follows:$${r}_{i,m,y}^{{{{{{rm{GCM}}}}}}}=frac{{widetilde{p}}_{l,m,y}^{{{{{{rm{GCM}}}}}}}}{{bar{p}}_{m,y}^{{{{{{rm{GCM}}}}}}}},$$
(12)
where ({widetilde{p}}_{l,m,y}^{{{{{{rm{GCM}}}}}}}) is the precipitation data that are bias-corrected using the method described previously49. Bias-corrected data were compared with the data without bias correction (Supplementary Figs. 8–11).Statistical analyses and reproducibilityWe adopted previously presented models in which environmental triggers for floral induction accumulate for n1 days prior to the onset of floral induction21 (Supplementary Fig. 2). Flowers then develop for n2 days before opening (Supplementary Fig. 2). The model assumption of the time lag between floral induction and anthesis, which is denoted as n2, was validated by a previous finding in which the expression peaks of flowering-time genes, which are used as molecular markers of floral induction, were shown to occur at least one month before anthesis in Shorea curtisii19. S. curtissi is included in our data set. The CU at time t, ({{{{{rm{CU}}}}}}left(t|{theta }^{C}right)), is calculated as follows:$${{{{{rm{CU}}}}}}left(t|{theta }^{C}right)=mathop{sum }limits_{n={n}_{2}}^{{n}_{2}+{n}_{1}-1}{{{{{rm{max }}}}}}{bar{C}-xleft(t-nright),0},$$
(13)
where ({theta }^{C}=left{{n}_{1},{n}_{2},bar{C}right}) is the set of parameters and x(t) is the temperature at time t. Here, (bar{C}) indicates the threshold temperature. The term max{x1, x2} is a function that returns a larger value for the two arguments. Similarly, given ({theta }^{D}={{n}_{1},{n}_{2},bar{D}},) the DU at time t, ({{{{{rm{DU}}}}}}left(t|{theta }^{D}right)), is defined as the difference between the mean daily accumulation of rainfall over n1 days and a threshold rainfall level ((bar{D})):$${{{{{rm{DU}}}}}}left(t|{theta }^{D}right)={{{{{rm{max }}}}}}left{bar{D}-mathop{sum }limits_{n={n}_{2}}^{{n}_{2}+{n}_{1}-1}yleft(t-nright)/{n}_{1},0right},$$
(14)
where y(t) is the rainfall value at time t. The term max{x1, x2} is defined similarly as in Eq. 13.Logistic regression was performed using only the DU and using the product of CU and DU (CU × DU) as the explanatory variables and using the presence or absence of a first flowering event as the dependent variable for each phenological cluster. Because the number of phenological clusters is unknown, we performed forward selection on the cluster number based on the AIC. Let m be the number of phenological clusters based on the dendrogram drawn from the time-series clustering explained above (Supplementary Fig. 5). Given m phenological clusters, let ({G}_{k}^{m}) be the kth set of clusters in which the DU model is adopted for model fitting. Here, ({G}_{k}^{m}) indicates the set of cluster IDs, and k ranges from 0 to m(m+1)/2. For example, when m = 2 (i.e., there are two clusters, clusters 1 and 2), there are four cluster sets, calculated as follows:$${G}_{0}^{m=2}={},{G}_{1}^{m=2}={1},{G}_{2}^{m=2}={2},{G}_{3}^{m=2}={1,2},$$
(15)
where the element in the bracket indicates the ID of the cluster in which the DU model is adopted for model fitting. When k = 0, the DU model is not used; instead, the CU × DU model is adopted for model fitting for both clusters 1 and 2. Let i be the ith element of the vector E, which is defined as follows:$${{{{{bf{E}}}}}}={{t}_{1}^{1},,{t}_{2}^{1},…,,,{t}_{n}^{1},,…,,,{t}_{1}^{m},,{t}_{2}^{m},…,,{t}_{n}^{m}},$$
(16)
where n is the length of the time-series data for each cluster. Notably, n = 223 is the same for all species and clusters. The term ({t}_{1}^{m}) in the above equation denotes the first time point of the time series of length n for the species included in cluster m. Given m and k, let ({p}^{(m,k)}(i)) be the flowering probability of element i of vector E. The term ({p}^{(m,k)}(i)) is expressed as follows:$${{log }}left[frac{{p}^{left(m,kright)}left(iright)}{1-{p}^{left(m,kright)}left(iright)}right]= mathop{sum }limits_{j=1}^{m}{alpha }_{m,j}cdot {Z}_{m,j}left(iright)+mathop{sum }limits_{jin {G}_{k}^{left(mright)}}^{m}{beta }_{m,j}cdot {Z}_{m,j}left(iright)cdot {{{{{{rm{DU}}}}}}}_{m,j}left(i|{theta }_{j}^{D}right)\ +mathop{sum }limits_{jnotin {G}_{k}^{left(mright)}}^{m}{beta }_{m,j}cdot {Z}_{m,j}left(iright)cdot {{{{{rm{CU}}}}}}left(i|{theta }_{j}^{C}right)times {{{{{{rm{DU}}}}}}}_{m,j}left(i|{theta }_{j}^{D}right),$$
(17)
where ({Z}_{m,j}(i)) is the dummy variable indicating a cluster for i; ({Z}_{m,j}(i)) equals 1 if the ith element of E belongs to the jth cluster, otherwise it is zero, and ({alpha }_{m,j}) and ({beta }_{m,j}) in Eq. (5) are regression coefficients for the jth cluster when the species are grouped into m clusters. We estimate the parameters and the number of clusters based on a finite number of observations. Given the number of clusters m, for each of m clusters, the parameters were estimated by maximizing the loglikelihood value calculated for all combinations of potential parameter values for ({n}_{1},{n}_{2},bar{C},) and (bar{D}) within the ranges of [1 (min), 50 (max)] for n1, [1,50] for n2, [19,25] for (bar{C}), and [1,9] for (bar{D}). We varied the days (n1 and n2) by integers, temperature ((bar{C})) by tenths of a degree C, and daily precipitation ((bar{D})) by tenths of a mm. Regression coefficients (({alpha }_{m,j}), ({beta }_{m,j})) for all j values under a given m value and associated likelihoods were determined using generalized linear models with binomial error structures.With the results of the parameter estimations, we determined the number of clusters in two steps. For the first step, for a given m, we obtained (hat{k}(m)) according to the following equation:$$hat{k}(m)={arg }mathop{{min }}limits_{k}{{{{{{rm{AIC}}}}}}{m,k(m)},,k(m),=,0,,…,{2}^{m}}.$$
(18)
For the second step, with the results of (hat{k}) obtained from the first step, we obtained the estimate of the number of clusters according to forward selection by searching for the (hat{m}) value that satisfies the following inequalities:$${{{{{rm{AIC}}}}}}(hat{m},,hat{k}(hat{m})), < ,{{{{{rm{AIC}}}}}}(hat{m}+1,,hat{k}(hat{m}+1))cap {{{{{rm{AIC}}}}}}(hat{m},,hat{k}(hat{m})), < ,{{{{{rm{AIC}}}}}}(hat{m}-1,,hat{k}(hat{m}-1)).$$ (19) For model fitting, the first flowering month was extracted from the flowering phenology data. When flowering lasted more than 1 month, the month after the first flowering month was replaced by a value of zero (absence of flowering). If the month before the first flowering month was a missing value, the first flowering month was treated as a missing value and was not used for further analyses. We assumed that phenology monitoring was performed on the first date of each month.Projections of 21st-century changes in flowering phenologyWe used two scenarios (RCP2.6 and RCP8.5) to forecast future reproductive phenology in dipterocarp species for each of the three GCMs (GFDL–ESM2M, IPSL–CM5A-LR, and MIROC5). We predicted the flowering probability per month for each phenological cluster during the periods from 1 May 1976–31 March 1996 and from 1 January 2050–31 December 2099 based on the best model (Supplementary Table 2). The predicted flowering probability during the 2050–2099 period was normalized to that during the 1976–1996 period for each climate scenario and for each of three GCMs. To compare the seasonal patterns between 1976–1996 and 2050–2099, the predicted flowering probability was averaged for each month from January to December and plotted for each month in Fig. 6. R version 3.6.3 was used for all analyses.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

Beerda, B., Schilder, M. B. H., Van Hooff, J. A., De Vries, H. W. & Mol, J. A. Chronic stress in dogs subjected to social and spatial restriction I. Behavioral responses. Physiol. Behav. 66, 233–242 (1999).CAS 
PubMed 
Article 

Google Scholar 
Rooney, N. J., Gaines, S. A. & Bradshaw, J. W. Behavioural and glucocorticoid responses of dogs (Canis familiaris) to kennelling: investigating mitigation of stress by prior habituation. Physiol. Behav. 92, 847–854. https://doi.org/10.1016/j.physbeh.2007.06.011 (2007).CAS 
PubMed 
Article 

Google Scholar 
Stephen, J. M. & Ledger, R. A. A longitudinal evaluation of urinary cortisol in kennelled dogs Canis familiaris. Physiol. Behav. 87, 911–916. https://doi.org/10.1016/j.physbeh.2006.02.015 (2006).CAS 
PubMed 
Article 

Google Scholar 
Mills, D., Karagiannis, C., Zulch, H. Stress its effects on health and behavior. Vet. Clin. North Am. Small Anim. Pract. 44, 525–541 (2014).Mormède, P. et al. Exploration of the hypothalamic–pituitary–adrenal function as a tool to evaluate animal welfare. Physiol. Behav. 92, 317–339 (2007).PubMed 
Article 

Google Scholar 
Hennessy, M. B. Using hypothalamic–pituitary–adrenal measures for assessing and reducing the stress of dogs in shelters: A review. Appl. Anim. Behav. Sci. 149, 1–12 (2013).Article 

Google Scholar 
Cobb, M. L., Iskandarani, K., Chinchilli, V. M. & Dreschel, N. A. A systematic review and meta-analysis of salivary cortisol measurement in domestic canines. Domest. Anim. Endocrinol. 57, 31–42 (2016).CAS 
PubMed 
Article 

Google Scholar 
Wester, V. L. & van Rossum, E. F. Clinical applications of cortisol measurements in hair. Eur. J. Endocrinol. 173, M1–M10 (2015).CAS 
PubMed 
Article 

Google Scholar 
Heimbürge, S., Kanitz, E. & Otten, W. The use of hair cortisol for the assessment of stress in animals. Gen. Comp. Endocrinol. 270, 10–17 (2019).PubMed 
Article 

Google Scholar 
Meyer, J. S. & Novak, M. A. Minireview: hair cortisol: A novel biomarker of hypothalamic-pituitary-adrenocortical activity. Endocrinology 153, 4120–4127 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Khoury, J. E., Bosquet Enlow, M., Plamondon, A. & Lyons-Ruth, K. The association between adversity and hair cortisol levels in humans: A meta-analysis. Psychoneuroendocrinology 103, 104–117 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Davenport, M. D., Tiefenbacher, S., Lutz, C. K., Novak, M. A. & Meyer, J. S. Analysis of endogenous cortisol concentrations in the hair of rhesus macaques. Gen. Comp. Endocrinol. 147, 255–261 (2006).CAS 
PubMed 
Article 

Google Scholar 
Greff, M. J. E. et al. Hair cortisol analysis: An update on methodological considerations and clinical applications. Clin. Biochem. 63, 1–9 (2019).CAS 
PubMed 
Article 

Google Scholar 
del Rosario, G. et al. Effects of adrenocorticotropic hormone challenge and age on hair cortisol concentrations in dairy cattle. Can. J. Vet. Res. 75, 216–221 (2011).
Google Scholar 
Macbeth, B. J., Cattet, M., Stenhouse, G. B., Gibeau, M. L. & Janz, D. M. Hair cortisol concentration as a noninvasive measure of long-term stress in free-ranging grizzly bears (Ursus arctos): considerations with implications for other wildlife. Can. J. Zool. 88, 935–949 (2010).CAS 
Article 

Google Scholar 
Accorsi, P. A. et al. Cortisol determination in hair and faeces from domestic cats and dogs. Gen. Comp. Endocrinol. 155, 398–402 (2008).CAS 
PubMed 
Article 

Google Scholar 
Bennett, A. & Hayssen, V. Measuring cortisol in hair and saliva from dogs: coat color and pigment differences. Domest. Anim. Endocrinol. 39, 171–180 (2010).CAS 
PubMed 
Article 

Google Scholar 
Bryan, H. M., Adams, A. G., Invik, R. M., Wynne-Edwards, K. E. & Smits, J. E. Hair as a meaningful measure of baseline cortisol levels over time in dogs. J. Am. Assoc. Lab. Anim. Sci. 52, 189–196 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Siniscalchi, M., McFarlane, J. R., Kauter, K. G., Quaranta, A. & Rogers, L. J. Cortisol levels in hair reflect behavioural reactivity of dogs to acoustic stimuli. Res. Vet. Sci. 94, 49–54 (2013).CAS 
PubMed 
Article 

Google Scholar 
Stella, J., Shreyer, T., Ha, J. & Croney, C. Improving canine welfare in commercial breeding (CB) operations: Evaluating rehoming candidates. Appl. Anim. Behav. Sci. 220, 104861. https://doi.org/10.1016/j.applanim.2019.104861 (2019).Article 

Google Scholar 
Nicholson, S. L. & Meredith, J. E. Should stress management be part of the clinical care provided to chronically ill dogs?. J. Vet. Behav. 10, 489–495 (2015).Article 

Google Scholar 
Maxwell, N., Buchanan, C. & Evans, N. Hair cortisol concentrations, as a measure of chronic activity within the hypothalamic-pituitary-adrenal axis, is elevated in dogs farmed for meat, relative to pet dogs South Korea. Anim. Welf. 28, 389–395 (2019).Article 

Google Scholar 
Roth, L. S., Faresjö, Å, Theodorsson, E., Jensen, P. Hair cortisol varies with season and lifestyle and relates to human interactions in German shepherd dogs. Sci. Rep. 6, 19631; https://doi.org/10.1038/srep19631 (2016).Packer, R. M. et al. What can we learn from the hair of the dog? Complex effects of endogenous and exogenous stressors on canine hair cortisol. PLoS ONE 14, e0216000. https://doi.org/10.1371/journal.pone.0216000 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Sundman, A. et al. Long-term stress levels are synchronized in dogs and their owners. Sci. Rep. 9, 7391; https://doi.org/10.1038/s41598-019-43851-x (2019).Höglin, A. et al. Long-term stress in dogs is related to the human-dog relationship and personality traits. Sci. Rep. 11, 8612; https://doi.org/10.1038/s41598-021-88201-y (2021).Bowland, G. B. et al. Fur color and nutritional status predict hair cortisol concentrations of dogs in Nicaragua. Front. Vet. Sci. 7, 565346. https://doi.org/10.3389/fvets.2020.565346 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Veronesi, M. C. et al. Coat and claws as new matrices for noninvasive long-term cortisol assessment in dogs from birth up to 30 days of age. Theriogenology 84, 791–796 (2015).CAS 
PubMed 
Article 

Google Scholar 
Davenport, M. D., Lutz, C. K., Tiefenbacher, S., Novak, M. A. & Meyer, J. S. A rhesus monkey model of self-injury: Effects of relocation stress on behavior and neuroendocrine function. Biol. Psychiatry 63, 990–996 (2008).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
van der Laan, J. E., Vinke, C. M., van der Borg, J. A. M. & Arndt, S. S. Restless nights? Nocturnal activity as a useful indicator of adaptability of shelter housed dogs. Appl. Anim. Behav. Sci. 241, 105377. https://doi.org/10.1016/j.applanim.2021.105377 (2021).Article 

Google Scholar 
Pollinger, J. P. et al. Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication. Nature 464, 898–902 (2010).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Voith, V. L., Ingram, E., Mitsouras, K. & Irizarry, K. Comparison of adoption agency breed identification and DNA breed identification of dogs. J. Appl. Anim. Welf. Sci. 12, 253–262 (2009).CAS 
PubMed 
Article 

Google Scholar 
Gunter, L. M., Barber, R. T. & Wynne, C. D. L. A canine identity crisis: Genetic breed heritage testing of shelter dogs. PLoS ONE 13, e0202633. https://doi.org/10.1371/journal.pone.0202633 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pinheiro, J., Bates, D., DebRoy, S. & Sarkar, D., R Core Team. Nlme: linear and nonlinear mixed effects models. R package version 3. 1–148 (2020).Protopopova, A. & Gunter, L. Adoption and relinquishment interventions at the animal shelter: a review. Anim. Welf. 26, 35–48 (2017).Article 

Google Scholar 
Müntener, T., Doherr, M. G., Guscetti, F., Suter, M. M. & Welle, M. M. The canine hair cycle – a guide for the assessment of morphological and immunohistochemical criteria. Vet. Dermatol. 22, 383–395 (2011).PubMed 
Article 

Google Scholar 
Wennig, R. Potential problems with the interpretation of hair analysis results. Forensic Sci. Int. 107, 5–12 (2000).CAS 
PubMed 
Article 

Google Scholar 
Heimbürge, S., Kanitz, E., Tuchscherer, A. & Otten, W. Within a hair’s breadth – Factors influencing hair cortisol levels in pigs and cattle. Gen. Comp. Endocrinol. 288, 113359. https://doi.org/10.1016/j.ygcen.2019.113359 (2020).CAS 
PubMed 
Article 

Google Scholar 
Diaz, S. F., Torres, S. M., Dunstan, R. W. & Lekcharoensuk, C. An analysis of canine hair re-growth after clipping for a surgical procedure. Vet. Dermatol. 15, 25–30 (2004).PubMed 
Article 

Google Scholar 
Zeugswetter, F., Bydzovsky, N., Kampner, D. & Schwendenwein, I. Tailored reference limits for urine corticoid:creatinine ratio in dogs to answer distinct clinical questions. Vet. Rec. 167, 997–1001 (2010).CAS 
PubMed 
Article 

Google Scholar 
Jones, S. et al. Use of accelerometers to measure stress levels in shelter dogs. J. Appl. Anim. Welf. Sci. 17, 18–28 (2014).CAS 
PubMed 
Article 

Google Scholar 
Gunter, L. M., Feuerbacher, E. N., Gilchrist, R. J. & Wynne, C. D. Evaluating the effects of a temporary fostering program on shelter dog welfare. PeerJ 7, e6620. https://doi.org/10.7717/peerj.6620 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Van den Brom, W. E. & Biewenga, W. J. Assessment of glomerular filtration rate in normal dogs: analysis of the 51Cr-EDTA clearance and its relation to several endogenous parameters of glomerular filtration. Res. Vet. Sci. 30, 152–157 (1981).PubMed 
Article 

Google Scholar 
Sandri, M., Colussi, A., Perrotta, M. G. & Stefanon, B. Salivary cortisol concentration in healthy dogs is affected by size, sex, and housing context. J. Vet. Behav. 10, 302–306 (2015).Article 

Google Scholar 
Haase, C. G., Long, A. K. & Gillooly, J. F. Energetics of stress: linking plasma cortisol levels to metabolic rate in mammals. Biol. Lett. 12, 20150867. https://doi.org/10.1098/rsbl.2015.0867 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Garnier, F., Benoit, E., Virat, M., Ochoa, R. & Delatour, P. Adrenal cortical response in clinically normal dogs before and after adaptation to a housing environment. Lab. Anim. 24, 40–43 (1990).CAS 
PubMed 
Article 

Google Scholar 
Beerda, B. et al. Chronic stress in dogs subjected to social and spatial restriction. II. Hormonal and immunological responses. Physiol. Behav. 66, 243–254 (1999).Rincón-Cortés, M., Herman, J. P., Lupien, S., Maguire, J. & Shansky, R. M. Stress: Influence of sex, reproductive status and gender. Neurobiol. Stress 10, 100155. https://doi.org/10.1016/j.ynstr.2019.100155 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Oyola, M. G. & Handa, R. J. Hypothalamic–pituitary–adrenal and hypothalamic–pituitary–gonadal axes: sex differences in regulation of stress responsivity. Stress 20, 476–494 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Willen, R. M., Mutwill, A., MacDonald, L. J., Schiml, P. A. & Hennessy, M. B. Factors determining the effects of human interaction on the cortisol levels of shelter dogs. Appl. Anim. Behav. Sci. 186, 41–48 (2017).Article 

Google Scholar 
Protopopova, A. Effects of sheltering on physiology, immune function, behavior, and the welfare of dogs. Physiol. Behav. 159, 95–103 (2016).CAS 
PubMed 
Article 

Google Scholar 
Mesarcova, L., Kottferova, J., Skurkova, L., Leskova, L. & Kmecova, N. Analysis of cortisol in dog hair-a potential biomarker of chronic stress: a review. Vet. Med. (Praha) 62, 363–376 (2017).CAS 
Article 

Google Scholar 
Neumann, A. et al. Predicting hair cortisol levels with hair pigmentation genes: a possible hair pigmentation bias. Sci. Rep. 7, 8529 (2017).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Romero, L. M. & Beattie, U. K. Common myths of glucocorticoid function in ecology and conservation. J. Exp. Zool. A. Ecol. Integr. Physiol. https://doi.org/10.1002/jez.2459 (2021).PubMed 
Article 

Google Scholar 
Heimbürge, S., Kanitz, E., Tuchscherer, A. & Otten, W. Is it getting in the hair? – Cortisol concentrations in native, regrown and segmented hairs of cattle and pigs after repeated ACTH administrations. Gen. Comp. Endocrinol. 295, 113534. https://doi.org/10.1016/j.ygcen.2020.113534 (2020).CAS 
PubMed 
Article 

Google Scholar 
Van Ockenburg, S. L. et al. The relationship between 63 days of 24-h urinary free cortisol and hair cortisol levels in 10 healthy individuals. Psychoneuroendocrinology 73, 142–147 (2016).PubMed 
Article 

Google Scholar 
Short, S. J. et al. Correspondence between hair cortisol concentrations and 30-day integrated daily salivary and weekly urinary cortisol measures. Psychoneuroendocrinology 71, 12–18 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Mack, Z. & Fokidis, H. B. A novel method for assessing chronic cortisol concentrations in dogs using the nail as a source. Domest. Anim. Endocrinol. 59, 53–57 (2017).CAS 
PubMed 
Article 

Google Scholar  More

Bach, L. T. et al. Testing the climate intervention potential of ocean afforestation using the Great Atlantic Sargassum Belt. Nat. Commun. 12, 2556 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
N‘Yeurt, A. D. R., Chynoweth, D. P., Capron, M. E., Stewart, J. R. & Hasan, M. A. Negative carbon via ocean afforestation. Process Saf. Environ. Prot. 90, 467–474 (2012).Article 
CAS 

Google Scholar 
Duarte, C. M., Bruhn, A. & Krause-Jensen, D. A seaweed aquaculture imperative to meet global sustainability targets. Nat. Sustain. 5, 185–193 (2022).Article 

Google Scholar 
Woody, T. Seaweed ‘forests’ can help fight climate change. National Geographic https://www.nationalgeographic.co.uk/environment-and-conservation/2019/08/seaweed-forests-can-help-fight-climate-change (2019).Godin, M. The ocean farmers trying to save the world with seaweed. Time https://time.com/5848994/seaweed-climate-change-solution/ (2020).Marshall, M. Kelp is coming: how seaweed could prevent catastrophic climate change. New Scientist https://www.newscientist.com/article/mg24632821-100-kelp-is-coming-how-seaweed-could-prevent-catastrophic-climate-change/ (2020).Bever, F. ‘Run the oil industry in reverse’: fighting climate change by farming kelp. NPR https://www.npr.org/2021/03/01/970670565/run-the-oil-industry-in-reverse-fighting-climate-change-by-farming-kelp (2021).Running Tide. https://www.runningtide.com/ (2022).IPCC: Summary for Policymakers. In Global Warming of 1.5 °C (eds Masson-Delmotte, V. et al.) (WMO, 2018).IPCC: Summary for Policymakers. In Climate Change 2021: The Physical Science Basis (eds Masson-Delmotte, V. et al.) (Cambridge Univ. Press) (in the press).GESAMP. High Level Review of a Wide Range of Proposed Marine Geoengineering Techniques (eds Boyd, P. W. & Vivian, C. M. G.) GESAMP Working Group 41 (International Maritime Organization, 2019).Boyd, P. & Vivian, C. Should we fertilize oceans or seed clouds? No one knows. Nature 570, 155–157 (2019).CAS 
PubMed 
Article 

Google Scholar 
Law, C. S. Predicting and monitoring the impact of large-scale iron fertilisation on marine trace gas emissions. Mar. Ecol. Prog. Ser. 364, 283–288 (2008).CAS 
Article 

Google Scholar 
Russell, L. M. et al. Ecosystem impacts of geoengineering: a review for developing a science plan. Ambio 41, 350–369 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Costello, C., Fries, L. & Gaines, S. Transformational opportunities in ocean-based food & nutrition. Zenodo https://zenodo.org/record/4646319#.YkBFxhPMLAw (2021).Jouffray, J.-B., Blasiak, R., Norström, A. V., Österblom, H. & Nyström, M. The blue acceleration: the trajectory of human expansion into the ocean. One Earth 2, 43–54 (2020).Article 

Google Scholar 
Cullen, J. J. & Boyd, P. W. Predicting and verifying the intended and uninterested consequence of large-scale iron fertilization. Mar. Ecol. Prog. Ser. 364, 295–301 (2008).CAS 
Article 

Google Scholar 
Bach, L. T., Gill, S. J., Rickaby, R. E. M., Gore, S. & Renforth, P. CO2 removal with enhanced weathering and ocean alkalinity enhancement: potential risks and co-benefits for marine pelagic ecosystems. Front. Clim. https://doi.org/10.3389/fclim.2019.00007 (2019).Moore, C. M. et al. Processes and patterns of oceanic nutrient limitation. Nat. Geosci. 6, 701–710 (2013).CAS 
Article 

Google Scholar 
Suchet, P. A., Probst, J.-L. & Ludwig, L. Worldwide distribution of continental rock lithology: implications for the atmospheric/soil CO2 uptake by continental weathering and alkalinity river transport to the oceans. Glob. Biogeochem. Cycles 17, 1038 (2003).
Google Scholar 
Macreadie, P. I. et al. The future of blue carbon science. Nat. Commun. 10, 3998 (2019).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Fraser, C. I., Nikula, R. & Waters, J. M. Oceanic rafting by a coastal community. Proc. Biol. Sci. 278, 649–655 (2011).PubMed 

Google Scholar 
Fraser, C. I., Davies, I. D., Bryant, D. & Waters, J. M. How disturbance and dispersal influence intraspecific structure. J. Ecol. 106, 1298–1306 (2018).Article 

Google Scholar 
Fraser, C. I. et al. Antarctica’s ecological isolation will be broken by storm-driven dispersal and warming. Nat. Clim. Change 8, 704–708 (2018).Article 

Google Scholar 
Chung, I. K., Beardall, J., Mehta, S., Sahoo, D. & Stojkovic, S. Using marine macroalgae for carbon sequestration: a critical appraisal. J. Appl. Phycol. 23, 877–886 (2011).CAS 
Article 

Google Scholar 
Krause-Jensen, D. & Duarte, C. M. Substantial role of macroalgae in marine carbon sequestration. Nat. Geosci. 9, 737–742 (2016).CAS 
Article 

Google Scholar 
Hurd, C. L. et al. Forensic carbon accounting: assessing the role of seaweeds for carbon sequestration. J. Phycol., https://doi.org/10.1111/jpy.13249 (2022).Stripe commits $8M to six new carbon removal companies. Stripe https://stripe.com/newsroom/news/spring-21-carbon-removal-purchases (2021).General application. Stripe https://github.com/stripe/carbon-removal-source-materials/blob/master/Project%20Applications/Spring2021/Running%20Tide%20-%20Stripe%20Spring21%20CDR%20Purchase%20Application.pdf (2021).Coston-Clements, L. Utilization of the Sargassum Habitat by Marine Invertebrates and Vertebrates: a Review. NOAA Technical Memorandum NMFS-SEFSC, 296 (U.S. Department of Commerce, National Oceanic and Atmospheric Administration, National Marine Fisheries Service, Southeast Fisheries Science Center & Beaufort Laboratory, 1991).Egan, S. et al. The seaweed holobiont: understanding seaweed–bacteria interactions. FEMS Microbiol. Rev. 37, 462–476 (2013).CAS 
PubMed 
Article 

Google Scholar 
Califano, G., Kwantes, M., Abreu, M. H., Costa, R. & Wichard, T. Cultivating the macroalgal holobiont: effects of integrated multi-trophic aquaculture on the microbiome of Ulva rigida (Chlorophyta)Front. Mar. Sci. 7, 52 (2020).Article 

Google Scholar 
Selvarajan, R. et al. Distribution, interaction and functional profiles of epiphytic bacterial communities from the rocky intertidal seaweeds, South Africa. Sci. Rep. 9, 19835 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bonthond, G. et al. The role of host promiscuity in the invasion process of a seaweed holobiont. ISME J. 15, 1668–1679 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Wang, M. et al. The great Atlantic Sargassum belt. Science 365, 83–87 (2019).CAS 
PubMed 
Article 

Google Scholar 
Johns, E. M. et al. The establishment of a pelagic Sargassum population in the tropical Atlantic: biological consequences of a basin-scale long distance dispersal event. Prog. Oceanogr. 182, 102269 (2020).Article 

Google Scholar 
Martiny, A. C. et al. Biogeochemical controls of surface ocean phosphate. Sci. Adv. 5, eaax0341 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zehr, J. P. & Capone, D. G. Changing perspectives in marine nitrogen fixation. Science 368, eaay9514 (2020).CAS 
PubMed 
Article 

Google Scholar 
Harrison, P. J., Druehl, L. D., Lloyd, K. E. & Thompson, P. A. Nitrogen uptake kinetics in three year-classes of Laminaria groenlandica (Laminariales: Phaeophyta). Mar. Biol. 93, 29–35 (1986).CAS 
Article 

Google Scholar 
Hurd, C. L. & Dring, M. L. Phosphate uptake by intertidal algae in relation to zonation and season. Mar. Biol. 107, 281–289 (1990).Article 

Google Scholar 
Ohtake, M. et al. Growth and nutrient uptake characteristics of Sargassum macrocarpum cultivated with phosphorus-replete wastewater. Aquat. Bot. 163, 103208 (2020).Article 

Google Scholar 
MacFarlane, J. J. & Raven, J. A. C, N and P nutrition of Lemanea mamillosa Kütz. (Batrachospermales, Rhodophyta) in the Dighty Burn, Angus, U.K. Plant Cell Environ. 13, 1–13 (1990).CAS 
Article 

Google Scholar 
Wu, J., Keller, D. P. & Oschlies, A. Carbon dioxide removal via macroalgae open-ocean mariculture and sinking: an Earth system modeling study. Preprint at Earth System Dynamics Discuss https://doi.org/10.5194/esd-2021-104 (2022).Kwiatkowski, L. et al. Twenty-first century ocean warming, acidification, deoxygenation, and upper-ocean nutrient and primary production decline from CMIP6 model projections. Biogeosciences 17, 3439–3470 (2020).CAS 
Article 

Google Scholar 
Chapman, A. R. O. & Craigie, J. S. Seasonal growth in Laminaria longicruris: relations with dissolved inorganic nutrients and internal reserves of nitrogen. Mar. Biol. 40, 197–205 (1977).CAS 
Article 

Google Scholar 
Dutkiewicz, S., Scott, J. R. & Follows, M. J. Winners and losers: ecological and biogeochemical changes in a warming ocean. Glob. Biogeochem. Cycles 27, 463–477 (2013).CAS 
Article 

Google Scholar 
Thomas, M. K. et al. Temperature–nutrient interactions exacerbate sensitivity to warming in phytoplankton. Glob. Change Biol. 2, 3269–3280 (2017).Article 

Google Scholar 
Lapointe, B. E. et al. Nutrient content and stoichiometry of pelagic Sargassum reflects increasing nitrogen availability in the Atlantic Basin. Nat. Commun. 12, 3060 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Fan, W. et al. A sea trial of enhancing carbon removal from Chinese coastal waters by stimulating seaweed cultivation through artificial upwelling. Appl. Ocean Res. 101, 102260 (2020).Article 

Google Scholar 
Karl, D. M. & Letelier, R. M. Nitrogen fixation-enhanced carbon sequestration in low nitrate, low chlorophyll seascapes. Mar. Ecol. Prog. Ser. 364, 257–268 (2008).CAS 
Article 

Google Scholar 
Oschlies, A. S., Pahlow, M., Yool, A. & Matear, R. Climate engineering by artificial ocean upwelling: channelling the sorcerer’s apprentice. Geophys. Res. Lett. 37, L04701 (2010).Article 
CAS 

Google Scholar 
Thornton, D. C. O. Dissolved organic matter (DOM) release by phytoplankton in the contemporary and future ocean. Eur. J. Phycol. 49, 20–46 (2014).CAS 
Article 

Google Scholar 
Morán, X. A. G., Sebastián, M., Pedrós-Alió, C. & Estrada, M. Response of Southern Ocean phytoplankton and bacterioplankton production to short-term experimental warming. Limnol. Oceanogr. 51, 1791–1800 (2006).Article 

Google Scholar 
Marañón, E., Cermeño, P., Fernández, E., Rodríguez, J. & Zabala, L. Significance and mechanisms of photosynthetic production of dissolved organic carbon in a coastal eutrophic ecosystem. Limnol. Oceanogr. 49, 1652–1666 (2004).Article 

Google Scholar 
Paine, E. R., Schmid, M., Boyd, P. W., Diaz-Pulido, G. & Hurd, C. L. Rate and fate of dissolved organic carbon release by seaweeds: a missing link in the coastal ocean carbon cycle. J. Phycol. 57, 1375–1391 (2021).CAS 
PubMed 
Article 

Google Scholar 
Brylinsky, M. Release of dissolved organic matter by some marine macrophytes. Mar. Biol. 39, 213–220 (1977).Article 

Google Scholar 
Sieburth, J. M. Studies on algal substances in the sea. III. The production of extracellular organic matter by littoral marine algae. J. Exp. Mar. Biol. Ecol. 3, 290–309 (1969).CAS 
Article 

Google Scholar 
Hanson, R. B. Pelagic Sargassum community metabolism: carbon and nitrogen. J. Exp. Mar. Biol. Ecol. 29, 107–118 (1977).CAS 
Article 

Google Scholar 
Zark, M., Riebesell, U. & Dittmar, T. Effects of ocean acidification on marine dissolved organic matter are not detectable over the succession of phytoplankton blooms. Sci. Adv. 1, e1500531 (2015).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Zhang, Y., Liu, X., Wang, M. & Qin, B. Compositional differences of chromophoric dissolved organic matter derived from phytoplankton and macrophytes. Org. Geochem. 55, 26–37 (2013).Article 
CAS 

Google Scholar 
Hulatt, C. J., Thomas, D. N., Bowers, D. G., Norman, L. & Zhang, C. Exudation and decomposition of chromophoric dissolved organic matter (CDOM) from some temperate macroalgae. Estuar. Coast. Shelf Sci. 84, 147–153 (2009).CAS 
Article 

Google Scholar 
Liu, S., Trevathan-Tackett, S. M., Ewers Lewis, C. J., Huang, X. & Macreadie, P. I. Macroalgal blooms trigger the breakdown of seagrass blue carbon. Environ. Sci. Technol. 54, 14750–14760 (2020).CAS 
PubMed 
Article 

Google Scholar 
Vieira, H. C. et al. Ocean warming may enhance biochemical alterations induced by an invasive seaweed exudate in the mussel Mytilus galloprovincialis. Toxics 9, 121 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Brooks, S. D. & Thornton, D. C. O. Marine aerosols and clouds. Ann. Rev. Mar. Sci. 10, 289–313 (2018).PubMed 
Article 

Google Scholar 
Lewis, M. R., Carr, M.-E., Feldman, G. C., Esaias, W. & McClain, C. Influence of penetrating solar radiation on the heat budget of the equatorial Pacific Ocean. Nature 347, 543–545 (1990).Article 

Google Scholar 
Morel, A. Optical modeling of the upper ocean in relation to its biogenous matter content (case-I waters). J. Geophys. Res. 93, 10749–10768 (1988).Article 

Google Scholar 
Park, J.-Y., Kug, J.-S., Bader, J., Rolph, R. & Kwon, M. Amplified Arctic warming by phytoplankton under greenhouse warming. Proc. Natl Acad. Sci. USA 112, 5921–5926 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Denaro, G. et al. Dynamics of two picophytoplankton groups in Mediterranean Sea: analysis of the deep chlorophyll maximum by a stochastic advection-reaction-diffusion model. PLoS ONE 8, e66765 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kavanaugh, M. T. et al. Experimental assessment of the effects of shade on an intertidal kelp: do phytoplankton blooms inhibit growth of open-coast macroalgae? Limnol. Oceanogr. 54, 276–288 (2009).Article 

Google Scholar 
Omand, M. M., Steinberg, D. K. & Stamies, K. Cloud shadows drive vertical migrations of deep-dwelling marine life. Proc. Natl Acad. Sci. USA 118, e2022977118 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bach, L. T. & Boyd, P. W. Seeking natural analogs to fast-forward the assessment of marine CO2 removal. Proc. Natl Acad. Sci. USA 118, e2106147118 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
van Donk, E. & van de Bund, W. J. Impact of submerged macrophytes including charophytes on phyto- and zooplankton communities: allelopathy versus other mechanisms. Aquat. Bot. 72, 261–274 (2002).Article 

Google Scholar 
Jin, Q., Dong, S. & Wang, C. Allelopathic growth inhibition of Prorocentrum micans (Dinophyta) by Ulva pertusa and Ulva linza (Chlorophyta) in laboratory cultures. Eur. J. Phycol. 40, 31–37 (2005).Article 

Google Scholar 
Wallace, R. B. & Gobler, C. J.Factors controlling blooms of microalgae and macroalgae (Ulva rigida) in a eutrophic, urban estuary: Jamaica Bay, NY, USA. Estuaries Coast 38, 519–533 (2015).CAS 
Article 

Google Scholar 
Tang, Y. Z. & Gobler, C. J. The green macroalga, Ulva lactuca, inhibits the growth of seven common harmful algal bloom species via allelopathy. Harmful Algae 10, 480–488 (2011).Article 

Google Scholar 
Cagle, S. E., Roelke, D. L. & Muhl, R. W. Allelopathy and micropredation paradigms reconcile with system stoichiometry. Ecosphere 12, e03372 (2021).Article 

Google Scholar 
Hein, M., Pedersen, M. F. & Sand-Jensen, K. Size-dependent nitrogen uptake in micro- and macroalgae. Mar. Ecol. Prog. Ser. 118, 247–253 (1995).Article 

Google Scholar 
Stevens, C. L., Hurd, C. L. & Smith, M. J. Water motion relative to subtidal kelp fronds. Limnol. Oceanogr. 46, 668–678 (2001).Article 

Google Scholar 
Raut, Y., Morando, M. & Capone, D. G. Diazotrophic macroalgal associations with living and decomposing Sargassum. Front. Microbiol. 9, 3127 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Villareal, T. A., Woods, S., Moore, J. K. & CulverRymsza, K. Vertical migration of Rhizosolenia mats and their significance to NO3− fluxes in the central North Pacific gyre. J. Plankton Res. 18, 1103–1121 (1996).Article 

Google Scholar 
Gachon, C. M. M., Sime-Ngando, T., Strittmatter, M., Chambouvet, A. & Kim, G. H. Algal diseases: spotlight on a black box. Trends Plant Sci. 15, 633–640 (2010).CAS 
PubMed 
Article 

Google Scholar 
Sánchez-Baracaldo, P., Bianchini, G., Wilson, J. D. & Knoll, A. H. Cyanobacteria and biogeochemical cycles through Earth history. Trends Microbiol. 30, 143–157 (2022).PubMed 
Article 
CAS 

Google Scholar 
Thiel, M. & Gutow, L. in Oceanography and Marine Biology: an Annual Review Vol. 43 (eds Gibson, R. et al.) 279–418 (Taylor & Francis, 2005).Rech, S., Borrell Pichs, Y. J. & García-Vazquez, E. Anthropogenic marine litter composition in coastal areas may be a predictor of potentially invasive rafting fauna. PLoS ONE 13, e0191859 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Food and Agriculture Organization (FAO) of the United Nations. The State of World Fisheries and Aquaculture 2020: Sustainability in Action (FAO, 2020).Schell, J. M., Goodwin, D. S. & Siuda, A. N. S. Recent Sargassum inundation events in the Caribbean: shipboard observations reveal dominance of a previously rare form. Oceanography 28, 8–10 (2015).Article 

Google Scholar 
Rodríguez-Martínez, R. E. et al. Element concentrations in pelagic Sargassum along the Mexican Caribbean coast in 2018–2019. Peer J. 8, e8667 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Flannery, T. How farming giant seaweed can feed fish and fix the climate. The Conversation Trust https://theconversation.com/how-farming-giant-seaweed-can-feed-fish-and-fix-the-climate-81761 (2017).GESAMP. Methodology for the Evaluation of Ballast Water Management Systems Using Active Substances. GESAMP No. 101 (eds Linders, J. & Dock, A.) (International Maritime Organization, 2019).Lenton, A., Boyd, P. W., Thatcher, M. & Emmerson, K. M. Foresight must guide geoengineering research and development. Nat. Clim. Change 9, 342 (2019).Article 

Google Scholar 
Sumaila, U. R. Financing a sustainable ocean economy. Nat. Commun. 12, 3259 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Rockström, J. et al. Planetary boundaries: exploring the safe operating space for humanity. Ecol. Soc. 14, 32 (2009).Article 

Google Scholar 
Rech, S., Salmina, S., Borrell Pichs, Y. J. & García-Vazquez, E. Dispersal of alien invasive species on anthropogenic litter from European mariculture areas. Mar. Pollut. Bull. 131, 10–16 (2018).CAS 
PubMed 
Article 

Google Scholar 
Therriault, T. W. et al. The invasion risk of species associated with Japanese tsunami marine debris in Pacific North America and Hawaii. Mar. Pollut. Bull. 132, 82–89 (2018).CAS 
PubMed 
Article 

Google Scholar 
Miller, J. A., Carlton, J. T., Chapman, J. W., Geller, J. B. & Ruiz, G. M. Transoceanic dispersal of the mussel Mytilus galloprovincialis on Japanese tsunami marine debris: an approach for evaluating rafting of a coastal species at sea. Mar. Pollut. Bull. 132, 60–69 (2018).CAS 
PubMed 
Article 

Google Scholar 
Carlton, J. T. et al. Tsunami-driven rafting: transoceanic species dispersal and implications for marine biogeography. Science 357, 1402–1406 (2017).CAS 
PubMed 
Article 

Google Scholar 
Hunt, G. L. Jr et al. Advection in polar and sub-polar environments: impacts on high latitude marine ecosystems. Prog. Oceanogr. 149, 40–81 (2016).Article 

Google Scholar 
Hallegraeff, G. M. & Bolch, C. J. Transport of dinoflagellate cysts in ship’s ballast water: implications for plankton biogeography and aquaculture. J. Plankton Res. 14, 1067–1084 (1992).Article 

Google Scholar 
Russell, L. K., Hepburn, C. D., Hurd, C. L. & Stuart, M. D. The expanding range of Undaria pinnatifida in southern New Zealand: distribution, dispersal mechanisms and the invasion of wave-exposed environments. Biol. Invasions 10, 103–115 (2008).Article 

Google Scholar 
Uwai, S. et al. Genetic diversity in Undaria pinnatifida (Laminariales, Phaeophyceae) deduced from mitochondria genes—origins and succession of introduced populations. Phycologia 45, 687–695 (2006).Article 

Google Scholar  More

Notomi, T. et al. Loop-mediated isothermal amplification of DNA. Nucleic Acids Res. https://doi.org/10.1093/nar/28.12.e63 (2000).Article 

Google Scholar 
Nagamine, K., Hase, T. & Notomi, T. Accelerated reaction by loop-mediated isothermal amplification using loop primers. Mol. Cell. Probes 16, 223–229. https://doi.org/10.1006/mcpr.2002.0415 (2002).CAS 
Article 

Google Scholar 
Nagamine, K., Watanabe, K., Ohtsuka, K., Hase, T. & Notomi, T. Loop-mediated isothermal amplification reaction using a nondenatured template. Clin. Chem. 47, 1742–1743 (2001).CAS 
Article 

Google Scholar 
Thai, H. T. C. et al. Development and evaluation of a novel loop-mediated isothermal amplification method for rapid detection of severe acute respiratory syndrome coronavirus. J. Clin. Microbiol. 42, 1956–1961. https://doi.org/10.1128/jcm.42.5.1956-1961.2004 (2004).CAS 
Article 

Google Scholar 
Geojith, G., Dhanasekaran, S., Chandran, S. P. & Kenneth, J. Efficacy of loop mediated isothermal amplification (LAMP) assay for the laboratory identification of Mycobacterium tuberculosis isolates in a resource limited setting. J. Microbiol. Methods 84, 71–73. https://doi.org/10.1016/j.mimet.2010.10.015 (2011).CAS 
Article 

Google Scholar 
Saengsawang, N. et al. Development of a fluorescent distance-based paper device using loop-mediated isothermal amplification to detect Escherichia coli in urine. Analyst 145, 8077–8086. https://doi.org/10.1039/d0an01306d (2020).CAS 
Article 

Google Scholar 
Yoshikawa, R. et al. Development and evaluation of a rapid and simple diagnostic assay for COVID-19 based on loop-mediated isothermal amplification. Plos Neglect. Trop. Dis. 14, 14. https://doi.org/10.1371/journal.pntd.000885 (2021).Article 

Google Scholar 
Kim, J. et al. Development and evaluation of a multiplex loop-mediated isothermal amplification (LAMP) assay for differentiation of Mycobacterium tuberculosis and non-tuberculosis mycobacterium in clinical samples. PLoS ONE 16, 11. https://doi.org/10.1371/journal.pone.0244753 (2021).CAS 
Article 

Google Scholar 
Hongjaisee, S. et al. Rapid visual detection of hepatitis C virus using a reverse transcription loop-mediated isothermal ampli fi cation assay. Int. J. Infect. Dis. 102, 440–445. https://doi.org/10.1016/j.ijid.2020.10.082 (2021).CAS 
Article 

Google Scholar 
Niessen, L. & Vogel, R. F. Detection of Fusarium graminearum DNA using a loop-mediated isothermal amplification (LAMP) assay. Int. J. Food Microbiol. 140, 183–191. https://doi.org/10.1016/j.ijfoodmicro.2010.03.036 (2010).CAS 
Article 

Google Scholar 
Ren, W. C., Liu, N. & Li, B. H. Development and application of a LAMP method for rapid detection of apple blotch caused by Marssonina coronaria. Crop Prot. 141, 6. https://doi.org/10.1016/j.cropro.2020.105452 (2021).CAS 
Article 

Google Scholar 
Kong, G. H. et al. Detection of Peronophythora litchii on lychee by loop-mediated isothermal amplification assay. Crop Prot. 139, 6. https://doi.org/10.1016/j.cropro.2020.105370 (2021).CAS 
Article 

Google Scholar 
Zhou, Q. J. et al. Simultaneous detection of multiple bacterial and viral aquatic pathogens using a fluorogenic loop-mediated isothermal amplification-based dual-sample microfluidic chip. J. Fish Dis. https://doi.org/10.1111/jfd.13325 (2020).Article 

Google Scholar 
Huang, H. L. et al. Molecular method for rapid detection of the red tide dinoflagellate Karenia mikimotoi in the coastal region of Xiangshan Bay, China. J. Microbiol. Methods 168, 7. https://doi.org/10.1016/j.mimet.2019.105801 (2020).CAS 
Article 

Google Scholar 
Sridapan, T. et al. Rapid detection of Clostridium perfringens in food by loop-mediated isothermal amplification combined with a lateral flow biosensor. PLoS ONE 16, 14. https://doi.org/10.1371/journal.pone.0245144 (2021).CAS 
Article 

Google Scholar 
Xiong, X. et al. Using real time fluorescence loop-mediated isothermal amplification for rapid species authentication of Atlantic salmon (Salmo salar). J. Food Compos. Anal. 95, 7. https://doi.org/10.1016/j.jfca.2020.103659 (2021).CAS 
Article 

Google Scholar 
Huang, C. G., Hsu, J. C., Haymer, D. S., Lin, G. C. & Wu, W. J. Rapid identification of the Mediterranean fruit fly (Diptera: Tephritidae) by loop-mediated isothermal amplification. J. Econ. Entomol. 102, 1239–1246 (2009).CAS 
Article 

Google Scholar 
Ide, T., Kanzaki, N., Ohmura, W. & Okabe, K. Molecular identification of an invasive wood-boring insect Lyctus brunneus (Coleoptera: Bostrichidae: Lyctinae) using frass by loop-mediated isothermal amplification and nested PCR assays. J. Econ. Entomol. 109, 1410–1414. https://doi.org/10.1093/jee/tow030 (2016).CAS 
Article 

Google Scholar 
Stainton, K., Hall, J., Budge, G. E., Boonham, N. & Hodgetts, J. Rapid molecular methods for in-field and laboratory identification of the yellow-legged Asian hornet (Vespa velutina nigrithorax). J. Appl. Entomol. 142, 610–616. https://doi.org/10.1111/jen.12506 (2018).CAS 
Article 

Google Scholar 
Agarwal, A., Cunningham, J. P., Valenzuela, I. & Blacket, M. J. A diagnostic LAMP assay for the destructive grapevine insect pest, phylloxera (Daktulosphaira vitifoliae). Sci. Rep. 10, 10. https://doi.org/10.1038/s41598-020-77928-9 (2020).CAS 
Article 

Google Scholar 
Rizzo, D. et al. Molecular identification of Anoplophora glabripennis (Coleoptera: Cerambycidae) from frass by loop-mediated isothermal amplification. J. Econ. Entomol. 113, 2911–2919. https://doi.org/10.1093/jee/toaa206 (2020).CAS 
Article 

Google Scholar 
Hsieh, C. H., Wang, H. Y., Chen, Y. F. & Ko, C. C. Loop-mediated isothermal amplification for rapid identification of biotypes B and Q of the globally invasive pest Bemisia tabaci, and studying population dynamics. Pest Manag. Sci. 68, 1206–1213. https://doi.org/10.1002/ps.3298 (2012).CAS 
Article 

Google Scholar 
Williams, M. R. et al. Isothermal amplification of environmental DNA (eDNA) for direct field-based monitoring and laboratory confirmation of Dreissena sp. PLoS ONE 12, 18. https://doi.org/10.1371/journal.pone.0186462 (2017).CAS 
Article 

Google Scholar 
Ponting, S., Tomkies, V. & Stainton, K. Rapid identification of the invasive small hive beetle (Aethina tumida) using LAMP. Pest Manag. Sci. 77, 1476–1481. https://doi.org/10.1002/ps.6168 (2020).CAS 
Article 

Google Scholar 
Davis, C. N. et al. Rapid detection of Galba truncatula in water sources on pasture-land using loop-mediated isothermal amplification for control of trematode infections. Parasites Vectors 13, 11. https://doi.org/10.1186/s13071-020-04371-0 (2020).CAS 
Article 

Google Scholar 
Carvalho, J. et al. Faster monitoring of the invasive alien species (IAS) Dreissena polymorpha in river basins through isothermal amplification. Sci. Rep. 11, 10. https://doi.org/10.1038/s41598-021-89574-w (2021).CAS 
Article 

Google Scholar 
Treguier, A. et al. Environmental DNA surveillance for invertebrate species: Advantages and technical limitations to detect invasive crayfish Procambarus clarkii in freshwater ponds. J. Appl. Ecol. 51, 871–879. https://doi.org/10.1111/1365-2664.12262 (2014).CAS 
Article 

Google Scholar 
Cai, W. et al. Using eDNA to detect the distribution and density of invasive crayfish in the Honghe-Hani rice terrace World Heritage site. PLoS ONE https://doi.org/10.1371/journal.pone.0177724 (2017).Article 

Google Scholar 
Wilcox, T. M. et al. Understanding environmental DNA detection probabilities: A case study using a stream-dwelling char Salvelinus fontinalis. Biol. Conserv. 194, 209–216. https://doi.org/10.1016/j.biocon.2015.12.023 (2016).Article 

Google Scholar 
Hunter, M. E., Ferrante, J. A., Meigs-Friend, G. & Ulmer, A. Improving eDNA yield and inhibitor reduction through increased water volumes and multi-filter isolation techniques. Sci. Rep. https://doi.org/10.1038/s41598-019-40977-w (2019).Article 

Google Scholar 
Twardochleb, L. A., Olden, J. D. & Larson, E. R. A global meta-analysis of the ecological impacts of nonnative crayfish. Freshw. Sci. 32, 1367–1382. https://doi.org/10.1899/12-203.1 (2013).Article 

Google Scholar 
Andruszkiewicz, A. E., Zhang, W. G. & Govindarajan, A. F. Environmental DNA shedding and decay rates from diverse animal forms and thermal regimes. Environ. DNA 3, 492–514. https://doi.org/10.1002/edn3.141 (2021).Article 

Google Scholar 
Stedtfeld, R. D. et al. Static self-directed sample dispensing into a series of reaction wells on a microfluidic card for parallel genetic detection of microbial pathogens. Biomed. Microdev. 17, 89. https://doi.org/10.1007/s10544-015-9994-1 (2015).CAS 
Article 

Google Scholar 
Koloren, Z., Sotiriadou, I. & Karanis, P. Investigations and comparative detection of Cryptosporidium species by microscopy, nested PCR and LAMP in water supplies of Ordu, Middle Black Sea, Turkey. Ann. Trop. Med. Parasitol. 105, 607–615. https://doi.org/10.1179/2047773211y.0000000011 (2011).CAS 
Article 

Google Scholar 
Sabike, I. I. et al. Use of direct LAMP screening of broiler fecal samples for Campylobacter jejuni and Campylobacter coli in the positive flock identification strategy. Front. Microbiol. 7, 1582. https://doi.org/10.3389/fmicb.2016.01582 (2016).Article 

Google Scholar 
Gahlawat, S. K., Ellis, A. E. & Collet, B. A sensitive loop-mediated isothermal amplification (LAMP) method for detection of Renibacterium salmoninarum, causative agent of bacterial kidney disease in salmonids. J. Fish Dis. 32, 491–497. https://doi.org/10.1111/j.1365-2761.2009.01005.x (2009).CAS 
Article 

Google Scholar 
Levy, J. et al. Methods for rapid and effective PCR-based detection of ‘Candidatus Liberibacter solanacearum’ from the insect vector Bactericera cockerelli: Streamlining the DNA extraction/purification process. J. Econ. Entomol. 106, 1440–1445. https://doi.org/10.1603/ec12419 (2013).CAS 
Article 

Google Scholar 
Kaneko, H., Kawana, T., Fukushima, E. & Suzutani, T. Tolerance of loop-mediated isothermal amplification to a culture medium and biological substances. J. Biochem. Biophys. Methods 70, 499–501. https://doi.org/10.1016/j.jbbm.2006.08.008 (2007).CAS 
Article 

Google Scholar 
Curtis, A. N., Tiemann, J. S., Douglass, S. A., Davis, M. A. & Larson, E. R. High stream flows dilute environmental DNA (eDNA) concentrations and reduce detectability. Divers. Distrib. 27, 1918–1931. https://doi.org/10.1111/ddi.13196 (2020).Article 

Google Scholar 
Mauvisseau, Q. et al. Environmental DNA as an efficient tool for detecting invasive crayfishes in freshwater ponds. Hydrobiologia 805, 163–175. https://doi.org/10.1007/s10750-017-3288-y (2018).CAS 
Article 

Google Scholar 
RStudioTeam. Boston (ed. PBC) (2020).Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer, 2016).Book 

Google Scholar  More