Ecology

1.Nemergut DR, Schmidt SK, Fukami T, O’Neill SP, Bilinski TM, Stanish LF, et al. Patterns and processes of microbial community assembly. Microbiol Mol Biol Rev. 2013;77:342–56.PubMed 
PubMed Central 
Article 

Google Scholar 
2.Jousset A, Bienhold C, Chatzinotas A, Gallien L, Gobet A, Kurm V, et al. Where less may be more: how the rare biosphere pulls ecosystems strings. ISME J. 2017;11:853–62.PubMed 
PubMed Central 
Article 

Google Scholar 
3.Rivett DW, Bell T. Abundance determines the functional role of bacterial phylotypes in complex communities. Nat Microbiol. 2018;3:767–72.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
4.Bell T, Newman JA, Silverman BW, Turner SL, Lilley AK. The contribution of species richness and composition to bacterial services. Nature. 2005;436:1157–60.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
5.Starke R, Capek P, Morais D, Callister SJ, Jehmlich N. The total microbiome functions in bacteria and fungi. J Proteom. 2020;213:1–5.Article 
CAS 

Google Scholar 
6.Saleem M, Hu J, Jousset A. More than the sum of its parts: microbiome biodiversity as a driver of plant growth and soil health. Annu Rev Ecol Evol Syst. 2019;50:145–68.Article 

Google Scholar 
7.Wagg C, Schlaeppi K, Banerjee S, Kuramae EE, Heijden van der MGA. Fungal-bacterial diversity and microbiome complexity predict ecosystem functioning. Nat Commun. 2019;10:1–10.CAS 
Article 

Google Scholar 
8.Delgado-Baquerizo M, Maestre FT, Reich PB, Jeffries TC, Gaitan JJ, Encinar D, et al. Microbial diversity drives multifunctionality in terrestrial ecosystems. Nat Commun. 2016;7:1–8.Article 
CAS 

Google Scholar 
9.Delgado-Baquerizo M, Reich PB, Trivedi C, Eldridge DJ, Abades S, Alfaro FD, et al. Multiple elements of soil biodiversity drive ecosystem functions across biomes. Nat Ecol Evol. 2020;4:210–20.PubMed 
Article 
PubMed Central 

Google Scholar 
10.Aanderud ZT, Jones SE, Fierer N, Lennon JT. Resuscitation of the rare biosphere contributes to pulses of ecosystem activity. Front Microbiol. 2015;6:1–11.Article 

Google Scholar 
11.Song H-K, Song W, Kim M, Tripathi BM, Kim H, Jablonski P, et al. Bacterial strategies along nutrient and time gradients, revealed by metagenomic analysis of laboratory microcosms. FEMS Microbiol Ecol. 2017;93:1–13.Article 
CAS 

Google Scholar 
12.Jiao S, Chen W, Wei G. Biogeography and ecological diversity patterns of rare and abundant bacteria in oil-contaminated soils. Mol Ecol. 2017;26:5305–5317.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
13.Delgado-Baquerizo M, Oliverio AM, Brewer TE, Benavent-González A, Eldridge DJ, Bardgett RD, et al. A global atlas of the dominant bacteria found in soil. Science. 2018;359:320–5.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
14.Yu X, Polz MF, Alm EJ. Interactions in self-assembled microbial communities saturate with diversity. ISME J. 2019;13:1602–17.PubMed 
PubMed Central 
Article 

Google Scholar 
15.Li P, Liu J, Jiang C, Wu M, Liu M, Li Z. Distinct successions of common and rare bacteria in soil under humic acid amendment—a microcosm study. Front Microbiol. 2019;10:1–14.Article 

Google Scholar 
16.Nemergut DR, Costello EK, Hamady M, Lozupone C, Jiang L, Schmidt SK, et al. Global patterns in the biogeography of bacterial taxa. Environ Microbiol. 2011;13:135–44.PubMed 
PubMed Central 
Article 

Google Scholar 
17.Bickel S, Chen X, Papritz A, Or D. A hierarchy of environmental covariates control the global biogeography of soil bacterial richness. Sci Rep. 2019;9:1–10.CAS 
Article 

Google Scholar 
18.Clarke RT, Murphy JF. Effects of locally rare taxa on the precision and sensitivity of RIVPACS bioassessment of freshwaters. Freshw Biol. 2006;51:1924–40.Article 

Google Scholar 
19.Kurm V, Putten WH, van der, Boer W, de, Naus‐Wiezer S, Hol WHG. Low abundant soil bacteria can be metabolically versatile and fast growing. Ecology. 2017;98:555–64.PubMed 
Article 
PubMed Central 

Google Scholar 
20.Kurm V, Putten WH, van der, Hol WHG. Cultivation-success of rare soil bacteria is not influenced by incubation time and growth medium. PLoS ONE. 2019;14:1–14.Article 
CAS 

Google Scholar 
21.Meyer KM, Memiaghe H, Korte L, Kenfack D, Alonso A, Bohannan BJM. Why do microbes exhibit weak biogeographic patterns? ISME J. 2018;12:1404–13.PubMed 
PubMed Central 
Article 

Google Scholar 
22.Escalas A, Hale L, Voordeckers JW, Yang Y, Firestone MK, Alvarez‐Cohen L, et al. Microbial functional diversity: from concepts to applications. Ecol Evol. 2019;9:12000–16.PubMed 
PubMed Central 
Article 

Google Scholar 
23.Barberán A, Ramirez KS, Leff JW, Bradford MA, Wall DH, Fierer N. Why are some microbes more ubiquitous than others? Predicting the habitat breadth of soil bacteria. Ecol Lett. 2014;17:794–802.PubMed 
Article 
PubMed Central 

Google Scholar 
24.Dee LE, Cowles J, Isbell F, Pau S, Gaines SD, Reich PB. When do ecosystem services depend on rare species? Trends Ecol Evol. 2019;34:746–58.PubMed 
Article 
PubMed Central 

Google Scholar 
25.Pueyo S, He F, Zillio T. The maximum entropy formalism and the idiosyncratic theory of biodiversity. Ecol Lett. 2007;10:1017–28.PubMed 
PubMed Central 
Article 

Google Scholar 
26.Bahram M, Hildebrand F, Forslund SK, Anderson JL, Soudzilovskaia NA, Bodegom PM, et al. Structure and function of the global topsoil microbiome. Nature. 2018;560:233–7.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
27.Zhou J, Deng Y, Shen L, Wen C, Yan Q, Ning D, et al. Temperature mediates continental-scale diversity of microbes in forest soils. Nat Commun. 2016;7:1–10.
Google Scholar 
28.Thompson LR, Jex AR, Campbell AH, Linz AM, Berry A, Williams AE, et al. A communal catalogue reveals Earth’s multiscale microbial diversity. Nature. 2017;551:457–63.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
29.Bickel S, Or D. Soil bacterial diversity mediated by microscale aqueous-phase processes across biomes. Nat Commun. 2020;11:1–9.
Google Scholar 
30.Xu X, Thornton PE, Post WM. A global analysis of soil microbial biomass carbon, nitrogen and phosphorus in terrestrial ecosystems. Glob Ecol Biogeogr. 2013;22:737–49.Article 

Google Scholar 
31.Serna-Chavez HM, Fierer N, van Bodegom PM. Global drivers and patterns of microbial abundance in soil: global patterns of soil microbial biomass. Glob Ecol Biogeogr. 2013;22:1162–72.Article 

Google Scholar 
32.Wang G, Or D. A hydration-based biophysical index for the onset of soil microbial coexistence. Sci Rep. 2012;2:1–5.
Google Scholar 
33.Li CH, Lee CK. Minimum cross entropy thresholding. Pattern Recognit. 1993;26:617–625.Article 

Google Scholar 
34.Walt S, van der, Schönberger JL, Nunez-Iglesias J, Boulogne F, Warner JD, Yager N, et al. scikit-image: image processing in Python. PeerJ. 2014;2:1–18.
Google Scholar 
35.Homem-de-Mello T, Rubinstein RY. Estimation of rare event probabilities using cross-entropy. Proc Winter Simul Conf. 2002;1:310–19.Article 

Google Scholar 
36.Murali A, Bhargava A, Wright ES. IDTAXA: a novel approach for accurate taxonomic classification of microbiome sequences. Microbiome. 2018;6:1–14.Article 

Google Scholar 
37.Šťovíček A, Kim M, Or D, Gillor O. Microbial community response to hydration-desiccation cycles in desert soil. Sci Rep. 2017;7:1–9.Article 
CAS 

Google Scholar 
38.Zhao M, Heinsch FA, Nemani RR, Running SW. Improvements of the MODIS terrestrial gross and net primary production global data set. Remote Sens Environ. 2005;95:164–76.Article 

Google Scholar 
39.Fick SE, Hijmans RJ. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas: new climate surfaces for global land areas. Int J Climatol. 2017;37:4302–15.Article 

Google Scholar 
40.Schoolfield RM, Sharpe PJH, Magnuson CE. Non-linear regression of biological temperature-dependent rate models based on absolute reaction-rate theory. J Theor Biol. 1981;88:719–31.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
41.Beck HE, Wood EF, Pan M, Fisher CK, Miralles DG, van Dijk AIJM, et al. MSWEP V2 Global 3-hourly 0.1° precipitation: methodology and quantitative assessment. Bull Am Meteorol Soc. 2019;100:473–500.Article 

Google Scholar 
42.Wang G, Or D. Hydration dynamics promote bacterial coexistence on rough surfaces. ISME J. 2013;7:395–404.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
43.Kim M, Or D. Individual-based model of microbial life on hydrated rough soil surfaces. PLoS ONE. 2016;11:1–31.
Google Scholar 
44.Hermsen R, Okano H, You C, Werner N, Hwa T. A growth-rate composition formula for the growth of E.coli on co-utilized carbon substrates. Mol Syst Biol. 2015;11:1–6.Article 
CAS 

Google Scholar 
45.García FC, Bestion E, Warfield R, Yvon-Durocher G. Changes in temperature alter the relationship between biodiversity and ecosystem functioning. Proc Natl Acad Sci. 2018;115:10989–94.PubMed 
Article 
CAS 
PubMed Central 

Google Scholar 
46.Slessarev EW, Lin Y, Bingham NL, Johnson JE, Dai Y, Schimel JP, et al. Water balance creates a threshold in soil pH at the global scale. Nature. 2016;540:567–9.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
47.Treves DS, Xia B, Zhou J, Tiedje JM. A two-species test of the hypothesis that spatial isolation influences microbial diversity in soil. Micro Ecol. 2003;45:20–8.CAS 
Article 

Google Scholar 
48.Campbell BJ, Yu L, Heidelberg JF, Kirchman DL. Activity of abundant and rare bacteria in a coastal ocean. Proc Natl Acad Sci. 2011;108:12776–81.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
49.Stauffer D. Scaling theory of percolation clusters. Phys Rep. 1979;54:1–74.Article 

Google Scholar 
50.Scher H, Zallen R. Critical density in percolation processes. J Chem Phys. 1970;53:3759–61.CAS 
Article 

Google Scholar 
51.Hengl T, de Jesus JM, Heuvelink GB, Gonzalez MR, Kilibarda M, Blagotić A, et al. SoilGrids250m: global gridded soil information based on machine learning. PloS ONE. 2017;12:1–40.Article 
CAS 

Google Scholar 
52.Chase AB, Arevalo P, Brodie EL, Polz MF, Karaoz U, Martiny JBH. Maintenance of sympatric and allopatric populations in free-living terrestrial bacteria. mBio. 2019;10:1–11.Article 

Google Scholar 
53.Fisher CK, Mehta P. The transition between the niche and neutral regimes in ecology. Proc Natl Acad Sci. 2014;111:13111–6.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
54.Ratzke C, Barrere J, Gore J. Strength of species interactions determines biodiversity and stability in microbial communities. Nat Ecol Evol. 2020;4:376–83.PubMed 
Article 
PubMed Central 

Google Scholar 
55.Doud DFR, Bowers RM, Schulz F, Raad MD, Deng K, Tarver A, et al. Function-driven single-cell genomics uncovers cellulose-degrading bacteria from the rare biosphere. ISME J. 2020;14:659–75.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
56.Shade A, Jones SE, Caporaso JG, Handelsman J, Knight R, Fierer N, et al. Conditionally rare taxa disproportionately contribute to temporal changes in microbial diversity. mBio. 2014;5:1–9.Article 
CAS 

Google Scholar 
57.Kaminsky R, Morales SE. Conditionally rare taxa contribute but do not account for changes in soil prokaryotic community structure. Front Microbiol. 2018;9:1–6.Article 

Google Scholar 
58.Price PB, Sowers T. Temperature dependence of metabolic rates for microbial growth, maintenance, and survival. Proc Natl Acad Sci U S A. 2004;101:4631–6.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar  More

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

IACUC approval: all work described below was approved by the Boise State Institutional Animal Care and Use Committee: AC15-021.Site layoutWe selected 20 sites, across five drainages, within the Pioneer Mountains of Idaho—matched for elevation and riparian habitat. We split these 20 sites into 10 noise playback sites, and 10 control sites (Fig. 1A; S1). The control sites ranged from quiet, slow-moving streams to relatively loud whitewater torrents. Noise playback sites, on the other hand, were relatively quiet (not whitewater) sites, where we broadcast loud whitewater river recordings with speaker arrays hung from towers (Fig. S1; S2; S3; S4; see supplementary information for more details on noise file creation, playback equipment, and experimental setup). At five of the noise playback sites we broadcast normal river noise (hereafter referred to as ‘river noise’ sites), and at the other five noise sites we broadcast spectrally-altered river recordings (hereafter referred to as “shifted noise” sites).Our field sites were oriented along the riparian zone, with data collection occurring at three primary locations within each site (Fig. S1): (1) roughly in the middle of the speaker tower systems, (2) at a shorter distance from the middle location (mean: 198.2 ± 54.5 m SD; range: 117.6–384.5 m), and (3) and a longer distance from the middle location (in the opposite direction from the nearer location; mean: 312.7 ± 64.7 m SD; range: 249.1–479.6 m). Thus, sites were approximately 510.9 ± 98.3 m long (range: 374.7–850.6 m), along the riparian corridor. All control sites were, at minimum, 1 km apart along the riparian corridor from any noise site, to maintain acoustic independence (see Fig. 1A; S1).Data collectionBirds
We conducted three-minute avian point counts between one half hour before sunrise and 6 h after sunrise (roughly 0530–1130 h). During the project, we conducted 1330 point-counts from 28 May to 20 July 2017 and 1639 point-count events occurred from 7 May to 24 July in 2018.
Caterpillar deploymentWe deployed a total of 720 clay caterpillars throughout the 2018 breeding season. Forty caterpillars were glued to stems and branches of trees between 1 and 2.5 m high at each site (Fig. S8). Twenty caterpillars surrounded the middle point count location at each site (a set of 10 were placed upstream, and another set of 10 were placed downstream starting from the middle ARU location), while the other twenty were at upstream and downstream sampling locations (10 each at upstream and downstream locations). We placed each caterpillar along the riparian corridor, at least 1 m apart from each other30. See Supplementary information for details on caterpillar predation scoring.Bird trait analysisWe performed a trait-based analysis to understand the mechanistic patterns of bird distributions in our study paradigm. Avian vocal frequencies and body mass were collected from Hu and Cardoso 2009, Cardoso 2014, and Francis 201516,31,32. When multiple sources contained data, the values were averaged. There were a few cases where none of those sources contained a vocal frequency or mass measurement for species of interest. Thus, representative songs were downloaded from the Macaulay Library of the Cornell Lab of Ornithology based on recording quality and geographical relevance (MacGillivray’s warbler: ML42249; dusky flycatcher: ML534684; red-naped sapsuckers: ML6956), and analyzed with Avisoft SASLab Pro to obtain a peak frequency measure. Mass measurements for these ‘missing’ birds were taken from the ‘All about birds’ webpage of the Cornell Lab of Ornithology.BatsMeasuring and identifying bat callsWe measured bat activity using Song Meter 3 (hereafter “SM3”) recording units (Wildlife Acoustics Inc., Massachusetts, USA) equipped with a single SMU (Wildlife Acoustics Inc.) ultrasonic microphone. One recording unit was used at each site and we pseudo-randomly rotated the unit between the three point-count locations so that each location was monitored for at least 21 days. We mounted microphones on metal conduit at a height of ~3 m, oriented perpendicular to the ground and facing away from the stream to optimize recording conditions (Fig. S9; S10; see Supplementary information for more information).Robotic insectsWe used a modified version of Lazure and Fenton’s26 apparatus to present bats with a fluttering target (Fig. S12). This consisted of a 3 cm2 piece of masking tape affixed to a metal rod [30.48 cm length × 3.25 mm diameter], which itself was connected to a 12-volt brushed DC motor (AndyMark 9015 12 V, AndyMark Inc., Kokomo, IN, USA). The no-load revolution speed of these motors (267 Hz) falls within the range of wingbeat frequency measured in Chironomidae27,33, a group that is an important food source for many North American bat species34.We attached each motor to a tripod made of PVC piping and positioned the tripod such that the target was approximately 1.2 m above the ground. Each motor was powered by a 12 V battery (35Ah AGM; DURA12-35C, Duracell) which was controlled by a programmable 12 V timer (CN101, FAVOLCANO) to automatically start and stop the motor each night. The rotors were powered for 2 h following sunset.Prey-sound speaker playbackWe created a playlist composed of several insect acoustic cues to present gleaning bats: a beetle (Tenebrio molitor) walking on dried grass, a cricket (Acheta domesticus) walking on leaves, mealworm larvae (Tenebrio molitor) on leaves, fall field cricket (Gryllus pennsylvanicus) calls, and fork-tailed bush katydid (Scudderia furcata) calls. The cricket and katydid calls were sourced from the Macaulay Library (ML527360 and ML107505, respectively).Experimental setup for bat foraging testsMost sites received two rotors (Fig. S12) and two speakers (Fig. S13): one of each at the center of the site, and one of each at approximately 125 m from the center of the site (in opposite directions in order to have tests in a range of acoustic environments), placed roughly 10 m from the edge of the riparian zone. Rotors and speakers at the center locations were separated by at least 50 m. The exception to this setup were the four positive control (loud whitewater river) sites, which only received a single rotor and speaker separated by 50 m because of logistical difficulties of accessing those sites. We paired each rotor and speaker with an SM2BAT + bat detector equipped with an SMX-US microphone (Wildlife Acoustics Inc.)35, using tripods to elevate the microphones approximately 1 m off the ground and ~1 m from the speaker/rotor. We programmed the bat detectors with a gain of 36 dB and a trigger level of 18 dB to limit recordings to bats that were passing within the immediate vicinity. To allow for a comparison of activity between speakers and rotors, bat activity was only considered for the first two hours following sunset.Bat trait analysisWe collected bat foraging behavior and peak echolocation frequency information to use as predictors in a phylogenetically controlled trait analysis (Tables S8; S13). We based our behavioral foraging classifications on the categories of Ratcliffe et al.36 and followed the classifications of Gordon et al.37 where possible, and others38,39,40,41,42,43 where necessary. We extracted peak echolocation frequency from the 2017 and 2018 SM3 field recordings and employed two controls to decrease variability in call parameters potentially introduced via this method. First, we selected only recordings made on control sites in 2017 and 2018 (n = 740,848 calls), as echolocation call characteristics may be affected by local acoustic environments (e.g., Bunkley et al.)22. Secondly, we averaged all call parameters per species per hour at each site to decrease the possible effects of few individuals driving measurements. This resulted in 9538 species-hours of recordings, which themselves were averaged per species (Table S13).Quantifying environmental variablesWe used long-term monitoring of the acoustic environment (via Roland R05 recorders) to calculate daily sound pressure level (L50 dBA) and median frequency (kHz) values for each location (see supplementary information for details on quantification of all predictor variables).Sound pressure level (SPL)We converted 106,769 h of long-term ARU recordings into daily-averaged median sound pressure levels (L50; measured as dBA rel. 20 µPa) see refs. 13,44 using custom software ‘AUDIO2NVSPL’ and ‘Acoustic Monitoring Toolbox’ (Damon Joyce, Natural Sounds and Night Skies Division, National Park Service).Acoustic environment spectrumWe used custom software45 in the programming language R and the package ‘FFmpeg’ in command prompt to convert 106,769 h of long-term recordings into 71,282 individual 3-minute files starting each hour of the day (Fig. S5). Thus 24, 3-min files were created per acoustic recording location per day (one for every hour). We then used the packages “tuneR” and “seewave” to read in and measure the median frequency of sound files, respectively45,46,47. These hourly metrics were then averaged by date to create a daily metric.StatisticsAll models of abundance, activity, and foraging transects were generalized linear mixed effects models (glmm) in R48 using the package ‘lme4’49,50 or ‘glmmTMB’51. All distribution families were selected based on theoretical sampling processes of the data, models were checked for collinearity (VIF scores)52, and model fits were visually checked with residual plots (see supplemental R code)53.Bird abundance and bat activity
Model predictors and covariates
Both bird and bat models had the following variables in a glmm: site and bird/bat species were random effects terms and sound pressure level (dBA L50), sound spectrum (median frequency), the interaction between sound pressure level and spectrum, elevation, percent riparian vegetation, ordinal date (and a quadratic version of this), and year as fixed effects. While year is sometimes used as a random-effect term, it is suggested to be used as a fixed effect if fewer than five levels exist for that factor, as variance estimates become imprecise54,55. Additionally, moon phase was a fixed effect in the bat models56, while spectral overlap (the absolute difference between sound spectrum and bird species vocalization frequencies) and the interaction between sound pressure level and spectral overlap were fixed effects in bird models.
We attempted to fit both sound pressure level and spectrum as having random slopes for each species, yet both bat and bird models would not converge with such complex model structure. Thus, we followed group models with individual species models (see Supplementary information).

Model family distribution and link function
For both bird and bat counts, we used a negative binomial distribution with a log link, rather than a Poisson distribution, because data were over-dispersed. We plotted variance-mean relationships and residuals of multiple models to select the appropriate variance structure, and compared these with AIC to select the best-fitting distribution (see R script for further justification of these methods)54.

Individual species models
Individual species models were parameterized the same as above (except without the species term). All 12 bat species (see Tables S6; S10) and 26 of the most common birds (see Tables S2; S9) were modeled individually to be able to interpret model parameter estimates, with complex interactions, for each species.
Clay caterpillar predationWe modeled caterpillar predation with a glmm (binomial family; logit link function), using the number of individual scorers as weights in the model. Like the bird abundance model, we used site as a random effect and sound pressure level (dBA L50), spectral frequency (median), elevation, percent riparian vegetation, ordinal date, and year as fixed effects (Table S4). Additionally, the predicted number of birds at a site were modeled as fixed effects to control for varying amounts of foraging birds on the landscape.Robotic moths and prey-sound speakersRobotic moth and prey-sound speaker models were parameterized exactly the same as the overall bat activity model. That is, the model was fit with a negative binomial family (log link) with site and species as random effects and sound pressure level (dBA L50), sound spectrum (median frequency), the interaction between sound pressure level and spectrum, moon phase, elevation, percent riparian vegetation, ordinal date (and a quadratic version of this), and year as fixed effects. Additionally, the predicted number of bats at a site were modeled as fixed effects to control for varying amounts of foraging bats on the landscape.Trait analysesWe performed trait analyses with phylogenetic generalized least squares (PGLS) to control for relatedness while predicting species responses to noise12. We performed PGLS analyses with the gls function in the R package nlme57, and accounted for error in the response variable with a fixed-variance weighting function of one divided by the square root of the standard error of the response estimate58,59. We accounted for phylogenetic structure by estimating Pagel’s λ60. When λ estimates fell outside of the zero to 1 range, we fixed λ at the nearest boundary. For bird models, we used a pruned consensus tree from a recent class-wide phylogeny61. For bats, we used a pruned mammalian tree62. We used initial global models with all traits as variables that explained the responses to sound pressure level (SPL; birds and bats), spectral overlap with birdsong (birds), background frequency (bats), and the interaction between SPL and each measure of frequency (birds and bats). We then used AIC model selection63 to choose top models in explaining these patterns. Models with dAIC ≤4 are included in Table S3 (birds) and Table S8 (bats), and the top model is interpreted in the main text.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. 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