Ecology

Zimmer, H. D. et al. Memory for Action: A Distinct Form of Episodic Memory? (Oxford University Press, 2001).
Google Scholar 
Goswami, U. The Wiley-Blackwell Handbook of Childhood Cognitive Development (Wiley, 2013).
Google Scholar 
Fujita, K., Morisaki, A., Takaoka, A., Maeda, T. & Hori, Y. Incidental memory in dogs (Canis familiaris): Adaptive behavioral solution at an unexpected memory test. Anim. Cogn. 15, 1055–1063 (2012).Article 
PubMed 

Google Scholar 
Lind, J., Enquist, M. & Ghirlanda, S. Animal memory: A review of delayed matching-to-sample data. Behav. Processes 117, 52–58 (2015).Article 
PubMed 

Google Scholar 
Kuczaj, S. A. II. & Eskelinen, H. C. (2014) The “creative dolphin” revisited: What do dolphins do when asked to vary their behavior. Anim. Behav. Cogn. 1, 66–77 (2014).Article 

Google Scholar 
Tulving, E. Episodic and semantic memory. Organ. Mem. 1, 381–403 (1972).
Google Scholar 
Tulving, E. How many memory systems are there?. Am. Psychol. 40, 385 (1985).Article 

Google Scholar 
Fugazza, C., Pongrácz, P., Pogány, Á., Lenkei, R. & Miklósi, Á. Mental representation and episodic-like memory of own actions in dogs. Sci. Rep. 10, 1–8 (2020).Article 

Google Scholar 
Hauser, M. D., Chomsky, N. & Fitch, W. T. The faculty of language: What is it, who has it, and how did it evolve?. Science 298, 1569–1579 (2002).Article 
PubMed 

Google Scholar 
Conway, M. A. Memory and the self. J. Mem. Lang. 53, 594–628 (2005).Article 

Google Scholar 
Scagel, A. & Mercado, E. III. Do that again! Memory for self-performed actions in dogs (Canis familiaris). J. Comp. Psychol. 20, 25 (2022).
Google Scholar 
Mercado, E., Murray, S. O., Uyeyama, R. K., Pack, A. A. & Herman, L. M. Memory for recent actions in the bottlenosed dolphin (Tursiops truncatus): Repetition of arbitrary behaviors using an abstract rule. Learn. Behav. 26, 210–218 (1998).Article 

Google Scholar 
Paukner, A., Anderson, J. R., Donaldson, D. I. & Ferrari, P. F. Cued repetition of self-directed behaviors in macaques (Macaca nemestrina). J. Exp. Psychol. Anim. Behav. Process. 33, 139 (2007).Article 
PubMed 

Google Scholar 
Smeele, S. Q. et al. Memory for own behaviour in pinnipeds. Anim. Cogn. 20, 1–12 (2019).
Google Scholar 
Clayton, N. S. Episodic-like memory and mental time travel in animals. (2017).Clayton, N. S., Griffiths, D. P. & Dickinson, A. Declarative and episodic-like memory in animals: Personal musings of a Scrub Jay (2000).Clayton, N. S. & Dickinson, A. Episodic-like memory during cache recovery by scrub jays. Nature 395, 272–274 (1998).Article 
PubMed 

Google Scholar 
Tulving, E. Episodic memory and autonoesis: Uniquely human. Missing Link Cogn. Orig. Self-Reflect. Conscious 20, 3–56 (2005).
Google Scholar 
Suddendorf, T. & Corballis, M. C. Mental time travel and the evolution of the human mind. Genet. Soc. Gen. Psychol. Monogr. 123, 133–167 (1997).PubMed 

Google Scholar 
Suddendorf, T. & Corballis, M. C. The evolution of foresight: What is mental time travel, and is it unique to humans?. Behav. Brain Sci. 30, 299–313 (2007).Article 
PubMed 

Google Scholar 
Crystal, J. D. Evaluating evidence from animal models of episodic memory. J. Exp. Psychol. Anim. Learn. Cogn. 47, 337 (2021).Article 
PubMed 

Google Scholar 
Mercado, E. III., Uyeyama, R. K., Pack, A. A. & Herman, L. M. Memory for action events in the bottlenosed dolphin. Anim. Cogn. 2, 17–25 (1999).Article 

Google Scholar 
Zentall, T. R. Coding of stimuli by animals: Retrospection, prospection, episodic memory and future planning. Learn. Motiv. 41, 225–240 (2010).Article 
PubMed 
PubMed Central 

Google Scholar 
Fugazza, C., Pogány, Á. & Miklósi, Á. Recall of others’ actions after incidental encoding reveals episodic-like memory in dogs. Curr. Biol. 26, 3209–3213 (2016).Article 
PubMed 

Google Scholar 
Lambert, M. L., Jacobs, I., Osvath, M. & von Bayern, A. M. Birds of a feather? Parrot and corvid cognition compared. Behaviour 156, 505–594 (2019).Article 

Google Scholar 
Emery, N. J. Cognitive ornithology: The evolution of avian intelligence. Philos. Trans. R. Soc. B Biol. Sci. 361, 23–43 (2006).Article 

Google Scholar 
Olkowicz, S. et al. Birds have primate-like numbers of neurons in the forebrain. Proc. Natl. Acad. Sci. 113, 7255–7260 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Emery, N. J. & Clayton, N. S. Evolution of the avian brain and intelligence. Curr. Biol. 15, R946–R950 (2005).Article 
PubMed 

Google Scholar 
Bradbury, J. W. & Balsby, T. J. The functions of vocal learning in parrots. Behav. Ecol. Sociobiol. 70, 293–312 (2016).Article 

Google Scholar 
Baciadonna, L., Cornero, F. M., Emery, N. J. & Clayton, N. S. Convergent evolution of complex cognition: Insights from the field of avian cognition into the study of self-awareness. Learn. Behav. 49, 9–22 (2021).Article 
PubMed 

Google Scholar 
Osvath, M., Kabadayi, C. & Jacobs, I. Independent evolution of similar complex cognitive skills (2014).Zentall, T. R., Clement, T. S., Bhatt, R. S. & Allen, J. Episodic-like memory in pigeons. Psychon. Bull. Rev. 8, 685–690 (2001).Article 
PubMed 

Google Scholar 
Zentall, T. R., Singer, R. A. & Stagner, J. P. Episodic-like memory: Pigeons can report location pecked when unexpectedly asked. Behav. Processes 79, 93–98 (2008).Article 
PubMed 

Google Scholar 
Healy, S. D. & Hurly, T. A. Spatial learning and memory in birds. Brain. Behav. Evol. 63, 211–220 (2004).Article 
PubMed 

Google Scholar 
Taylor, A. H. Corvid cognition. Wiley Interdiscip. Rev. Cogn. Sci. 5, 361–372 (2014).Article 
PubMed 

Google Scholar 
Boeckle, M. & Bugnyar, T. Long-term memory for affiliates in ravens. Curr. Biol. 22, 801–806 (2012).Article 
PubMed 
PubMed Central 

Google Scholar 
Marzluff, J. M., Walls, J., Cornell, H. N., Withey, J. C. & Craig, D. P. Lasting recognition of threatening people by wild American crows. Anim. Behav. 79, 699–707 (2010).Article 

Google Scholar 
Pepperberg, I. M. & Pepperberg, I. M. The Alex Studies: Cognitive and Communicative Abilities of Grey Parrots (Harvard University Press, 2009).Book 

Google Scholar 
Emery, N. J. & Clayton, N. S. Effects of experience and social context on prospective caching strategies by scrub jays. Nature 414, 443–446 (2001).Article 
PubMed 

Google Scholar 
Herzog, S. K. et al. First systematic sampling approach to estimating the global population size of the Critically Endangered Blue-throated Macaw Ara glaucogularis. Bird Conserv. Int. 31, 293–311 (2021).Article 

Google Scholar 
Auersperg, A. M. & von Bayern, A. M. Who’sa clever bird—now? A brief history of parrot cognition. Behaviour 156, 391–407 (2019).Article 

Google Scholar 
Tassin de Montaigu, C., Durdevic, K., Brucks, D., Krasheninnikova, A. & von Bayern, A. Blue-throated macaws (Ara glaucogularis) succeed in a cooperative task without coordinating their actions. Ethology 126, 267–277 (2020).Article 

Google Scholar 
Auersperg, A. M. et al. Social transmission of tool use and tool manufacture in Goffin cockatoos (Cacatua goffini). Proc. R. Soc. B Biol. Sci. 281, 20140972 (2014).Article 

Google Scholar 
Brucks, D. & von Bayern, A. M. Parrots voluntarily help each other to obtain food rewards. Curr. Biol. 30, 292–297 (2020).Article 
PubMed 

Google Scholar 
Krasheninnikova, A., Höner, F., O’Neill, L., Penna, E. & von Bayern, A. M. Economic decision-making in parrots. Sci. Rep. 8, 1–10 (2018).Article 

Google Scholar 
Jarvis, E. D. et al. Avian brains and a new understanding of vertebrate brain evolution. Nat. Rev. Neurosci. 6, 151–159 (2005).Article 
PubMed 

Google Scholar 
Gutiérrez-Ibáñez, C., Iwaniuk, A. N. & Wylie, D. R. Parrots have evolved a primate-like telencephalic-midbrain-cerebellar circuit. Sci. Rep. 8, 1–11 (2018).Article 

Google Scholar 
Smeele, S. Q. et al. Coevolution of relative brain size and life expectancy in parrots. Proc. R. Soc. B 289, 20212397 (2022).Article 
PubMed 
PubMed Central 

Google Scholar 
Kirsch, J. A., Güntürkün, O. & Rose, J. Insight without cortex: Lessons from the avian brain. Conscious. Cogn. 17, 475–483 (2008).Article 
PubMed 

Google Scholar 
Dunbar, R. I. & Shultz, S. Evolution in the social brain. Science 317, 1344–1347 (2007).Article 
PubMed 

Google Scholar 
Wright, A. A. & Katz, J. S. Mechanisms of same/different concept learning in primates and avians. Behav. Processes 72, 234–254 (2006).Article 
PubMed 

Google Scholar 
Smirnova, A. A., Obozova, T. A., Zorina, Z. A. & Wasserman, E. A. How do crows and parrots come to spontaneously perceive relations-between-relations?. Curr. Opin. Behav. Sci. 37, 109–117 (2021).Article 

Google Scholar 
Schusterman, R. J. & Kastak, D. A California sea lion (Zalophus californianus) is capable of forming equivalence relations. Psychol. Rec. 43, 823–839 (1993).Article 

Google Scholar 
Kastak, D. & Schusterman, R. J. Transfer of visual identity matching-to-sample in two California sea lions (Zalophus californianus). Anim. Learn. Behav. 22, 427–435 (1994).Article 

Google Scholar 
Zentall, T. R., Wasserman, E. A., Lazareva, O. F., Thompson, R. K. & Rattermann, M. J. Concept learning in animals. Comp. Cogn. Behav. Rev. 20, 25 (2008).
Google Scholar 
Marino, L. Convergence of complex cognitive abilities in cetaceans and primates. Brain. Behav. Evol. 59, 21–32 (2002).Article 
PubMed 

Google Scholar 
Huber, L., Range, F. & Virányi, Z. Dog imitation and its possible origins. In Domestic dog Cognition and Behavior 79–100 (Springer, 2014).Chapter 

Google Scholar 
Schmidjell, T., Range, F., Huber, L. & Virányi, Z. Do owners have a Clever Hans effect on dogs? Results of a pointing study. Front. Psychol. 3, 558 (2012).Article 
PubMed 
PubMed Central 

Google Scholar 
Hare, B., Brown, M., Williamson, C. & Tomasello, M. The domestication of social cognition in dogs. Science 298, 1634–1636 (2002).Article 
PubMed 

Google Scholar 
Lindqvist, C. & Jensen, P. Domestication and stress effects on contrafreeloading and spatial learning performance in red jungle fowl (Gallus gallus) and White Leghorn layers. Behav. Processes 81, 80–84 (2009).Article 
PubMed 

Google Scholar 
Pack, A. A., Herman, L. M. & Roitblat, H. L. Generalization of visual matching and delayed matching by a California sea lion (Zalophus californianus). Anim. Learn. Behav. 19, 37–48 (1991).Article 

Google Scholar 
Bennett, M. S. Five breakthroughs: A first approximation of brain evolution from early bilaterians to humans. Front. Neuroanat. 15, 25 (2021).Article 

Google Scholar 
Cisek, P. Resynthesizing behavior through phylogenetic refinement. Atten. Percept. Psychophys. 81, 2265–2287 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Toft, C. A. & Wright, T. F. Parrots of the wild. Nat. Hist. World’s Most Captiv. Birds 20, 25 (2015).
Google Scholar 
Merkle, J. A., Sigaud, M. & Fortin, D. To follow or not? How animals in fusion–fission societies handle conflicting information during group decision-making. Ecol. Lett. 18, 799–806 (2015).Article 
PubMed 

Google Scholar 
Stevens, J. R. & Gilby, I. C. A conceptual framework for nonkin food sharing: Timing and currency of benefits. Anim. Behav. 67, 603–614 (2004).Article 

Google Scholar 
Kamil, A. C. & Roitblat, H. L. The ecology of foraging behavior—Implications for animal learning and memory. Annu. Rev. Psychol. 36, 141–169 (1985).Article 
PubMed 

Google Scholar 
Ortiz, S. T., Castro, A. C., Balsby, T. J. S. & Larsen, O. N. Problem-solving in a cooperative task in peach-fronted conures (Eupsittula aurea). Anim. Cogn. 23, 265–275 (2020).Article 

Google Scholar 
Krasheninnikova, A., Brucks, D., Blanc, S. & von Bayern, A. M. Assessing African grey parrots’ prosocial tendencies in a token choice paradigm. R. Soc. Open Sci. 6, 190696 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Krasheninnikova, A. et al. Parrots do not show inequity aversion. Sci. Rep. 9, 1–12 (2019).Article 

Google Scholar 
Clayton, N. S., Griffiths, D. P., Emery, N. J. & Dickinson, A. Elements of episodic–like memory in animals. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 356, 1483–1491 (2001).Article 
PubMed 
PubMed Central 

Google Scholar 
McElreath, R. Statistical Rethinking: A Bayesian Course with Examples in R and Stan (Chapman and Hall, 2020).Book 

Google Scholar  More

Aristegui, J., Gasol, J. M., Duarte, C. M. & Herndl, G. J. Microbial oceanography of the dark ocean’s pelagic realm. Limnol. Oceanogr. 54, 1501–1529 (2009).Article 

Google Scholar 
Jannasch, H. W., Eimhjellen, K., Wirsen, C. O. & Farmanfarmaian, A. Microbial degradation of organic matter in the deep sea. Science 171, 672–675 (1971).Article 

Google Scholar 
Tamburini, C., Boutrif, M., Garel, M., Colwell, R. R. & Deming, J. W. Prokaryotic responses to hydrostatic pressure in the ocean – a review. Environ. Microbiol. 15, 1262–1274 (2013).Article 

Google Scholar 
Yayanos, A. A. Microbiology to 10,500 meters in the deep-sea. Annu. Rev. Microb. 49, 777–805 (1995).Article 

Google Scholar 
Jebbar, M., Franzetti, B., Girard, E. & Oger, P. Microbial diversity and adaptation to high hydrostatic pressure in deep-sea hydrothermal vents prokaryotes. Extremophiles 19, 721–740 (2015).Article 

Google Scholar 
Yayanos, A. A. Evolutional and ecological implications of the properties of deep-sea barophilic bacteria. Proc. Natl Acad. Sci. USA 83, 9542–9546 (1986).Article 

Google Scholar 
Nagata, T. et al. Emerging concepts on microbial processes in the bathypelagic ocean – ecology, biogeochemistry, and genomics. Deep-Sea Res. II 57, 1519–1536 (2010).Article 

Google Scholar 
Picard, A. & Daniel, I. Pressure as an environmental parameter for microbial life – a review. Biophys. Chem. 183, 30–41 (2013).Article 

Google Scholar 
Herndl, G. J. & Reinthaler, T. Microbial control of the dark end of the biological pump. Nat. Geosci. 6, 718–724 (2013).Article 

Google Scholar 
Marietou, A. & Bartlett, D. H. Effects of high hydrostatic pressure on coastal bacterial community abundance and diversity. Appl. Environ. Microbiol. 80, 5992–6003 (2014).Article 

Google Scholar 
Lauro, F. M. & Bartlett, D. H. Prokaryotic lifestyles in deep sea habitats. Extremophiles 12, 15–25 (2008).Article 

Google Scholar 
Peoples, L. M. et al. Distinctive gene and protein characteristics of extremely piezophilic Colwellia. BMC Genom. 21, 692 (2020).Article 

Google Scholar 
Reinthaler, T. et al. Prokaryotic respiration and production in the meso- and bathypelagic realm of the eastern and western North Atlantic basin. Limnol. Oceanogr. 51, 1262–1273 (2006).Article 

Google Scholar 
Steinberg, D. K. et al. Bacterial vs. zooplankton control of sinking particle flux in the ocean’s twilight zone. Limnol. Oceanogr. 53, 1327–1338 (2008).Article 

Google Scholar 
Burd, A. B. et al. Assessing the apparent imbalance between geochemical and biochemical indicators of meso- and bathypelagic biological activity: what the @$#! is wrong with present calculations of carbon budgets? Deep-Sea Res. II 57, 1557–1571 (2010).Article 

Google Scholar 
Boyd, P. W., Claustre, H., Levy, M., Siegel, D. A. & Weber, T. Multi-faceted particle pumps drive carbon sequestration in the ocean. Nature 568, 327–335 (2019).Article 

Google Scholar 
Kirchman, D., Knees, E. & Hodson, R. Leucine incorporation and its potential as a measure of protein-synthesis by bacteria in natural aquatic systems. Appl. Environ. Microbiol. 49, 599–607 (1985).Article 

Google Scholar 
Nielsen, J. L., Christensen, D., Kloppenborg, M. & Nielsen, P. H. Quantification of cell-specific substrate uptake by probe-defined bacteria under in situ conditions by microautoradiography and fluorescence in situ hybridization. Environ. Microbiol. 5, 202–211 (2003).Article 

Google Scholar 
Sintes, E. & Herndl, G. J. Quantifying substrate uptake by individual cells of marine bacterioplankton by catalyzed reporter deposition fluorescence in situ hybridization combined with micro autoradiography. Appl. Environ. Microbiol. 72, 7022–7028 (2006).Article 

Google Scholar 
Garel, M. et al. Pressure-retaining sampler and high-pressure systems to study deep-sea microbes under in situ conditions. Front. Microbiol 10, 453 (2019).Article 

Google Scholar 
Peoples, L. M. et al. A full-ocean-depth rated modular lander and pressure-retaining sampler capable of collecting hadal-endemic microbes under in situ conditions. Deep-Sea Res. I 143, 50–57 (2019).Article 

Google Scholar 
Gross, M. & Jaenicke, R. Proteins under pressure – the influence of high hydrostatic pressure on structure, function and assembly of proteins and protein complexes. Eur. J. Biochem. 221, 617–630 (1994).Article 

Google Scholar 
Kirchman, D. L. Growth rates of microbes in the oceans. Annu. Rev. Mar. Sci. 8, 285–309 (2016).Article 

Google Scholar 
Ashburner, M. et al. Gene ontology: tool for the unification of biology. Nat. Genet. 25, 25–29 (2000).Article 

Google Scholar 
Xie, Z., Jian, H., Jin, Z. & Xiao, X. Enhancing the adaptability of the deep-sea bacterium Shewanella piezotolerans WP3 to high pressure and low temperature by experimental evolution under H2O2 stress. Appl. Environ. Microbiol. 84, e02342–02317 (2018).Article 

Google Scholar 
Tamburini, C. et al. Effects of hydrostatic pressure on microbial alteration of sinking fecal pellets. Deep-Sea Res. II 56, 1533–1546 (2009).Article 

Google Scholar 
Ivars-Martinez, E. et al. Comparative genomics of two ecotypes of the marine planktonic copiotroph Alteromonas macleodii suggests alternative lifestyles associated with different kinds of particulate organic matter. ISME J. 2, 1194–1212 (2008).Article 

Google Scholar 
Zhao, Z., Baltar, F. & Herndl, G. J. Linking extracellular enzymes to phylogeny indicates a predominantly particle-associated lifestyle of deep-sea prokaryotes. Sci. Adv. 6, eaaz4354 (2020).Article 

Google Scholar 
Bochdansky, A. B., van Aken, H. M. & Herndl, G. J. Role of macroscopic particles in deep-sea oxygen consumption. Proc. Natl Acad. Sci. USA 107, 8287–8291 (2010).Article 

Google Scholar 
Chikuma, S., Kasahara, R., Kato, C. & Tamegai, H. Bacterial adaptation to high pressure: a respiratory system in the deep-sea bacterium Shewanella violacea DSS12. FEMS Microbiol. Lett. 267, 108–112 (2007).Article 

Google Scholar 
Qin, Q. L. et al. Oxidation of trimethylamine to trimethylamine N-oxide facilitates high hydrostatic pressure tolerance in a generalist bacterial lineage. Sci. Adv. 7, eabf9941 (2021).Article 

Google Scholar 
Mestre, M. et al. Sinking particles promote vertical connectivity in the ocean microbiome. Proc. Natl Acad. Sci. USA 115, E6799–E6807 (2018).Article 

Google Scholar 
Thiele, S., Fuchs, B. M., Amann, R. & Iversen, M. H. Colonization in the photic zone and subsequent changes during sinking determine bacterial community composition in marine snow. Appl. Environ. Microbiol. 81, 1463–1471 (2015).Article 

Google Scholar 
Tada, Y. et al. Differing growth responses of major phylogenetic groups of marine bacteria to natural phytoplankton blooms in the western North Pacific Ocean. Appl. Environ. Microbiol. 77, 4055–4065 (2011).Article 

Google Scholar 
Cottrell, M. T. & Kirchman, D. L. Natural assemblages of marine proteobacteria and members of the Cytophaga-Flavobacter cluster consuming low- and high-molecular-weight dissolved organic matter. Appl. Environ. Microbiol. 66, 1692–1697 (2000).Article 

Google Scholar 
Poff, K. E., Leu, A. O., Eppley, J. M., Karl, D. M. & DeLong, E. F. Microbial dynamics of elevated carbon flux in the open ocean’s abyss. Proc. Natl Acad. Sci. USA 118, e2018269118 (2021).Article 

Google Scholar 
Ducklow, H. in Microbial Ecology of the Oceans (ed. Kirchman, D. L.) Ch. 4, 85–120 (Wiley-Liss, 2000).Herndl, G. J. et al. Contribution of archaea to total prokaryotic production in the deep Atlantic Ocean. Appl. Environ. Microbiol. 71, 2303–2309 (2005).Article 

Google Scholar 
Baltar, F., Aristegui, J., Gasol, J. M. & Herndl, G. J. Prokaryotic carbon utilization in the dark ocean: growth efficiency, leucine-to-carbon conversion factors, and their relation. Aquat. Microb. Ecol. 60, 227–232 (2010).Article 

Google Scholar 
Edgcomb, V. P. et al. Comparison of Niskin vs. in situ approaches for analysis of gene expression in deep Mediterranean Sea water samples. Deep-Sea Res. II 129, 213–222 (2016).Article 

Google Scholar 
Cario, A., Oliver, G. C. & Rogers, K. L. Exploring the deep marine biosphere: challenges, innovations, and opportunities. Front. Earth Sci. 7, 225 (2019).Article 

Google Scholar 
Giering, S. L. C. et al. Reconciliation of the carbon budget in the ocean’s twilight zone. Nature 507, 480–483 (2014).Article 

Google Scholar 
Simon, M. & Azam, F. Protein content and protein synthesis rates of planktonic marine bacteria. Mar. Ecol. Prog. Ser. 51, 201–213 (1989).Article 

Google Scholar 
Gasol, J. M. et al. Mesopelagic prokaryotic bulk and single-cell heterotrophic activity and community composition in the NW Africa-Canary Islands coastal-transition zone. Prog. Oceanogr. 83, 189–196 (2009).Article 

Google Scholar 
DeLong, E. F. et al. Community genomics among stratified microbial assemblages in the ocean’s interior. Science 311, 496–503 (2006).Article 

Google Scholar 
Teira, E., Reinthaler, T., Pernthaler, A., Pernthaler, J. & Herndl, G. J. Combining catalyzed reporter deposition-fluorescence in situ hybridization and microautoradiography to detect substrate utilization by bacteria and archaea in the deep ocean. Appl. Environ. Microbiol. 70, 4411–4414 (2004).Article 

Google Scholar 
Woebken, D., Fuchs, B. M., Kuypers, M. M. M. & Amann, R. Potential interactions of particle-associated anammox bacteria with bacterial and archaeal partners in the Namibian upwelling system. Appl. Environ. Microbiol. 73, 4648–4657 (2007).Article 

Google Scholar 
Wand, M. P. Data-based choice of histogram bin width. Am. Stat. 51, 59–64 (1997).
Google Scholar 
Acinas, S. G. et al. Deep ocean metagenomes provide insight into the metabolic architecture of bathypelagic microbial communities. Commun. Biol. 4, 604 (2021).Article 

Google Scholar 
Sunagawa, S. et al. Structure and function of the global ocean microbiome. Science 348, 1261359 (2015).Article 

Google Scholar 
Delmont, T. O. et al. Nitrogen-fixing populations of Planctomycetes and Proteobacteria are abundant in surface ocean metagenomes. Nat. Microbiol. 3, 804–813 (2018).Article 

Google Scholar 
Li, D., Liu, C. M., Luo, R., Sadakane, K. & Lam, T. W. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31, 1674–1676 (2015).Article 

Google Scholar 
Wu, Y. W., Tang, Y. H., Tringe, S. G., Simmons, B. A. & Singer, S. W. MaxBin: an automated binning method to recover individual genomes from metagenomes using an expectation-maximization algorithm. Microbiome 2, 26 (2014).Article 

Google Scholar 
Kang, D. D. et al. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. Peerj 7, e7359 (2019).Article 

Google Scholar 
Olm, M. R., Brown, C. T., Brooks, B. & Banfield, J. F. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. ISME J. 11, 2864–2868 (2017).Article 

Google Scholar 
Chaumeil, P. A., Mussig, A. J., Hugenholtz, P. & Parks, D. H. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics 36, 1925–1927 (2020).
Google Scholar 
Hyatt, D. et al. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinf. 11, 119 (2010).Article 

Google Scholar 
Li, W. & Godzik, A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658–1659 (2006).Article 

Google Scholar 
Eng, J. K., McCormack, A. L. & Yates, J. R. An approach to correlate tandem mass spectral data of peptides with amino acid sequences in a protein database. J. Am. Soc. Mass. Spectrom. 5, 976–989 (1994).Article 

Google Scholar 
Elias, J. E. & Gygi, S. P. Target-decoy search strategy for increased confidence in large-scale protein identifications by mass spectrometry. Nat. Methods 4, 207–214 (2007).Article 

Google Scholar 
Riffle, M. et al. MetaGOmics: a web-based tool for peptide-centric functional and taxonomic analysis of metaproteomics data. Proteomes 6, 2 (2017).Article 

Google Scholar 
Reinthaler, T., van Aken, H. M. & Herndl, G. J. Major contribution of autotrophy to microbial carbon cycling in the deep North Atlantic’s interior. Deep-Sea Res. II 57, 1572–1580 (2010).Article 

Google Scholar 
Yokokawa, T., Yang, Y. H., Motegi, C. & Nagata, T. Large-scale geographical variation in prokaryotic abundance and production in meso- and bathypelagic zones of the central Pacific and Southern Ocean. Limnol. Oceanogr. 58, 61–73 (2013).Article 

Google Scholar 
Frank, A. H., Garcia, J. A., Herndl, G. J. & Reinthaler, T. Connectivity between surface and deep waters determines prokaryotic diversity in the North Atlantic Deep Water. Environ. Microbiol. 18, 2052–2063 (2016).Article 

Google Scholar 
Herndl, G. J., Bayer, B., Baltar, F. & Reinthaler, T. Prokaryotic life in the deep ocean’s water column. Annu. Rev. Mar. Sci. (in the press).Uchimiya, M., Ogawa, H. & Nagata, T. Effects of temperature elevation and glucose addition on prokaryotic production and respiration in the mesopelagic layer of the western North Pacific. J. Oceanogr. 72, 419–426 (2016).Article 

Google Scholar 
Antia, A. N. et al. Basin-wide particulate carbon flux in the Atlantic Ocean: regional export patterns and potential for atmospheric CO2 sequestration. Glob. Biogeochem. Cycles 15, 845–862 (2001).Article 

Google Scholar 
Behrenfeld, M. J. & Falkowski, P. G. Photosynthetic rates derived from satellite-based chlorophyll concentration. Limnol. Oceanogr. 42, 1–20 (1997).Article 

Google Scholar  More

Austin, K. G. et al. Land Use Policy 69, 41–48 (2017).Article 

Google Scholar 
Busch, J. et al. Environ. Res. Lett. 17, 014035 (2022).Article 
CAS 

Google Scholar 
Fleiss, S. et al. Nat. Ecol. Evol. https://doi.org/10.1038/s41559-022-01941-6 (2022).Qaim, M. et al. Annu. Rev. Resour. Econ. 12, 321–344 (2020).Article 

Google Scholar 
Haupt, F. et al. Progress on Corporate Commitments and their Implementation (Tropical Forest Alliance, 2018).Brooks, T. et al. Nat. Ecol. Evol. 1, 0099 (2017).Article 

Google Scholar 
Buisson, E. et al. Biol. Rev. 94, 590–609 (2019).Article 
PubMed 

Google Scholar 
López-Ricaurte, L. et al. Biol. Conserv. 213, 225–233 (2017).Article 

Google Scholar 
Furumo, P. R. & Aide, T. M. Environ. Res. Lett. 12, 024008 (2017).Article 

Google Scholar 
RTRS Standard for Responsible Soy Production Version 3.1 (RTRS, 2017). More

United Nations. [World population prospects 2019]. United Nations. Department of Economic and Social Affairs. World Population Prospects 2019. (2019).Consortium, I. & Commission, E. The circular Bio-society in 2050. (2018).Ramaswami, A., Russell, A. G., Culligan, P. J., Rahul Sharma, K. & Kumar, E. Meta-principles for developing smart, sustainable, and healthy cities. Science (1979) 352, 940–943 (2016).CAS 

Google Scholar 
Cooper, C. M., Troutman, J. P., Awal, R., Habibi, H. & Fares, A. Climate change-induced variations in blue and green water usage in U.S. urban agriculture. J. Clean. Prod. 348, 567–579 (2022).Article 

Google Scholar 
Crippa, M. et al. Food systems are responsible for a third of global anthropogenic GHG emissions. Nat. Food 2, 198–209 (2021).Article 
CAS 

Google Scholar 
Paul, S., Dutta, A., Defersha, F. & Dubey, B. Municipal food waste to biomethane and biofertilizer: A circular economy concept. Waste Biomass Valorizat. 9, 601–611 (2018).Article 
CAS 

Google Scholar 
Zhang, X. et al. Managing nitrogen for sustainable development. Nature 528, 51–59 (2015).Article 
CAS 
PubMed 

Google Scholar 
Bergstrand, K. J. Organic fertilizers in greenhouse production systems—A review. Sci. Hortic. 295, 1–8 (2022).Article 

Google Scholar 
Chiaregato, C. G., França, D., Messa, L. L., dos Santos Pereira, T. & Faez, R. A review of advances over 20 years on polysaccharide-based polymers applied as enhanced efficiency fertilizers. Carbohydr. Polym. 279, 1–10 (2022).Article 

Google Scholar 
Timilsena, Y. P. et al. Enhanced efficiency fertilisers: A review of formulation and nutrient release patterns. J. Sci. Food Agric. 95, 1131–1142 (2015).Article 
CAS 
PubMed 

Google Scholar 
Chen, J. et al. Environmentally friendly fertilizers: A review of materials used and their effects on the environment. Sci. Total Environ. 613–614, 829–839 (2018).Article 
PubMed 

Google Scholar 
Aguilera, E., Lassaletta, L., Sanz-Cobena, A., Garnier, J. & Vallejo, A. The potential of organic fertilizers and water management to reduce N2O emissions in Mediterranean climate cropping systems. A review. Agric. Ecosyst. Environ. 164, 32–52 (2013).Article 
CAS 

Google Scholar 
Lv, G. et al. Biochar-based fertilizer enhanced Cd immobilization and soil quality in soil-rice system. Ecol. Eng. 171, 1–12 (2021).Article 

Google Scholar 
Clark, M. J. & Zheng, Y. Fertilizer rate influences production scheduling of sedum-vegetated green roof mats. Ecol. Eng. 71, 644–650 (2014).Article 

Google Scholar 
Samoraj, M. et al. Biochar in environmental friendly fertilizers—Prospects of development products and technologies. Chemosphere 296, 1–7 (2022).Article 

Google Scholar 
Dimkpa, C. O., Fugice, J., Singh, U. & Lewis, T. D. Development of fertilizers for enhanced nitrogen use efficiency—Trends and perspectives. Sci. Total Environ. 731, 1–9 (2020).Article 

Google Scholar 
Fertahi, S., Ilsouk, M., Zeroual, Y., Oukarroum, A. & Barakat, A. Recent trends in organic coating based on biopolymers and biomass for controlled and slow release fertilizers. J. Control. Release 330, 341–361 (2021).Article 
CAS 
PubMed 

Google Scholar 
García-Garizábal, I., Causapé, J. & Abrahao, R. Nitrate contamination and its relationship with flood irrigation management. J. Hydrol. (AMST) 442–443, 15–22 (2012).Article 

Google Scholar 
Adu-Poku, D., Ackerson, N. O. B., Devine, R. N. O. A. & Addo, A. G. Climate mitigation efficiency of nitrification and urease inhibitors: Impact on N2O emission—A review. Sci. Afr. 16, 1–7 (2022).
Google Scholar 
Ding, W., Qin, H., Yu, S. & Yu, S. L. The overall and phased nitrogen leaching from a field bioretention during rainfall runoff events. Ecol. Eng. 179, 1–9 (2022).Article 

Google Scholar 
Li, X. et al. Loss of nitrogen and phosphorus from farmland runoff and the interception effect of an ecological drainage ditch in the North China Plain—A field study in a modern agricultural park. Ecol. Eng. 169, 1–10 (2021).Article 

Google Scholar 
Michalsky, R. & Pfromm, P. H. Thermodynamics of metal reactants for ammonia synthesis from steam, nitrogen and biomass at atmospheric pressure. AIChE J. 58, 3203–3213 (2012).Article 
CAS 

Google Scholar 
Pleissner, D. Decentralized utilization of wasted organic material in urban areas: A case study in Hong Kong. Ecol. Eng. 86, 120–125 (2016).Article 

Google Scholar 
Masullo, A. Organic wastes management in a circular economy approach: Rebuilding the link between urban and rural areas. Ecol. Eng. 101, 84–90 (2017).Article 

Google Scholar 
Zeng, Y., de Guardia, A., Ziebal, C., de Macedo, F. J. & Dabert, P. Nitrogen dynamic and microbiological evolution during aerobic treatment of digested sludge. Waste Biomass Valorizat. 5, 441–450 (2014).CAS 

Google Scholar 
Nagarajan, S., Eswaran, P., Masilamani, R. P. & Natarajan, H. Chicken feather compost to promote the plant growth activity by using Keratinolytic Bacteria. Waste Biomass Valorizat. 9, 531–538 (2018).Article 
CAS 

Google Scholar 
Bhat, S. A., Singh, J. & Vig, A. P. Earthworms as organic waste managers and biofertilizer producers. Waste Biomass Valorizat. 9, 1073–1086 (2018).Article 
CAS 

Google Scholar 
Mekki, A., Arous, F., Aloui, F. & Sayadi, S. Treatment and valorization of agro-wastes as biofertilizers. Waste Biomass Valorizat. 8, 611–619 (2017).Article 
CAS 

Google Scholar 
Liu, T. et al. Black soldier fly larvae for organic manure recycling and its potential for a circular bioeconomy: A review. Sci. Total Environ. 833, 1–10 (2022).Article 

Google Scholar 
Siddiqui, S. A. et al. Black soldier fly larvae (BSFL) and their affinity for organic waste processing. Waste Manag. 140, 1–13 (2022).Article 
PubMed 

Google Scholar 
Bortolini, S. et al. Hermetia illucens (L.) larvae as chicken manure management tool for circular economy. J. Clean. Prod. 262, 1–10 (2020).Article 

Google Scholar 
Diener, S., Studt Solano, N. M., Roa Gutiérrez, F., Zurbrügg, C. & Tockner, K. Biological treatment of municipal organic waste using black soldier fly larvae. Waste Biomass Valorizat. 2, 357–363 (2011).Article 
CAS 

Google Scholar 
Cai, M. et al. Rapidly mitigating antibiotic resistant risks in chicken manure by Hermetia illucens bioconversion with intestinal microflora. Environ. Microbiol. 20, 4051–4062 (2018).Article 
CAS 
PubMed 

Google Scholar 
Yang, C. et al. Characteristics and mechanisms of ciprofloxacin degradation by black soldier fly larvae combined with associated intestinal microorganisms. Sci. Total Environ. 811, 1–8 (2022).Article 

Google Scholar 
Pang, W. et al. The influence on carbon, nitrogen recycling, and greenhouse gas emissions under different C/N ratios by black soldier fly. Environ. Sci. Pollut. Res. 27, 42767–42777 (2020).Article 
CAS 

Google Scholar 
Beskin, K. v. et al. Larval digestion of different manure types by the black soldier fly (Diptera: Stratiomyidae) impacts associated volatile emissions. Waste Manag. 74, 213–220 (2018).Gligorescu, A. et al. Pilot scale production of Hermetia illucens (L.) larvae and frass using former foodstuffs. Clean Eng. Technol. 10, 1–10 (2022).Rosa, R. et al. Life cycle assessment of chemical vs enzymatic-assisted extraction of proteins from black soldier fly prepupae for the preparation of biomaterials for potential agricultural use. ACS Sustain. Chem. Eng. 8, 14752–14764 (2020).Article 
CAS 

Google Scholar 
Surendra, K. C. et al. Rethinking organic wastes bioconversion: Evaluating the potential of the black soldier fly (Hermetia illucens (L.)) (Diptera: Stratiomyidae) (BSF). Waste Manag. 117, 58–80 (2020).Hasnol, S. et al. A review on insights for green production of unconventional protein and energy sources derived from the larval biomass of black soldier fly. Processes 8, 1–13 (2020).Article 

Google Scholar 
Wong, C. Y. et al. Rhizopus oligosporus-assisted valorization of coconut endosperm waste by black soldier fly larvae for simultaneous protein and lipid to biodiesel production. Processes 9, 1–14 (2021).Article 

Google Scholar 
Raksasat, R. et al. Blended sewage sludge–palm kernel expeller to enhance the palatability of black soldier fly larvae for biodiesel production. Processes 9, 1–13 (2021).Article 

Google Scholar 
Dortmans B.M.A., Diener S. & Verstappen B.M. Black Soldier Fly Biowaste Processing A Step-by-Step Guide. (2017).European Parliament. Regulation (EC) No 767/2009 of the European Parliament and of the council. (2009).Italian Government. Norme in materia ambientale. (Dlgs, 2006).European Parliament. Regulation (EC) No 178/2002 of the European Parliament and of the Council. Official Journal of the European Communities (2002).Palma, L., Fernandez-Bayo, J., Niemeier, D., Pitesky, M. & VanderGheynst, J. S. Managing high fiber food waste for the cultivation of black soldier fly larvae. NPJ Sci. Food 3, 1–7 (2019).Article 

Google Scholar 
Righi, C. et al. Suitability of porous inorganic materials from industrial residues and bioproducts for use in horticulture: A multidisciplinary approach. Appl. Sci. 12, 5437 (2022).Article 
CAS 

Google Scholar 
Barbi, S. et al. Preliminary study on sustainable NPK slow-release fertilizers based on byproducts and leftovers: A design-of-experiment approach. ACS Omega 5, 27154–27163 (2020).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Macavei, L. I., Benassi, G., Stoian, V. & Maistrello, L. Optimization of Hermetia illucens (L.) egg laying under different nutrition and light conditions. PLoS ONE 15, 1–12 (2020).Article 

Google Scholar 
Leni, G., Maistrello, L., Pinotti, G., Sforza, S. & Caligiani, A. Production of carotenoid-rich Hermetia illucens larvae using specific agri-food by-products. J. Insects Food Feed 1, 1–12 (2022).
Google Scholar 
Caligiani, A. et al. Composition of black soldier fly prepupae and systematic approaches for extraction and fractionation of proteins, lipids and chitin. Food Res. Int. 105, 812–820 (2018).Article 
CAS 
PubMed 

Google Scholar 
Montgomery, D. C. Design and Analysis of Experiments Eighth Edition. Design vol. 2 (2012).Barbi, S., Messori, M., Manfredini, T., Pini, M. & Montorsi, M. Rational design and characterization of bioplastics from Hermetia illucens prepupae proteins. Biopolymers 110–118, (2019).Eriksson, L., Johansson, E., Kettaneh-Wold, N., WikstrÄom, C. & Wold, S. Design of Experiments: Principles and Applications. (2008).Morris, P. & John, P. W. M. Statistical Design and Analysis of Experiments. Math. Gaz. 83, 189–200 (1999).Article 

Google Scholar 
Kros, J. F. & Mastrangelo, C. M. Comparing multi-response design methods with mixed responses. Qual Reliab Eng Int 20, 527–539 (2004).Article 

Google Scholar 
Fernandez Pulido, C. R., Caballero, J., Bruns, M. A. & Brennan, R. A. Recovery of waste nutrients by duckweed for reuse in sustainable agriculture: Second-year results of a field pilot study with sorghum. Ecol Eng 168, 1–8 (2021).Kaya, M. et al. Biological, mechanical, optical and physicochemical properties of natural chitin films obtained from the dorsal pronotum and the wing of cockroach. Carbohydr. Polym. 163, 162–169 (2017).Article 
CAS 
PubMed 

Google Scholar 
Kaya, M. et al. On chemistry of γ-chitin. Carbohydr. Polym. 176, 177–186 (2017).Article 
CAS 
PubMed 

Google Scholar 
Poerio, A. et al. Extraction and physicochemical characterization of chitin from cicada orni sloughs of the south-eastern French mediterranean basin. Molecules 25, 1–12 (2020).Article 

Google Scholar 
Sagheer, F. A. A., Al-Sughayer, M. A., Muslim, S. & Elsabee, M. Z. Extraction and characterization of chitin and chitosan from marine sources in Arabian Gulf. Carbohydr. Polym. 77, 410–419 (2009).Article 

Google Scholar 
Waśko, A. et al. The first report of the physicochemical structure of chitin isolated from Hermetia illucens. Int. J. Biol. Macromol. 92, 316–320 (2016).Article 
PubMed 

Google Scholar 
Wang, K. et al. Preparation of bacterial cellulose/silk fibroin double-network hydrogel with high mechanical strength and biocompatibility for artificial cartilage. Cellulose 27, 1845–1852 (2020).Article 
CAS 

Google Scholar 
Morin, A. & Dufresne, A. Nanocomposites of Chitin Whiskers from Riftia Tubes and Poly(caprolactone). Macromolecules 35, 2190–2199 (2002).Article 
CAS 

Google Scholar 
George Socrates. Infrared and Raman Characteristic Group Frequencies: Tables and Charts. (John Wiley & Sons, 2004).Chen, P. & Zhang, L. New evidences of glass transitions and microstructures of soy protein plasticized with glycerol. Macromol. Biosci. 5, 237–245 (2005).Article 
CAS 
PubMed 

Google Scholar 
Robertson, N.-L.M., Nychka, J. A., Alemaskin, K. & Wolodko, J. D. Mechanical performance and moisture absorption of various natural fiber reinforced thermoplastic composites. J. Appl. Polym. Sci. 130, 969–980 (2013).Article 
CAS 

Google Scholar 
Chavez, M. The sustainability of industrial insect mass rearing for food and feed production: Zero waste goals through by-product utilization. Curr. Opin. Insect. Sci. 48, 44–49 (2021).Article 
PubMed 

Google Scholar 
Fisher, H. J. et al. Black soldier fly larvae meal as a protein source in low fish meal diets for Atlantic salmon (Salmo salar). Aquaculture 521, 1–12 (2020).Article 

Google Scholar 
Figueiredo, L. R. F., Nepomuceno, N. C., Melo, J. D. D. & Medeiros, E. S. Glycerol-based polymer adhesives reinforced with cellulose nanocrystals. Int. J. Adhes. Adhes. 110, (2021). More

Simulation experimental set upLaboratory monitoring of leakage migration process can provide an important basis for field tests. The designed and improved ERT device can better describe the migration range of leakage in soil41. In this experiment, a parallel potential monitoring device was used to improve the monitoring of leakage fluid migration. The simulation experiment in the laboratory is carried out in a (100 cm*100 cm*50 cm) plexiglass tank. Sand and clay shall be screened with a 2.36 mm square sieve, watered and compacted with a board to ensure that the soil layer is in close contact with the measuring electrode.Electrode arrangementThe ground wire of high-density electrical method instrument is connected to the electrodes arranged around the bottom of the tank as the power electrode C2, as shown in Fig. 2a. The host is connected to the electrode system. The electrode system consists of 47 electrode grids with a spacing of 0.08 m. The measuring electrode P1 is connected to the mainframe through a wire 0.05 m below the grid center. The geomembrane is located 0.03 m above the measuring electrode P1. The collection device is used as a monitoring system for various leachate. The arrangement of electrodes is shown in Fig. 2b. The power supply electrode C1 is placed at a certain depth in the middle of the saturated sand to provide a constant current. The location of electrode C1 and leakage point is shown in Fig. 2c. The layers from the bottom of the tank are silty clay, geomembrane, silty clay and saturated sand, as shown in Fig. 2d.Figure 2Set-up of leachate migration simulation experiment: (a) Schematic diagram of electrode C2 layout; (b) Schematic diagram of electrical system laying; (c) Position of electrode C1 and leakage point; (d) Schematic diagram of simulated experimental soil layer.Full size imageComposition of monitoring systemThe electrode system is used to monitor the background electric field and artificial electric field of the landfill site. In the experiment, the electrode system is laid in the clay layer under the geomembrane. It is composed of detection electrodes distributed in a grid at a certain distance.The electrical signal conversion system adjusts the measurement mode, sampling accuracy, acquisition frequency and other parameters of the electrode in the field according to the instructions of the mainframe, and transmits the collected electrical signal to the mainframe.The mainframe can control the operation of the monitoring system. The possible leachate points and their pollution range are determined by collecting data. The system mainly includes mainframe and its software system, power supply, etc., as shown in Fig. 3.Figure 3Se2432 parallel electric method instrument.Full size imageLeachate devicePlace 4 leakage bottles above the tank. No.1 and No.4 bottled water are used to simulate the leakage liquid formed by the direct infiltration of rainwater in slag through geomembrane and as a reference. Because Cl-1 is a typical pollutant in the landfill. No. 2 bottle containing 20 g/L NaCl solution is used to simulate inorganic salt leakage in urban life. No. 3 bottle containing 20 ml/L ethanol solution is used to simulate the leakage liquid containing a large amount of organic matter in municipal solid waste. The characteristics of leachate have been summarized in Table1.Table 1 The characteristics of leachate.Full size tableBefore the experiment, configure four solutions, close the injection, use an electric meter to check the conductivity of each measuring point. After each measuring point has no open circuit, supply power to the soil layer through the mainframe to measure the background electric field of the soil. Then open the injection, adjust the flow rate, release the solution at a fixed flow rate, record the soil electric field in the process of leakage every half an hour, collect the potential values of each measuring point, process the data through the potentiometry and potential difference method, and form the relevant potential horizontal profile and longitudinal section of the soil.Principle of potentiometric detection technologyWhen there are leakage points in the landfill, power is supplied to the landfill, and the current forms a current loop through the geomembrane. If there are n (n = 1,2,3…) leakage points in the geomembrane, the power supply current is I, and the artificial electric field will form a leakage electric field at the leakage point, which can be used as a point power supply.$$I = int dI = int j cdot dS$$
(1)
where I is the current intensity, j is the current density vector, and S is the area passing through the current.When there are n leakage points, I will be shunted. If a leakage point is regarded as a finite surface, the current intensity I as:$$I = {I_1} + {I_2} + cdot cdot cdot + {I_{text{n}}} = sumlimits_{i = 1}^n {int_{S_i} {jdS} }$$
(2)
Generally, the power supply current field of landfill site will be affected by the formation medium structure. It is assumed that the formation medium structure is composed of three layers, each layer has uniform properties and stable conductivity, and the layers from top to bottom are: landfill layer, with resistivity of ρ1. The saturated leakage liquid layer above the geomembrane has a resistivity of ρ2. The clay layer under the geomembrane has a resistivity of ρ3. The electrode C1 is arranged in the garbage layer for power supply, and the electrode C2 is arranged at the lower part of the geomembrane away from the electrode system area. The electrode C2 can be regarded as a far pole.Because of the ρ1  > ρ2, the conductivity of the saturated leakage liquid layer at the upper part of the geomembrane is better than that of the landfill layer, so that there is almost no reflected current between the ρ1 layer and the ρ2 layer, that is, the current generated by the power supply electrode C1 is all transmitted to the ρ2 layer. Because of the ρ3  > ρ2, it can be considered that the interface between ρ2 layer and ρ3 layer has both a reflection current, and a transmission current through the leakage point. The potential generated at the detection electrode P1 under the geomembrane is formed by the action of transmission current. The total potential of point P1 is obtained by the superposition of the potential of point power supply passing through n leakage points at P1.$${U_{P1}} = sumlimits_{i = 1}^n {frac{{{I_i}{rho_3}}}{{2pi {{text{r}}_{iP1}}}}}$$
(3)
Parallel potential difference methodThe test adopts pole–pole arrangement, and the calculation formula of apparent resistivity is as follows:$$rho = 2pi {text{aR}}$$
(4)
where ρ is apparent resistivity; a is the distance between electrodes C1 and P1; R is measuring resistivity.When there are loopholes in the geomembrane of the landfill, the leakage liquid will gradually penetrate into the soil layer under the geomembrane through the loopholes, resulting in the change of the conductivity of the soil layer under the geomembrane. The pole-pole acquisition mode of Se2432 parallel electrical instrument is used to obtain the original data (potential difference) of each measuring point on the grid. After current normalization, the apparent resistivity of the soil layer is obtained. The electrical properties of different depths of the soil layer can be obtained by inversion of the apparent resistivity data of the soil layer, so as to determine the occurrence point and distribution range of leakage.The monitoring grid is 5 × 5. The spacing between measuring points is 0.08 m. The measurement method adopted by Se2432 parallel electric method instrument is cross diagonal measurement method. Figure 4 shows that it only needs to measure the potential values on the measuring points on the horizontal, vertical and 45° diagonal lines.Figure 4Schematic diagram of cross-diagonal measurement method.Full size imageTheoretical calculation of test modelTheoretical results of 10 × 10 grid monitoringAccording to the experimental model and statistical data, the resistivity of the clay layer under the geomembrane is assumed ρ = 10 Ω· m, the resistivity ratio of tap water, NaCl solution and ethanol solution after penetrating into the soil layer ρNo.1:ρNo.2:ρNo.3 = 5:3:10. If the four leakage points set by the model are regarded as four conductive resistors, the ratio of the current passing through the four leakage points is INo. 1:INo. 2:INo. 3:INo. 4 = 6:10:3:6.The calculation model is 10 × 10 grid, and the spacing of measuring points is 0.05 m. The potential value on each measuring point is calculated according to Eq. 3, and the obtained data is processed with surfer software to obtain the potential contour map, as shown in Fig. 5. Among them, points 1, 2, 3 and 4 are the leakage positions of water, NaCl solution, ethanol solution and water respectively, and the spacing between leakage points is 0.15 m.Figure 510 × 10 Grid theory detection potential contour map.Full size imageFigure 5 shows that the leakage fields formed by the four kinds of leaking liquids interfere with each other from the theoretical calculation results. The leachate current at point 2 is larger, the high potential closed loop is obvious, and its center corresponds to the leakage center. The reason for this is that the NaCl solution contains conductive particles that increase the conductivity of the leak point. Point 1 and 4 are the same as water, and the leakage electric field is almost the same. Its closed loop is obvious, and the high potential center also corresponds to their leakage position. There is almost no closed loop effect at point 3 under the influence of 1, 2 and 4. The results show that the leakage field formed by high resistance leakage liquid is not easy to be detected by potentiometric detection, and low resistance leakage is suitable to be detected by potentiometric detection.Theoretical results of 12 × 12 grid monitoringThe resistivity of the clay layer under the geomembrane is assumed ρ = 10Ω·m. In consideration of the mutual influence between the leachate and appropriately reduce its influence effect, the resistivity ratio of water, NaCl solution, and ethanol solution after penetrating into the soil layer is set as ρNo.1:ρNo.2:ρNo.3 = 20:15:24, the ratio of the current passing through the four leakage points is INo.1:INo.2:INo.3:INo.4 = 6:8:5:6. And adjust the distance between the two points to 0.28 m. 12 × 12 grid was used for detection, and the spacing of detection points is 0.04 m. Calculate the potential value of each detection point according to Eq. 3, and use Surfer to obtain the detection contour map of four kinds of leakage, as shown in Fig. 6.Figure 612 × 12 Grid theory detection potential contour map.Full size imageTheoretical calculation results show that when the distance between the leakage points is large and the distance between the detection points is small, the potentiometric method can detect the leakage position of various leachates well. At the same time, the diffusion range of different leachates in the same plane is roughly the same, and they all gradually diffuse outward from the center of the leakage point, and the potential value gradually decrease. Point 2 has the largest potential closed loop range, which also has a certain impact on the leakage points of adjacent points 1 and 3. Point 1 and point 4 are water leakage. Affected by different leakage liquids, the leakage electric field of the two same leakage liquids is obviously different. The potential closed loop range of point 1 is larger than that of point 4. Point 3 is the leakage of ethanol solution. Because its resistance is the largest, the range of potential closed loop is the smallest.Figure 7 shows that the leakage fields around the leachates are funnel-shaped, and its size is related to the type of leachate. Therefore, different network density should be designed for different types of leakage liquid, so as to use the most economical scheme to detect the leakage point.Figure 712 × 12 Grid theory detects potential 3d view.Full size imageInterpretation and discussion of resultsLaboratory simulation experiment researchFigure 8a shows the background electric field potential of soil layer. The four injection pipes are opened at the same time and adjusted to the same flow rate. Under the condition of continuous leakage, monitor the leakage field potential at an interval of 1 h. Figure 8b shows the leakage electric field potential value for 1 h. Reduce the injection pipes flow rate to 1/2 of the initial value. Figure 8c shows the monitoring results of 2 h soil layer leakage field potential. Figure 8d shows the soil leakage field potential monitored after 30 min of sealing the injection pipes.Figure 8Leakage field potential diagram of soil layer: (a) Background electric field of soil layer; (b) Potential distribution of soil layer after 1 h of leakage; (c) Potential distribution of soil layer after 2 h of leakage; (d) Potential distribution of soil layer after closing the injection tube for 30 min.Full size imageFigure 8a shows that the background potential contour of the experimental soil layer is at a lower value. Few current lines pass through the monitoring area. A dense closed potential circle of high potential value is formed at point 2. The current flow at point 2 is greater than the other points 1, 3 and 4. The analysis result may be that in the process of watering and compaction, the clay layer under the geomembrane is not uniform, and the compaction degree of the soil layer is different, resulting in different potential values ​​obtained by monitoring. The permeability at point 2 is better than other points, so when the flow rate of the leakage liquid is large, the leakage liquid under the geomembrane gathers near point 2 and spreads out around. After the clay is watered and compacted, the soil compaction is small and the pore water content is large, resulting in a high potential abnormal area in the lower left corner of point 3.Point 2 forms a closed loop of anomaly potential contour much higher than the background electric field, while the value of potential contour coil at leakage point 3 is lower than the surrounding value. It can be analyzed that positions 2 and 3 are leakage points. The leachate at point 2 is a high concentration NaCl solution containing more conductive particles, which will reduce the resistivity of the soil layer under the geomembrane at point 2. Thus, the passing current is increased to form a high potential closed loop. The leachate at point 3 is ethanol solution, which will increase the resistivity of the soil layer under the geomembrane at point 3. So as to reduce the passing current and form a low potential closed loop. Figure 8b shows that the potential contour is consistent with the influence of NaCl solution and ethanol solution on the soil layer under the geomembrane. It can be concluded that point 2 and point 3 are leakage points. The electric field formed after water leakage at point 1 and point 4 cannot clearly distinguish the leakage points.During the monitoring process, the leachate was continuously released from the injection pipe, and the results reflected the dynamic characteristics. Figure 8b shows the phenomenon that the leachate from point 1 and point 4 aggregates around point 2, which is consistent with the inference of better permeability at point 2. Figure 8b,c show that when the flow rate of the leachate is changed and the flow rate of the injection pipe is reduced, the high-potential region of the entire electric field is reduced. Under the influence of gravity, the leachate will migrate longitudinally, and the closed-loop abnormally high-potential regions and abnormally low-potential regions at points 2 and 3 also decrease.Compared with the surrounding potential contours, the difference is more obvious. Figure 8d shows that when the injection pipe stops leaking for a period of time, the leachate migrates longitudinally along the leakage point. At this time, the electric field of the soil layer is similar to the original background electric field, but the potential value is higher than the background electric field, indicating that the leachate is stagnant in the pores of the soil layer, it is the result of changing the electrical properties of the soil layer. The parallel potential method can collect the potential value of each point in the field at one time, which provides a basis for real-time monitoring of landfill leachate.Figure 9 shows the inversion results of the horizontal section of the experimental model. The blue area corresponds to the distribution range of the low resistance anomaly. There are no jump or distortion points in the profile. The resistivity in the longitudinal direction basically shows a change from low to high. The upper layer seepage liquid migrates, and the bottom soil layer is characterized by low humidity and high resistivity. The low-resistance areas formed by the leakage of NaCl solution are widely distributed in the horizontal section. The distribution range is 0–0.28 m, and the migration scale of the leakage liquid can be clearly seen. The morphological characteristics of water leakage in different parts are basically the same. The distribution range is 0–0.18 m. The leakage of ethanol solution is only reflected at 0–0.06 m, and the distribution range is the smallest at the same depth. The ethanol solution also had the slowest migration rate.Figure 9Inversion map of plane section at different depths.Full size imageFigure 10 shows the inversion results of the X–Z longitudinal section of the test model. The two apparent resistivity profiles at Y = 0.24 m and Y = 0.32 m show that there is no low-resistance area in the shallow layer on the soil layer, indicating that the geomembrane in this area is not damaged. The low resistance zone in the middle is caused by the lateral migration of leakage fluid. The low-resistance anomaly area at the top of the profile can be judged as a leak point or formed by the migration of nearby leachate. Combined with the horizontal section, the leakage depth is similar, and the lateral migration speed of leachate is faster than the longitudinal migration speed. Four leak points can be distinguished, delineating the general location of the leak.Figure 10X–Z longitudinal section on different Y axes.Full size imagePhysical model experimentThe potential value of each electrode was monitored after 2 h of leakage, and the resistivity profiles at different positions were obtained by the potential difference method.It can be seen from Fig. 11 that the potential difference method can monitor the leakage of leachate in different directions. The morphological features of the plume formed by the downward migration of the leak point are approximately funnel-shaped in longitudinal section. The affected area of ​​the soil layer can be obtained in time. Figure 11b shows that the potential difference at the monitoring point is very different on both sides. After 2 h of leakage, a large amount of leakage liquid exists in the soil layer. When the water content in the soil layer increases, the diffusion rate of the ethanol solution increases, showing high resistance characteristics. At the same time, due to the action of gravity, there is a lot of vertical migration, and the potential value changes greatly. The profile clearly shows that the distribution area of ​​high potential difference is large, and the distribution of low potential is small. Figure 11c shows that since the migration rate of leachate in the horizontal direction is greater than that in the vertical direction, the potential difference of the monitoring point in the middle region is smaller, and a closed region of a high-potential circle is formed in the middle. The difference between the two results in a smaller potential difference area. Figure 11d shows that almost all the low-potential areas on the monitoring point are on the left side, because the leakage rate of NaCl solution in the horizontal direction is similar to that in the vertical direction under the condition of good soil compaction. At this time, a large number of conductive particles are contained, resulting in a large high-potential region. The difference between the two forms a large area of ​​low potential difference on the left. This is in good agreement with the lower resistance characteristics of the NaCl solution. Figure 11e shows that the two low-resistance regions correspond to the two leakage centers. The low potential difference region is formed by migration around the leak point. The migration speed in the horizontal direction is similar to that in the vertical direction, and the water migration speed on the left is slower than that of the sodium chloride solution on the right. Figure 11e,f show that the monitoring results are the same, but the resulting potential difference is also increased. This is affected by the distance between the monitoring point and the leak point. When the monitored point and the leakage point are located on the same section, the soil layer is the most severely affected area by leakage. Through the change of the potential difference, the leakage range and the location of the leakage point can be better judged.Figure 11Electrical resistivity tomograms of profile: (a) Resistivity of the slitting profile Y = 0; (b) Resistivity of the slitting profile Y = 0.08; (c) Resistivity of the slitting profile Y = 0.16; (d) Resistivity of the slitting profile Y = 0.24; (e) Resistivity of the slitting profile Y = 0.32; (f) Resistivity of the slitting profile Y = 0.4.Full size image More

Phakamani M’Afrika Xaba speaks at a botanical workshop.Credit: Nong Nooch/Tropical Botanical Garden

For Black communities in today’s South Africa, the legacies of colonialism and apartheid still prevail, shaping social structure and limiting access to opportunities. Colonialism displaced Black South Africans from the mid-seventeenth century, eroding cultural and social systems.From the 1950s, apartheid legalized systematic racial discrimination against disenfranchised, mainly Black, people. It limited their economic opportunities and social standing, prescribing an inferior education system to deliberately shape a poor quality of life. The policy fuelled systemic sexism, sexual-orientation discrimination, ageism, and the use of ethnicity as a divide-and-conquer strategy.Seventy years later, even after more than 25 years of democracy following the end of apartheid in 1994, schools and suburbs are still predominantly segregated, with government funding unevenly allocated in terms of facilities and quality of education.Former South African president Nelson Mandela once said, “In Africa there is a concept known as ubuntu — the profound sense that we are human only through the humanity of others; that if we are to accomplish anything in this world, it will in equal measure be due to the work and achievement of others.”As three past and present employees of the South African National Biodiversity Institute (SANBI), a conservation organization founded in 2004 to manage the country’s biodiversity resources, we have been advocating for a culture of treating others in the way we want to be treated: by applying universal shared human values, redefining institutional culture and systems to be inclusive, and opening safe spaces for a diversity of ideas. We have proposed a ground-up approach that aims to focus on holistic transformation at different levels in our organization.Our approach was to initiate a platform to identify inclusivity challenges, foster awareness and collaboration among staff and collectively develop innovative ideas and solutions. These would be aligned to existing organizational values, such as ubuntu, growth, respect and tolerance, excellence, accountability and togetherness. We strive to bring about institutional cultural change through facilitated, constructive conversations, by strengthening connections and cohesion among staff and through creative and proactive problem-solving across our institution.Mentorship that thrivesInstitutional culture needs to enable successful mentoring by creating a safe space. For example, SANBI’s mentoring programme for interns, students and early-career scientists involves quarterly meetings between them and dedicated human-resources staff — check-ins that provide a space to engage with programme coordinators without early-career researchers’ supervisors being present. In addition to sharing feedback on institutional policies and procedures, early-career scientists have the opportunity to discuss challenges they might face because of their supervisor or work placement. When issues are identified early, transfers or exchanges between work programmes can be arranged.Every year, we each sign up to mentor junior researchers to provide a supportive environment for guidance, counselling and the transfer of skills. We develop structured workplans with specific goals and outputs, and we discuss expectations with our protégés. In addition, we offer shared workspaces for interns and encourage peer learning, so that interns can form a peer support network. In these relationships, trust is crucial and can be a gateway to broader professional and personal networks.

Early-career researchers doing fieldwork training at the Stellenbosch University Experimental Farms in South Africa.Credit: Tlou Masehela

Institutions should recruit outside of their walls, if necessary, to ensure that appropriately skilled mentors are paired with early-career researchers. For mentorship to thrive, institutions must also create an enabling environment. In non-supportive environments, staff — particularly those from under-represented groups — who remain inadequately skilled and work without guidance become frustrated. Some can even feel they don’t belong because they see themselves as lagging behind their peers.Institutions often focus too strongly on outputs — such as publications, products or technologies — at the expense of reflecting on the values that uphold the institution. These values might be outdated and out of touch with those of staff, or with those held by partners, stakeholders or society at large. If staff cannot relate to the institutional culture and systems, job satisfaction and retention rates can suffer.Until a few years ago, for example, venues at our organization were named after former staff, as a way of acknowledging their contributions. Inevitably, these were mostly white, male, senior staff, such as Harold Pearson, the first director of Kirstenbosch National Botanical Garden, and Brian Rycroft, who served as director in the 1950s. But the contributions of staff who were employed in junior positions for 20–30 years also needed to be acknowledged. After an outcry around 2014, then-chief-executive Tanya Abrahamse, the first Black woman to hold the post, decided to acknowledge contributions of staff no matter their position. As a result, we now have Richard Crowie Hall, an exhibition space named after one of SANBI’s longest-serving staff members.The protracted legacy of apartheid in South Africa means that if institutional implicit biases are left unaddressed, they can create a fertile ground for racial, ethnic, tribal, financial and gender tensions. We urge more institutional recognition of the contributions of all.Fostering safe spacesThrough our engagements with each other, we have discovered a set of shared values, aligned with those of our institution, and have set out to establish a space to build our vision of a supportive, safe environment based on these values. Safe spaces are required for expressing controversial or uncomfortable views and to do the hard work of finding solutions to inequities. Confidentiality and trust cultivate such safe spaces, which can be created initially in small groups, then expanded to constructive formal or informal spaces. The conversations and suggestions of informal discussion groups about staff development and transformation can be elevated to management for implementation.
Decolonizing science toolkit
Safe spaces are a necessity for institutions that wish to truly address their legacies of racism and colonialism. Policies alone will not create these spaces — they require dedicated staff, too. Such spaces should ensure that those who speak up can do so without fear of being labelled as troublemakers or victimized.As a first step in pursuing this vision, we met with the senior teams at our organization to share ideas around the need for and benefits of an inclusive culture. We highlighted that inclusivity improves work–life balance, productivity and mental well-being for all employees.Any change, transformative or otherwise, is a process that takes perseverance, patience and determination. For any individual scientist to grow and flourish, they need a supportive environment, rich mentorship, a safe space and an enabling culture. It’s time for those factors to apply to all scientists equitably, no matter their gender, race, ethnicity or tribe. By fostering this mindset, we aim to reframe the narrative of our history and, in doing so, give all South African scientists their chance to thrive. More

Variation in environmental parameters across space and time in Bohai BayThe environmental parameters of samples collected near the Tianjin coastal area from different stations and seasons exhibited high temporal and spatial heterogeneity. The seawater temperature was 28.09 ± 0.53 °C in Aug, 17.48 ± 2.36 °C in May, and 19.55 ± 1.26 °C in Oct (Table 1). The seasonal variation in seawater temperature corresponded to the meteorological characteristics in Bohai Bay, with warm seawater in summer and relatively cool seawater in spring. The salinity was 29.69 ± 2.71‰ in Aug, 33.19 ± 0.33‰ in May, and 30.15 ± 1.63‰ in Oct. Seasonal variations in salinity may be mainly related to freshwater loading. According to the precipitation observed data of Bohai Bay in previous years, the rainfall amount and days in summer are the most19, which may lead to the increase in runoff and the relatively low salinity in summer. Chlorophyll a (Chl a) was highest in May, with lower levels in Aug and Oct. The dissolved inorganic nitrogen (DIN) was significantly higher in May and Aug than in Oct. The higher level of DIN in May and Aug may be related to terrestrial input and supply for the demand of phytoplankton growth. In October, the temperature and DIN content were both not suitable for phytoplankton growth, and Chl_a showed the lowest value. Spatially, the DIN distribution across the three seasons was rather similar, with high values observed in nearshore waters and low values in offshore waters (Dataset S1 & Fig. S1), which suggested that terrestrial input was an important source of DIN. The pH, soluble reactive phosphate (SRP) and chemical oxygen demand (COD) showed relatively higher values in October than in August and May, which may be caused by the dead phytoplankton release and terrestrial loadings through coasts and rivers. The dissolved oxygen (DO), conductivity and depth did not show significant variation among sampling times (Table 1), while the conductivity and depth had relatively higher values at offshore stations (Dataset S1) since the more remote the sampling water was, the greater the depth was in Bohai Bay and the closer it was to the open sea with higher salinity and conductivity. The ordination plot showed distinct partitioning of samples from nearshore and offshore sites along principal component axis 1 (PC1) (Fig. 1). The ordination plot could explain 73.49% of the total variation in the geo-physical–chemical parameters and revealed a linear positive correlation between different parameters (Fig. 1). AN, DIN, nitrate and Chl_a were most crucial in the partitioning of samples from May and the other 2 months; salinity, longitude, depth and conductivity were crucial for the partitioning of samples from offshore and nearshore stations; pH, COD, SRP, nitrite and temperature were crucial for the partitioning of samples from nearshore stations in August and October and samples from offshore stations. Overall, the principal component analysis (PCA) plot clearly showed both the temporal and spatial variation of the measured environmental parameters, indicating that complex biogeochemical processes and hydrodynamic conditions lead to the variation among sites and seasons.Table 1 The independent-samples t test of environmental variables and α-diversity among different months.Full size tableFigure 1Biplot of the principal component analysis (PCA) for environmental parameters in the seawater samples of the Bohai Bay coastal area across different seasons and sites. The two principal components (PC1 and PC2) explained 73.49% of the total variation in the environmental data and showed clear partitioning of offshore samples (in blue font) from other nearshore samples along PC1 and partitioning of May samples from August and October along PC2. The variables transparency and latitude were strongly correlated with PC1, and the variables ammonia nitrogen (AN), COD, pH, soluble reactive phosphate (SRP), and nitrite were strongly correlated with PC2. Chlorophyll a (Chl_a), dissolved inorganic nitrogen (DIN), nitrate and DO were mainly positively correlated with samples from May, while salinity, longitude, depth and conductivity were mainly positively correlated with offshore samples. Blue arrows represent environmental parameters, and circles in color represent sampling points.Full size imageProkaryotic α/β-diversity variationMeasures of α-diversity showed significant differences in shannon, evenness, faith_pd and OTU richness between samples from May/Aug and Oct (Fig. 2, Table 1). Principal coordinates analysis (PCoAs) based on weighted UniFrac (WUF) distance and unweighted UniFrac (UUF) distance showed that the PCC from different sampling months separated across the first and second principal coordinates (Fig. 3A-B). Both the analysis of similarity (ANOSIM) and permutational multivariate analysis of variance (PERMANOVA/ADONIS) results indicated that the prokaryotic communities varied significantly across different sampling months when using a WUF distance metric (ANOSIM, r = 0.709, P = 0.001; ADONIS, R2 = 40.0%, P = 0.001) and UUF distance metric (ANOSIM, r = 0.934, P = 0.001; ADONIS, R2 = 38.7%, P = 0.001). At the same time, the prokaryotic α– and β-diversity both showed high within-month variability in Aug (Figs. 2, 3C–D), which indicated that the community varied greatly among different sites in Aug.Figure 2Alpha diversity of shannon, eveness, faith_pd (phylogenetic diversity) and OTU richness value of the prokaryotic community of all the samples from different stations at different sampling times.Full size imageFigure 3Principal coordinate analysis (PCoA) based on unweighted (A) and weighted (B) UniFrac distances for prokaryotic communities in the surface waters; box plots showing the unweighted (C) and weighted (D) UniFrac distances among each station at different sampling times.Full size imageCorrelation between prokaryotic α/β-diversity and physical, chemical and geographic factorsThe α-diversity measurements exhibited significant positive correlations with temperature, pH, SRP, AN and un_ionN (Dataset S2). The correlation between α-diversity indexes and geo factors (longitude and latitude) was not strong or significant both in samples across the three sampling times or from each sampling time (Dataset S2).The environmental variation significantly correlated with β-diversity among the three seasons (r_weighted = 0.4558, r_unweighted = 0.4631, P = 0.001, Table 2), with pH, AN, temperature, un_ionN, COD, nitrite, SRP, salinity, DO and DIN as the main individual determinants. However, it did not show significant correlations with β-diversity at any sampling time except in Oct (Table S1).Table 2 Spearman’s rank correlation between environmental/spatial variability (Euclidean distance) and prokaryotic β-diversity (weighted/unweighted UniFrac distance) among all samples from different season.Full size tableThe geographic distance was not correlated with prokaryotic β-diversity (variation in community composition; r  0.05; Table 2) among the three sampling times. However, samples from Aug and Oct exhibited a significant correlation between β-diversity and geographic distance (Table S1).Factors driving the PCC variationPERMANOVA using the UUF/WUF distance indicated that temperature variation explained the largest part of community variation among the investigated factors (34.90%/19.83%, P = 0.001, Dataset S3), with AN (31.84%/13.56%, P = 0.001) and salinity (12.91%/6.21%, P = 0.001) as the second and third most significant sources of variation.The variance partitioning analysis (VPA) conducted on both UUF/WUF distances showed that almost 100% percent of the variation in PCC among all three sampling times was explained by the detected environmental factors. In May, no environmental or spatial factors could be selected as significantly explain the PCC variation; in Aug, the joint effects of environmental and spatial factors could explain 49.5% of the variation; in Oct, based on WUF distance, the spatial factors could purely explain 10.5%, environmental factors could purely explain 38.8%, their joint effects could explain 28.2%, and based on UUF distance, the joint effects of environmental factors and trend could explain 13.7% of the PCC variation. These results indicated dramatic shifts in the spatial or environmental factor effects on the PCC variation at different sampling times in Bohai Bay (Table 3).Table 3 Variance partitioning analysis of prokaryotic community in Bohai Bay according to seawater environmental factors and geospatial factors. The spatial factors including linear trend and PCNM variables. Forward selection procedures were used to select the best subset of environmental, trend, and PCNM variables explaining community variation, respectively. The community variation was calculated on the weighted and unweighted UniFrac distance matrix, respectively. Monte Carlo permutation test was performed on each set without the effect of the other by permuting samples freely (999 permutations).Full size tableDistinct seasonal features at the phylum and OTU levelsThere were notable differences in the proportions of various phyla among different seasons (sampling month). In May, there was a greater proportion of Alphaproteobacteria (41.41%), Planctomycetes (6.42%), Actinobacteria (3.86%), Firmicutes (1.48%), Acidobacteria (0.45%), TM7 (0.16%), Tenericutes (0.16%), OD1 (0.13%), and WPS-2 (0.09%) than in Aug and Oct, whereas Gammaproteobacteria (44.23%), GN02 (0.08%) and SAR406 (0.04%) were depleted in May and Aug but enriched in Oct. In Aug, Bacteroidetes (13.98%), Deltaproteobacteria (6.93%), Verrucomicrobia (4.5%), Chloroflexi (0.36%), Lentisphaerae (0.97%), TM6 (0.25%), Nitrospirae (0.08%), Chlamydiae (0.07%), Chlorobi (0.07%), Spirochaetes (0.04%) and OP8 (0.03%) were significantly enriched than in the other two sampling times (Duncan test; Table S2).At the OTU level, OTUs with relative abundance  > 0.01% (1040 OTUs) were used to perform the difference analysis, and 175 OTUs in May, 281 OTUs in Aug, and 210 OTUs in Oct were identified as seasonal specific OTUs (ssOTUs). The cooccurrence network showed that the ssOTUs were clustered separately (Fig. 4A). Furthermore, the separation of the three modules contained most of the ssOTUs specific to different seasons (Fig. 4A-B). All the ssOTUs of different seasons comprised a taxonomically broad set of prokaryotes at the phylum (phylum Proteobacteria is grouped at the class level) level (Fig. 4C) belonging to various phyla with different proportions. Betaproteobacteria, Verrucomicrobia, Gemmatimonadetes, Epsilonproteobacteria, PAUC34f., and Euryarchaeota did not show significant differences among the three sampling times at the phylum level, but features belonging to these phyla showed differences at the OTU level (Fig. 4C, Dataset S4). In addition, the phylum ssOTUs belonging to, such as Alphaproteobacteria, Gammaproteobacteria, Bacteroidetes, Actinobacteria, and Deltaproteobacteria, were not only enriched at one sampling time (Dataset S4) but also enriched at the other two sampling times (Fig. 4C, Dataset S4). These results revealed that different seasons do not strictly select specific microbial lineages at the phylum level, but a finer level analysis could more strictly reflect the seasonal variation.Figure 4Co-occurrence patterns of seasonal sensitive OTUs (A). Co-occurrence network visualizing significant correlations (ρ  > 0.7, P  0.01%. Different colors represent ssOTUs in May (green), Aug (red) and Oct (blue). Cumulative relative abundance (as counts per million, CPM; y-axis in × 1000) of all the sensitive modules in the networks (B). The phylum attribution of ssOTUs in each season (C). The y-axis is the percentage of the number of OTUs that belong to a particular phylum that accounts for the total number of all the OTUs.Full size imageRegression analysis between the relative abundance of modules to which the ssOTUs belonged and the environmental factors was also conducted, and module 1 abundance, to which the Aug-ssOTUs belonged, showed a significant positive correlation with temperature (R2 = 0.77, P = 6.609e−62), AN (R2 = 0.43, P = 7.416e−25), and un_ionN (R2 = 0.75, P = 1.366e−58) and a negative correlation with SRP (R2 = 0.81, P = 6.762e-17). This may be caused by the functional role of the microbes in Aug. In the Aug-ssOTUs, Deltaproteobacteria showed a higher ratio than in the other 2 months (Fig. 4c), and in the following functional analysis, Deltaproteobacteria contributed to the genes related to nitrogen fixation, which may help to explain why there was a positive correlation of Aug-ssOTUs to AN and un_ionN. The module 2 abundance to which the May-ssOTUs belonged showed a significant negative correlation with pH (R2 = 0.65, P = 4.026e−44), temperature (R2 = 0.19, P = 2.325e−10), un_ionN (R2 = 0.025, P = 0.01779), and SRP (R2 = 0.12, P = 4.104e−07) and a positive correlation with AN (R2 = 0.26, P = 5.174e−14). In the May-ssOTUs, the ratio of Alphaproteobacteria was the highest, and Alphaproteobacteria were reported to be pH-sensitive groups in marine environments20, which prefer neutral pH environments21. In this study, the pH in May was 8.04 ± 0.07, in Aug was 8.39 ± 0.09, in Oct was 8.38 ± 0.07, and the pH in May was the closest to neutral, and the ratio decreased with increasing pH in Oct and Aug. The abundance of module 3, to which the Oct-ssOTUs belonged, showed a significant positive correlation with SRP (R2 = 0.81, P = 0.16e-10) and pH (R2 = 0.054, P = 0.00075) and a negative correlation with temperature (R2 = 0.44, P = 2.276e−25), AN (R2 = 0.75, P = 4.51e−58), and un_ionN (R2 = 0.6, P = 3.995e-39) (Fig. S2). Phosphate has been identified to limit primary productivity22, which is of great importance in the structure of dominant bacterial taxa in marine environments23. In the Oct-ssOTUs, the ratio of Gammaproteobacteria was the highest, as reported. Gammaproteobacteria was significantly explained by SRP during the seasonal variation in the Western English Channel, with Rho equal to 0.7523, which suggested the sensitivity of it to SRP, and in that study, it also showed a negative correlation between temperature and Gammaproteobacteria and a positive correlation between SRP and Gammaproteobacteria. Although the correlation was not significant, the variation trend was consistent, which indicates that the phenomenon observed in this study was not unexpected. In addition, most ammonia-oxidizing bacteria belong to the Betaproteobacteria and Gammaproteobacteria classes are chemolithoautotrophs that oxidize ammonia to nitrite24. Gammaproteobacteria and Betaproteobacteria both had higher ratios in Oct-ssOTUs, and the functional prediction results also showed that pmoA/amoA and pmoB/amoB, which encode ammonia monooxygenase, were mainly contributed by OTUs from Gammaproteobacteria and Betaproteobacteria (Dataset S10). The utilization of ammonia may explain the negative correlation between the Oct-ssOTUs and AN.Community assembly processes across different sampling months and sitesBased on the analysis of phylogenetic turnover, unweighted βNTI mostly ranged from -2 to 2 across different sites at a single sampling time in May, Aug and Oct, revealing that PCC variations across different sampling sites at a single time were mostly affected by stochastic processes. The unweighted βNTI was greater than 2 across May–Aug, May–Oct and Aug-Oct (Fig. 5A), which revealed that the variations in PCC across different sampling times were mostly affected by deterministic processes. The RCbray values across any two sampling times were equal to 1, and in each sampling time, the RCbray values ranged from − 1 to 1 (Fig. 5B). Combining the βNTI and RCbray values, the community assembly processes were quantified at each sampling time and at any two sampling times. As shown in Fig. 5C, turning over of the community during different sampling times was mainly governed by selection; among the different sites in May and Oct, it was mainly governed by “undominated” processes; community turn over in Aug was mainly governed by the influence of “Dispersal Limitation”. These results indicated that the shifts in the assembly of prokaryotic communities during different sampling times were caused by strong “heterogeneous selection” (βNTI  > 2), and the community variation at each sampling time was mainly caused by stochastic processes.Figure 5Patterns of distribution of unweighted βNTI (A) and RCbray (B) values across different sampling times. Quantification of the features that impose community assembly processes in and among different sampling times. (C) Pie charts give the percent of turnover in community composition governed primarily by Selection acting alone (white fill), Dispersal Limitation (green line fill), Homogenizing Dispersal (blue line fill) and undominated process (cyan fill).Full size imagePrediction of the metabolic potential at different sampling timesThe NSTI scores of each sample ranged from 0.033 to 0.096, with a mean of 0.058 (Dataset S5). Microbial functions were detected in all the samples from the three sampling times, and it was found that the relative abundances of 242 pathways were significantly changed between samples from May and samples from Aug (Dataset S6). The relative abundances of 321 pathways were significantly changed between samples from May and Oct (Dataset S7), and the relative abundances of 370 pathways were significantly changed between samples from Aug and Oct (Dataset S8).Genes related to energy metabolism were given more attention. For nitrogen metabolism genes relevant with nitrogen fixation (nifD, nifK) were detected only enriched in Aug, while genes relevant with nitrate reduction and denitrification (narG, narZ, nxrA, narH, narY, nxrB, narI, narV, nirD, nasA, nasB) were detected enriched in May, genes related with ammonia oxidation were both detected enriched in Oct and Aug. For sulfur metabolism, genes relevant with thiosulfate oxidation (soxA, soxB, soxC, soxX, soxY and soxZ) were only enriched in Aug, while genes relevant with sulfate and sulfite reduction (cysNC, aprA, aprB, cysJ, cysI, cysK, dsrA) were detected enriched in May and Oct (Fig. 6).Figure 6The LEfSe analysis indicated significantly differential abundances of PICRUSt predicted genes relevant to energy metabolism in different months of samples.Full size imageProkaryotic taxa contributed to the significantly varied functional genes related to nitrogen and sulfur metabolism at different sampling times. At the species level, the taxa contributing to nifK and nifD mainly belonged to Deltaproteobacteria and Firmicutes, and the taxa contributing to the sox-series genes mainly belonged to Alphaproteobacteria and Gammaproteobacteria (Fig. S3). The denitrification-related functional genes that were enriched in May were mainly contributed by members from Alphaproteobacteria and Gammaproteobacteria. The taxa contributing to dsrA, aprA and aprB were mainly from Deltaproteobacteria, including members of Desulfarculaceae, Desulfobacteraceae, Desulfobulbaceae, Desulfovibrionaceae and Syntrophobacteraceae (Fig. S4). More

Study sites and coloniesAll the experiments were carried out during the apple and pear blooming seasons of 2007, 2008, 2011, 2013 and 2014 in different locations of the province of Rio Negro, Argentina, while some laboratory experiments performed in the city of Buenos Aires. We used individual foragers of Apis mellifera L. and their colonies containing a mated queen, brood, and food reserves in ten-frame Langstroth hives. All beehives used had similar sizes and the same management history from the beekeeper. The honey bees studied belonged to commercial Langstroth-type hives rented to pollinate these plots. Each hive had a fertilized queen, 3 or 4 capped brood frames, reserves and approximately 15,000 individuals56.Testing generalization of memories from pear mimic odours to pear and apple natural floral scentsThe absolute conditioning assays were performed in the laboratories of the School of Exacts and Natural Sciences of the University of Buenos Aires (34° 32′ S, 58° 26′ W), Buenos Aires, Argentina. We used honey bee foragers collected at the entrance of the hives settle in the experimental field of the School of Exacts and Natural Sciences. The apple (‘Granny Smith’ and ‘Red Delicious’ varieties) and pear (‘Packham’ and ‘D’anjou’ varieties) bud samples that we used as conditioned stimuli (CS) during the conditioning were collected at the end of the blossom of 2011 in Ingeniero Huergo (39° 03′ 27.5″ S; 67° 13′ 53.5″ W), province of Río Negro, Argentina, and taken to the laboratory in the city of Buenos Aires, Argentina, to be used within the following 2 days.We first developed the three different synthetic mixtures (PM, PMI and PMII) that could be generalized to the fragrance of the pear flower by foraging bees. The pear synthetic mixtures were formulated considering the previously reported volatile profile of pear blossoms57. Then, we chose the synthetic mixture most perceptually similar to the pear flower fragrance and measured its generalisation response to the apple flower fragrance to test the compounds’ specificity. The chemical compounds used to prepare the different synthetic mixtures for the behavioural assays were obtained from Sigma-Aldrich, Steinheim, Germany. The compounds used for the three pear mixtures (PM, PMI and PMII) were composed by alpha-pinene, 2-ethyl-hexanol, (R)-(+)-limonene, and (±)-linalool. For details of the PM and mixture proportions see Patent PCT/IB2018/05555058.To test generalization, we took advantages of the fact that honey bees reflexively extend their proboscises when sugar solution is applied to their antennae59. The proboscis extension reflex (PER) can be used to condition bees to an odour if a neutral olfactory stimulus (CS) is paired with a sucrose reward as unconditioned stimulus, US60. Conditioned honey bees extend their proboscises towards the odour alone, a response that indicates that this stimulus has been learned and predicts the oncoming food reward. Conditioned bees can generalize such a learned response to a novel odour if it is perceived like the conditioned one (CS). Then we performed three absolute PER conditionings where we paired each of the three PMs with a sucrose-water solution (30%) reward along three learning trials (exp. 4.2a). Afterwards, pear floral scent was presented as novel odour to test generalization. Based on the generalization level to the pear odour, we chose the synthetic mixture that showed the highest generalisation towards pear flower fragrance, and we used it in all the experiments that follow. In an additional 3-trial PER conditioning with the chosen mixture, we quantified generalisation towards both the pear and apple fragrances as novel stimuli (exp. 4.2b).The experimental bees were all foragers, captured from colonies that had no access to any pear and/or apple tree, hence completely naïve for the CSs. Immediately after capture, bees were anaesthetized at 4 °C and harnessed in metal tubes so that they could only move their mouthparts and antennae60. They were fed 30% weight/weight unscented sucrose solution for about three seconds and kept in a dark incubator (30 °C, 55% relative humidity) for about two hours. Only those bees that showed the unconditioned response (the reflexive extension of the proboscis after applying a 30% w/w sucrose solution to the antennae) and did not respond to the mechanical air flow stimulus were used. Trials lasted 46 s and presented three steps: 20 s of clean air, 6 s of odour presentation (CS) and the last 20 s of clean air. During rewarded trials (CS), the reward (US, a drop of 30% w/w sucrose solution) was delivered upon the last 3 s of CS presentation. The synthetic mixtures (PM) were delivered in a constant air flow (15 ml/s) that passed through a 1 ml syringe containing 4 µl of the synthetic mixture on a small strip of filter paper. On the other hand, pear and apple floral volatiles were swept from a 100 g of fresh pear buds (var. ‘D’Anjou’ and ‘Packham’) or apple buds (var. ‘Granny Smith’, ‘Gala’ and ‘Red Delicious’) inside a kitasato by means of an air flow (54 ml/s).Testing discrimination between mimics and natural floral scentsThe differential conditioning assays were performed in a field laboratory in Ingeniero Huergo, province of Río Negro, Argentina. Conditioning trials with AM as CS were carried out in September 2007 and 2008, prior to the beginning of flowering of the fruit trees. Conditioning trials with PM as CS were carried out in September 2011 in the same area (Ingeniero Huergo, province of Río Negro, Argentina). Apple and pear bud samples used as CS were collected in plots that start blooming located around Ingeniero Huergo, but distant (more than 1 km) from the plot where we collected the bees. The bud samples presented the following varieties: M. domesticus sp., ‘Granny Smith’, ‘Gala’, and ‘Red Delicious’; P. communis sp., ‘Packham’ and ‘D’Anjou’.With the aim to develop a synthetic mixture that presents difficult to discriminate with the fragrance of the apple flower by foraging bees, an apple synthetic mixture (AM) was formulated considering the previously reported volatile profile of apple blossoms61. The chemical compounds used to prepare the apple synthetic mixtures for the behavioural assays were obtained from Sigma-Aldrich, Steinheim, Germany. Apple mimic (AM) was composed by benzaldehyde, limonene and citral. For details of the AM proportions see Patent AR2011010244162. Jasmine mimic (JM) was a commercial extract obtained from Firmenich S.A.I.C. y F, Argentina.If the synthetic mixture chosen were perceptually similar to the apple flower fragrance, experimental bees should have difficult to discriminate to the apple flower fragrance to test the compounds’ specificity. Thus, we performed differential PER conditioning between synthetic mixtures (AM and Jasmine mimic, JM) or between synthetic mixtures (AM or JM) and the apple natural fragrance. We followed a differential PER conditioning34 to assess to what extent the bees were able to discriminate the synthetic mimics from their natural flower scents. PER differential conditioning consisted of four pairs of trials, four rewarded trials (CS+) and four non-rewarded trials (CS−) that were presented in a pseudo-randomized manner. Conditionings were performed using the synthetic mixtures PM and AM and the natural floral scents, pear and apple, either as CS+ and CS−. We followed the same procedure that in 3.3 to capture the bees and to present the stimuli during trials.Feeding protocolWe used the offering of scented sucrose solution in the hive as a standardized procedure to establish long-term olfactory memory in honey bees23,24,24,26,63. Scented sucrose solution was obtained by diluting 50 µl of PM or AM per litre of sucrose solution (50% weight/weight, henceforth: w/w). For the ‘apple’ series, colonies were fed 1500 ml of sugar solution offered in an internal plastic feeder for 2 days, about 3 days before the apple trees began to bloom. For the ‘pear’ series, hives were fed 500 ml of sugar solution that we spread over the top of the central frames. Both feeding procedures have been found to be functional for establishing olfactory in-hive memories26. Depending on the pear varieties, the scented sucrose solution was offered when the pear trees were 10–40% in bloom.Colony activityThe effects of the AM-treatment on colony nest entrance activity were studied in 18 colonies located in an agricultural setting of apple and pear trees in Ingeniero Huergo, on an 8-ha plot, half of which was planted with apple trees (varieties: ‘Granny Smith’, ‘Gala’ and ‘Red Delicious’) and the other 4 ha with pear trees (varieties: ‘Packham’ and ‘D’anjou’). The effect of the PM-treatment on colony activity was studied in 14 colonies located in three adjoining pear plots (total surface: 8 ha) in Otto Krause (39° 06′ 22″ S 66° 59′ 46″ O, Supplementary Fig. S5), province of Río Negro, Argentina. The varieties of these plots corresponded to ‘Packham’ and ‘Williams’. Pollen collection (exp. 4.5.2) was also studied in colonies located in these plots.We focused on the nest entrance activity since once the first successful foragers return to the hive and display dances and/or unload the food collected, it promotes the activation or reactivation of inactive foragers and, in a minor proportion, those hive mates ready to initiate foraging tasks39,65,66,67,67. Then, we choose number of incoming bees as an indicator of colony foraging activity, since most of these bees are expected to return from foraging sites33. Thus, we compared the activity level at the nest entrance between 7 SS + PM-treated colonies and 7 SS-treated colonies. We also compared the nest entrance activity level between 5 colonies treated with SS + AM and 5 colonies fed with SS. This activity value was estimated by the number of incoming foragers at the entrance of the hive for one minute, every morning at the same time (10:30 a.m.) during the entire experiment (9 consecutive days). A first measurement was done one day before feeding the colonies (used as covariate) and 7 measurements afterwards.We measured the amount of pollen loads collected by two colonies: one fed with SS + PM and one fed with SS. Pollen loads were collected using conventional pollen traps (frontal-entrance trap), consisting of a wooden structure with a removable metal mesh inside. Pollen samples were collected for 3 days, two hours per day during the late morning, 3, 7 and 8 days after the offering of SS + PM or SS. Pollen pellets identified based on pollen colour as coming from the pear flower or from other species were separated and counted. In addition, we estimated the weight of pear pollen loads during a 5 days period, from 6 to 10 days after the offering of scented or unscented sucrose solution. To reduce measurement error, pollen loads were weighed in groups of 10.Crop yieldPear crop yield was studied in pear plots in General Roca (39° 02′ 00″ S; 67° 35′ 00″ O, Supplementary Fig. S4, Supplementary Table S3), province of Río Negro, Argentina. In an area of 15.2 ha (4 plots of 3.8 ha each), 45 beehives were equidistantly located in groups. We measured the number of fruits per tree set of 30 trees in the surrounding areas of the PM-treated colonies (2 groups of 8 hives) and control colonies (2 groups of 8 hives). A third group category contained 13 untreated colonies. The varieties of the pear trees were ‘D’Anjou’ and ‘Packham’.Apple crop yield estimated by means of number of fruits per plant was studied in General Roca (Supplementary Fig. S2, Supplementary Table S1), province of Río Negro, Argentina. We measured fruit set in the two plots that covered a surface of 3.8 ha and contained a total of 74 colonies distributed in groups (the control plot, 39 SS-treated-colonies treated with SS; and the treated plot, 35 SS + AM-treated-colonies treated with SS + AM). The varieties of the apple trees were ‘Red Delicious’ (clone 1), ‘Royal Gala’ and ‘Granny Smith’.A second studied on apple fruit yield by means of kg of fruits per hectare was performed in Coronel Belisle (39° 11′ 00″ S 65° 59′ 00″ O, Supplementary Fig. S3, Supplementary Table S2), province of Río Negro, Argentina. Four apple plots with ‘Granny Smith’, ‘Hi Early’ and ‘Red Delicious’, clone 1 varieties of 15.4 ha each were randomly assigned to different treatments (treated plot 1, 40 SS + AM-treated-hives treated with SS + AM; treated plot 2, 40 SS + AM-treated-hives treated with SS + AM; control plot 1, 40 SS-treated-hives treated with SS; control plot 2, 40 SS-treated-hives treated with SS).During the fruit harvest, the fruit yield was estimated in the surroundings (150 m around) of two groups of 8 colonies each. We fed one group SS + PM and the other unscented sucrose solution (SS). Yield was estimated as the number of fruits per trees in 30 randomly selected trees within each area, alternating the counts between the North and South faces of the plots. Following the same procedure, we also estimated the number of fruits per trees in the surroundings of two groups of 14 colonies each that pollinated apple crops. Again, we fed one group SS + AM and the other SS. Additionally, a total of 218 colonies in General Roca and 180 colonies in Coronel Belisle have been separated in the two experimental groups, in which yield had been provided by the producer and expressed in kg of fruits per ha. It is worth remarking that in some plots the distance between treated and control beehive groups was around 300 m, suggesting that might have been overlapping flying areas between treated and control hives. Additionally, the apple fields studied in the surrounding of Coronel Belisle, presented many trees without flowers. It was considered that the absence of flowers in numerous trees would bias the counts performed in those fields. Then, to quantify this situation, which might be associated with the masting phenomenon68, samples with the proportions of trees without flowers for every 20 trees in each plot was done. Trees that had between 80 and 100% of their surface devoid of flowers were considered “without flowers” trees, and “trees with available flowers” those that had more than 20% of their surface covered with flowers. An average of 30% of the trees within these plots were devoid of flowers. Thus, a correction factor was considered to evaluate the yield data provided by the grower per plot analysed (Supplementary Table S4).StatisticsAll statistical analyses were performed with R Core Team 201969. For Experiment 4.2 and 4.3, we analysed PER proportion by means of a binomial multiplicative generalized linear mixed model using the “glmer” function of the ‘lme4’ package70.For experiment 4.2a we considered the pear mimics (three-level factor corresponding to PM, PMI and PMII) and the event (two-level factor corresponding to 3rd trial and test) as fixed factors and each “bee” as a random factor.For experiment 4.2b we considered the tested odours (three-level factor corresponding to Apple, Pear and PM) as fixed factors.For experiment 4.3 we considered the tested odours (two-level factor corresponding to CS+ and CS−) as fixed factors. Post hoc contrasts were conducted on models to assess effects and significance between fixed factors using the “emmeans” function of the ‘emmeans’ package version 1.7.071 with a significance level of 0.05.For experiment 4.5.1 we analysed “rate of incoming bees” using a generalized linear mixed model. As Poisson model for incoming bees was overdispersed72, we used a negative binomial distribution using the ‘glmmTMB’ package (function ‘glmmTMB’73. We considered “treatment” [two-level factor corresponding to SS + AM (or SS + PM) and SS], “days” (7-level factor corresponding to the date after treatment), the rate of incoming bees before the offering of food (to control for pre-existing colony differences) as covariate (a quantitative fixed effects variable), and “colony” as a random factor.For experiment 4.6, we analysed fruits per trees by means of a negative binomial multiplicative generalized linear mixed model using the “log” function of the ‘ml’ package70. Post hoc contrasts were conducted on models to assess effects and significance between fixed factors using the “emmeans” function of the ‘emmeans’ package version 1.8.071 with a significance level of 0.05. For experiment 4.6b we analysed “yield” (as weight of fruits per unit area) using a general linear mixed model. We checked homoscedasticity and normality assumptions (Levene and Shapiro–Wilk tests, respectively). We considered “treatment” (two-level factor corresponding to SS + AM and SS) and “apple varieties” (3-level factor corresponding to Hi Early, Granny Smith and Chañar 28) as fixed factors and “location” (2-level factor corresponding to General Roca and Coronel Belisle) as random factors. More