Ecology

1.Addo-Bediako, A., Chown, S. L. & Gaston, K. J. Thermal tolerance, climatic variability and latitude. Proc. R. Soc. B. 267, 739–745 (2000).CAS 
PubMed 
Article 

Google Scholar 
2.Andersen, J. L. et al. How to assess Drosophila cold tolerance: chill coma temperature and lower lethal temperature are the best predictors of cold distribution limits. Funct. Ecol. 29, 55–65 (2015).Article 

Google Scholar 
3.Kimura, M. T. Cold and heat tolerance of drosophilid flies with reference to their latitudinal distributions. Oecologia 140, 442–449 (2004).ADS 
PubMed 
Article 
PubMed Central 

Google Scholar 
4.Gaston, K. J. & Chown, S. L. Elevation and climatic tolerance: A test using dung beetles. Oikos 86, 584–590 (1999).Article 

Google Scholar 
5.MacMillan, H. A. Dissecting cause from consequence: a systematic approach to thermal limits. J. Exp. Biol. 222, jeb191593 (2019).6.Overgaard, J. & MacMillan, H. A. The integrative physiology of insect chill tolerance. Annu. Rev. Physiol. 79, 187–208 (2017).CAS 
PubMed 
Article 

Google Scholar 
7.Armstrong, G. A. B., Rodríguez, E. C. & Meldrum Robertson, R. Cold hardening modulates K+ homeostasis in the brain of Drosophila melanogaster during chill coma. J. Insect Physiol. 58, 1511–1516 (2012).8.Rodgers, C. I., Armstrong, G. A. B. & Robertson, R. M. Coma in response to environmental stress in the locust: a model for cortical spreading depression. J. Insect Physiol. 56, 980–990 (2010).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
9.Andersen, M. K. & Overgaard, J. The central nervous system and muscular system play different roles for chill coma onset and recovery in insects. Comp. Biochem. Physiol. Part A Mol. Integr. Physiol. 233, 10–16 (2019).10.Koštál, V., Vambera, J. & Bastl, J. On the nature of pre-freeze mortality in insects: Water balance, ion homeostasis and energy charge in the adults of Pyrrhocoris apterus. J. Exp. Biol. 207, 1509–1521 (2004).PubMed 
Article 

Google Scholar 
11.Zachariassen, K. E., Kristiansen, E. & Pedersen, S. A. Inorganic ions in cold-hardiness. Cryobiology 48, 126–133 (2004).CAS 
PubMed 
Article 

Google Scholar 
12.MacMillan, H. A. & Sinclair, B. J. The role of the gut in insect chilling injury: Cold-induced disruption of osmoregulation in the fall field cricket, Gryllus pennsylvanicus. J. Exp. Biol. 214, 726–734 (2011).CAS 
PubMed 
Article 

Google Scholar 
13.MacMillan, H. A., Williams, C. M., Staples, J. F. & Sinclair, B. J. Reestablishment of ion homeostasis during chill-coma recovery in the cricket Gryllus pennsylvanicus. PNAS 109, 20750–20755 (2012).ADS 
CAS 
PubMed 
Article 

Google Scholar 
14.MacMillan, H. A., Findsen, A., Pedersen, T. H. & Overgaard, J. Cold-induced depolarization of insect muscle: Differing roles of extracellular K+ during acute and chronic chilling. J. Exp. Biol. 217, 2930–2938 (2014).PubMed 
PubMed Central 

Google Scholar 
15.Bayley, J. S., Sørensen, J. G., Moos, M., Koštál, V. & Overgaard, J. Cold-acclimation increases depolarization resistance and tolerance in muscle fibers from a chill-susceptible insect, Locusta migratoria. Am. J. Physiol. Regul. Integr. Comp. Physiol. 319, R439–R447 (2020).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
16.Bayley, J. S. et al. Cold exposure causes cell death by depolarization-mediated Ca2+ overload in a chill-susceptible insect. PNAS 115, E9737–E9744 (2018).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
17.Carrington, J., Andersen, M. K., Brzezinski, K. & MacMillan, H. A. Hyperkalemia, not apoptosis, accurately predicts chilling injury in individual locusts. Proc. R. Soc. B. (in press).18.Koštál, V., Yanagimoto, M. & Bastl, J. Chilling-injury and disturbance of ion homeostasis in the coxal muscle of the tropical cockroach (Nauphoeta cinerea). Comp. Biochem. Physiol. B. Biochem. Mol. Biol. 143, 171–179 (2006).19.MacMillan, H. A., Baatrup, E. & Overgaard, J. Concurrent effects of cold and hyperkalaemia cause insect chilling injury. Proc. R. Soc. B. 282 (2015).20.Garcia, M. J., Littler, A. S., Sriram, A. & Teets, N. M. Distinct cold tolerance traits independently vary across genotypes in Drosophila melanogaster. Evolution 74, 1437–1450 (2020).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
21.Gerken, A. R., Mackay, T. F. C. & Morgan, T. J. Artificial selection on chill-coma recovery time in Drosophila melanogaster: Direct and correlated responses to selection. J. Therm. Biol. 59, 77–85 (2016).PubMed 
Article 
PubMed Central 

Google Scholar 
22.Colinet, H. & Hoffmann, A. A. Comparing phenotypic effects and molecular correlates of developmental, gradual and rapid cold acclimation responses in Drosophila melanogaster. Funct. Ecol. 26, 84–93 (2012).Article 

Google Scholar 
23.MacMillan, H. A., Andersen, J. L., Loeschcke, V. & Overgaard, J. Sodium distribution predicts the chill tolerance of Drosophila melanogaster raised in different thermal conditions. Am. J. Physiol. Regul. Integr. Comp. Physiol. 308, R823–R831 (2015).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
24.Ransberry, V. E., MacMillan, H. A. & Sinclair, B. J. The relationship between chill-coma onset and recovery at the extremes of the thermal window of Drosophila melanogaster. Physiol. Biochem. Zool. 84, 553–559 (2011).PubMed 
Article 
PubMed Central 

Google Scholar 
25.Sunday, J. M., Bates, A. E. & Dulvy, N. K. Global analysis of thermal tolerance and latitude in ectotherms. Proc. R. Soc. B. 278, 1823–1830 (2011).PubMed 
Article 
PubMed Central 

Google Scholar 
26.Hoffmann, A. A., Anderson, A. & Hallas, R. Opposing clines for high and low temperature resistance in Drosophila melanogaster. Ecol. Lett. 5, 614–618 (2002).Article 

Google Scholar 
27.Hoffmann, A. A., Shirriffs, J. & Scott, M. Relative importance of plastic vs genetic factors in adaptive differentiation: Geographical variation for stress resistance in Drosophila melanogaster from eastern Australia. Funct. Ecol. 19, 222–227 (2005).Article 

Google Scholar 
28.Overgaard, J., Hoffmann, A. A. & Kristensen, T. N. Assessing population and environmental effects on thermal resistance in Drosophila melanogaster using ecologically relevant assays. J. Therm. Biol. 36, 409–416 (2011).Article 

Google Scholar 
29.Ayrinhac, A. et al. Cold adaptation in geographical populations of Drosophila melanogaster: Phenotypic plasticity is more important than genetic variability. Funct. Ecol. 18, 700–706 (2004).Article 

Google Scholar 
30.Gibert, P. & Huey, R. B. Chill-coma temperature in Drosophila: Effects of developmental temperature, latitude, and phylogeny. Physiol. Biochem. Zool. 74, 429–434 (2001).CAS 
PubMed 
Article 

Google Scholar 
31.Hori, Y. & Kimura, M. T. Relationship between cold stupor and cold tolerance in Drosophila (Diptera: Drosophilidae). Environ. Entomol. 27, 1297–1302 (1998).Article 

Google Scholar 
32.Teets, N. M. & Hahn, D. A. Genetic variation in the shape of cold-survival curves in a single fly population suggests potential for selection from climate variability. J. Evol. Biol. 31, 543–555 (2018).CAS 
PubMed 
Article 

Google Scholar 
33.Kellermann, V. et al. Phylogenetic constraints in key functional traits behind species’ climate niches: Patterns of desiccation and cold resistance across 95 Drosophila species. Evolution 66, 3377–3389 (2012).PubMed 
Article 

Google Scholar 
34.Pool, J. E., Braun, D. T. & Lack, J. B. Parallel evolution of cold tolerance within Drosophila melanogaster. Mol. Biol. Evol. 34, 349–360 (2017).CAS 
PubMed 

Google Scholar 
35.Mansourian, S. et al. Wild African Drosophila melanogaster are seasonal specialists on marula fruit. Curr. Biol. 28, 3960-3968.e3 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
36.Pool, J. E. et al. Population genomics of Sub-Saharan Drosophila melanogaster: African diversity and non-African admixture. PLoS Genetics 8, e1003080 (2012).37.Baudry, E., Viginier, B. & Veuille, M. Non-African populations of Drosophila melanogaster have a unique origin. Mol. Biol. Evol. 21, 1482–1491 (2004).CAS 
PubMed 
Article 

Google Scholar 
38.MacMillan, H. A., Andersen, J. L., Davies, S. A. & Overgaard, J. The capacity to maintain ion and water homeostasis underlies interspecific variation in Drosophila cold tolerance. Sci. Rep. 5, 18607 (2015).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
39.Chen, C.-P. & Walker, V. K. Cold-shock and chilling tolerance in Drosophila. J. Insect Physiol. 40, 661–669 (1994).Article 

Google Scholar 
40.Hoffmann, A. A. & Watson, M. Geographical variation in the acclimation responses of Drosophila to temperature extremes. Am. Nat. 142, S93–S113 (1993).PubMed 
Article 
PubMed Central 

Google Scholar 
41.Ørsted, M., Hoffmann, A. A., Rohde, P. D., Sørensen, P. & Kristensen, T. N. Strong impact of thermal environment on the quantitative genetic basis of a key stress tolerance trait. Heredity 122, 315–325 (2019).PubMed 
Article 
PubMed Central 

Google Scholar 
42.Gerken, A. R., Eller, O. C., Hahn, D. A. & Morgan, T. J. Constraints, independence, and evolution of thermal plasticity: probing genetic architecture of long- and short-term thermal acclimation. PNAS 112, 4399–4404 (2015).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
43.Nyamukondiwa, C., Terblanche, J. S., Marshall, K. E. & Sinclair, B. J. Basal cold but not heat tolerance constrains plasticity among Drosophila species (Diptera: Drosophilidae). J. Evol. Biol. 24, 1927–1938 (2011).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
44.van Heerwaarden, B. & Kellermann, V. Does plasticity trade off with basal heat tolerance?. Trends Ecol. Evol. 35, 874–885 (2020).PubMed 
Article 
PubMed Central 

Google Scholar 
45.Gilchrist, G. W., Huey, R. B. & Partridge, L. Thermal sensitivity of Drosophila melanogaster: evolutionary responses of adults and eggs to laboratory natural selection at different temperatures. Physiol. Zool. 70, 403–414 (1997).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
46.Maclean, H. J., Kristensen, T. N., Sørensen, J. G. & Overgaard, J. Laboratory maintenance does not alter ecological and physiological patterns among species: A Drosophila case study. J. Evol. Biol. 31, 530–542 (2018).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
47.Henry, Y., Renault, D. & Colinet, H. Hormesis-like effect of mild larval crowding on thermotolerance in Drosophila flies. J. Exp. Biol. 221, jeb169342 (2018).48.Nilson, T. L., Sinclair, B. J. & Roberts, S. P. The effects of carbon dioxide anesthesia and anoxia on rapid cold-hardening and chill coma recovery in Drosophila melanogaster. J. Insect Physiol. 52, 1027–1033 (2006).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
49.Hazell, S. P. & Bale, J. S. Low temperature thresholds: are chill coma and CTmin synonymous?. J. Insect Physiol. 57, 1085–1089 (2011).CAS 
PubMed 
Article 

Google Scholar 
50.Bertram, G. C. L. The low temperature limit of activity of arctic insects. J. Anim. Ecol. 4, 35–42 (1935).Article 

Google Scholar 
51.Sinclair, B. J., Coello Alvarado, L. E. & Ferguson, L. V. An invitation to measure insect cold tolerance: Methods, approaches, and workflow. J. Therm. Biol. 53, 180–197 (2015).52.MacMillan, H. A. et al. Anti-diuretic activity of a CAPA neuropeptide can compromise Drosophila chill tolerance. J. Exp. Biol. 221, jeb185884 (2018).53.R Core Team. R: A Language and Environment for Statistical Computing. (R Foundation for Statistical Computing, 2020). More

1.Hoegh-Guldberg, O. & Bruno, J. F. The impact of climate change on the world’s marine ecosystems. Science 328, 1523–1528 (2010).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
2.Rosenberg, E., Kellogg, C. A. & Rohwer, F. L. Coral microbiology. Oceanography 20, 146 (2007).Article 

Google Scholar 
3.Rohwer, F. & Youle, M. Coral Reefs in the Microbial Seas (Plaid Press, 2010).4.Haas, A. F. et al. Global microbialization of coral reefs. Nat. Microbiol. 1, 16042 (2016).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
5.Roth, E., Jeon, K. & Stacey, G. Homology in endosymbiotic systems: the term ‘symbiosome’. (1988).6.Zvuloni, A. et al. Spatio-temporal transmission patterns of black-band disease in a coral community. PLoS ONE 4, e4993 (2009).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
7.Peters, E. C. in Coral Reefs in the Anthropocene 147–178 (Springer Press, 2015).8.Kushmaro, A., Rosenberg, E., Fine, M., Ben Haim, Y. & Loya, Y. Effect of temperature on bleaching of the coral Oculina patagonica by Vibrio AK-1. Mar. Ecol. Prog. Ser. 171, 131–137 (1998).Article 

Google Scholar 
9.Ben-Haim, Y. et al. Vibrio coralliilyticus sp. nov., a temperature-dependent pathogen of the coral Pocillopora damicornis. Int. J. Syst. Evol. Microbiol. 53, 309–315 (2003).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
10.Rosenberg, E., Kushmaro, A., Kramarsky-Winter, E., Banin, E. & Yossi, L. The role of microorganisms in coral bleaching. ISME J. 3, 139–146 (2009).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
11.Kramarsky-Winter, E., Downs, C., Downs, A. & Loya, Y. Cellular responses in the coral Stylophora pistillata exposed to eutrophication from fish mariculture. Evol. Ecol. Res. 11, 381–401 (2009).
Google Scholar 
12.Santos, E. et al. Genomic and proteomic analyses of the coral pathogen Vibrio coralliilyticus reveal a diverse virulence repertoire. ISME J. 5, 1471–1483 (2011).Article 

Google Scholar 
13.Ben Haim, Y. & Rosenberg, E. A novel Vibrio sp. pathogen of the coral Pocillopora damicornis. Mar. Biol. 141, 47–55 (2002).14.Vidal-Dupiol, J. et al. Coral bleaching under thermal stress: putative involvement of host/symbiont recognition mechanisms. BMC Physiol. 9, 14 (2009).15.Zvuloni, A., Artzy-Randrup, Y., Katriel, G., Loya, Y. & Stone, L. Modeling the impact of white-plague coral disease in climate change scenarios. PLoS Comput. Biol. 11, e1004151 (2015).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
16.Wright, R. M. et al. Intraspecific differences in molecular stress responses and coral pathobiome contribute to mortality under bacterial challenge in Acropora millepora. Sci. Rep. 7, 2609-017-02685-1 (2017).
Google Scholar 
17.Garren, M., Son, K., Tout, J., Seymour, J. R. & Stocker, R. Temperature-induced behavioral switches in a bacterial coral pathogen. ISME J. 10, 1363–1372 (2016).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
18.Sussman, M. et al. Vibrio zinc-metalloprotease causes photoinactivation of coral endosymbionts and coral tissue lesions. PLoS ONE 4, e4511 (2009).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
19.Rosenberg, E. & Kushmaro, A. in Coral Reefs: An Ecosystem in Transition 451–464 (Springer, 2011).20.Bourne, D. et al. Microbial disease and the coral holobiont. Trends Microbiol. 17, 554–562 (2009).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
21.Pollock, F. J., Morris, P. J., Willis, B. L. & Bourne, D. G. The urgent need for robust coral disease diagnostics. PLoS Pathog. 7, e1002183 (2011).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
22.Weis, V., Davy, S., Hoegh-Guldberg, O., Rodriguez-Lanetty, M. & Pringle, J. Cell biology in model systems as the key to understanding corals. Trends Ecol. Evol. 23, 369–376 (2008).PubMed 
Article 
PubMed Central 

Google Scholar 
23.Work, T. & Meteyer, C. To understand coral disease, look at coral cells. EcoHealth 11, 610–618 (2014).24.Shapiro, O. H., Kramarsky-Winter, E., Gavish, A. R., Stocker, R. & Vardi, A. A coral-on-a-chip microfluidic platform enabling live-imaging microscopy of reef-building corals. Nat. Commun. 7, 10860 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
25.Dunn, A. K., Millikan, D. S., Adin, D. M., Bose, J. L. & Stabb, E. V. New rfp- and pES213-derived tools for analyzing symbiotic Vibrio fischeri reveal patterns of infection and lux expression in situ. Appl. Environ. Microbiol. 72, 802–810 (2006).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
26.Sammarco, P. W. Polyp bail-out: an escape response to environmental stress and a new means of reproduction in corals. Mar. Ecol. Prog. Ser. Oldendorf 10, 57–65 (1982).Article 

Google Scholar 
27.Work, T. M. & Aeby, G. S. Pathology of tissue loss (white syndrome) in Acropora sp. corals from the Central Pacific. J. Invertebr. Pathol. 107, 127–131 (2011).PubMed 
Article 
PubMed Central 

Google Scholar 
28.Work, T. M. & Aeby, G. S. Systematically describing gross lesions in corals. Dis. Aquat. Org. 70, 155–160 (2006).Article 

Google Scholar 
29.Ainsworth, T., Fine, M., Roff, G. & Hoegh-Guldberg, O. Bacteria are not the primary cause of bleaching in the Mediterranean coral Oculina patagonica. ISME J. 2, 67–73 (2008).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
30.Ainsworth, T. D., Fine, M., Blackall, L. L. & Hoegh-Guldberg, O. Fluorescence in situ hybridization and spectral imaging of coral-associated bacterial communities. Appl. Environ. Microbiol. 72, 3016–3020 (2006).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
31.Boyett, H. V., Bourne, D. G. & Willis, B. L. Elevated temperature and light enhance progression and spread of black band disease on staghorn corals of the Great Barrier Reef. Mar. Biol. 151, 1711–1720 (2007).Article 

Google Scholar 
32.Gignoux-Wolfsohn, S., Marks, C. J. & Vollmer, S. V. White band disease transmission in the threatened coral, Acropora cervicornis. Sci. Rep. 2, 804 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
33.Kaczmarsky, L. T. Coral disease dynamics in the central Philippines. Dis. Aquat. Org. 69, 9–21 (2006).Article 

Google Scholar 
34.Horridge, G. A. The co-ordination of the protective retraction of coral polyps. Philos. Trans. R. Soc. Lond. B Biol. Sci. 240, 495–528 (1957).Article 

Google Scholar 
35.Katz, S. M., Pollock, F. J., Bourne, D. G. & Willis, B. L. Crown-of-thorns starfish predation and physical injuries promote brown band disease on corals. Coral Reefs 33, 705–716 (2014).Article 

Google Scholar 
36.Shapiro, O. H. et al. Vortical ciliary flows actively enhance mass transport in reef corals. Proc. Natl Acad. Sci. USA 111, 13391–13396 (2014).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
37.Gladfelter, E. Circulation of fluids in the gastrovascular system of the reef coral Acropora cervicornis. Biol. Bull. 165, 619–625 (1983).PubMed 
Article 
PubMed Central 

Google Scholar 
38.Patterson, M. R. A chemical engineering view of cnidarian symbioses. Am. Zool. 32, 566–582 (1992).Article 

Google Scholar 
39.Lewis, J. & Price, W. Feeding mechanisms and feeding strategies of Atlantic reef corals. J. Zool. 176, 527–544 (1975).Article 

Google Scholar 
40.Brown, B. E. & Bythell, J. C. Perspectives on mucus secretion in reef corals. Mar. Ecol. Prog. Ser. 296, 291–309 (2005).CAS 
Article 

Google Scholar 
41.Zetsche, E., Baussant, T., Meysman, F. J. & van Oevelen, D. Direct visualization of mucus production by the cold-water coral Lophelia pertusa with digital holographic microscopy. PLoS ONE 11, e0146766 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
42.Lam, E., Kato, N. & Lawton, M. Programmed cell death, mitochondria and the plant hypersensitive response. Nature 411, 848–853 (2001).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
43.Chuang, P. S. & Mitarai, S. Signaling pathways in the coral polyp bailout response. Coral Reefs 39, 532–534 (2020).Article 

Google Scholar 
44.Kvitt, H. et al. Breakdown of coral colonial form under reduced pH conditions is initiated in polyps and mediated through apoptosis. Proc. Natl Acad. Sci. USA 112, 2082–2086 (2015).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
45.Garren, M. et al. A bacterial pathogen uses dimethylsulfoniopropionate as a cue to target heat-stressed corals. ISME J. 8, 999–1007 (2014).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
46.Garren, M. & Azam, F. Corals shed bacteria as a potential mechanism of resilience to organic matter enrichment. ISME J. 6, 1159–1165 (2012).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
47.Banin, E., Israely, T., Fine, M., Loya, Y. & Rosenberg, E. Role of endosymbiotic zooxanthellae and coral mucus in the adhesion of the coral-bleaching pathogen Vibrio shiloi to its host. FEMS Microbiol. Lett. 199, 33–37 (2001).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
48.Meron, D. et al. Role of flagella in virulence of the coral pathogen Vibrio coralliilyticus. Appl. Environ. Microbiol. 75, 5704–5707 (2009).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
49.Certner, R. H., Dwyer, A. M., Patterson, M. R. & Vollmer, S. V. Zooplankton as a potential vector for white band disease transmission in the endangered coral, Acropora cervicornis. PeerJ 5, e3502 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
50.Ushijima, B. et al. Vibrio coralliilyticus strain OCN008 is an etiological agent of acute Montipora white syndrome. Appl. Environ. Microbiol. 80, 2102–2109 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
51.Butler, S. M. & Camilli, A. Going against the grain: chemotaxis and infection in Vibrio cholerae. Nat. Rev. Microbiol. 3, 611 (2005).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
52.Stocker, R. The 100 µm length scale in the microbial ocean. Aquat. Microb. Ecol. 76, 189–194 (2015).Article 

Google Scholar 
53.Csaszar, N. B. M., Seneca, F. O. & van Oppen, M. J. H. Variation in antioxidant gene expression in the scleractinian coral Acropora millepora under laboratory thermal stress. Mar. Ecol. Prog. Ser. 392, 93–102 (2009).CAS 
Article 

Google Scholar 
54.Agostini, S. et al. Biological and chemical characteristics of the coral gastric cavity. Coral Reefs 31, 147–156 (2012).Article 

Google Scholar 
55.Houlbrèque, F., Rodolfo‐Metalpa, R. & Ferrier‐Pagès, C. Heterotrophic nutrition of tropical, temperate and deep‐sea corals. Dis. Coral 150–163 (2015).56.Sorokin, Y. I. Trophical role of bacteria in the ecosystem of the coral reef. Nature 242, 415–417 (1973).Article 

Google Scholar  More

The contexts of FinBIF’s emergence and developmentBuild-up of impetus and consensusA gradually strengthening drive from different sectors of society led to the eventual realisation of FinBIF as a comprehensively integrated BDDI. The earliest push came from a working group convened by the Ministry of Education and commissioned to analyse and plan the future of natural history museums (NHMs) in Finland. It identified as a goal to establish a national central NHM, one task of which would be to “act as an information centre maintaining an Automated Data Processing (ADP) based national register of collections”26. A subsequent working group commissioned to plan the central NHM concluded that one of the most urgent of the new nation-wide functions was “to create an ADP-based central registry to serve all Finnish natural history museums”27.While the central NHM was established in 1988, the proposed inclusive, national database did not materialise, but separate actors (institutions, projects, research groups) developed a variety of information systems independently. When Finland joined GBIF in 2001, an increased enthusiasm on BD informatics emerged, but significant national funding for GBIF-related activities was not mobilised. However, the demand for national progress grew rapidly with the international development of BD informatics. A project on developing the efficiency of nature conservation made recommendations on IT development, and on collating and opening data through a national solution28. The National Red List29 included recommendations on developing databases for species observations. An action plan on developing species conservation recommended a separate project for overall development of biodiversity data management, which was already close to what FinBIF later came to be30.In 2012, the Finnish Museum of natural History ‘Luomus’ successfully lobbied for a national, cross-sectoral meeting to discuss the way forward in building a national BDDI. This kick-off meeting was convened by the Ministry of the Environment on 7 September 2012. All major organisations that use, hold and/or produce biodiversity data, or fund such activities, were represented (Table 3). Subsequently, the urgency of establishing the national BDDI was re-emphasised in Finland’s Biodiversity Strategy and Action Plan31.Table 3 The organisations represented at the meeting where the establishment of a national BDDI was agreed upon on 7 September 2012.Full size tableCross-sectoral co-creation, collaboration and fundingThe host institution of FinBIF is the Finnish Museum of Natural History’Luomus’, which has led its development and the acquisition of funding, and coordinated the national and international cooperation. However, a wide national collaborative network has been an integral part of FinBIF from the start (Fig. 2), and many organisations from all sectors of society have participated in the development of FinBIF’s services and in data mobilisation. Luomus has convened seven advisory and co-creation groups in which c. 140 specialists have shared their expertise on a voluntary basis to advise policy and service development (Table 4). The wide collaborative network and the all-inclusive business model of FinBIF has helped in attracting funding for its development from many different sources with different funding criteria, such as national and EU-level research infrastructure funding, national funds for developing governance and administration, EU-funding for nature conservation actions, funds for Nordic collaboration in e-infrastructure development, private foundation grants to research and citizen science, and collaborative projects with private companies.Table 4 The expert advisory groups convened by FinBIF for co-creation of content, services, and policies. CMS = collection management system; OMS = observations management system; API = application programming interface; GIS = geographical information system; INSPIRE = infrastructure for spatial information in Europe.Full size tableIT ArchitectureFinBIF’s architecture has four key goals:

(1)

To harvest occurrence data from Finland and surrounding areas and share them as one interoperable data mass

(2)

To provide a customisable system for recording species observations

(3)

To provide a collection management system, and digitisation technologies and workflows for Finnish natural history collections

(4)

To compile, maintain, and share master data, including a taxonomic backbone of Finnish species.

These goals are accomplished by using Service-Oriented Architecture (SOA)32 and by providing over twenty individual background services that interact with each other, mostly using RESTful HTTP APIs (for explanation of REST, see33). These include: (1) A triplestore RDF database34,35 and an API for master data (taxonomy, schema vocabularies, people, collection metadata, image metadata, data download metadata) and collection specimen data; (2) a JavaScript Object Notation (JSON36) storage database and an API of primary occurrence data; (3) a Data Warehouse database, and an extract-transform-load (ETL) process and an API built on top of an HP Vertica database; (4) an Elasticsearch search engine (see https://www.elastic.co/what-is/elasticsearch) and an API for performance-critical access of primary occurrence data and taxonomy data; (5) multimedia services (images, audio storage and conversions); and (6) map services.FinBIF is provided as a service to fulfil the needs of the Finnish biodiversity data community and to share Finnish data to all interested parties. The architecture and complete software of FinBIF cannot be easily adopted by other countries or institutions. However, FinBIF maintains several software libraries that are specifically designed to be reusable by anyone; such as the specimen label generator library, and the powerful form generation tools of the Notebook OMS. Data and services can be used by third party applications via FinBIF’s public API (https://api.laji.fi).Occurrence data architectureFinBIF maintains systems for both primary and secondary occurrence data storage (Fig. 4). A primary data storage is a database or other type of storage where the data are maintained and used for the original purposes and use cases. A secondary data storage holds a copy of the primary data for distribution, analysis, and other derived purposes.Fig. 4Services, processes and data flows of FinBIF (a). The relations of services and processes to data life cycle phases as depicted in Fig. 1. Dashed lines denote planned services (b). Occurrence data flow from FinBIF’s and external actors’ primary sources through FinBIF’s data warehouse to FinBIF’s portal. CMS = Collection Management System, OMS = Observations Management System, DW = data warehouse, ETL = Extract-Transform-Load, API = Application Programming Interface, GIS = Geographical Information System, IAS = Invasive Alien Species, AI = Artificial Intelligence.Full size imageFinBIF has two IT systems for storing and maintaining primary occurrence data: the Kotka CMS for specimen data and the Notebook OMS for observational data (see below for descriptions). FinBIF also co-operates with iNaturalist and maintains iNaturalist Finland.FinBIF’s Data Warehouse harvests occurrence data from primary data sources and stores them as secondary data. The data warehouse transforms primary occurrence data into a single interoperable data mass through various ETL processes. Currently, live updates are received from 34 different data sources using eight different ETL processes. Many more datasets have been loaded using a one-time process. For primary data sources that are not databases or IT systems (for example Excel, MapInfo files), a tool is provided for harmonizing the data into FinBIF format and uploading them to the data warehouse as secondary data. The primary data remain under the ownership and management of the source, where all changes and updates are made. However, FinBIF maintains an annotation system, which stores added information about the occurrence entries (as primary data), for example identifications and quality markings by taxon experts. These annotations are stamped on top of the original data in the data warehouse to provide an enriched version of the data.Each dataset must have metadata that describe, e.g., the name, type and owner of the dataset. In the Data Warehouse, primary occurrence data, annotations, dataset metadata, the taxonomic backbone, and information on locations and people are linked to provide an enriched query service. FinBIF is a visible platform to make data FAIR16, thus increasing the prestige of the shared dataset, but the data sharer also gains through receiving quality feedback about the occurrences from the many taxon experts and volunteers that annotate the data in FinBIF. Based on the annotations, the data can be improved in the primary data source. To make data more interoperable, FinBIF encourages primary data source owners to harmonise their data with FinBIF-maintained master data, by using the same taxonomy and schema vocabularies.The Data Warehouse has two sides: public and private. The public side contains open-access data, some of which have been coarsened because of the sensitivity of the species in question or because of a research embargo, as specified in FinBIF’s data policy (https://laji.fi/en/about/2982; https://laji.fi/en/about/875). In addition, individual observers may hide, e.g., the exact location of their occurrence or their name. The restricted-access side of the Data Warehouse contains uncoarsened data. Some data sources provide a limited version to the public side and a full version to the restricted side. Government officials can access the restricted-use, uncoarsened data via a separate authorities’ portal. Researchers and other users can do so by issuing a data request (see below). Maintaining two versions of the data has required developing unique designs. They do not reveal information about sensitive species, which might endanger their preservation, but still allow joining annotations made on the restricted side so these can be shown also on the public side. Occurrences that need to be coarsened are detached from their original concept and uploaded with a random delay to the public side to make it difficult to discover the exact location from accompanying occurrences. Apart from this delay, the private and public sides are automatically fully synchronised.API and GIS services have been built on top of the Data Warehouse to allow open-access use of the data, both for the public and the restricted side. Most of the occurrence data from Finland go to GBIF via FinBIF. This data transfer is still being developed, and all datasets are not yet automatically copied to GBIF.Master data management architectureA triplestore is used in FinBIF for all small datasets, including taxon data. More specifically, the data are stored according to the RDF35 specification. An RDF Schema defines the allowed properties for each class. FinBIF’s triplestore34 implementation is an Oracle relational database with two tables (resource and statement), which provides the ability to do Structured Query Language (SQL) queries and updates. Doing small, atomic updates is easy, as only a small subset of the triplets can be updated instead of the entire data entity. Maintaining a complete record of history comes without much effort, as it can be done on an individual triplet level.As an example, the FinBIF taxon data model – including adjacent classes such as publication, person, image, and threat assessments – consists of 260 properties. If the data model were stored in a normalized relational database, there would be an estimated 56 tables, which could be difficult to maintain. Thus, in FinBIF, non-relational database solutions are preferred.IdentifiersFinBIF uses a persistent HTTP-URI identifier for all types of real-life and digital objects (specimens, occurrences, taxa, metadata, persons, organisations, information systems, etc.), as recommended by the World Wide Web Consortium37. The identifier takes the user to an ID redirect service, which redirects the user to a page that shows information about the object in human-readable format. For example, specimen identifiers redirect to information about the specimen and taxon identifiers to a page describing the taxon.The redirect service can also provide machine-readable data about the object, if the user (client software) requests that using Accept headers. Supported formats vary based on data types, e.g., for specimens, the system can offer data in RDF + XML38 and JSON-LD (see https://json-ld.org/) formats using CETAF compliant vocabulary (CETAF Specimen Preview Profile CSPP39). This is also compatible with MIDS (Minimum Information about a Digital Specimen40).If partner organisations do not provide HTTP-URI identifiers for their occurrences, FinBIF will use the persistent internal IDs of the data source to generate globally unique URI identifiers. DOI identifiers for data downloads and dataset metadata will be created in the near future.IT solutions in key services and processesFinBIF runs numerous processes to provide a rich set of services (Fig. 4). The IT solutions applied in building key services and in enabling central processes are described below. The order of the descriptions follows the data life cycle phases identified in Fig. 1.Kotka CMSKotka is one of the two primary data management systems of FinBIF. It is designed to fit the needs of different types of collections and can be further adapted when new needs arise.Kotka differs in many ways from traditional CMS solutions. It applies simple and pragmatic approaches. This has helped it grow into a nationally used system despite limited development resources – on average less than one full-time equivalent developer. The aim is to improve collection management efficiency by providing practical tools. Kotka emphasises the quantity of digitised specimens over completeness of the data. It harmonises practices by bringing all types of collections under one system; the types currently covered include zoological, botanical, mycological and palaeontological museum collections, tissue and DNA samples, and botanic garden and microbial living collections.Kotka stores data mostly in a denormalised free text format using a triplestore and a simple hierarchical data model. This allows greater flexibility of use and faster development compared to a normalized relational database. New data fields and structures can be added easily as needs arise. Kotka does some data validation, but quality control is seen as a continuous process and is mostly done after the data have been recorded into the system. The data model is loosely based on the Access to Biological Collection Data (ABCD) standard41, but has been adapted for practical needs.Kotka is a web application. Data can be entered, edited, searched, and exported through a browser-based user interface (UI). However, most users prefer to enter new data in customizable MS-Excel templates, which support the hierarchical data model, and upload these to Kotka. Batch updates can also be done using Excel. Kotka stores all revisions of the data to avoid any data loss due to technical or human error.Kotka supports designing and printing specimen labels9, annotations by external users, and handling accessions, loan transactions, and the Nagoya protocol10.Notebook OMSNotebook is the other primary data management system of FinBIF. It is a web solution for recording opportunistic as well as sampling-event-based species observations. It is being used for systematic monitoring schemes, various citizen science projects, and platforms for species enthusiasts.Notebook’s main software component is LajiForm, which is the engine that renders a given JSON Schema into a web form. LajiForm is a separate, reusable module that is fully independent from other FinBIF systems. Notebook as a whole includes other features embedded in FinBIF, such as adding complex geographical shapes to observation documents, importing data from spreadsheets, and form templates.All Notebook forms use FinBIF’s ontological schema in the JSON36 Schema format. Rendering user-friendly web forms based on a single schema is difficult, because the web form should be asking meaningful questions, instead of just rendering the schema fields according to the form description. Questions should be presented in an interactive manner. For instance, after drawing a geographical location on a map for a potential flying squirrel nesting tree, one would ask “did you see droppings at the nest?”, and answering “yes” would update the document to include a flying squirrel taxon identification with fields “breeding” and “record basis” filled in but not rendered to the form. A simpler form engine without a user interface (UI) customisation layer would just render the fields “taxon”, “breeding” and “record basis”, and the user would have no understanding why there are so many fields to fill in and how they relate to their work or study.Some Notebook forms are complex, e.g., for experienced biology enthusiasts who need a form that is advanced, customisable, and compact. Some forms are simple, e.g., for elementary school children. To tackle this, LajiForm uses a separate schema for the UI that allows everything from simple customisation, such as defining widgets for fields, changing field order or customising field labels, to more complex customisation like transforming the schema object structure, defining conditions when certain fields are shown, or if updating a field should have an effect on other fields. All the functionality is split into a loosely coupled collection of components, which can be either used as standalone components or composed together in order to achieve more advanced customisation. The programming philosophy has drawn inspiration from functional programming, which has been helpful in writing isolated, composable functionality.LajiForm is written with the JavaScript framework ‘React’. LajiForm is built on top of react-jsonschema-form (RJSF), which is an open source JSON schema web form library founded by Mozilla (see https://react-jsonschema-form.readthedocs.io/en/latest/). RJSF handles only simple customization, but it is very flexible in design and allows building extensions with features that are more powerful. Some features and design proposals were submitted to Mozilla – FinBIF is the largest outsider code contributor to RJSF, with a dozen pull requests merged.iNaturalistFiiNaturalist (https://www.inaturalist.org/) is an international observation and citizen science application and platform. The iNaturalist Network is a collection of websites that are localized to national use in c. 10 different countries. FinBIF supports the Finnish network site iNaturalist Finland through translations, instructions, communication, moderation, and user support.Finnish iNaturalist data are automatically synchronized weekly to FinBIF’s data warehouse, where they are available for local use. iNaturalist observations are linked to the observer’s own account in the FinBIF portal, if they have linked their iNaturalist and FinBIF accounts. Both features are important in encouraging observers to use iNaturalist, and to allow it to work seamlessly with other FinBIF services.Taxon EditorFinBIF has developed its own taxon database, ‘Taxon Editor’. It allows taxon specialists to maintain their own, expert-validated view of Finnish species. The aggregation of these is used as a backbone taxonomy for all FinBIF services, and the national checklist of Finnish taxa is extracted from it (see https://laji.fi/en/theme/checklist). Each taxon is given a globally unique persistent HTTP-URI identifier, which refers to the taxon concept, not to the name. The identifier does not change if the taxon concept remains unchanged. Compatibility with checklists from other countries is sought by linking taxon concepts as Linked Data.The taxon specialists (currently c. 60) maintain the taxon data using a web application. All changes made go live every night. The nightly update interval allows the specialists a grace period to make their changes. To maintain the integrity of critical data, such as lists of protected species, limitations to what the specialists can do have been imposed. Changes to critical data are carried out by an administrator.Taxon Editor has special features for linking observations to the taxonomy. These include hidden species aggregates and tools to override how a certain name used in observations is linked to the taxonomy. Misapplied names, however, remain an unresolved problem. Most observations are still recorded using plain names, but it is possible for the observer to pick a taxon concept instead, which is the most precise way. When data are published through the FinBIF portal from other information systems, the data providers can link their observations to FinBIF’s taxon concepts by providing the concept’s identifier. The ability to use taxon concepts as a basis of observations means the concepts have to be maintained over time, a task that may become arduous in the future. For further description of the functionalities of Taxon Editor, see7.Taxon Editor is also used in Red List assessments21. The threat assessment is carried out using the criteria of the International Union for Conservation of Nature (IUCN). FinBIF offers a documentation tool and an archive for the assessment, which is based on the national checklist of Finnish taxa. Information about previous assessments is available in the tool, and the assessor can copy and confirm, e.g., area of occupancy, extent of occurrence, generation length and habitat preferences of a species from the previous assessment. The service offers the possibility to add notes to most of the fields separately and commenting on the assessments by other authorised users. In line with the IUCN instructions, the tool automatically chooses the criteria leading to the highest possible threat category of criteria filled out for the species, although the assessor confirms the final evaluation. In several fields, the tool automatically checks the validity of values entered, e.g., criteria, threat category, length of observation period, causes of threat, and current threat factors. The tool includes necessary fields for backcasting the categories of previous assessments to count the Red List Index42. There is also a possibility to do regional threat assessments. Data are stored to the triplestore34, which archives the history of all changes.Media serviceThe FinBIF media service currently supports receiving, transforming, storing and serving images and audio. For images, the original media is stored, and the service generates a smaller JPEG version to be used in the web. The service also generates different sized thumbnails. Only a handpicked set of Exchangeable image file format (Exif43) metadata is kept, so that the metadata would not leak location information about sensitive species occurrences. For audio, wav and mp3 formats are supported. The original file is processed (cleaned) to prevent any malicious content. Mp3, wav and a spectrogram are generated and stored. Support for the International Image Interoperability Framework (IIIF; https://iiif.io/about/) standard is under planning.Data upload serviceMany occurrence datasets are not yet maintained in modern IT systems or databases, which could use an API to transfer data. FinBIF is able to receive data from Excel and GIS systems as secondary data to its data warehouse. First, metadata are generated about the dataset, and selected users are given access to upload data to that dataset. Then the owner of the data must transform the data to a row/column-based table, i.e., MS-Excel or tab-separated value (TSV) file. Each row has one occurrence and must have an ID that is unique to that dataset. The dataset owner then proceeds to upload the row/column-based file using a Web UI, in which the owner maps the fields and values of their data to FinBIF schema fields and values. Updates are done by re-uploading the entire dataset or only a part of it. Deletions are handled by a specific column that annotates that occurrence as deleted.Taxonomy backendThe taxonomy backend transfers the data created in Taxon Editor to an Elasticsearch search engine on a nightly basis. FinBIF’s API is built on top of the data in Elasticsearch. GBIF’s taxonomic backbone is currently being integrated into FinBIF’s taxonomy for taxa with occurrences but no taxonomy in FinBIF. This allows FinBIF users to browse a taxonomy that is a combination of the FinBIF and GBIF taxonomies. The taxonomy backend could harvest descriptions and images from other sources, but currently these functionalities are disabled because of data quality problems.Collection metadata backendCollection and dataset metadata are maintained primarily using the Kotka CMS, but metadata are also harvested from partner organisations’ metadatabases. The metadata of a collection contain the taxonomic, geographic and temporal coverage of the dataset, as well as information about its quality using a three-level grading: (1) Professional; (2) Expert hobbyist / expert curated; (3) Citizen science / mostly non-curated. Occurrence data can be filtered in the FinBIF data warehouse based on these levels. Each occurrence also has its occurrence-specific quality grading, which is different from the dataset grading. The collection metadata define, e.g., the people who handle data requests, which are done in FinBIF’s restricted data request service.Taxonid.orgFinland and Sweden are piloting, with a subset of taxonomic groups, to connect national checklists using Linked Open Data standards44 and agreed vocabularies. By using HTTP-URI as globally unique, persistent identifiers for taxon concepts45, the service provides both human-readable (Hypertext Markup Language, HTML) and machine-readable (Extensible Markup Language, XML) responses for client requests via a central checklist (http://taxonid.org/). Future steps include linking the national lists with the global reference checklist developed by Catalogue of Life (https://www.catalogueoflife.org/). Currently the service includes only taxonomic information, but the eventual goal is to share information also on genetics, images, and traits, as well as on conservation status and observations, in a standardised way. The work was part of the DeepDive project which was funded by the Nordic e-Infrastructure Collaboration (https://neic.no/deepdive/). The vision is to establish a regional infrastructure network consisting of Nordic and Baltic data centres and information systems, and to provide seamlessly operating regional data services, tools, and virtual laboratories.Application programming interface (API)A fundamental goal of FinBIF is to provide all data in machine-readable formats, so that the data and services can be used by third party applications. This can be accomplished using FinBIF’s public API (https://api.laji.fi). It provides access to all data available in FinBIF. The FinBIF portal is built using solely the public API, which should ensure the API is robust, of high performance, reliable and easy to use.GIS servicesFinBIF aims to share occurrence data in GIS formats (WFS, WMS; see https://www.ogc.org/standards/). These services are currently under construction.Restricted data request serviceThe restricted-use data that FinBIF harbours play a crucial role in, e.g., land-use decisions and conservation. To make these data findable, and conditionally accessible, through the public portal, a Restricted-use Data Request Service (RDRS) has been employed. It allows any user with a valid justification to request access to restricted-use datasets. After selecting the required compilation of data, the user submits a standardised form with the required information. The often multiple data owners receive a notification and a request to log in to the data owners’ section of the portal to scrutinise the request. The owners may discuss the request privately among themselves to facilitate decisions. In case of discrepancy by separate owners, the user who submitted the request may still download the released part of the data.Occurrence annotationFinBIF maintains an annotation system, which allows adding information to all occurrence records in the data warehouse. Any registered user can add comments and mark occurrences as needing verification. Trusted users can be given an expert annotator status, which allows them to have more effect on how the occurrences are shown. Experts can change the identification (taxon) that is displayed by default, grade occurrences based on their quality (verified, unreliable, erroneous), and override other annotations. However, annotations never change the original occurrence, which is always kept available. Data owners can be notified about annotations regarding their data, so that they can check and correct possible errors in the primary data source.R-libraryFinBIF provides an R programming language interface to the FinBIF API. The FinBIF R package makes the publicly available data in FinBIF accessible from within R. Biodiversity information is available on taxonomy and taxon occurrence. Occurrence data can be filtered by taxon, time, location and other variables. The data accessed are conveniently preformatted for subsequent analyses. Documentation and download can be found at https://luomus.github.io/finbif/.Map servicesFinBIF has built a JavaScript map service on top of a popular Leaflet library. It provides import and export in various formats, support for the Finnish national coordinate systems and a legacy coordinate grid layout, which is still widely used as a basis of monitoring schemes, though no longer officially supported by Finland’s geographical authorities. The tool provides a rich collection of national and international map layers that are useful in reporting and evaluating occurrence data. The tool also allows calculating lengths and areas as well as drawing complex geographical shapes (e.g. survey polygons and polylines with buffer areas).AI-based species identificationFinBIF cooperates with projects that utilise machine learning in species identification and data classification. They involve creating (semi-)automated identification pipelines for Finnish fungus species from images, and bat and bird species from audio recordings. FinBIF has provided images of fungi, is building a crowd-sourcing platform on which bird experts can produce training material, and aims to implement the Animal Sound Identifier software46 for building on-line identification services.User authenticationFinBIF provides an authentication and authorisation service that can be used also by third party applications and websites. Separate authentication flows exist for websites and native applications. Users can create one FinBIF account per e-mail address and associate multiple authentication methods and other user identities to their account. For example, users can login with Google or Facebook credentials and associate their iNaturalist, Finnish Wildlife Agency and Finnish Bird Ringing user identities to their FinBIF account. This enables users to see the occurrence data they have entered into various systems as their own occurrences in the FinBIF portal.Servers, service continuity, backup and disaster recoveryFinBIF uses services provided by CSC, the Finnish IT centre for science (https://www.csc.fi) and the University of Helsinki IT department. Both provide FinBIF several OpenStack (https://www.openstack.org/) based virtual servers and an OpenShift (https://www.openshift.com/) cloud platform. FinBIF is not committed to any certain Service Level Agreement (SLA), but its availability is above 95%. All primary data gathered by FinBIF are professionally backed up by the University of Helsinki. CSC IDA (https://www.fairdata.fi/en/ida/) is used to archive larger datasets, such as specimen images. The CSC digital preservation service (https://www.fairdata.fi/en/fairdata-pas/) will be used for long-term archiving.Development methodologiesAgile developmentThe FinBIF in-house ICT team uses a three-level development process. The first level provides capacity to do long-term planning. Epic level requests (entirely new parts of the infrastructure, inclusion of new data sources, new monitoring schemes) are compiled, prioritised and to some degree scheduled using a Trello board. The aim is to open this board for all stakeholders and the public, so that interested parties can track progress and focus of development. This board is maintained in biweekly sessions.The second level is used to direct current focus. A combination of Kanban (see https://www.atlassian.com/agile/kanban) and Scrum47 methodologies is used. Every two weeks the status of ongoing epics is checked, and the epics for the next two weeks are decided. Unlike in Scrum, there is no attempt to estimate or set goals for the sprints. Rather, a Kanban-like free-flowing system is used, where things take as long as they take.The third level is the implementation process for individual epics. Like in Scrum, each epic has a named product owner (PO). The PO communicates with the stakeholders and defines user stories. The development team turns the user stories into smaller tasks, which are maintained in a Pivotal tracker backlog (www.pivotaltracker.com).It is hoped that, in the future, developers from partner organisations could be attracted to participate in the development process.Data acquisition for this paperComparative data for other BDDIsThe data shown in Table 1 and supporting the drawing of Fig. 1 were acquired by searching for information on the public internet portals of the infrastructures between July 2019 and March 2020. The data were then sent to the representatives of the infrastructures for checking and possible corrections. Altogether 27 infrastructures (eight global and 19 national or institutional) were contacted, and answers were received from 23 (six global and 16 national or institutional) infrastructures. The following information was gathered.
What is:

The taxonomic coverage of your infra?

Current number of records?

Number of species you have information on?

Does your infrastructure share the following types of data and information:

Natural history collection data?

Opportunistic observation data?

Observation data collected in systematic monitoring schemes?

Verbal descriptions of species?

Standardized taxonomic core data?

DNA-barcodes (or links to BOLD database)?

Does your infrastructure provide the following kinds of user services:

Possibility to enter event-based data, meaning several species observations linked together?

A Collection Management System?

IUCN red-listing tools?

Red-list classifications or administrative statuses of species?

e-lab services? What kind of?

Is the infrastructure designed to enable:

Data generation, such as digitization of natural history specimens?

Collation / aggregation of data?

Entering observations by citizen science users?

Annotation of records/identifications?

FinBIF’s user dataTotal numbers of users and use sessions were acquired from Google Analytics (https://analytics.google.com) for the website laji.fi in its entirety. To determine which traffic belongs to which user, Google Analytics sends a unique identifier associated with each user with each hit. This is accomplished via a Client ID, a unique, randomly generated string that gets stored in the browser’s cookies, so subsequent visits to the same site can be associated with the same user. Using cookies allows identifying unique users across browsing sessions, but it cannot identify unique users across different browsers or devices. Hence, the user numbers presented here are slight overestimates of true numbers of different people using the service. The use sessions are total number of sessions within the given range of dates. A session is the period during which a user is actively engaged with the website. All usage is associated with a session.The number of registered users comes from FinBIF’s own registry. Until the end of 2019, users had to register at the service to be able to download data from laji.fi to their own device. Subsequently, the requirement has been relaxed so that registration is required only to obtain a citable persistent unique HTTP-URI identifier for a downloaded data batch, but non-citeable so-called light downloads can be done without registration. Registration is also required to be able to record data to FinBIF’s primary data management systems and to annotate records at https://laji.fi/.Numbers of data downloads and downloaded data points are logged by FinBIF itself. As of the beginning of 2020, the numbers reported in Fig. 3b include light downloads. More

‘Foraging habitat suitability’ is a reference to the favorability of a patch of land for day-to-day subsistence. Here, suitability is an index value ascribed to each potential foraging patch (grid cell) captured in a raster image, based on terrain movement costs and the proximity of each patch to water and green vegetation. We constructed our foraging habitat suitability model using satellite-derived environmental data, digital terrain information and anthropological field data on foraging range (Fig. 3). The model’s environmental foundation is based on more than two decades of continuous near bi-weekly Landsat-5 satellite observations, allowing for the systematic detection and measurement of water recurrence and vegetation condition for every 30-x-30 m image pixel. This period of observation is long enough to observe multiple fluctuations in this highly variable environment and to not be restricted to a single short-term climatic state, such as a bushfire or drought. Therefore the time frame provides a reliable observation and measurement of maximum vegetation greenness, regardless of temporary drops in NDVI. Similarly, maximum extent and occurrence of surface water is systematically measured through long-term satellite observations, avoiding measurements only of phases of drought or irregular rainfall. For this reason our model focuses on maximal values to represent the best environmental conditions that would have been available for past foraging activities since the last glacial, based on the contemporary climatic regime.Figure 3A satellite derived model of foraging habitat suitability for the Australian Western Desert. Foraging habitat suitability is highly variable within IBRA boundaries and throughout the Western Desert. Several massive areas of low-ranked foraging habitats are evident throughout the region. IBRA codes and excavated rockshelter sites (lime green- numbered) are defined in the Fig. 1 caption. Map created in ESRI ArcGIS Desktop 10.5.1 (https://desktop.arcgis.com), linear stretch (1.0%) visualization. See “Methods” section for source raster information.Full size imageThe model also uses the ALOS World 3D 30 m digital elevation data product to quantify terrain ruggedness across the study area46. Terrain ruggedness is a geomorphometric measure of land surface ruggedness, where elevation variability is used to infer ease of traversal when walking between locations in the landscape. Terrain ruggedness is suggestive of potential energy expenditure, assuming that increasingly rugged terrains necessitate higher levels of physical activity and caloric intake. Here, we integrated measures of ruggedness with environmental satellite data, providing an indication of which patches of vegetation and water are most easily accessed in regards to minimum changes in elevation.Walking time to observed surface water is the final spatial parameter incorporated with the model. It is calculated using Tobler’s47 hiking algorithm and information on daily foraging practices. Historic anthropological data indicates Western Desert foraging activities typically operated for 4–6 h each day1,48, with foragers moving up to a day from ephemeral water sources in their food quest1. In accordance with these ethnographic statements we spatially delineated land areas where regular foraging activities may have occurred by first calculating the walking time from water, then weighting all areas that were less than 8 h walk from water more heavily in the input which went into our final suitability model. Since resources are said to be permanent in uplands4,5, we assume mountainous refugia were always suitable foraging habitats, so these refuge areas have been masked and removed from consideration (see mountain ranges in Fig. 3).Appropriate elements from all of the aforementioned satellite datasets were combined to produce our foraging habitat suitability model (Fig. 3). The ~ 30 m spatial resolution of the data facilitates the construction of a spatially-explicit, geographically broad, yet fine-grained ecological model to visually observe and critically appraise foraging habitat suitability at a variety of scales, offering new perspectives on regional human behavioral ecology. The model provides a continuous ranking of the relative foraging value for each landscape patch (or 30 m grid cell in this instance). Interpretation of patch values is based on the proposition that foragers know the conditions in all parts of the landscape they visit, and they organized their daily foraging movements in accordance with the factors outlined above.Our habitat suitability model illustrates the highly varied favorability of foraging patches across the Western Desert (Fig. 3), as calculated from data on natural resource distribution, terrain attributes, and daily foraging range. The model is conceptual, based on quantitative environmental variables that have been well documented to influence desert foraging activity. In regards to the model’s robustness, the input variables are equally weighted and statistically independent (see “Methods” section). The equal weighting reflects the concepts and assumptions of earlier research, particularly of existing landscape mapping, offering a coherent and consistent modelling approach. Advanced mathematical modelling, incorporating sensitivity analysis49,50, could be used to modify the weighted contribution of each variable, and such modelling will be the subject of future papers. Until more detailed knowledge of past forager land use and contemporary resources becomes available there is little benefit in arbitrarily substituting other input values in our model.The model comprises a matrix of nearly 1.3 billion data cells, each of which has been individually analyzed and ascribed each foraging patch a value indicative of potential habitat suitability. The computational power required to statistically analyze the dataset is massive, so to simplify computing and broadly characterize intraregional variation, we scaled up using nationally defined IBRA subregions. We used IBRA boundaries to group and rank the patch values into low, moderate, and high foraging habitat suitability classes and then calculated the land area occupied by each class (Fig. 4 and Table S1). Higher-ranked localities are well positioned in relation to suitable resources and easily traversed terrains. Lower-ranked patches are considered poorly-suited habitats due to their considerable distance from water and plant resources, and they are in comparatively rugged terrains. Areas deemed to have moderate foraging suitability have mixed accessibility to resources and variable terrain ruggedness.Figure 4Percentage of land area (km2) occupied by low, moderate, and high-ranked habitat suitability patches for the eleven largest IBRA bioregions of the Western Desert (Table S1). The histogram is ordered left to right based of the percentage of high-ranked foraging habitat within each bioregion. The percentages for the entire Western Desert are presented on the far right.Full size imageThe results show that during times of maximum water abundance and vegetation greenness, 36.6% of the Western Desert has high-ranked habitat suitability (Fig. 4 and Table S1). Moderately suitable areas constitute 48.9% and low-ranked patches encompass 13.1%. Breaking these findings down further, we calculated the ranked land areas for the eleven largest IBRA subregions ( > 10,000 km2) of the Western Desert (Fig. 4 and Table S1). At a broad bioregional level, intra-upland zones (Fig. 4; CER01) and desert plains (Fig. 4; GAS02, GVD01, and NUL01) offer a greater percentage of high-ranked foraging habitats. Bioregions dominated by dunefields have considerably less high ranked land areas compared to uplands, plains, and areas of low relief, although it is important to note that there is also considerable patchiness amongst suitable foraging areas in sandridge desert regions (Fig. 3). For instance, the centrally located Gibson Desert dunefield area (Fig. 4; GID02) has very little area of high-ranked habitat (10.8%), which is far less than other sandridge desert bioregions (Fig. 4; GSD02, LSD02, GVD02, GVD03, and GVD04) where high-ranked suitability areas range between 22.3 and 39.9%. Similarly, the Gibson Desert stony desert bioregion (Fig. 4; GID01), which is dominated by lateritic surface gravels, records only 25.5% high-ranked habitat areas. Thus, at a coarse-grained scale, it seems that some central core regions of the Western Desert are more environmentally hostile and offer less high-ranked foraging opportunities compared to more peripheral bioregions. This generality does not imply such areas were unutilized by desert peoples, but rather some areas were on average volatile and had low productivity.Foraging potential is highly varied amongst bioregions and land systemsWhen viewed at a fine-grained scale, our model clearly shows that there is an uneven gradient of suitable foraging habitats across the Western Desert, and foraging suitability trends are not pervasive throughout particular bioregions or land systems (Fig. 3). Away from montane uplands, water permanence is always temporary, and land systems with low topography, such as plains, stony plains, and sandridge desert, have highly varied foraging suitability, even when characterized in the best environmental conditions.The implications of this variation are important to understanding human ecology of the ethnohistoric period and the late Holocene archaeological record of the past 2000 years, when climatic conditions and landscapes were much like the present day36,40,42. Many scholars have noted that the historic desert peoples were familiar with the distribution of regional natural resources1,5,7. It has been argued that resource knowledge was articulated with socioeconomic strategies, and that groups routinely utilized all areas of the Western Desert during times of good rainfall and resource abundance. However, our suitability model reveals that there are large, expansive areas of the desert landscape that would have presented substantial challenges for survival, even in the best environmental circumstances (Fig. 3).Our model further suggests that low-ranked locations of foraging suitability were always below average productivity and were always comparatively unsuited as foraging habitats. To carry out that measure, we needed an independent indicator of land productivity, NDVI. We used satellite observations of maximum vegetation greenness to quantify how land productivity differs amongst low, moderate, and high ranked foraging habitats (Table S2). Variation in mean (µ) NDVI for each habitat class illustrates how land productivity differs within and amongst the most prominent Western Desert bioregions (Fig. 5). Given the below average NDVI of all low-ranked desert lowlands, we hypothesize that broad clusters of extremely unsuitable localities would be unlikely to provide adequate returns (Fig. 5), even when foragers were pursuing low-variance or lower quality resources51. Based on the distribution of low-ranked patches (Fig. 3), we agree with earlier research that the entire desert region was not equally economically viable for foraging, and that substantial tracts of land were not economically attractive to resident populations4,5,14,32. We also recognise that the distribution of massive-sized sub-optimal patches may be an important factor shaping the patterns of movement through the landscape, with foragers potentially preferring movement along high suitability corridors. However, unlike earlier research, our suitability model shows that unfavorable foraging areas are not correlated with large units of biogeography alone. Our model depicts the environmental variability of the Western Desert at a much higher resolution than its predecessors, revealing several massive land tracts where unfavorable foraging conditions occur (Fig. 3). If ethnographic patterns of land use were in place, we predict that many of these large areas would have been rarely utilized or perhaps some were purposefully avoided due to known deficiencies in the resource energy base12. This proposition is readily testable because it predicts that archaeological sites with poorly sorted, low densities of artefacts will be found in these places12. Defining the appropriate scale will be the key to testing our model, since we have demonstrated that broad biogeographic units are heterogenous and yet at a fine-grained scale, small areas of low suitability, which are often a local geographic feature (e.g., sand dune, bare rock outcrop, or erosional area), need not have been obstacles. Model testing will need an intermediate scale commensurate with daily foraging radii.Figure 5Boxplot of mean NDVI values and one standard deviation for low, moderate, and high-ranked suitability classes for the eleven largest IBRA subregions (a–k) of the Western Desert (l). Mean NDVI for individual bioregions and the spatial bounds of the Western Desert study area denoted as dashed black line and solid green line, respectively. IBRA subregion boxplot groups (a–k) are presented in order of increasing percentage of high-ranked foraging habitat, after Fig. 4. Table S2 offers the precise summarised NDVI values for each bioregion and suitability class.Full size imageAt present, the archaeological land use pattern of low-ranked foraging habitats is not something that is well-understood from the Western Desert, although periodic and short term use of impoverished, low productivity patches has been predicted12. Studies of contemporary Western Desert groups indicate that human-induced firing of the landscape enhances biodiversity and land productivity51,52,53,54, so it is possible that low productivity patches may have occasionally benefited from anthropogenic burning, especially in the past 1500 years51. However, research also suggests that cultural burning practices did not have widespread regional impacts51,52,53,54,55. Human influence on landscape modification is localized within day-range foraging areas around residential camps and frequently traversed pathways51,52,53,54. Low productivity patches away from residential camps were probably unlikely targets for either anthropogenic burning or foraging if higher-ranked patches were closer.Elsewhere, in the eastern Australian arid zone, periodic use of climatically harsh desert localities is known from archaeological sequences. While in some cases preservation may explain chronological discontinuities56, there is compelling evidence for irregular occupation in several desert areas10,57,58,59,60. For instance, in the western Strzelecki Desert broad portions of dunefield landscapes were periodically abandoned for centuries or even millennia57,60, while in semi-arid portions of southeastern Australia sequences of occupation were separated by decades or centuries of local/regional abandonment58,59. Fluctuations in local foraging suitability may well be a factor producing discontinuous land use across the Australian arid lands, and we suggest that in the Western Desert there were patches with chronologically varying foraging potential. The key test of this prediction would be to investigate whether archaeological sites in locations of fluctuating habitat suitability over time also display histories of discontinuous visitation. Such sites could be identified through local palaeoenvironmental records but we suggest that selections based on time-series analysis of vegetation greenness from the past few decades would be more readily used to establish samples and would facilitate comparison of archaeological sites in terms of local foraging suitability and NDVI values, as well as archaeological records of continuous or discontinuous visitation.Satellite data reveals a more nuanced understanding of land useAustralian archaeological research has relied heavily on biogeographic principles to distinguish the ‘barriers and boundaries’ of Aboriginal subsistence and settlement in the arid zone4,5,61. While equating particular land use practices with specific bioregional areas was initially useful for generalized conceptualizations of traditional foraging behaviors, the coarse analytical scale of earlier approaches is now problematic. Subsequent research has shown the dynamics of Aboriginal occupation and land usage in the Western Desert to be more complex and variable across spatial and temporal scales than originally conceived9,24,30,33. To gain a more nuanced understanding of past land use and foraging patterns, finer-scale methods of analysis are required.We used satellite imagery to tackle the issue of scale, allowing for a sharper and more spatially explicit examination of desert environments and landscapes. For example, as we focus at higher resolution on various areas of the Western Desert, our model clearly shows that foraging suitability is highly varied across all desert lowlands (Figs. 3, 4 and 6). In sandridge desert areas, proposed to have been a barrier at times in the past4, the model shows there are many well-watered and amply vegetated localities where good foraging is possible when rainfall is high and surface water is abundant (Fig. 6a). In this context, interdunal swales are hardly barriers to occupation because they can be lush with water, plant, and wildlife resources after local rain, and the energy expenditure required to walk along interdunal swales is low in comparison to the requirements needed to repeatedly scramble across a sea of loose sands and undulating dunes. Thus, it seems entirely plausible that resident groups could navigate and forage in many dunefield areas by following a well-resourced network of swales during times of good environmental conditions. The fine-grained nature of this observation opens up the possibility that many sandridge deserts were not necessarily broad barriers to occupation and that precontact land use behaviors varied in different dunefield contexts.Figure 6High resolution perspectives of various Western Desert landforms (e.g., sandridge, stony plain and sandy plain contexts) with generally higher-ranked and lower-ranked areas of foraging suitability. This figure illustrates the fine-grained scale of our habitat suitability model (Fig. 3), which has implications for better understanding localized land use behaviors. Juxtaposed areas, as mapped in Fig. 3 are: (a) Higher-ranked sandridge habitats vs. (b) lower-ranked sandridge land system. (c) Higher-ranked stony desert habitat vs. (d) lower-ranked stony desert areas. (e) Higher-ranked sandplain land systems vs. (f) lower-ranked plain habitats. Maps created in ESRI ArcGIS Desktop 10.5.1 (https://desktop.arcgis.com), linear stretch (1.0%) visualization. See “Methods” section for source raster information.Full size imageWe also highlight that the resource-rich swale pattern is not found in all dune systems (Fig. 6b), and it is plausible that some of these areas were periodic barriers to occupation, as previously suggested in more generalized ecological models4. There are substantial areas of sandridge desert, especially within central areas of the Western Desert (e.g., GID02), where survival would have always been extremely difficult, even during times of abundance (Fig. 6b). This variability is also expressed in stony desert contexts, where southern areas of the lateritic Gibson Desert (GID01) offer better habitat suitability (Fig. 6c) than the northern areas (Fig. 6d). On a fine scale, plain land systems also exhibit a wide range of habitat suitability, where high-ranked habitat suitability appears fairly widespread in some areas of the Nullabor Plain (NUL01; Fig. 6e), yet other areas of the plain were poorly-suited for foraging (Fig. 6f).In previous ecological models, stony desert and plain land systems are considered more favorable than sandridge desert4; however, as shown above, the modelled data clearly illustrate that there are substantial areas of plains and stony desert landscapes that vary considerably between high and low-ranked habitat suitability (Figs. 3, 4 and 6). The fine-grained scale of our model adds to a growing body of research5,9,24,30,33 that demonstrates how previous pan-continental characterizations of deserts as ‘corridors’ and ‘barriers’ for foragers oversimplify the link between human behavior and biogeography. When scrutinized at high resolution, extremely unsuitable foraging and very well-suited foraging areas can potentially occur in any area of the Western Desert, regardless of the biogeography or other physical characteristics. Thus, fine-grained ecological models allow for a more nuanced and spatially-explicit understanding of the past land use behaviors that led to the formation the desert archaeological record.Using environmental remote sensing to infer LGM habitat suitabilityThere is no doubt that the Western Desert environment has changed and evolved over time, through both natural and human-induced processes8,36,37,51. The region has undergone considerable environmental fluctuations over time, resulting in landform transformations (dune aggradation, in particular) and changes in vegetation cover. The long-term physical impact of these environmental changes clearly place limitations on how modern satellite data can be used to interpret deep-time patterns of occupation and land use. However, our model shows the likely distribution of low-ranked foraging habitats when climatic conditions were much drier than present.Ethnoarchaeological accounts depict resident populations as being low density and highly mobile, frequently moving and foraging across vast expanses of territory, thereby necessitating intermittent patterns of settlement1,2,3,4. Such a mobile strategy means that large swathes of the desert could not have been continuously occupied. Our habitat suitability model (Fig. 3) makes sense of the impermanent and mobile land use strategy seen historically. For example, we document several massive areas of the Western Desert where, in combination, surface terrains are physically challenging, the nearest proximity to surface water is greater than 2 days walk, and vegetation cover, density and condition is substandard, even in the best of documented environmental conditions! These exceptionally large areas were poorly-suited to foraging, and could not have been permanently occupied in the historical period. They would only have been visited rarely, perhaps only in atypical short-term climatic events, and may even have constrained forger movement between more favorable parts of their territories. Given current palaeoclimatic evidence we infer that in the Pleistocene these low ranked habitats would have been even more inhospitable to foragers than in recent times. During the LGM, the resource yields in such areas would have been more diminished than present, making conditions for survival even more difficult than today. Consequently, we predict that unless radically different economic strategies were being employed in the Pleistocene those areas would have been only rarely visited since the peak of the last glacial cycle, ~ 24–18 ka, even though adjacent desert areas may have supported regular or at least sporadic visitation.Our hypothesis is clear, detailed, and framed to be testable by archaeological fieldwork. The number of Western Desert sites with old archaeological sequences is growing, but the sample is small, site distribution is widely scattered, and none are located in the harsher core areas identified in this study (Figs. 1 and 3). Thus, it is evident that archaeological fieldwork in those impoverished landscapes as well as environmentally richer and more reliable landscapes is necessary to understand historical land use patterns and to make statements about earlier phases of regional occupation. Our work highlights how future models of forager land use across Australia’s desert regions can comprehend the environmental complexity and fine scale of resource variability in these vast, remote and diverse places. More

1.He, T., Belcher, C. M., Lamont, B. B. & Lim, S. L. A 350-million-year legacy of fire adaptation among conifers. J. Ecol. 104, 352–363 (2016).Article 

Google Scholar 
2.Doerr, S. H. & Santín, C. Global trends in wildfire and its impacts: Perceptions versus realities in a changing world. Philos. Trans. R. Soc. B Biol. Sci. 371, 20150345 (2016).Article 

Google Scholar 
3.Hoegh-Guldberg, O. et al. Impacts of 1.5°C Global Warming on Natural and Human Systems. in Global Warming of 1.5°C. An IPCC Special Report on the impacts of global warming of 1.5°C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, (ed. Masson-Delmotte, V., P. Zhai, H.-O. Pörtner, D. Roberts, J. Skea, P.R. Shukla, A. Pirani, W. Moufouma-Okia, C. Péan, R. Pidcock, S. Connors, J.B.R. Matthews, Y. Chen, X. Zhou, M.I. Gomis, E. Lonnoy, T. Maycock, M. Tignor, and T. W.) 175–311 (2018).4.Dennison, P. E., Brewer, S. C., Arnold, J. D. & Moritz, M. A. Large wildfire trends in the western United States, 1984–2011. Geophys. Prospect. 41, 2928–2933 (2014).ADS 

Google Scholar 
5.Westerling, A. L. R. Increasing western US forest wildfire activity: Sensitivity to changes in the timing of spring. Philos. Trans. R. Soc. B Biol. Sci. 371, 20150178 (2016).Article 

Google Scholar 
6.Bailey, R. & Yeo, J. The Burning Issue (Marsh & McLennan Insights, 2019).
Google Scholar 
7.Province of British Columbia. 2018 Wildfire Season Summary. 2018 Wildfire Season Summary (2019). https://www2.gov.bc.ca/gov/content/safety/wildfire-status/about-bcws/wildfire-history/wildfire-season-summary?keyword=total&keyword=area&keyword=burned&keyword=by&keyword=wildfire&keyword=2018.8.Cal Fire. https://www.fire.ca.gov/incidents/2018/. https://www.fire.ca.gov/incidents/2018/ (2020). https://www.fire.ca.gov/incidents/2018/.9.McCullough, I. et al. Do lakes feel the burn? Ecological consequences of increasing exposure of lakes to fire in the continental US. Glob. Chang. Biol. https://doi.org/10.1111/gcb.14732 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
10.Westerling, A. L. et al. Climate change and growth scenarios for California wildfire. Clim. Change 109, 445–463 (2011).Article 

Google Scholar 
11.Nagy, C. R., Fusco, E., Bradley, B., Abatzoglou, J. T. & Balch, J. Human-related ignitions increase the number of large wildfires across U.S. Ecoregions. Fire 1, 1–14 (2018).
Google Scholar 
12.Balch, J. K. et al. Human-started wildfires expand the fire niche across the United States. Proc. Natl. Acad. Sci. U. S. A. 114, 2946–2951 (2017).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
13.Radeloff, V. C. et al. Rapid growth of the US wildland-urban interface raises wildfire risk. Proc. Natl. Acad. Sci. U. S. A. 115, 3314–3319 (2018).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
14.Wright, R. F. The Impact of Forest Fire on the Nutrient Influxes to Small Lakes in Northeastern Minnesota Author (s): Richard F . Wright Published by : Ecological Society of America Stable URL : http://www.jstor.org/stable/1936180 THE IMPACT OF FOREST FIRE ON THE NUT. 57, 649–663 (1976).15.Carignan, R., D’Arcy, P. & Lamontagne, S. Comparative impacts of fire and forest harvesting on water quality in Boreal Shield lakes. Can. J. Fish. Aquat. Sci. 57, 105–117 (2000).CAS 
Article 

Google Scholar 
16.Tecle, A. & Neary, D. Water quality impacts of forest fires. J. Pollut. Eff. Control 03, (2015).17.Abney, R. B., Sanderman, J., Johnson, D., Fogel, M. L. & Berhe, A. A. Post-wildfire Erosion in mountainous terrain leads to rapid and major redistribution of soil organic carbon. Front. Earth Sci. 5, 1–16 (2017).Article 

Google Scholar 
18.Williamson, C. E. et al. Sentinel responses to droughts, wildfires, and floods: Effects of UV radiation on lakes and their ecosystem services. Front. Ecol. Environ. 14, 102–109 (2016).Article 

Google Scholar 
19.Goldman, C. R., Jassby, A. D. & De Amezaga, E. Forest fires, atmospheric deposition and primary productivity at Lake Tahoe, California-Nevada. Int. Vereinigung Theor. Angew. Limnol. Verhandlungen 24, 499–503 (1990).
Google Scholar 
20.Allen, E. W., Prepas, E. E., Gabos, S., Strachan, W. & Chen, W. Surface water chemistry of burned and undisturbed watersheds on the Boreal Plain: An ecoregion approach. J. Environ. Eng. Sci. 2, S73–S86 (2003).CAS 
Article 

Google Scholar 
21.Earl, S. R. & Blinn, D. W. Effects of wildfire ash on water chemistry and biota in south-western U.S.A. streams. Freshw. Biol. 48, 1015–1030 (2003).CAS 
Article 

Google Scholar 
22.Overholt, E. P., Rose, K. C., Williamson, C. E., Fischer, J. M. & Cabrol, N. A. Behavioral responses of freshwater calanoid copepods to the presence of ultraviolet radiation: Avoidance and attraction. J. Plankton Res. 38, 16–26 (2015).Article 

Google Scholar 
23.Urmy, S. S. et al. Vertical redistribution of zooplankton in an oligotrophic lake associated with reduction in ultraviolet radiation by wildfire smoke. Geophys. Res. Lett. 43, 3746–3753 (2016).ADS 
Article 

Google Scholar 
24.Williamson, C. E. et al. The interactive effects of stratospheric ozone depletion, UV radiation, and climate change on aquatic ecosystems. Photochem. Photobiol. Sci. 18, 717–746 (2019).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
25.Aguilera, R., Gershunov, A., Ilango, S. D., Guzman-Morales, J. & Benmarhnia, T. Santa ana winds of Southern California impact PM2.5 with and without smoke from wildfires. GeoHealth 4, 1–9 (2020).Article 

Google Scholar 
26.Liu, J. C. et al. Wildfire-specific fine particulate matter and risk of hospital admissions in urban and rural counties. Epidemiology 28, 77–85 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
27.Environmental Protection Agency. Air Quality Index, A Guide to Air Quality and Your Health. Encyclopedia of Quality of Life and Well-Being Research (2014).28.Melack, J. M., Sadro, S., Sickman, S. & Dozier, J. Lakes and Watersheds in the Sierra Nevada of California: Responses to Environmental Change. (University of California Press, 2020). https://doi.org/10.2307/j.ctv17hm9sr29.Goldman, C. R., Jassby, A. & Powell, T. Interannual fluctuations in primary production: Meteorological forcing at two subalpine lakes. Limnol. Oceanogr. 34, 310–323 (1989).ADS 
CAS 
Article 

Google Scholar 
30.Jassby, A. D., Powell, T. M. & Goldman, C. R. Interannual fluctuations in primary production: Direct physical effects and the trophic cascade at Castle Lake, California. Limnol. Oceanogr. 35, 1021–1038 (1990).ADS 
CAS 
Article 

Google Scholar 
31.Park, S., Brett, M. T., Müller-Solger, A. & Goldman, C. R. Climatic forcing and primary productivity in a subalpine lake: Interannual variability as a natural experiment. Limnol. Oceanogr. 49, 614–619 (2004).ADS 
Article 

Google Scholar 
32.Winslow, L. et al. Package ‘ rLakeAnalyzer ’. Lake Physics Tools. (2019).33.Read, J. S. et al. Derivation of lake mixing and stratification indices from high-resolution lake buoy data. Environ. Model. Softw. 26, 1325–1336 (2011).Article 

Google Scholar 
34.Goldman, C. R. Primary productivity, nutrients, and transparency during the early onset of eutrophication in ultra-oligotrophic Lake Tahoe Califomia-Nevada. Limnol. Oceanogr. 33, 1321–1333 (1988).ADS 
CAS 
Article 

Google Scholar 
35.Marker, A. F. H. The use of acetone and methanol in the estimation of chlorophyll in the presence of phaeophytin. Freshw. Biol. 2, 361–385 (1972).Article 

Google Scholar 
36.Redfield, G. W. & Goldman, C. R. Diel vertical migration and dynamics of zooplankton biomass in the epilimnion of Castle Lake, California. Verhandlungen des Int. Verein Limnol. 20, 381–387 (1978).
Google Scholar 
37.Elser, J. J. et al. Factors associated with interannual and intraannual variation in nutrient limitation of phytoplankton growth in Castle Lake, California. Can. J. Fish. Aquat. Sci. 52, 93–104 (1995).Article 

Google Scholar 
38.Huovinen, P. S., Brett, M. T. & Goldman, C. R. Temporal and vertical dynamics of phytoplankton net growth in Castle Lake, California. J. Plankton Res. 21, 373–385 (1999).Article 

Google Scholar 
39.Maberly, S. C., King, L., Dent, M. M., Jones, R. I. & Gibson, C. E. Nutrient limitation of phytoplankton and periphyton growth in upland lakes. Freshw. Biol. 47, 2136–2152 (2002).Article 

Google Scholar 
40.R Core Team. A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. (2020).41.Zuur, A. F., Ieno, E. N., Walker, N., Saveliev, A. A. & Smith, G. M. Mixed effects models and extensions in ecology with R. (Springer, 2009).42.Lenth, R. V. emmeans: Estimated Marginal Means, aka Least-Squares Means. (2021). https://cran.r-project.org/package=emmeans.43.Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D. & R Core Team. nlme: Linear and Nonlinear Mixed Effects Models. (2020). https://cran.r-project.org/package=nlme.44.Anderson, M. J. Permutational Multivariate Analysis of Variance (PERMANOVA). Wiley StatsRef Stat. Ref. Online 1–15 (2017). https://doi.org/10.1002/9781118445112.stat0784145.Oksanen, J. F. et al. vegan: Community Ecology Package. (2019). https://cran.r-project.org/package=vegan%0A.46.Environmental Systems Research Institute. ArcGIS 10.8.1. (2020). https://www.esri.com/en-us/home.47.Inkscape Project. Inkscape. (2020). https://inkscape.org.48.Bachmann, R. W. & Goldman, C. R. Hypolimnetic heating in Castle Lake. California. Limnol. Oceanogr. 10, 233–239 (1965).ADS 
Article 

Google Scholar 
49.Kochanski, A. K. et al. Modeling wildfire smoke feedback mechanisms using a coupled fire-atmosphere model with a radiatively active aerosol scheme. J. Geophys. Res. Atmos. 124, 9099–9116 (2019).ADS 
Article 

Google Scholar 
50.David, A. T., Asarian, J. E. & Lake, F. K. Wildfire smoke cools summer river and stream water temperatures. Water Resour. Res. 54, 7273–7290 (2018).ADS 
Article 

Google Scholar 
51.Moeller, R. Contribution of ultraviolet radiation (UV-A, UV-B) to photoinhibition of epilimnetic phytoplankton in lakes of differing UV transparency. Arch. Hydrobiol. Beihefte Ergebnisse Limnol. 43, 157–170 (1994).
Google Scholar 
52.Morris, D. P. & Hargreaves, B. R. The role of photochemical degradation of dissolved organic carbon in regulating the UV transparency of three lakes on the Pocono Plateau. Limnol. Oceanogr. 42, 239–249 (1997).ADS 
CAS 
Article 

Google Scholar 
53.Meyers, P. A. & Lallier-Vergès, E. Lacustrine sedimentary organic matter records of Late Quaternary paleoclimates. J. Paleolimnol. 21, 345–372 (1999).ADS 
Article 

Google Scholar 
54.Lamb, A. L., Wilson, G. P. & Leng, M. J. A review of coastal palaeoclimate and relative sea-level reconstructions using d 13 C and C/N ratios in organic material. (2005). https://doi.org/10.1016/j.earscirev.2005.10.00355.Maxwell, T. M., Silva, L. C. R. & Horwath, W. R. Integrating effects of species composition and soil properties to predict shifts in montane forest carbon–water relations. Proc. Natl. Acad. Sci. U. S. A. 115, E4219–E4226 (2018).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
56.Bao, H., Niggemann, J., Luo, L., Dittmar, T. & Kao, S. J. Aerosols as a source of dissolved black carbon to the ocean. Nat. Commun. 8, 1–7 (2017).ADS 
Article 
CAS 

Google Scholar 
57.Zhang, Y. et al. Dissolved organic carbon in glaciers of the southeastern Tibetan Plateau: Insights into concentrations and possible sources. PLoS ONE 13, e0205414 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
58.Solomon, C. T. et al. Ecosystem consequences of changing inputs of terrestrial dissolved organic matter to lakes: Current knowledge and future challenges. Ecosystems 18, 376–389 (2015).Article 

Google Scholar 
59.Banse, K. Rates of growth, respiration and photosynthesis of unicellular algae as related to cell size—A review. J. Phycol. 12, 135–140 (1976).
Google Scholar 
60.Gao, K., Li, G., Helbling, E. W. & Villafañe, V. E. Variability of UVR effects on photosynthesis of summer phytoplankton assemblages from a tropical coastal area of the South China Sea. Photochem. Photobiol. 83, 802–809 (2007).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
61.Häder, D. P., Helbling, E. W., Williamson, C. E. & Worrest, R. C. Effects of UV radiation on aquatic ecosystems and interactions with climate change. Photochem. Photobiol. Sci. 10, 242–260 (2011).PubMed 
Article 
CAS 
PubMed Central 

Google Scholar 
62.Priscu, J. C. & Goldman, C. R. Seasonal dynamics of the deep-chlorophyll maximum in Castle Lake, California. Can. J. Fish. Aquat. Sci. 40, 208–214 (1983).CAS 
Article 

Google Scholar 
63.Leach, T. H. et al. Patterns and drivers of deep chlorophyll maxima structure in 100 lakes: The relative importance of light and thermal stratification. Limnol. Oceanogr. 63, 628–646 (2018).ADS 
CAS 
Article 

Google Scholar 
64.Priscu, J. C. & Goldman, C. R. The effect of temperature on photosynthetic and respiratory electron transport system activity in the shallow and deep-living phytoplankton of a subalpine lake. Freshw. Biol. 14, 143–155 (1984).Article 

Google Scholar 
65.Modenutti, B. E. et al. Effect of volcanic eruption on nutrients, light, and phytoplankton in oligotrophic lakes. Limnol. Oceanogr. 58, 1165–1175 (2013).ADS 
Article 

Google Scholar 
66.Horne, J. A. & Goldman, C. R. Zooplankton and zoobenthos. in Limnology 265–298 (McGraw-Hill Inc, 1994).67.Caldwell, T. J., Chandra, S., Feher, K., Simmons, J. B. & Hogan, Z. Ecosystem response to earlier ice break-up date: Climate-driven changes to water temperature, lake-habitat-specific production, and trout habitat and resource use. Glob. Chang. Biol. 26, 5475–5491 (2020).ADS 
PubMed 
Article 
PubMed Central 

Google Scholar 
68.Elser, J. J., Luecke, C., Brett, M. T. & Goldman, C. R. Effects of food web compensation after manipulation of rainbow trout in an oligotrophic lake. Ecology 76, 52–69 (1995).Article 

Google Scholar 
69.Cohen, J. H. & Forward Jr., R. B. Zooplankton diel vertical migration-a review of proximate control. in Oceanography and marine biology: An annual review (eds. Gibson, R. N., Atkinson, R. J. A. & Gordon, J. D. M.) 89–122 (Taylor & Francis, 2009).70.Williamson, C. E., Fischer, J. M., Bollens, S. M., Overholt, E. P. & Breckenridgec, J. K. Toward a more comprehensive theory of zooplankton diel vertical migration: Integrating ultraviolet radiation and water transparency into the biotic paradigm. Limnol. Oceanogr. 56, 1603–1623 (2011).ADS 
Article 

Google Scholar 
71.Storz, U. C. & Paul, R. J. Phototaxis in water fleas (Daphnia magna) is differently influenced by visible and UV light. J. Comp. Physiol. Sens. Neural Behav. Physiol. 183, 709–717 (1998).Article 

Google Scholar 
72.National Interagency Fire Center. Total Wildland Fires and Acres (1983–2020). (2021). https://www.nifc.gov/fire-information/statistics/wildfires.73.MTBS. https://www.mtbs.gov/. MTBS (2020). https://www.mtbs.gov/. More

1.Maslo, B. & Fefferman, N. H. A case study of bats and white-nose syndrome demonstrating how to model population viability with evolutionary effects. Conserv. Biol. 29, 1176–1185 (2015).PubMed 
Article 
PubMed Central 

Google Scholar 
2.Morris, W. F. & Doak, D. F. Quantitative Conservation Biology (Sinauer, Sunderland, 2002).
Google Scholar 
3.Liebhold, A. & Bascompte, J. The Allee effect, stochastic dynamics and the eradication of alien species. Ecol. Lett. 6, 133–140 (2003).Article 

Google Scholar 
4.Stephens, P. A. & Sutherland, W. J. Consequences of the Allee effect for behaviour, ecology and conservation. Trends Ecol. Evol. 14, 401–405 (1999).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
5.Nunney, L. & Elam, D. R. Estimating the effective population size of conserved populations. Conserv. Biol. 8, 175–184 (1994).Article 

Google Scholar 
6.Lande, R. & Barrowclough, G. Effective population size, genetic variation, and their use in population. Viable populations for conservation, 87 (1987).7.Frankham, R. Effective population size/adult population size ratios in wildlife: A review. Genet. Res. 66, 95–107 (1995).Article 

Google Scholar 
8.Tallmon, D. A., Luikart, G. & Waples, R. S. The alluring simplicity and complex reality of genetic rescue. Trends Ecol. Evol. 19, 489–496 (2004).PubMed 
Article 
PubMed Central 

Google Scholar 
9.Whiteley, A. R., Fitzpatrick, S. W., Funk, W. C. & Tallmon, D. A. Genetic rescue to the rescue. Trends Ecol. Evol. 30, 42–49 (2015).PubMed 
Article 
PubMed Central 

Google Scholar 
10.Jiao, J., Gilchrist, M. A. & Fefferman, N. H. The impact of host metapopulation structure on short-term evolutionary rescue in the face of a novel pathogenic threat. Glob. Ecol. Conserv. 23, 01174 (2020).
Google Scholar 
11.Hanski, I. Metapopulation Ecology (Oxford University Press, Oxford, 1999).
Google Scholar 
12.Mortier, F., Jacob, S., Vandegehuchte, M. L. & Bonte, D. Habitat choice stabilizes metapopulation dynamics by enabling ecological specialization. Oikos 128, 529–539 (2019).Article 

Google Scholar 
13.Jiao, J., Riotte-Lambert, L., Pilyugin, S. S., Gil, M. A. & Osenberg, C. W. Mobility and its sensitivity to fitness differences determine consumer–resource distributions. R. Soc. Open Sci. 7, 200247 (2020).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
14.Anderson, S. C., Moore, J. W., McClure, M. M., Dulvy, N. K. & Cooper, A. B. Portfolio conservation of metapopulations under climate change. Ecol. Appl. 25, 559–572 (2015).PubMed 
Article 
PubMed Central 

Google Scholar 
15.Case, T. J. Invasion resistance, species build-up and community collapse in metapopulation models with interspecies competition. Biol. J. Lin. Soc. 42, 239–266 (1991).Article 

Google Scholar 
16.Gyllenberg, M. & Hanski, I. Habitat deterioration, habitat destruction, and metapopulation persistence in a heterogenous landscape. Theor. Popul. Biol. 52, 198–215 (1997).CAS 
PubMed 
MATH 
Article 
PubMed Central 

Google Scholar 
17.Jiao, J., Pilyugin, S. S. & Osenberg, C. W. Random movement of predators can eliminate trophic cascades in marine protected areas. Ecosphere 7, e01421 (2016).Article 

Google Scholar 
18.Nee, S. & May, R. M. Dynamics of metapopulations: Habitat destruction and competitive coexistence. J. Anim. Ecol. 61, 37–40 (1992).Article 

Google Scholar 
19.Ying, Y., Chen, Y., Lin, L. & Gao, T. Risks of ignoring fish population spatial structure in fisheries management. Can. J. Fish. Aquat. Sci. 68, 2101–2120 (2011).Article 

Google Scholar 
20.Hess, G. Disease in metapopulation models: Implications for conservation. Ecology 77, 1617–1632 (1996).Article 

Google Scholar 
21.Daszak, P., Cunningham, A. A. & Hyatt, A. D. Emerging infectious diseases of wildlife: Threats to biodiversity and human health. Science 287, 443–449 (2000).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
22.Harding, K. C., Begon, M., Eriksson, A. & Wennberg, B. Increased migration in host–pathogen metapopulations can cause host extinction. J. Theor. Biol. 298, 1–7 (2012).MathSciNet 
CAS 
PubMed 
MATH 
Article 
PubMed Central 

Google Scholar 
23.Dowling, A. J., Hill, G. E. & Bonneaud, C. Multiple differences in pathogen-host cell interactions following a bacterial host shift. Sci. Rep. 10, 1–12 (2020).Article 
CAS 

Google Scholar 
24.Kuzmin, I. V. et al. Molecular inferences suggest multiple host shifts of rabies viruses from bats to mesocarnivores in Arizona during 2001–2009. PLoS Pathog 8, e1002786 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
25.Levine, R. S. et al. Supersuppression: Reservoir competency and timing of mosquito host shifts combine to reduce spillover of West Nile virus. Am. J. Trop. Med. Hyg. 95, 1174–1184 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
26.Langwig, K. E. et al. Context-dependent conservation responses to emerging wildlife diseases. Front. Ecol. Environ. 13, 195–202 (2015).Article 

Google Scholar 
27.Smith, K. F., Acevedo-Whitehouse, K. & Pedersen, A. B. The role of infectious diseases in biological conservation. Anim. Conserv. 12, 1–12 (2009).Article 

Google Scholar 
28.Xiao, Y., Tang, B., Wu, J., Cheke, R. A. & Tang, S. Linking key intervention timing to rapid decline of the COVID-19 effective reproductive number to quantify lessons from mainland China. Int. J. Infect. Dis. 97, 296–298 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
29.Cintrón-Arias, A., Castillo-Chávez, C., Betencourt, L., Lloyd, A. L. & Banks, H. T. The Estimation of the Effective Reproductive Number from Disease Outbreak Data. (North Carolina State University, Center for Research in Scientific Computation, 2008).30.Salpeter, E. E. & Salpeter, S. R. Mathematical model for the epidemiology of tuberculosis, with estimates of the reproductive number and infection-delay function. Am. J. Epidemiol. 147, 398–406 (1998).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
31.Grenfell, B. & Harwood, J. (Meta) population dynamics of infectious diseases. Trends Ecol. Evol. 12, 395–399 (1997).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
32.Millard, A. R., Roberts, C. A. & Hughes, S. S. Isotopic evidence for migration in Medieval England: The potential for tracking the introduction of disease. Soc. Biol. Human Affairs. 70, 9–13 (2005).
Google Scholar 
33.Chen, M. et al. The introduction of population migration to SEIAR for COVID-19 epidemic modeling with an efficient intervention strategy. Inf. Fusion 64, 252–258 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
34.Reed, K. D., Meece, J. K., Henkel, J. S. & Shukla, S. K. Birds, migration and emerging zoonoses: West Nile virus, Lyme disease, influenza A and enteropathogens. Clin. Med. Res. 1, 5–12 (2003).PubMed 
PubMed Central 
Article 

Google Scholar 
35.Roy, B. & Kirchner, J. Evolutionary dynamics of pathogen resistance and tolerance. Evolution 54, 51–63 (2000).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
36.Bliven, K. A. & Maurelli, A. T. Antivirulence genes: Insights into pathogen evolution through gene loss. Infect. Immun. 80, 4061–4070 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
37.Wild, G., Gardner, A. & West, S. A. Adaptation and the evolution of parasite virulence in a connected world. Nature 459, 983–986 (2009).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
38.Laine, A. L. Resistance variation within and among host populations in a plant–pathogen metapopulation: Implications for regional pathogen dynamics. J. Ecol. 92, 990–1000 (2004).Article 

Google Scholar 
39.Thrall, P. H. et al. Rapid genetic change underpins antagonistic coevolution in a natural host-pathogen metapopulation. Ecol. Lett. 15, 425–435 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
40.Juhas, M. Horizontal gene transfer in human pathogens. Crit. Rev. Microbiol. 41, 101–108 (2015).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
41.Soanes, D. & Richards, T. A. Horizontal gene transfer in eukaryotic plant pathogens. Annu. Rev. Phytopathol. 52, 583–614 (2014).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
42.Brunham, R. C., Plummer, F. A. & Stephens, R. S. Bacterial antigenic variation, host immune response, and pathogen-host coevolution. Infect. Immun. 61, 2273 (1993).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
43.Sasaki, A. Evolution of antigen drift/switching: Continuously evading pathogens. J. Theor. Biol. 168, 291–308 (1994).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
44.Lange, A. & Ferguson, N. M. Antigenic diversity, transmission mechanisms, and the evolution of pathogens. PLoS Comput. Biol. 5, 1000536 (2009).ADS 
MathSciNet 
Article 
CAS 

Google Scholar 
45.Alizon, S., Hurford, A., Mideo, N. & Van Baalen, M. Virulence evolution and the trade-off hypothesis: History, current state of affairs and the future. J. Evol. Biol. 22, 245–259 (2009).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
46.Messenger, S. L., Molineux, I. J. & Bull, J. Virulence evolution in a virus obeys a trade off. Proc. R. Soc. Lond. B 266, 397–404 (1999).CAS 
Article 

Google Scholar 
47.Alizon, S., de Roode, J. C. & Michalakis, Y. Multiple infections and the evolution of virulence. Ecol. Lett. 16, 556–567 (2013).PubMed 
Article 
PubMed Central 

Google Scholar 
48.Bull, J. J. & Lauring, A. S. Theory and empiricism in virulence evolution. PLoS Pathog 10, e1004387 (2014).49.Gray, M. J. & Chinchar, V. G. Ranaviruses: Lethal Pathogens of Ectothermic Vertebrates (Springer Science+ Business Media, New York, 2015).Book 

Google Scholar 
50.Dobbelaere, T., Muller, E. M., Gramer, L. J., Holstein, D. M. & Hanert, E. Coupled epidemio-hydrodynamic modeling to understand the spread of a deadly coral disease in Florida. Front. Mar. Sci. 7, 1016 (2020).Article 

Google Scholar 
51.Stoddard, S. T. et al. House-to-house human movement drives dengue virus transmission. Proc. Natl. Acad. Sci. USA 110, 994–999 (2013).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
52.Diekmann, O., Heesterbeek, J. & Roberts, M. G. The construction of next-generation matrices for compartmental epidemic models. J. R. Soc. Interface 7, 873–885 (2010).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
53.Smith, K. F., Sax, D. F. & Lafferty, K. D. Evidence for the role of infectious disease in species extinction and endangerment. Conserv. Biol. 20, 1349–1357 (2006).PubMed 
Article 
PubMed Central 

Google Scholar 
54.Anderson, R. M., Anderson, B. & May, R. M. Infectious Diseases of Humans: Dynamics and Control (Oxford University Press, Oxford, 1992).
Google Scholar 
55.O’Brien, S. J. et al. Genetic basis for species vulnerability in the cheetah. Science 227, 1428–1434 (1985).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
56.Ingvarsson, P. K. & Lundberg, S. The effect of a vector-borne disease on the dynamics of natural plant populations: A model for Ustilago violacea infection of Lychnis viscaria. J. Ecol. 81, 263–270 (1993).Article 

Google Scholar 
57.Carlson, S. M., Cunningham, C. J. & Westley, P. A. Evolutionary rescue in a changing world. Trends Ecol. Evol. 29, 521–530 (2014).PubMed 
Article 
PubMed Central 

Google Scholar 
58.Gonzalez, A., Ronce, O., Ferriere, R. & Hochberg, M. E. (The Royal Society, 2013).59.Fine, P. E. Herd immunity: History, theory, practice. Epidemiol. Rev. 15, 265–302 (1993).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
60.Fontanet, A. & Cauchemez, S. COVID-19 herd immunity: Where are we?. Nat. Rev. Immunol. 20, 583–584 (2020).PubMed 
Article 
CAS 
PubMed Central 

Google Scholar 
61.Fine, P., Eames, K. & Heymann, D. L. “Herd immunity”: A rough guide. Clin. Infect. Dis. 52, 911–916 (2011).PubMed 
Article 
PubMed Central 

Google Scholar 
62.Barbarossa, M. V. & Röst, G. Immuno-epidemiology of a population structured by immune status: A mathematical study of waning immunity and immune system boosting. J. Math. Biol. 71, 1737–1770 (2015).MathSciNet 
CAS 
PubMed 
MATH 
Article 
PubMed Central 

Google Scholar 
63.Hamami, D., Cameron, R., Pollock, K. G. & Shankland, C. Waning immunity is associated with periodic large outbreaks of mumps: A mathematical modeling study of Scottish data. Front. Physiol. 8, 233 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
64.Klepac, P. & Caswell, H. The stage-structured epidemic: Linking disease and demography with a multi-state matrix approach model. Thyroid Res. 4, 301–319 (2011).
Google Scholar 
65.Anderson, R. M. & May, R. M. Population biology of infectious diseases: Part I. Nature 280, 361–367 (1979).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
66.Colizza, V. & Vespignani, A. Epidemic modeling in metapopulation systems with heterogeneous coupling pattern: Theory and simulations. J. Theor. Biol. 251, 450–467 (2008).MathSciNet 
PubMed 
MATH 
Article 
PubMed Central 

Google Scholar 
67.Stone, C. M., Schwab, S. R., Fonseca, D. M. & Fefferman, N. H. Human movement, cooperation and the effectiveness of coordinated vector control strategies. J. R. Soc. Interface 14, 20170336 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
68.Jiao, J., Pilyugin, S. S., Riotte-Lambert, L. & Osenberg, C. W. Habitat-dependent movement rate can determine the efficacy of marine protected areas. Ecology 99, 2485–2495 (2018).PubMed 
Article 
PubMed Central 

Google Scholar 
69.Keeling, M. J., Rohani, P. & Grenfell, B. T. Seasonally forced disease dynamics explored as switching between attractors. Physica D 148, 317–335 (2001).ADS 
MATH 
Article 

Google Scholar  More

1.Malthus, T. An Essay on the Principle of Population (J. Johnson Publisher, London, 1798).
Google Scholar 
2.Nicholson, A. J. The balance of animal populations. J. Anim. Ecol. 2(1), 132–588 (1933).Article 

Google Scholar 
3.Andrewartha, H. G. & Birch, L. C. The distribution and abundance of animals (University of Chicago Press, 1954).
Google Scholar 
4.Turchin, P. Complex Population Dynamics: A Theoretical/Empirical Synthesis. Monographs in Population Biology Vol. 35 (Princeton University Press, 2003).MATH 

Google Scholar 
5.Bellows, T. S. & Hassell, M. P. Theories and mechanisms of natural population regulation. In Handbook of Biological Control (eds Bellows, T. S. & Fisher, T. W.) 17–44 (Academic Press, 1999).
Google Scholar 
6.Cock, M. J. W. et al. Do new access and benefit sharing procedures under the convention on biological diversity threaten the future of biological control?. Biocontrol 55, 199–218. https://doi.org/10.1007/s10526-009-9234-9 (2010).Article 

Google Scholar 
7.Biondi, A., Guedes, R. N. C., Wan, F. H. & Desneux, N. Ecology, worldwide spread, and management of the invasive South American tomato pinworm, Tuta absoluta: past, present, and future. Annu. Rev. Entomol. 63, 239–258. https://doi.org/10.1146/annurev-ento-031616-034933 (2018).CAS 
Article 
PubMed 

Google Scholar 
8.Ferracini, C. et al. Natural enemies of Tuta absoluta in the Mediterranean basin, Europe and South America. Biocontrol Sci. Technol. 29, 578–609. https://doi.org/10.1080/09583157.2019.1572711 (2019).Article 

Google Scholar 
9.Guedes, N.C., Picanco M. Tuta absoluta in South America: pest status, management & insecticide resistance. Proceedings of the EPPO/IOBC/FAO/NEPPO Joint International Symposium on Management of Tuta absoluta (tomato borer). Agadir, Marocco, Nov. 16–18, 2011, 15–16 (2011).10.Urbaneja, A., Monton, H. & Mollá, O. Suitability of the tomato borer Tuta absoluta as prey for Macrolophus pygmaeus and Nesidiocoris tenuis. J. Appl. Entomol. 133, 292–296 (2009).Article 

Google Scholar 
11.Parra, J. R. P. & Zucchi, R. A. Trichogramma in Brazil: feasibility of use after 20 years of research. Neotrop. Entomol. 33, 271–281 (2004).Article 

Google Scholar 
12.Pérez-Hedo, M. & Urbaneja, A. The zoophytophagous predator Nesidiocoris tenuis: a successful but controversial biocontrol agent in tomato crops. In Advances in Insect Control and Resistance Management (eds Horowitz, A. R. & Ishaaya, I.) 121–138 (Springer, Dordrecht, 2016).
Google Scholar 
13.Mollá, O., Biondi, A., Alonso-Valiente, M. & Urbaneja, A. A comparative life history study of two mirid bugs preying on Tuta absoluta and Ephestia kuehniella eggs on tomato crops: implications for biological control. Biocontrol 59, 175–183. https://doi.org/10.1007/s10526-013-9553-8 (2014).Article 

Google Scholar 
14.Bajonero, J.G. Tuta absoluta (Meyrick, 1917) (Lepidoptera: Gelechiidae): adequação de uma dieta artificial e avaliação do seu controle biológico com Trichogramma pretiosum Riley em tomateiro. PhD Thesis, Universidade de São Paulo Escola Superior de Agricultura “Luiz de Queiroz”, Piracicaba, Sao Paolo, Brazil, p. 87 (2016).15.Calvo, F. J., Lorente, M. J., Stansly, P. A. & Belda, J. E. Pre-plant release of Nesidiocoris tenuis and supplementary tactics for control of Tuta absoluta and Bemisia tabaci in greenhouse tomato. Entomol. Exp. Appl. 143, 111–119 (2012).Article 

Google Scholar 
16.van Lenteren, J. C. et al. Pest kill rate as aggregate evaluation criterion to rank biological control agents: a case study with Neotropical predators of Tuta absoluta on tomato. Bull. Entomol. Res. 109, 812–820. https://doi.org/10.1017/S0007485319000130 (2019).CAS 
Article 
PubMed 

Google Scholar 
17.Tommasini, M. G., van Lenteren, J. C. & Burgio, G. Biological traits and predation capacity of four Orius species on two prey species. Bull. Insectol. 57, 79–94 (2004).
Google Scholar 
18.van Lenteren, J. C. Ecology: Cool Science, But Does It Help? 44 (Wageningen University, Wageningen, 2010).
Google Scholar 
19.Biondi, A., Desneux, N., Amiens-Desneux, E., Siscaro, G. & Zappalà, L. Biology and developmental strategies of the Palaearctic parasitoid Bracon nigricans (Hymenoptera: Braconidae) on the Neotropical moth Tuta absoluta (Lepidoptera: Gelechiidae). J. Econ. Entomol. 106, 1638–1647 (2013).Article 

Google Scholar 
20.Bin, F., Vinson, S. B. Efficacy assessment in egg parasitoids (Hymenoptera): proposal for a unified terminology. In Trichogramma and other egg Parasitoids (eds. Wajnberg E. & Vinson S.B.). Proceedings 3rd International Symposium, San Antonio, Texas, pp. 175–179 (1990).21.Abram, P. K., Brodeur, J., Urbaneja, A. & Tena, A. Nonreproductive effects of insect parasitoids on their hosts. Annu. Rev. Entomol. 64, 259–276. https://doi.org/10.1146/annurev-ento-011118-111753 (2019).CAS 
Article 
PubMed 

Google Scholar 
22.Chailleux, A., Droui, A., Bearez, P. & Desneux, N. Survival of a specialist natural enemy experiencing resource competition with an omnivorous predator when sharing the invasive prey Tuta absoluta. Ecol. Evol. 7, 8329–8337 (2017).Article 

Google Scholar 
23.Han, P. et al. Bottom-up efects of irrigation, fertilization and plant resistance on Tuta absoluta: implications for Integrated Pest Management. J. Pest Sci. 92, 1359–1370. https://doi.org/10.1007/s10340-018-1066-x (2019).Article 

Google Scholar 
24.Calvo, F. J., Bolckmans, K. & Belda, J. E. Release rate for a pre-plant application of Nesidiocoris tenuis for Bemisia tabaci control in tomato. Biocontrol 57, 809–817 (2012).Article 

Google Scholar 
25.van Lenteren, J. C., Bueno, V. H. P., Calvo, F. J., Calixto, A. M. & Montes, F. C. Comparative effectiveness and injury to tomato plants of three Neotropical mirid predators of Tuta absoluta (Lepidoptera: Gelechiidae). J. Econ. Entomol. 111, 1080–1086. https://doi.org/10.1093/jee/toy057 (2018).Article 
PubMed 

Google Scholar 
26.Calvo, F. J., Soriano, J. D., Stansly, P. A. & Belda, J. E. Can the parasitoid Necremnus tutae (Hymenoptera: Eulophidae) improve existing biological control of the tomato leafminer Tuta aboluta (Lepidoptera: Gelechiidae). Bull. Entomol. Res. 406, 502–511. https://doi.org/10.1017/S0007485316000183 (2016).CAS 
Article 

Google Scholar 
27.Crisol-Marinez, E. & van der Blom, J. Necremnus tutae (Hymenoptera, Eulophidae) is widespread and efficiently controls Tuta absoluta in tomato greenhouses in SE Spain. IOBC/WPRS Bull. 147, 22–29 (2019).
Google Scholar 
28.Castañé, C., van der Blom, J. & Nicot, P. C. Tomatoes. In Integrated Pest And Disease Management In Greenhouse Crops (eds Gullino, M. L. et al.) 487–511 (Springer, Switzerland, 2020). https://doi.org/10.1007/978-3-030-22304-5_17.
Google Scholar 
29.Knapp, M., Palevsky, E. & Rapisarda, C. Insect and mite pests. In Integrated Pest and Disease Management in Greenhouse Crops (eds Gullino, M. L. et al.) 101–144 (Springer, Switzerland, 2020). https://doi.org/10.1007/978-3-030-22304-5_17.
Google Scholar 
30.van Lenteren, J. C., Alomar, O., Ravensberg, W. J. & Urbaneja, A. Biological control agents for control of pests in greenhouses. In Integrated Pest And Disease Management In Greenhouse Crops (eds Gullino, M. L. et al.) 409–440 (Springer, Switzerland, 2020).
Google Scholar 
31.Pratissoli, D. & de Carvalho, J.D. Guia de Campo: Pragas da Cultura do Tomateiro. Alegre, ES: NUDEMAFI, Centro de Ciências Agrárias, UFES, 35pp. (Série Técnica/NUDEMAFI, ISSN 2359-4179; 1) (2015).32.Bodino, N., Ferracini, C. & Tavella, L. Functional response and age-specific foraging behaviour of Necremnus tutae and N. cosmopterix, native natural enemies of the invasive pest Tuta absoluta in Mediterranean area. J. Pest Sci. 92, 1467–1478. https://doi.org/10.1007/s10340-018-1025-6 (2019).Article 

Google Scholar 
33.Pérez-Hedo, M., Riahi, C. & Urbaneja, A. Use of zoophytophagous mirid bugs in horticultural crops: current challenges and future perspectives. Pest Manag. Sci. 77, 33–42 (2021).Article 

Google Scholar 
34.Wheeler, A. G. Jr. & Krimmel, B. A. Mirid (Hemiptera: Heteroptera) specialists of sticky plants: adaptations, interactions and ecological implications. Annu. Rev. Entomol. 60, 393–414 (2015).CAS 
Article 

Google Scholar 
35.Bueno, V. H. P., Lins, J. C., Silva, D. B. & van Lenteren, J. C. Is predation of Tuta absoluta by three Neotropical mirid predators affected by tomato lines with different densities in glandular trichomes?. Arthropod-Plant Int. 13, 41–48. https://doi.org/10.1007/s11829-018-9658-1 (2019).Article 

Google Scholar 
36.Pérez-Hedo, M., Arias-Sanguino, A. M. & Urbaneja, A. Induced tomato plant resistance against Tetranychus urticae triggered by the phytophagy of Nesidiocoris tenuis. Front. Plant Sci. 9, 1419. https://doi.org/10.3389/fpls.2018.01419 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
37.Barra-Bucarei, L., Devotto Moreno, L. & Iglesis, A. F. Biological control in Chile. In Biological Control in Latin America and the Caribbean: Its Rich History and Bright Future (eds van Lenteren, J. C. et al.) 108–123 (CAB International, Wallingford, 2020).
Google Scholar 
38.López, S. N., OrozcoMuñoz, A., Andorno, A. V., Cuello, E. M. & Cagnotti, C. L. Predatory capacity of Tupiocoris cucurbitaceus (Hemiptera Miridae) on several pests of tomato. Bull. Insectol. 72, 201–205 (2019).
Google Scholar 
39.Fauvel, G., Malausa, J. C. & Kaspar, B. Etude en laboratoire des principales caracteristiques biologiques de Macrolophus caliginosus (Heteroptera: Miridae). Entomophaga 32, 529–543 (1987).Article 

Google Scholar 
40.Mollá, O., Montón, H., Vanaclocha, P., Beitia, F. & Urbaneja, A. Predation by the mirids Nesidiocoris tenuis and Macrolophus pygmaeus on the tomato borer Tuta absoluta. IOBC/WPRS Bull. 49, 203–208 (2009).
Google Scholar 
41.Sánchez, J. A., Lacasa, A., Arnó, J., Castañé, C. & Alomar, O. Life history parameters for Nesidiocoris tenuis (Reuter) (Het., Miridae) under different temperature regimes. J. Appl. Entomol. 133, 125–132. https://doi.org/10.1111/j.1439-0418.2008.01342.x (2009).Article 

Google Scholar 
42.Sánchez, J. A., La-Spina, M. & Lacasa, A. Numerical response of Nesidiocoris tenuis (Hemiptera: Miridae) preying on Tuta absoluta (Lepidoptera: Gelechiidae) in tomato crops. Eur. J. Entomol. 11, 387–395 (2014).Article 

Google Scholar 
43.Arnó, J. et al. Tuta absoluta, a new pest in IPM tomatoes in the northeast of Spain. IOBC/WPRS Bull. 49, 203–208 (2009).
Google Scholar 
44.Mahdavi, T. S., Madadi, H. & Biondi, A. Predation and reproduction of the generalist predator Nabis pseudoferus preying on Tuta absoluta. Entomol. Exp. Appl. 168, 732–741. https://doi.org/10.1111/eea.12975 (2020).Article 

Google Scholar 
45.Salas Gervassio, N. G., Aquino, D., Vallina, C., Biondi, A. & Luna, M. G. A re-examination of Tuta absoluta parasitoids in South America for optimized biological control. J. Pest Sci. 92, 1343–1357. https://doi.org/10.1007/s10340-018-01078-1 (2019).Article 

Google Scholar 
46.Idriss, G. E. A., Mohamed, S. A., Khamis, F., Du Plessis, H. & Ekesi, S. Biology and performance of two indigenous larval parasitoids on Tuta absoluta (Lepidoptera: Gelechiidae) in Sudan. Biocontrol Sci. Technol. 28, 614–628. https://doi.org/10.1080/09583157.2018.1477117 (2018).Article 

Google Scholar 
47.Cagnotti, C. L., Riquelme Virgala, M., Botto, E. N. & López, S. N. Dispersion and persistence of Trichogrammatoidea bactrae (Nagaraja) over Tuta absoluta (Meyrick), in tomato greenhouses. Neotrop. Entomol. 47, 553–559. https://doi.org/10.1007/s13744-017-0573-4 (2018).CAS 
Article 
PubMed 

Google Scholar 
48.Pratissoli, D. Bioecologia de Trichogramma pretiosum Riley, 1879, nas traças, Scrobipalpuloides absoluta (Meyrick, 1917) e Phthorimaea operculella (Zeller, 1873), em tomateiro. Piracicaba: Doutorado – Escola Superior de Agricultura “Luiz de Queiroz”/USP, p. 153 (1995).49.Pratissoli, D. & Parra, J. R. P. Fertility life table of Trichogramma pretiosum (Hym., Trichogrammatidae) in eggs of Tuta absoluta and Phthorimaea operculella (Lep., Gelechiidae) at different temperatures. J. Appl. Entomol. 124, 339–342 (2000).Article 

Google Scholar 
50.Riquelme Virgala, M. B. & Botto, E. N. Estudios biológicos de Trichogrammatoidea bactrae Nagaraja (Hymenoptera: Trichogrammatidae), parasitoide de huevos de Tuta absoluta Meyrick (Lepidoptera: Gelechiidae). Neotrop. Entomol. 36, 612–617 (2010).Article 

Google Scholar 
51.Aigbedion-Atalor, P.O, Abuelgasim Mohamed, S., Hill, M.P., Zalucki, M.P., Azrag, A.G.A., Srinivasan, R. & Ekesi, S. Host stage preference and performance of Dolichogenidea gelechiidivoris (Hymenoptera: Braconidae), a candidate for classical biological control of Tuta absoluta in Africa. Biol. Control. Preprint at https://doi.org/https://doi.org/10.1016/j.biocontrol.2020.104215 (2020).52.Guleria, P., Sharma, P. L., Verma, S. C. & Chandel, R. S. Life history traits and host-killing rate of Neochrysocharis formosa on Tuta absoluta. Biocontrol 65, 401–411. https://doi.org/10.1007/s10526-020-10016-z (2020).CAS 
Article 

Google Scholar 
53.Calvo, F. J., Soriano, J. D., Bolckmans, K. & Belda, J. E. Host instar suitability and life-history parameters under different temperature regimes of Necremnus artynes on Tuta absoluta. Biocontrol Sci. Technol. 23, 803–815 (2013).Article 

Google Scholar 
54.Nieves, E. L., Pereyra, P. C., Luna, M. G., Medone, P. & Sánchez, N. E. Laboratory population parameters and field impact of the larval endoparasitoid Pseudapanteles dignus (Hymenoptera: Braconidae) on its host Tuta absoluta (Lepidoptera: Gelechiidae) in tomato crops in Argentina. J. Econ. Entomol. 108, 1553–1559 (2015).Article 

Google Scholar 
55.Luna, M. G., Wada, V. & Sánchez, N. E. Biology of Dineulophus phtorimaeae (Hymenoptera: Eulophidae), and field interaction with Pseudapanteles dignus (Hymenoptera: Braconidae), larval parasitoids of Tuta absoluta (Lepidoptera: Gelechiidae) in tomato. Ann. Entomol. Soc. Am. 106, 936–942 (2010).Article 

Google Scholar 
56.Savino, V., Coviella, C. E. & Luna, M. G. Reproductive biology and functional response of Dineulophus phtorimaeae a natural enemy of the tomato moth Tuta absoluta. J. Insect Sci. 12, 1–14 (2012).Article 

Google Scholar 
57.Birch, L. C. The intrinsic rate of natural increase of an insect population. J. Anim. Ecol. 17, 15–26 (1948).Article 

Google Scholar 
58.Lotka, A. J. Relation between birth rates and death rates. Science 26, 21–22 (1907).ADS 
CAS 
Article 

Google Scholar 
59.Lotka, A. J. & Sharpe, F. R. A problem in age distribution. Philos. Mag. 6(21), 339–345 (1911).MATH 

Google Scholar 
60.Dalgaard, P. Introductory Statistics With R 2nd edn. (Springer, New York, 2008).Book 

Google Scholar  More

1.Ryther, J. H. The Sargasso Sea. Sci. Am. 194, 98–108 (1956).ADS 
Article 

Google Scholar 
2.Littler, D. S. & Littler, M. M. Caribbean Reef Plants (Offshore Graphics, 2000).3.Winge, O. The Sargasso Sea, Its Boundaries and Vegetation In Report of the Danish Oceanographic Expedition, Vol. III, 1908–1910, (Copenhagen: Andr. Fred. Hòst and Sòn) 34 pp. Miscellaneous Paper Number 2. (1923).4.Parr, A. E. Quantitative observations on the pelagic Sargassum vegetation of the western North Atlantic. Bull. Bingham Oceanogr. Collect. 6, 1–94 (1939).
Google Scholar 
5.Lapointe, B. E. A comparison of nutrient-limited productivity in Sargassum natans from neritic vs. oceanic waters of the western North Atlantic Ocean. Limnol. Oceanogr. 40, 625–633 (1995).ADS 
CAS 
Article 

Google Scholar 
6.Lapointe, B. E., West, L. E., Sutton, T. T. & Hu, C. Ryther revisited: nutrient excretions by fishes enhance productivity of pelagic Sargassum in the western North Atlantic Ocean. J. Exp. Mar. Biol. Ecol. 458, 46–56 (2014).CAS 
Article 

Google Scholar 
7.Gower, J., Hu, C., Borstad, G. & King, S. Ocean color satellites show extensive lines of floating Sargassum in the Gulf of Mexico. IEEE Trans. Geosci. Remote Sens. 44, 3619–3625 (2006).ADS 
Article 

Google Scholar 
8.Williams, A., Feagin, R. & Stafford, A. W. Environmental impacts of beach raking of Sargassum spp. on Galveston Island, TX. Shore Beach 76, 63–69 (2008).
Google Scholar 
9.Moritsugu, K. Tampa Bay Times (Times Publishing Company, 1991).10.Turner, R. E. & Rabalais, N. N. Coastal eutrophication near the Mississippi river delta. Nature 368, 619–621 (1994).ADS 
Article 

Google Scholar 
11.Gower, J. F. R. & King, S. A. Distribution of floating Sargassum in the Gulf of Mexico and the Atlantic Ocean mapped using MERIS. Int. J. Remote Sens. 32, 1917–1929 (2011).ADS 
Article 

Google Scholar 
12.Johnson, D. R., Ko, D. S., Franks, J. S., Moreno, P. & Sanchez-Rubio, G. The Sargassum invasion of the Eastern Caribbean and dynamics of the Equatorial North Atlantic. In Proceedings of the 65th Annual Gulf and Caribbean Fisheries Institute Conference pp. 102–103 (2013). http://aquaticcommons.org/21444/1/GCFI_65-17.pdf.13.Gower, J., Young, E. & King, S. Satellite images suggest a new Sargassum source region in 2011. Remote Sens. Lett. 4, 764–773 (2013).Article 

Google Scholar 
14.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–102269 (2020).Article 

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

Google Scholar 
16.Djakouré, S., Araujo, M., Hounsou-Gbo, A., Noriega, C. & Bourlès, B. On the potential causes of the recent Pelagic Sargassum blooms events in the tropical North Atlantic Ocean. Biogeosci. Discuss. https://doi.org/10.5194/bg-2017-346 (2017).17.Oviatt, C. A., Huizenga, K., Rogers, C. S. & Miller, W. J. What nutrient sources support anomalous growth and the recent Sargassum mass stranding on Caribbean beaches? A review. Mar. Pollut. Bull. 145, 517–525 (2019).CAS 
PubMed 
Article 

Google Scholar 
18.McGillicuddy, D. J., Jr, Anderson, L. A., Doney, S. C. & Maltrud, M. E. Eddy‐driven sources and sinks of nutrients in the upper ocean: results from a 0.1 resolution model of the North Atlantic. Global Biogeochem. Cycles 17, 1035 (2003).19.Barkley, A. E. et al. African biomass burning is a substantial source of phosphorus deposition to the Amazon, Tropical Atlantic Ocean, and Southern Ocean. Proc. Natl Acad. Sci. USA 116, 16216–16221 (2019).ADS 
CAS 
PubMed 
Article 

Google Scholar 
20.Qi, L., Hu, C., Xing, Q. & Shang, S. Long-term trend of Ulva prolifera blooms in the western Yellow Sea. Harmful Algae 58, 35–44 (2016).PubMed 
Article 

Google Scholar 
21.Qi, L., Hu, C., Wang, M., Shang, S. & Wilson, C. Floating algae blooms in the East China Sea. Geophys. Res. Lett. 44, 501–511,509 (2017).Article 
CAS 

Google Scholar 
22.Smetacek, V. & Zingone, A. Green and golden seaweed tides on the rise. Nature 504, 84–88 (2013).ADS 
CAS 
PubMed 
Article 

Google Scholar 
23.Van Tussenbroek, B. I. et al. Severe impacts of brown tides caused by Sargassum spp. on near-shore Caribbean seagrass communities. Mar. Pollut. Bull. 122, 272–281 (2017).PubMed 
Article 
CAS 

Google Scholar 
24.Alvarez-Filip, L., Estrada-Saldívar, N., Pérez-Cervantes, E., Molina-Hernández, A. & González-Barrios, F. J. A rapid spread of the stony coral tissue loss disease outbreak in the Mexican Caribbean. PeerJ 7, e8069–e8069 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
25.Cabanillas-Terán, N., Hernández-Arana, H. A., Ruiz-Zárate, M.-Á., Vega-Zepeda, A. & Sanchez-Gonzalez, A. Sargassum blooms in the Caribbean alter the trophic structure of the sea urchin Diadema antillarum. PeerJ 7, e7589–e7589 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
26.Maurer, A. S., De Neef, E. & Stapleton, S. Sargassum accumulation may spell trouble for nesting sea turtles. Front. Ecol. Environ. 13, 394–395 (2015).Article 

Google Scholar 
27.Webster, R. K. & Linton, T. Development and implementation of Sargassum early advisory system (SEAS). Shore Beach 81, 1–1 (2013).
Google Scholar 
28.Resiere, D. et al. Sargassum seaweed on Caribbean islands: an international public health concern. Lancet 392, 2691–2691 (2018).Article 

Google Scholar 
29.Glibert, P. et al. The role of in the global proliferation of harmful algal blooms: new perspectives and approaches. Oceanography 18, 196–207 (2005).
Google Scholar 
30.Glibert, P. M. Eutrophication, harmful algae and biodiversity — Challenging paradigms in a world of complex nutrient changes. Mar. Pollut. Bull. 124, 591–606 (2017).CAS 
PubMed 
Article 

Google Scholar 
31.Steffen, W. et al. Planetary boundaries: guiding human development on a changing planet. Science 347, 6223 https://doi.org/10.1126/science.1259855 (2015).32.Ryther, J. H. The ecology of phytoplankton blooms in Moriches bay and Great South bay, Long Island, New York. Biol. Bull. 106, 198–209 (1954).Article 

Google Scholar 
33.Ryther, J. H. & Dunstan, W. M. Nitrogen, Phosphorus, and Eutrophication in the coastal marine environment. Science 171, 1008 LP-1013 (1971).Article 

Google Scholar 
34.Howarth, R. W. & Marino, R. Nitrogen as the limiting nutrient for eutrophication in coastal marine ecosystems: evolving views over three decades. Limnol. Oceanogr. 51, 364–376 (2006).ADS 
CAS 
Article 

Google Scholar 
35.Oelsner, G. P. & Stets, E. G. Recent trends in nutrient and sediment loading to coastal areas of the conterminous U.S.: insights and global context. Sci. Total Environ. 654, 1225–1240 (2019).ADS 
CAS 
PubMed 
Article 

Google Scholar 
36.Falkowski, P. G. Evolution of the nitrogen cycle and its influence on the biological sequestration of CO2 in the ocean. Nature 387, 272–275 (1997).ADS 
CAS 
Article 

Google Scholar 
37.Tyrrell, T. The relative influences of nitrogen and phosphorus on oceanic primary production. Nature 400, 525–531 (1999).ADS 
CAS 
Article 

Google Scholar 
38.Lapointe, B. E., Littler, M. M. & Littler, D. S. A comparison of nutrient-limited productivity in macroalgae from a Caribbean barrier reef and from a mangrove ecosystem. Aquat. Bot. 28, 243–255 (1987).Article 

Google Scholar 
39.Culliney, J. L. Measurements of reactive phosphorus associated with pelagic Sargassum in the Northwest Sargasso Sea1. Limnol. Oceanogr. 15, 304–305 (1970).ADS 
CAS 
Article 

Google Scholar 
40.Schaffelke, B. Particulate organic matter as an alternative nutrient source for tropical Sargassum species (Fucales, Phaeophyceae). J. Phycol. 35, 1150–1157 (1999).CAS 
Article 

Google Scholar 
41.Vonk, J. A., Middelburg, J. J., Stapel, J. & Bouma, T. J. Dissolved organic nitrogen uptake by seagrasses. Limnol. Oceanogr. 53, 542–548 (2008).ADS 
CAS 
Article 

Google Scholar 
42.Han, T., Qi, Z., Huang, H., Liao, X. & Zhang, W. Nitrogen uptake and growth responses of seedlings of the brown seaweed Sargassum hemiphyllum under controlled culture conditions. J. Appl. Phycol. 30, 507–515 (2018).CAS 
Article 

Google Scholar 
43.Fujita, R., Wheeler, P. & Edwards, R. Assessment of macroalgal nitrogen limitation in a seasonal upwelling region. Mar. Ecol. Prog. Ser. 53, 293–303 (1989).ADS 
Article 

Google Scholar 
44.Prospero, J. M. et al. in Nitrogen Cycling in the North Atlantic Ocean and its Watersheds (ed. Robert, W. H.) (Springer, 1996).45.Howarth, R. W. Coastal nitrogen pollution: a review of sources and trends globally and regionally. Harmful Algae 8, 14–20 (2008).CAS 
Article 

Google Scholar 
46.Rockström, J. & Karlberg, L. The quadruple squeeze: defining the safe operating space for freshwater use to achieve a triply green revolution in the Anthropocene. Ambio 39, 257–265 (2010).PubMed 
PubMed Central 
Article 

Google Scholar 
47.Hanisak, M. D. & Samuel, M. A. Twelfth International Seaweed Symposium (Springer, 1986).48.Rabalais, N. N. et al. Hypoxia in the northern Gulf of Mexico: does the science support the plan to reduce, mitigate, and control hypoxia? Estuar. Coasts 30, 753–772 (2007).CAS 
Article 

Google Scholar 
49.Tian, H. et al. Long-term trajectory of nitrogen loading and delivery from Mississippi river basin to the Gulf of Mexico. Glob. Biogeochem. Cycles 34, e2019GB006475–e002019GB006475 (2020).ADS 
CAS 
Article 

Google Scholar 
50.Lapointe, B. E., Brewton, R. A., Herren, L. W., Porter, J. W. & Hu, C. Nitrogen enrichment, altered stoichiometry, and coral reef decline at Looe Key, Florida Keys, USA: a 3-decade study. Mar. Biol. 166, 108–108 (2019).Article 
CAS 

Google Scholar 
51.Lapointe, B. E., Barile, P. J. & Littler, M. M. & Littler, D. S. Macroalgal blooms on southeast Florida coral reefs: II. Cross-shelf discrimination of nitrogen sources indicates widespread assimilation of sewage nitrogen. Harmful Algae 4, 1106–1122 (2005).CAS 
Article 

Google Scholar 
52.Dunn, D. E. Trends in Nutrient Inflows to the Gulf of Mexico from Streams Draining the Conterminous United States, 1972-93. Report No. 96-4113 (Austin, TX, 1996).53.Turner, R. E. & Rabalais, N. N. Changes in Mississippi River water quality this century: implications for coastal food webs. Bioscience 41, 140–147 (1991).Article 

Google Scholar 
54.Rabalais, N. N. et al. Nutrient changes in the Mississippi River and system responses on the adjacent continental shelf. Estuaries 19, 386–407 (1996).CAS 
Article 

Google Scholar 
55.Weber, S. C. et al. Amazon River influence on nitrogen fixation and export production in the western tropical North Atlantic. Limnol. Oceanogr. 62, 618–631 (2017).ADS 
Article 

Google Scholar 
56.Ryther, J. H., Menzel, D. W. & Corwin, N. Influence of Amazon River outflow on ecology of Western Tropical Atlantic. I. Hydrography and nutrient chemistry. J. Mar. Res. 25, 69–69 (1967).
Google Scholar 
57.Subramaniam, A. et al. Amazon River enhances diazotrophy and carbon sequestration in the tropical North Atlantic Ocean. Proc. Natl Acad. Sci.USA 105, 10460 LP–10410465 (2008).ADS 
Article 

Google Scholar 
58.Barichivich, J. et al. Recent intensification of Amazon flooding extremes driven by strengthened Walker circulation. Sci. Adv. 4, eaat8785–eaat8785 (2018).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
59.Howarth, R. W. et al. Regional nitrogen budgets and riverine N & P fluxes for the drainages to the North Atlantic Ocean: Natural and human influences. In Nitrogen Cycling in the North Atlantic Ocean and its Watersheds (ed. Robert, W. Howarth) (Springer, Dordrecht, 1996). https://doi.org/10.1007/978-94-009-1776-7_3.60.Galloway, J. N. et al. Regional nitrogen budgets and riverine N & P fluxes for the drainages to the North Atlantic Ocean: Natural and human influences. Biogeochemistry (ed. Robert, W. Howarth) 35, 181–226 (Springer, 1996).61.Gower, J. & King, S. Satellite images show the movement of floating Sargassum in the Gulf of Mexico and Atlantic Ocean. Nat. Preced. https://doi.org/10.1038/npre.2008.1894.1 (2008).62.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 
63.Zimmerman, R. C. & Kremer, J. N. Episodic nutrient supply to a kelp forest ecosystem in Southern California. J. Mar. Res. 42, 591–604 (1984).Article 

Google Scholar 
64.Kain, J. M. The seasons in the subtidal. Br. Phycol. J. 24, 203–215 (1989).Article 

Google Scholar 
65.Dorado, S., Rooker, J. R., Wissel, B. & Quigg, A. Isotope baseline shifts in pelagic food webs of the Gulf of Mexico. Mar. Ecol. Prog. Ser. 464, 37–49 (2012).ADS 
CAS 
Article 

Google Scholar 
66.Kendall, C., Elliott, E. M. & Wankel, S. D. Wiley Online Books 375-449 (2007).67.Altieri, K. E., Hastings, M. G., Peters, A. J., Oleynik, S. & Sigman, D. M. Isotopic evidence for a marine ammonium source in rainwater at Bermuda. Glob. Biogeochem. Cycles 28, 1066–1080 (2014).ADS 
CAS 
Article 

Google Scholar 
68.Bateman, A. S. & Kelly, S. D. Fertilizer nitrogen isotope signatures. Isotopes Environ. Health Stud. 43, 237–247 (2007).CAS 
PubMed 
Article 

Google Scholar 
69.Knapp, A. N., DiFiore, P. J., Deutsch, C., Sigman, D. M. & Lipschultz, F. Nitrate isotopic composition between Bermuda and Puerto Rico: implications for N2 fixation in the Atlantic Ocean. Global Biogeochem. Cycles 22, GB3014 https://doi.org/10.1029/2007GB003107 (2008).70.Knapp, A. N., Sigman, D. M. & Lipschultz, F. N isotopic composition of dissolved organic nitrogen and nitrate at the Bermuda Atlantic Time-series Study site. Global Biogeochem. Cycles 19, GB1018 https://doi.org/10.1029/2004GB002320 (2005).71.Montoya, J. P. Nitrogen stable isotopes in marine environments. Nitrogen Mar. Environ. 2, 1277–1302 (2008).Article 

Google Scholar 
72.Wissel, B. & Fry, B. Sources of particulate organic matter in the Mississippi River, USA. Large Rivers 15 105–118 (2003).73.Zaia Alves, G. H., Hoeinghaus, D. J., Manetta, G. I. & Benedito, E. Dry season limnological conditions and basin geology exhibit complex relationships with δ13C and δ15N of carbon sources in four Neotropical floodplains. PLoS ONE 12, e0174499 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
74.Smith, N. P. Upwelling in Atlantic shelf waters of South Florida. Florida Scientist 45, 125–138 (1982).75.Atkinson, L. P., O’Malley, P. G., Yoder, J. A. & Paffenhöfer, G. A. The effect of summertime shelf break upwelling on nutrient flux in southeastern United States continental shelf waters. J. Mar. Res. 42, 969–993 (1984).Article 

Google Scholar 
76.Subramaniam, A., Mahaffey, C., Johns, W. & Mahowald, N. Equatorial upwelling enhances nitrogen fixation in the Atlantic Ocean. Geophys. Res. Lett. 40, 1766–1771 (2013).ADS 
Article 

Google Scholar 
77.Carpenter, E. J. Nitrogen fixation by a blue-green epiphyte on Pelagic Sargassum. Science 178, 1207–1209 (1972).ADS 
CAS 
PubMed 
Article 

Google Scholar 
78.Phlips, E. J., Willis, M. & Verchick, A. Aspects of nitrogen fixation in Sargassum communities off the coast of Florida. J. Exp. Mar. Biol. Ecol. 102, 99–119 (1986).CAS 
Article 

Google Scholar 
79.Subramaniam, A., Montoya, J. P., Foster, R. A. & Capone, D. G. Nitrogen fixation in the eastern equatorial Atlantic: who and how much? European Geosciences Union General Assembly 11, 10156–10156 (2009).80.Carpenter, E. J. et al. Extensive bloom of a N2-fixing diatom/cyanobacterial association in the tropical Atlantic Ocean. Mar. Ecol. Prog. Ser. 185, 273–283 (1999).ADS 
CAS 
Article 

Google Scholar 
81.Zubkova, M., Boschetti, L., Abatzoglou, J. T. & Giglio, L. Changes in fire activity in Africa from 2002 to 2016 and their potential drivers. Geophys. Res. Lett. 46, 7643–7653 (2019).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
82.Baker, A. R., French, M. & Linge, K. L. Trends in aerosol nutrient solubility along a west–east transect of the Saharan dust plume. Geophys. Res. Lett. 33 L07805, https://doi.org/10.1029/2005GL024764 (2006).83.Baker, A. R., Jickells, T. D., Witt, M. & Linge, K. L. Trends in the solubility of iron, aluminium, manganese and phosphorus in aerosol collected over the Atlantic Ocean. Mar. Chem. 98, 43–58 (2006).CAS 
Article 

Google Scholar 
84.Shelley, R. U., Morton, P. L. & Landing, W. M. Elemental ratios and enrichment factors in aerosols from the US-GEOTRACES North Atlantic transects. Deep Sea Res. Part II Top. Stud. Oceanogr. 116, 262–272 (2015).ADS 
CAS 
Article 

Google Scholar 
85.Giglio, L., Schroeder, W. & Justice, C. O. The collection 6 MODIS active fire detection algorithm and fire products. Remote Sens. Environ. 178, 31–41 (2016).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
86.Giglio, L., Descloitres, J., Justice, C. O. & Kaufman, Y. J. An enhanced contextual fire detection algorithm for MODIS. Remote Sens. Environ. 87, 273–282 (2003).ADS 
Article 

Google Scholar 
87.Giglio, L., van der Werf, G. R., Randerson, J. T., Collatz, G. J. & Kasibhatla, P. Global estimation of burned area using MODIS active fire observations. Atmos. Chem. Phys. 6, 957–974 (2006).ADS 
CAS 
Article 

Google Scholar 
88.Roberts, G., Wooster, M. J. & Lagoudakis, E. Annual and diurnal african biomass burning temporal dynamics. Biogeosciences 6, 849–866 (2009).ADS 
Article 

Google Scholar 
89.Baker, A. R. & Jickells, T. D. Atmospheric deposition of soluble trace elements along the Atlantic Meridional Transect (AMT). Prog. Oceanogr. 158, 41–51 (2017).ADS 
Article 

Google Scholar 
90.Chance, R., Jickells, T. D. & Baker, A. R. Atmospheric trace metal concentrations, solubility and deposition fluxes in remote marine air over the south-east Atlantic. Mar. Chem. 177, 45–56 (2015).CAS 
Article 

Google Scholar 
91.Myriokefalitakis, S., Nenes, A., Baker, A. R., Mihalopoulos, N. & Kanakidou, M. Bioavailable atmospheric phosphorous supply to the global ocean: a 3-D global modeling study. Biogeosciences 13, 6519–6543 (2016).ADS 
CAS 
Article 

Google Scholar 
92.Kanakidou, M., Myriokefalitakis, S. & Tsigaridis, K. Aerosols in atmospheric chemistry and biogeochemical cycles of nutrients. Environ. Res. Lett. 13, 063004 (2018).93.Rosenzweig, M. L. Paradox of enrichment: destabilization of exploitation ecosystems in ecological time. Science 171, 385–387 (1971).ADS 
CAS 
PubMed 
Article 

Google Scholar 
94.McCann, K. S. et al. Landscape modification and nutrient‐driven instability at a distance. Ecol. Lett. 24, 398–414 (2021).PubMed 
Article 

Google Scholar 
95.Meybeck, M. Carbon, nitrogen, and phosphorus transport by world rivers. Am. J. Sci. 282, 401–450 (1982).ADS 
CAS 
Article 

Google Scholar 
96.Fanning, K. A. Nutrient provinces in the sea: concentration ratios, reaction rate ratios, and ideal covariation. J. Geophys. Res. Oceans 97, 5693–5712 (1992).ADS 
Article 

Google Scholar 
97.Ammerman, J. W., Hood, R. R., Case, D. A. & Cotner, J. B. Phosphorus deficiency in the Atlantic: an emerging paradigm in oceanography. Eos, Trans. Am. Geophys. Union 84, 165–170 (2003).ADS 
Article 

Google Scholar 
98.Lomas, M. W., Bonachela, J. A., Levin, S. A. & Martiny, A. C. Impact of ocean phytoplankton diversity on phosphate uptake. Proc. Natl Acad. Sci. USA 111, 17540–17545 (2014).ADS 
CAS 
PubMed 
Article 

Google Scholar 
99.Richey, J. E. et al. (ORNL Distributed Active Archive Center, 2008).100.Cochonneau, G. et al. The environmental observation and research project, ORE HYBAM, and the rivers of the Amazon basin. In Climate Variability and Change—Hydrological Impacts (eds Demuth, S. et al.) vol. 308, 44–50 (2006). More