Philippa Kaur

Flemming, H.-C. & Wuertz, S. Bacteria and archaea on Earth and their abundance in biofilms. Nat. Rev. Microbiol. 17, 247–260 (2019).CAS 
PubMed 
Article 

Google Scholar 
Inagaki, F. et al. Exploring deep microbial life in coal-bearing sediment down to ~2.5 km below the ocean floor. Science 349, 420–424 (2015).CAS 
PubMed 
Article 

Google Scholar 
Bar-On, Y. M., Phillips, R. & Milo, R. The biomass distribution on Earth. Proc. Natl Acad. Sci. USA 115, 6506–6511 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
D’Hondt, K. et al. Microbiome innovations for a sustainable future. Nat. Microbiol. 6, 138–142 (2021).PubMed 
Article 
CAS 

Google Scholar 
Lloyd, K. G., Steen, A. D., Ladau, J., Yin, J. & Crosby, L. Phylogenetically novel uncultured microbial cells dominate earth microbiomes. mSystems 3, e00055-18 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Lewis, W. H., Tahon, G., Geesink, P., Sousa, D. Z. & Ettema, T. J. G. Innovations to culturing the uncultured microbial majority. Nat. Rev. Microbiol. 19, 225–240 (2021).CAS 
PubMed 
Article 

Google Scholar 
Nayfach, S. et al. A genomic catalog of Earth’s microbiomes. Nat. Biotechnol. 39, 499–509 (2021).CAS 
PubMed 
Article 

Google Scholar 
Rinke, C. et al. Insights into the phylogeny and coding potential of microbial dark matter. Nature 499, 431–437 (2013).CAS 
PubMed 
Article 

Google Scholar 
Hugenholtz, P. Exploring prokaryotic diversity in the genomic era. Genome Biol. 3, REVIEWS0003 (2002).PubMed 
PubMed Central 
Article 

Google Scholar 
Engelen, B. & Imachi, H. Cultivation of subseafloor prokaryotic life in developments in marine geology. In Earth and Life Processes Discovered from Subseafloor Environment. Vol. 7 (eds. Stein, R., Blackman, D., Inagaki, F. & Larsen, H.-L.) 197–209 (Elsevier, 2014).Baker, B. J., Appler, K. E. & Gong, X. New microbial biodiversity in marine sediments. Annu. Rev. Mar. Sci. 13, 161–175 (2021).Article 

Google Scholar 
Hoshino, T. et al. Global diversity of microbial communities in marine sediment. Proc. Natl Acad. Sci. USA 117, 27587–27597 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tahon, G., Geesink, P. & Ettema, T. J. G. Expanding archaeal diversity and phylogeny: past, present, and future. Annu. Rev. Microbiol. 75, 359–381 (2021).PubMed 
Article 
CAS 

Google Scholar 
Bhattarai, S., Cassarini, C. & Lens, P. N. L. Physiology and distribution of archaeal methanotrophs that couple anaerobic oxidation of methane with sulfate reduction. Microbiol. Mol. Biol. Rev. 83, e00074-18 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Krukenberg, V. et al. Gene expression and ultrastructure of meso- and thermophilic methanotrophic consortia. Environ. Microbiol. 20, 1651–1666 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Wegener, G., Krukenberg, V., Ruff, S. E., Kellermann, M. Y. & Knittel, K. Metabolic capabilities of microorganisms involved in and associated with the anaerobic oxidation of methane. Front. Microbiol. 7, 46 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Agrawal, L. K. et al. Treatment of raw sewage in a temperate climate using a UASB reactor and the hanging sponge cubes process. Water Sci. Technol. 36, 433–440 (1997).CAS 
Article 

Google Scholar 
Tyagi, V. K. et al. Future perspectives of energy saving down-flow hanging sponge (DHS) technology for wastewater valorization—a review. Rev. Environ. Sci. Biotechnol. 20, 389–418 (2021).Article 

Google Scholar 
Namita Maharjan, N. et al. Downflow hanging sponge system: a self-sustaining option for wastewater treatment. In Wastewater Treatment (IntechOpen, London, UK, 2020) Available at https://www.intechopen.com/online-first/74120Nurmiyanto, A. & Ohashi, A. Downflow hanging sponge (DHS) reactor for wastewater treatment—a short review. MATEC Web Conf. 280, 05004 (2019).CAS 
Article 

Google Scholar 
Hatamoto, M., Okubo, T., Kubota, K. & Yamaguchi, T. Characterization of downflow hanging sponge reactors with regard to structure, process function, and microbial community compositions. Appl. Microbiol. Biotechnol. 102, 10345–10352 (2018).CAS 
PubMed 
Article 

Google Scholar 
Tandukar, M., Uemura, S., Ohashi, A. & Harada, H. Combining UASB and the ‘fourth generation’ down-flow hanging sponge reactor for municipal wastewater treatment. Water Sci. Technol. 53, 209–218 (2006).CAS 
PubMed 
Article 

Google Scholar 
Chuang, H.-P. et al. Microbial community that catalyzes partial nitrification at low oxygen atmosphere as revealed by 16S rRNA and amoA genes. J. Biosci. Bioeng. 104, 525–528 (2007).CAS 
PubMed 
Article 

Google Scholar 
Imachi, H. et al. Cultivation of methanogenic community from subseafloor sediments using a continuous-flow bioreactor. ISME J. 5, 1913–1925 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Aoki, M. et al. A long-term cultivation of an anaerobic methane-oxidizing microbial community from deep-sea methane-seep sediment using a continuous-flow bioreactor. PLoS ONE 9, e105356 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Kato, S. et al. Biotic manganese oxidation coupled with methane oxidation using a continuous-flow bioreactor system under marine conditions. Water Sci. Technol. 76, 1781–1795 (2017).CAS 
PubMed 
Article 

Google Scholar 
Imachi, H. et al. Cultivable microbial community in 2-km-deep, 20-million-year-old subseafloor coalbeds through ~1000 days anaerobic bioreactor cultivation. Sci. Rep. 9, 2305 (2019).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Imachi, H. et al. Pelolinea submarina gen. nov., sp. nov., an anaerobic, filamentous bacterium of the phylum Chloroflexi isolated from subseafloor sediment. Int. J. Syst. Evol. Microbiol. 64, 812–818 (2014).CAS 
PubMed 
Article 

Google Scholar 
Imachi, H. et al. Sedimentibacter acidaminivorans sp. nov., an anaerobic, amino-acid-utilizing bacterium isolated from marine subsurface sediment. Int. J. Syst. Evol. Microbiol. 66, 1293–1300 (2016).CAS 
PubMed 
Article 

Google Scholar 
Imachi, H. et al. Isolation of an archaeon at the prokaryote-eukaryote interface. Nature 577, 519–525 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Miyazaki, M. et al. Sphaerochaeta multiformis sp. nov., an anaerobic, psychrophilic bacterium isolated from subseafloor sediment, and emended description of the genus Sphaerochaeta. Int. J. Syst. Evol. Microbiol. 64, 4147–4154 (2014).PubMed 
Article 
CAS 

Google Scholar 
Miyazaki, M. et al. Spirochaeta psychrophila sp. nov., a psychrophilic spirochaete isolated from subseafloor sediment, and emended description of the genus Spirochaeta. Int. J. Syst. Evol. Microbiol. 64, 2798–2804 (2014).CAS 
PubMed 
Article 

Google Scholar 
Nakahara, N. et al. Aggregatilinea lenta gen. nov., sp. nov., a slow-growing, facultatively anaerobic bacterium isolated from subseafloor sediment, and proposal of the new order Aggregatilineales ord. nov. within the class Anaerolineae of the phylum Chloroflexi. Int. J. Syst. Evol. Microbiol. 69, 1185–1194 (2019).CAS 
PubMed 
Article 

Google Scholar 
Spang, A. et al. Complex archaea that bridge the gap between prokaryotes and eukaryotes. Nature 521, 173–179 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zaremba-Niedzwiedzka, K. et al. Asgard archaea illuminate the origin of eukaryotic cellular complexity. Nature 541, 353–358 (2017).CAS 
PubMed 
Article 

Google Scholar 
Liu, Y. et al. Expanded diversity of Asgard archaea and their relationships with eukaryotes. Nature 593, 553–557 (2021).CAS 
PubMed 
Article 

Google Scholar 
Ruff, S. E. et al. Global dispersion and local diversification of the methane seep microbiome. Proc. Natl Acad. Sci. USA 112, 4015–4020 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Chadwick, G. et al. Comparative genomics reveals electron transfer and syntrophic mechanisms differentiating methanotrophic and methanogenic archaea. PLoS Biol. 20, e3001508 (2022).PubMed 
PubMed Central 
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 
Jørgensen, B. B. & Boetius, A. Feast and famine—microbial life in the deep-sea bed. Nat. Rev. Microbiol. 5, 770–781 (2007).PubMed 
Article 
CAS 

Google Scholar 
Hoehler, T. M. & Jørgensen, B. B. Microbial life under extreme energy limitation. Nat. Rev. Microbiol. 11, 83–94 (2013).CAS 
PubMed 
Article 

Google Scholar 
Schink, B. & Stams, A. J. M. Syntrophism among prokaryotes. In The Prokaryotes: Prokaryotic Communities and Ecophysiology (eds. Rosenberg, E. et al.) 471–493 (Springer, 2013).de Bok, F. A. M. et al. The first true obligately syntrophic propionate-oxidizing bacterium, Pelotomaculum schinkii sp. nov., co-cultured with Methanospirillum hungatei, and emended description of the genus Pelotomaculum. Int. J. Syst. Evol. Microbiol. 55, 1697–1703 (2005).PubMed 
Article 
CAS 

Google Scholar 
Imachi, H. et al. Pelotomaculum propionicicum sp. nov., an anaerobic, mesophilic, obligately syntrophic, propionate-oxidizing bacterium. Int. J. Syst. Evol. Microbiol. 57, 1487–1492 (2007).PubMed 
Article 

Google Scholar 
Qiu, Y.-L. et al. Syntrophorhabdus aromaticivorans gen. nov., sp. nov., the first cultured anaerobe capable of degrading phenol to acetate in obligate syntrophic associations with a hydrogenotrophic methanogen. Appl. Environ. Microbiol. 74, 2051–2058 (2008).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Matsushita, S. et al. Anti-bacterial effects of MnO2 on the enrichment of manganese-oxidizing bacteria in downflow hanging sponge reactors. Microbes Environ. 35, ME20052 (2020).PubMed Central 

Google Scholar 
Momper, L. et al. Rectinema subterraneum sp. nov, a chemotrophic spirochaete isolated from the deep terrestrial subsurface. Int. J. Syst. Evol. Microbiol. 70, 4739–4747 (2020).CAS 
PubMed 
Article 

Google Scholar 
Chuang, H.-P., Ohashi, A., Imachi, H., Tandukar, M. & Harada, H. Effective partial nitrification to nitrite by down-flow hanging sponge reactor under limited oxygen condition. Water Res. 41, 295–302 (2007).CAS 
PubMed 
Article 

Google Scholar 
Hatamoto, M., Koshiyama, Y., Kindaichi, T., Ozaki, N. & Ohashi, A. Enrichment and identification of methane-oxidizing bacteria by using down-flow hanging sponge bioreactors under low methane concentration. Ann. Microbiol. 61, 683–687 (2010).Article 
CAS 

Google Scholar 
Hatamoto, M., Yamamoto, H., Kindaichi, T., Ozaki, N. & Ohashi, A. Biological oxidation of dissolved methane in effluents from anaerobic reactors using a down-flow hanging sponge reactor. Water Res. 44, 1409–1418 (2010).CAS 
PubMed 
Article 

Google Scholar 
Hatamoto, M., Miyauchi, T., Kindaichi, T., Ozaki, N. & Ohashi, A. Dissolved methane oxidation and competition for oxygen in down-flow hanging sponge reactor for post-treatment of anaerobic wastewater treatment. Bioresour. Technol. 102, 10299–10304 (2011).CAS 
PubMed 
Article 

Google Scholar 
Hatamoto, M. et al. Potential of nitrous oxide conversion in batch and down-flow hanging sponge bioreactor systems. Sustain. Environ. Res. 24, 117–128 (2014).
Google Scholar 
Cao, L. T. T. et al. Biological oxidation of Mn(II) coupled with nitrification for removal and recovery of minor metals by downflow hanging sponge reactor. Water Res. 68, 545–553 (2015).CAS 
PubMed 
Article 

Google Scholar 
Yamaguchi, T. et al. A novel approach for toluene gas treatment using a downflow hanging sponge reactor. Appl. Microbiol. Biotechnol. 102, 5625–5634 (2018).CAS 
PubMed 
Article 

Google Scholar 
Matsushita, S. et al. Production of biogenic manganese oxides coupled with methane oxidation in a bioreactor for removing metals from wastewater. Water Res. 130, 224–233 (2018).CAS 
PubMed 
Article 

Google Scholar 
Onodera, T. et al. Characterization of the retained sludge in a down-flow hanging sponge (DHS) reactor with emphasis on its low excess sludge production. Bioresour. Technol. 136, 169–175 (2013).CAS 
PubMed 
Article 

Google Scholar 
Miyaoka, Y., Hatamoto, M., Yamaguchi, T. & Syutsubo, K. Eukaryotic community shift in response to organic loading rate of an aerobic trickling filter (down-flow hanging sponge reactor) treating domestic sewage. Microb. Ecol. 73, 801–814 (2016).PubMed 
Article 
CAS 

Google Scholar 
Park, M.-O., Ikenaga, H. & Watanabe, K. Phage diversity in a methanogenic digester. Microb. Ecol. 53, 98–103 (2006).PubMed 
Article 
CAS 

Google Scholar 
Wu, Q. & Liu, W.-T. Determination of virus abundance, diversity and distribution in a municipal wastewater treatment plant. Water Res. 43, 1101–1109 (2009).CAS 
PubMed 
Article 

Google Scholar 
Chien, I.-C., Meschke, J. S., Gough, H. L. & Ferguson, J. F. Characterization of persistent virus-like particles in two acetate-fed methanogenic reactors. PLoS One 8, e81040 (2013).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Girguis, P. R., Orphan, V. J., Hallam, S. J. & DeLong, E. F. Growth and methane oxidation rates of anaerobic methanotrophic archaea in a continuous-flow bioreactor. Appl. Environ. Microbiol. 69, 5472–5482 (2003).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zhang, Y., Maignien, L., Zhao, X., Wang, F. & Boon, N. Enrichment of a microbial community performing anaerobic oxidation of methane in a continuous high-pressure bioreactor. BMC Microbiol. 11, 137 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Sauer, P., Glombitza, C. & Kallmeyer, J. A system for incubations at high gas partial pressure. Front. Microbiol. 3, 25 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
Bhattarai, S. et al. Enrichment of sulfate reducing anaerobic methane oxidizing community dominated by ANME-1 from Ginsburg Mud Volcano (Gulf of Cadiz) sediment in a biotrickling filter. Bioresour. Technol. 259, 433–441 (2018).CAS 
PubMed 
Article 

Google Scholar 
Cassarini, C. et al. Enrichment of anaerobic methanotrophs in biotrickling filters using different sulfur compounds as electron acceptor. Environ. Eng. Sci. 36, 431–443 (2018).Article 
CAS 

Google Scholar 
Cassarini, C. et al. Anaerobic methane oxidation coupled to sulfate reduction in a biotrickling filter: reactor performance and microbial community analysis. Chemosphere 236, 124290 (2019).CAS 
PubMed 
Article 

Google Scholar 
Adams, M. M., Hoarfrost, A. L., Bose, A., Joye, S. B. & Girguis, P. R. Anaerobic oxidation of short-chain alkanes in hydrothermal sediments: potential influences on sulfur cycling and microbial diversity. Front. Microbiol. 4, 110 (2013).PubMed 
PubMed Central 
Article 

Google Scholar 
Machdar, I., Sekiguchi, Y., Sumino, H., Ohashi, A. & Harada, H. Combination of a UASB reactor and a curtain type DHS (downflow hanging sponge) reactor as a cost-effective sewage treatment system for developing countries. Water Sci. Technol. 42, 83–88 (2000).CAS 
Article 

Google Scholar 
Judd, S. The status of membrane bioreactor technology. Trends Biotechnol. 26, 109–116 (2008).CAS 
PubMed 
Article 

Google Scholar 
Smith, A. L. et al. Perspectives on anaerobic membrane bioreactor treatment of domestic wastewater: a critical review. Bioresour. Technol. 122, 149–159 (2012).CAS 
PubMed 
Article 

Google Scholar 
Meulepas, R. J. W. et al. Enrichment of anaerobic methanotrophs in sulfate-reducing membrane bioreactors. Biotechnol. Bioeng. 104, 458–470 (2009).CAS 
PubMed 
Article 

Google Scholar 
Jagersma, G. C. et al. Microbial diversity and community structure of a highly active anaerobic methane-oxidizing sulfate-reducing enrichment. Environ. Microbiol. 11, 3223–3232 (2009).CAS 
PubMed 
Article 

Google Scholar 
Lettinga, G., van Velsen, A. F. M., Hobma, S. W., de Zeeuw, W. & Klapwijk, A. Use of the upflow sludge blanket (USB) reactor concept for biological wastewater treatment, especially for anaerobic treatment. Biotechnol. Bioeng. 22, 699–734 (1980).CAS 
Article 

Google Scholar 
Sekiguchi, Y., Kamagata, Y., Nakamura, K., Ohashi, A. & Harada, H. Fluorescence in situ hybridization using 16S rRNA-targeted oligonucleotides reveals localization of methanogens and selected uncultured bacteria in mesophilic and thermophilic sludge granules. Appl. Environ. Microbiol. 65, 1280–1288 (1999).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tawfik, A., Ohashi, A. & Harada, H. Sewage treatment in a combined up-flow anaerobic sludge blanket (UASB)-down-flow hanging sponge (DHS) system. Biochem. Eng. J. 29, 210–219 (2006).CAS 
Article 

Google Scholar 
Onodera, T. et al. Development of a sixth-generation down-flow hanging sponge (DHS) reactor using rigid sponge media for post-treatment of UASB treating municipal sewage. Bioresour. Technol. 152, 93–100 (2014).CAS 
PubMed 
Article 

Google Scholar 
Tandukar, M., Ohashi, A. & Harada, H. Performance comparison of a pilot-scale UASB and DHS system and activated sludge process for the treatment of municipal wastewater. Water Res. 41, 2697–2705 (2007).CAS 
PubMed 
Article 

Google Scholar 
Hallam, S. J. et al. Reverse methanogenesis: testing the hypothesis with environmental genomics. Science 305, 1457–1462 (2004).CAS 
PubMed 
Article 

Google Scholar 
Nauhaus, K., Albrecht, M., Elvert, M., Boetius, A. & Widdel, F. In vitro cell growth of marine archaeal-bacterial consortia during anaerobic oxidation of methane with sulfate. Environ. Microbiol. 9, 187–196 (2007).CAS 
PubMed 
Article 

Google Scholar 
Wegener, G., Niemann, H., Elvert, M., Hinrichs, K. & Boetius, A. Assimilation of methane and inorganic carbon by microbial communities mediating the anaerobic oxidation of methane. Environ. Microbiol. 10, 2287–2298 (2008).CAS 
PubMed 
Article 

Google Scholar 
Hu, H., Natarajan, V. P. & Wang, F. Towards enriching and isolation of uncultivated archaea from marine sediments using a refined combination of conventional microbial cultivation methods. Mar. Life Sci. Technol. 3, 231–242 (2021).Article 
CAS 

Google Scholar 
Balch, W. E., Fox, G. E., Magrum, L. J., Woese, C. R. & Wolfe, R. S. Methanogens: reevaluation of a unique biological group. Microbiol. Rev. 43, 260–296 (1979).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Yokooji, Y. et al. Pantoate kinase and phosphopantothenate synthetase, two novel enzymes necessary for CoA biosynthesis in the Archaea. J. Biol. Chem. 284, 28137–28145 (2009).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Moench, T. & Zeikus, J. G. An improved preparation method for a titanium (III) media reductant. J. Microbiol. Methods 1, 199–202 (1983).CAS 
Article 

Google Scholar 
Miyashita, A. et al. Development of 16S rRNA gene-targeted primers for detection of archaeal anaerobic methanotrophs (ANMEs). FEMS Microbiol. Lett. 297, 31–37 (2009).CAS 
PubMed 
Article 

Google Scholar 
Sekiguchi, Y. et al. Sequence-specific cleavage of 16S rRNA for rapid and quantitative detection of particular groups of anaerobes in bioreactors. Water Sci. Technol. 52, 107–113 (2005).CAS 
PubMed 
Article 

Google Scholar 
Meechan, P. J. & Wilson, C. Use of ultraviolet lights in biological safety cabinets: a contrarian view. Appl. Biosaf. 11, 222–227 (2006).Article 

Google Scholar 
Imachi, H., Sekiguchi, Y., Kamagata, Y., Ohashi, A. & Harada, H. Cultivation and in situ detection of a thermophilic bacterium capable of oxidizing propionate in syntrophic association with hydrogenotrophic methanogens in a thermophilic methanogenic granular sludge. Appl. Environ. Microbiol. 66, 3608–3615 (2000).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Stams, A. J., Grolle, K. C., Frijters, C. T. & van Lier, J. B. Enrichment of thermophilic propionate-oxidizing bacteria in syntrophy with Methanobacterium thermoautotrophicum or Methanobacterium thermoformicicum. Appl. Environ. Microbiol. 58, 346–352 (1992).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Uchino, Y. & Suzuki, K. A simple preparation of liquid media for the cultivation of strict anaerobes. J. Pet. Environ. Biotechnol. S3, 001 (2011).
Google Scholar 
Akinyemi, T. S., Shao, N. & Whitman, W. B. Genus Methanothrix. In Bergey’s Manual of Systematics of Archaea and Bacteria (John Wiley & Sons, 2020). More

Mack, R. N. et al. Biotic invasions: Causes, epidemiology, global consequences, and control. Ecol. Appl. 10, 689–710 (2000).
Google Scholar 
Turbelin, A. J., Malamud, B. D. & Francis, R. A. Mapping the global state of invasive alien species: Patterns of invasion and policy responses. Glob. Ecol. Biogeogr. 26, 78–92 (2017).
Google Scholar 
Jackson, M. C. Interactions among multiple invasive animals. Ecology 96, 2035–2041 (2015).CAS 
PubMed 

Google Scholar 
Rodriguez, L. F. Can invasive species facilitate native species? Evidence of how, when, and why these impacts occur. Biol. Invasions 8, 927–939 (2006).
Google Scholar 
Duenas, M. A. et al. The role played by invasive species in interactions with endangered and threatened species in the United States: A systematic review. Biodivers. Conserv. 27, 3171–3183 (2018).
Google Scholar 
Weidenhamer, J. D. & Callaway, R. M. Direct and indirect effects of invasive plants on soil chemistry and ecosystem function. J. Chem. Ecol. 36, 59–69 (2010).CAS 
PubMed 

Google Scholar 
Bajwa, A. A., Chauhan, B. S., Farooq, M., Shabbir, A. & Adkins, S. W. What do we really know about alien plant invasion? A review of the invasion mechanism of one of the world’s worst weeds. Planta 244, 39–57 (2016).CAS 
PubMed 

Google Scholar 
Tallamy, D. W., Narango, D. L. & Mitchell, A. B. Do non-native plants contribute to insect declines?. Ecol. Entomol. 46, 729–742. https://doi.org/10.1111/een.12973 (2021).Article 

Google Scholar 
Bezemer, T. M., Harvey, J. A. & Cronin, J. T. Response of native insect communities to invasive plants. Annu. Rev. Entomol. 59, 119 (2014).CAS 
PubMed 

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

Google Scholar 
Kalisz, S., Kivlin, S. N. & Bialic-Murphy, L. Allelopathy is pervasive in invasive plants. Biol. Invasions 23, 367–371 (2021).
Google Scholar 
Pyšek, P. et al. A global assessment of invasive plant impacts on resident species, communities and ecosystems: The interaction of impact measures, invading species’ traits and environment. Glob. Change Biol. 18, 1725–1737 (2012).ADS 

Google Scholar 
Zhang, P., Li, B., Wu, J. & Hu, S. Invasive plants differentially affect soil biota through litter and rhizosphere pathways: A meta-analysis. Ecol. Lett. 22, 200–210 (2019).ADS 
PubMed 

Google Scholar 
Dudareva, N., Klempien, A., Muhlemann, J. K. & Kaplan, I. Biosynthesis, function and metabolic engineering of plant volatile organic compounds. New Phytol. 198, 16–32 (2013).CAS 
PubMed 

Google Scholar 
Clavijo McCormick, A. Can plant–natural enemy communication withstand disruption by biotic and abiotic factors?. Ecol. Evol. 6, 8569–8582 (2016).PubMed 
PubMed Central 

Google Scholar 
Bruce, T. J., Wadhams, L. J. & Woodcock, C. M. Insect host location: A volatile situation. Trends Plant Sci. 10, 269–274 (2005).CAS 
PubMed 

Google Scholar 
Clavijo McCormick, A., Unsicker, S. B. & Gershenzon, J. The specificity of herbivore-induced plant volatiles in attracting herbivore enemies. Trends Plant Sci. 17, 303–310 (2012).CAS 
PubMed 

Google Scholar 
Baldwin, I. T., Halitschke, R., Paschold, A., Von Dahl, C. C. & Preston, C. A. Volatile signaling in plant–plant interactions: “Talking trees” in the genomics era. Science 311, 812–815 (2006).ADS 
CAS 
PubMed 

Google Scholar 
Kegge, W. & Pierik, R. Biogenic volatile organic compounds and plant competition. Trends Plant Sci. 15, 126–132 (2010).CAS 
PubMed 

Google Scholar 
Effah, E., Holopainen, J. K. & Clavijo McCormick, A. Potential roles of volatile organic compounds in plant competition. Perspect. Plant Ecol. Evol. Syst. 38, 58–63 (2019).
Google Scholar 
Kigathi, R. N., Weisser, W. W., Reichelt, M., Gershenzon, J. & Unsicker, S. B. Plant volatile emission depends on the species composition of the neighboring plant community. BMC Plant Biol. 19, 1–17 (2019).
Google Scholar 
Karban, R., Wetzel, W. C., Shiojiri, K., Pezzola, E. & Blande, J. D. Geographic dialects in volatile communication between sagebrush individuals. Ecology 97, 2917–2924 (2016).PubMed 

Google Scholar 
Wheeler, G. S., David, A. S. & Lake, E. C. Volatile chemistry, not phylogeny, predicts host range of a biological control agent of Old-World climbing fern. Biol. Control 159, 104636 (2021).CAS 

Google Scholar 
Li, N. et al. Manipulating two olfactory cues causes a biological control beetle to shift to non-target plant species. J. Ecol. 105, 1534–1546 (2017).CAS 

Google Scholar 
Buddenhagen, C. E. Broom Control Monitoring at Tongariro National Park (Department of Conservation Wellington, 2000).
Google Scholar 
Hayes, L. et al. Biocontrol of Weeds: Achievements to Date and Future Outlook. Ecosystem services in New Zealand-conditions and trends Vol. 2, 375–385 (Manaaki Whenua Press, 2013).
Google Scholar 
Bagnall, A. Heather at Tongariro. A study of a weed introduction. Tussock Grasslands Mt. Lands Inst. Rev 41, 17–21 (1982).
Google Scholar 
Chapman, H. M. & Bannister, P. The spread of heather, Calluna vulgaris (L.) Hull, into indigenous plant communities of Tongariro National Park. N. Z. J. Ecol. 7–16 (1990).Effah, E. et al. Effects of two invasive weeds on arthropod community structure on the Central Plateau of New Zealand. Plants 9, 919 (2020).CAS 
PubMed Central 

Google Scholar 
Keesing, V. F. Impacts of invasion on community structure: habitat and invertebrate assemblage responses to Calluna vulgaris (L.) Hull invasion, in Tongariro National Park, New Zealand, Massey University Palmerston North, New Zealand, (1995).Peterson, P. G., Fowler, S. V. & Barrett, P. Is the poor establishment and performance of heather beetle in Tongariro National Park due to the impact of parasitoids predators or disease. N. Z. Plant Prot. 57, 89–93. https://doi.org/10.30843/nzpp.2004.57.6977 (2004).Article 

Google Scholar 
Ajpark. The brands and the bees: trade marks and the mānuka challenge for honey businesses, https://www.ajpark.com/insights/the-brands-and-the-bees-trade-marks-and-the-manuka-challenge-for-honey-businesses/#:~:text=M%C4%81nuka%20is%20a%20taonga%20species,may%20be%20offensive%20to%20M%C4%81ori (2021).Effah, E. et al. Seasonal and environmental variation in volatile emissions of the New Zealand native plant Leptospermum scoparium in weed-invaded and non-invaded sites. Sci. Rep. 10, 1–11 (2020).
Google Scholar 
Effah, E., Min Tun, K., Rangiwananga, N. & Clavijo McCormick, A. Mānuka clones differ in their volatile profiles: Potential implications for plant defence, pollinator attraction and bee products. Agronomy 12, 169 (2022).CAS 

Google Scholar 
Effah, E. et al. Natural variation in volatile emissions of the invasive weed Calluna vulgaris in New Zealand. Plants 9, 283 (2020).CAS 
PubMed Central 

Google Scholar 
Team, R. C. R: A language and environment for statistical computing (R Foundation for Statistical Computing, 2021).Ripley, B. et al. Package ‘mass’. Cran r 538, 113–120 (2013).
Google Scholar 
Chen, B. M., Liao, H. X., Chen, W. B., Wei, H. J. & Peng, S. L. Role of allelopathy in plant invasion and control of invasive plants. Allelopathy J 41, 155–166 (2017).
Google Scholar 
Ninkovic, V., Markovic, D. & Rensing, M. Plant volatiles as cues and signals in plant communication. Plant Cell Environ. 44, 1030–1043 (2021).CAS 
PubMed 

Google Scholar 
Holopainen, J. K. Multiple functions of inducible plant volatiles. Trends Plant Sci. 9, 529–533 (2004).CAS 
PubMed 

Google Scholar 
Rhoades, D. F. Responses of alder and willow to attack by tent caterpillars and webworms: evidence for pheromonal sensitivity of willows. In Plant Resistance to Insects (ed. Hedin, P. A.) 55–68 (American Chemical Society, 1983).
Google Scholar 
Hedin, P. A. Plant Resistance to Insects (American Chemical Society, 1983).
Google Scholar 
Heil, M. & Karban, R. Explaining evolution of plant communication by airborne signals. Trends Ecol. Evol. 25, 137–144 (2010).PubMed 

Google Scholar 
Barbosa, P. et al. Associational resistance and associational susceptibility: Having right or wrong neighbors. Annu. Rev. Ecol. Evol. Syst. 40, 1 (2009).
Google Scholar 
Kigathi, R. N., Weisser, W. W., Veit, D., Gershenzon, J. & Unsicker, S. B. Plants suppress their emission of volatiles when growing with conspecifics. J. Chem. Ecol. 39, 537–545 (2013).CAS 
PubMed 

Google Scholar 
Peñuelas, J. & Llusià, J. Influence of intra-and inter-specific interference on terpene emission by Pinus halepensis and Quercus ilex seedlings. Biol. Plant. 41, 139–143 (1998).
Google Scholar 
Ormeno, E., Fernandez, C. & Mévy, J.-P. Plant coexistence alters terpene emission and content of Mediterranean species. Phytochemistry 68, 840–852 (2007).CAS 
PubMed 

Google Scholar 
Himanen, S. J. et al. Birch (Betula spp.) leaves adsorb and re-release volatiles specific to neighbouring plants—A mechanism for associational herbivore resistance?. New Phytol. 186, 722–732 (2010).CAS 
PubMed 

Google Scholar 
Kessler, A. & Kalske, A. Plant secondary metabolite diversity and species interactions. Annu. Rev. Ecol. Evol. Syst. 49, 115–138 (2018).
Google Scholar 
Quintana-Rodriguez, E. et al. Plant volatiles cause direct, induced and associational resistance in common bean to the fungal pathogen Colletotrichum lindemuthianum. J. Ecol. 103, 250–260 (2015).CAS 

Google Scholar 
Loreto, F. & D’Auria, S. How do plants sense volatiles sent by other plants? Trends Plant Sci. (2021).Giordano, D., Facchiano, A., D’Auria, S. & Loreto, F. A hypothesis on the capacity of plant odorant-binding proteins to bind volatile isoprenoids based on in silico evidences. Elife 10, e66741 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Ninkovic, V., Markovic, D. & Dahlin, I. Decoding neighbour volatiles in preparation for future competition and implications for tritrophic interactions. Perspect. Plant Ecol. Evol. Syst. 23, 11–17 (2016).
Google Scholar 
Kegge, W. et al. Red: far-red light conditions affect the emission of volatile organic compounds from barley (Hordeum vulgare), leading to altered biomass allocation in neighbouring plants. Ann. Bot. 115, 961–970 (2015).CAS 
PubMed 
PubMed Central 

Google Scholar 
Gershenzon, J. Metabolic costs of terpenoid accumulation in higher plants. J. Chem. Ecol. 20, 1281–1328 (1994).CAS 
PubMed 

Google Scholar 
Anderson, P., Sadek, M., Larsson, M., Hansson, B. & Thöming, G. Larval host plant experience modulates both mate finding and oviposition choice in a moth. Anim. Behav. 85, 1169–1175 (2013).
Google Scholar 
Cunningham, J. P., Moore, C. J., Zalucki, M. P. & West, S. A. Learning, odour preference and flower foraging in moths. J. Exp. Biol. 207, 87–94 (2004).PubMed 

Google Scholar 
McCormick, A. C., Reinecke, A., Gershenzon, J. & Unsicker, S. B. Feeding experience affects the behavioral response of polyphagous gypsy moth caterpillars to herbivore-induced poplar volatiles. J. Chem. Ecol. 42, 382–393 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Proffit, M., Khallaf, M. A., Carrasco, D., Larsson, M. C. & Anderson, P. ‘Do you remember the first time?’ Host plant preference in a moth is modulated by experiences during larval feeding and adult mating. Ecol. Lett. 18, 365–374 (2015).PubMed 

Google Scholar 
Mayhew, P. J. Herbivore host choice and optimal bad motherhood. Trends Ecol. Evol. 16, 165–167 (2001).PubMed 

Google Scholar 
Jackson, T. et al. Anticipating the unexpected–managing pasture pest outbreaks after large-scale land conversion (New Zealand Grassland Association, 2012).Townsend, R. J., Dunbar, J. E. & Jackson, T. A. Flight behaviour of the manuka chafers, Pyronota festiva (Fabricius) and Pyronota setosa (Given) (Coleoptera: Melolonthinae), on the flipped soils of Cape Foulwind on the West Coast of New Zealand. N. Z. Plant Prot. 71, 255–261. https://doi.org/10.30843/nzpp.2018.71.175 (2018).Article 

Google Scholar 
Ferguson, C. M. et al. Quantifying the economic cost of invertebrate pests to New Zealand’s pastoral industry. N. Z. J. Agric. Res. 62, 255–315 (2019).
Google Scholar 
Cunningham, J. Can mechanism help explain insect host choice?. J. Evol. Biol. 25, 244–251 (2012).CAS 
PubMed 

Google Scholar 
Syrett, P., Smith, L. A., Bourner, T. C., Fowler, S. V. & Wilcox, A. A European pest to control a New Zealand weed: Investigating the safety of heather beetle, Lochmaea suturalis (Coleoptera: Chrysomelidae) for biological control of heather, Calluna vulgaris. Bull. Entomol. Res. 90, 169–178. https://doi.org/10.1017/S0007485300000286 (2000).CAS 
Article 
PubMed 

Google Scholar 
Fowler, S., Harman, H., Memmott, J., Peterson, P. & Smith, L. In Proceedings of the XII International Symposium on Biological Control of Weeds (eds Julien, M. H. et al.) 495–502.Fowler, S. V. et al. Investigating the poor performance of heather beetle, Lochmaea suturalis (Thompson) (Coleoptera: Chrysomelidae), as a weed biocontrol agent in New Zealand: Has genetic bottlenecking resulted in small body size and poor winter survival?. Biol. Control 87, 32–38 (2015).
Google Scholar 
Effah, E. et al. Herbivory and attenuated UV radiation affect volatile emissions of the invasive weed Calluna vulgaris. Molecules 25, 3200 (2020).CAS 
PubMed Central 

Google Scholar 
Pearson, D. E. & Callaway, R. M. Indirect nontarget effects of host-specific biological control agents: Implications for biological control. Biol. Control 35, 288–298 (2005).
Google Scholar 
Rand, T. A. & Louda, S. M. Exotic weed invasion increases the susceptibility of native plants to attack by a biocontrol herbivore. Ecology 85, 1548–1554. https://doi.org/10.1890/03-3067 (2004).Article 

Google Scholar  More

MaterialsWeinan city is located in the middle reaches of the Yellow River and in the southern part of the Loess Plateau (34°13’–35°52’N, 108°58’–110°35’E) (Fig. 1). It has a temperate semihumid, semiarid climate. The modern MAT observations indicate a value of 13.8 °C, and MAP is 570 mm; these values were obtained from the China meteorological data network, comprising the meteorological data of 2000–2015 (http://data.cma.cn/). Weinan has four distinct seasons, with hot and rainy conditions occurring in the same season. Much of the annual precipitation falls from June to August. The Weinan profile contains 42.8 m of loess–paleosol sequences (LPSs), including five paleosol layers from S0–S4 and five loess layers from L1–L5 and covering five glacial–interglacial cycles. The sampling method involved collecting one sample every 10 cm without interruption. A total of 427 samples were collected from this profile.Modern brGDGTs dataset and MAP datasetPreviously published brGDGTs data from surface soil samples were extracted using an established brGDGT-MAP model. The surface soil samples contain various types of soil and cover nearly all climatic and latitudinal zones. These datasets contain 712 surface soil samples, which all have separated 5-methyl and 6-methyl brGDGTs isomers (Table 1). To reduce the errors in collecting data from different laboratories, the MAP datasets we entered into the brGDGT-MAP model were all published in their previous studies, and we calculated the fractional abundances of each brGDGTs compound for each sample (Table 1), although there were no data regarding changes in soil occurring based on the brGDGTs indices among various laboratories. To eliminate and test the error of the previous MAP dataset, in this study, we also extracted each soil site’s multiyear MAP (1990–2020) through TerraClimate, which is a dataset of high-spatial-resolution monthly climate for global terrestrial surfaces (1/24°, ∼4 km)48. TerraClimate datasets reveal significant advances in the overall mean absolute error and enhance spatial realism compared with coarser resolution gridded datasets. Supplementary Fig. 3 shows that the two MAP datasets have high correlations, with only a few sites exhibiting large deviations. In this study, we entered these two MAP datasets into the DLNN model to obtain the most suitable DLNN-MAP model.Grain-size and magnetic susceptibility measurementsSamples at 10-cm intervals were dried in an oven at 40 °C for 3 days. Then, 0.2 g of each sample was weighed using a clean beaker with an electronic balance. Then, 10 ml of 30% H2O2 and 10 ml 10% HCl were added to remove organic matter and carbonate, respectively. Before the grain-size measurement, 0.05 mol/L (NaPO3)6 was added, and the solutions were placed in an ultrasonic machine for 10 min. The magnetic susceptibility of the samples were measured with an MS-2B Bartington meter. The grain-size was measured using a Mastersizer 2000 produced by Marvern Company in the UK, with an error of less than 1%.ChronologyWe used the ages of LPS control points on the Loess Plateau to obtain the age of each sample in the Weinan profile40. We used the magnetic susceptibility as an indicator of the accumulation rate39 combined with the U–230Th-dated oxygen isotope records from Sanbao caves in central China14. Each sample’s magnetic susceptibility was analyzed at 10-cm intervals (Supplementary Fig. 7). The calculation was as follows:$${T}_{{{{{{rm{m}}}}}}}={T}_{1}+frac{left({sum }_{i=1}^{m}{a}_{i}{s}_{i}right)left({T}_{2}-{T}_{1}right)}{{sum }_{i=1}^{n}{a}_{i}{s}_{i}}$$
(1)
where T1 and T2 indicate the ages of the control points, ai indicates the thickness of the layer, and si indicates the magnetic susceptibility of the layer.GDGTs analysisLipids in a total of 238 LPS samples were extracted, including the 198 samples reported in ref. 49. Forty samples at depths from 34.9 m to 43.7 m were selected every 20 cm intervals from the Weinan profile, and dried in an oven at 40 °C for 3 days. Afterward, the loess and paleosol samples were ground into powder and passed through a 60-mesh sieve. Each sample was weighed and extracted with 80 ml of methanol: dichloromethane (DCM) (1:9, v/v) using accelerated solvent extractors (ASE 100 or 150, Dionex, USA). The temperature and pressure were set at 100 °C and 1400 psi, respectively. Then, the lipid extracts were condensed in a rotary evaporator at 40 °C and separated into apolar and polar fractions on a flash silica gel column (0.7 cm i.d. and 1.5 g activated silica gel) chromatography using n-hexane and methanol as eluents, respectively. All polar components were passed through a 0.45-µm PTFE syringe filter. All apolar and polar compositions were dried under a gentle stream of nitrogen gas.The GDGTs were analyzed by using an Agilent 1200 series liquid chromatography-atmospheric pressure chemical ionization-6460A triple quadrupole mass spectrometry (LC-APCI-MS/MS). Ten microlitres of C46 GTGTs (0.001157 μg/μl) were added to each polar fraction, and the samples were then dissolved in 300 μl of n-hexane: iso-isopropanol (IPA) (98.2:1.8, v/v)). Two silica gel columns in series (150 mm × 2.1 mm, 1.9 μm, Thermo Finnigan; USA) were used for the separation of 5-methyl and 6-methyl brGDGTs, with the column temperature kept at 40 °C. The mass spectrometry settings were as follows: the vaporizer pressure 60 psi, the vaporizer temperature 400 °C, the flow rate of dry gas (N2) 6 l/min, drying gas temperature 200 °C, the capillary voltage 3500 V, the corona current 5 μA (∼3200 V), and a single-ion monitoring mode (SIM) was used50, targeting the protonated molecular ions ([M + H]+) 1304, 1302, 1300, 1298, 1296, 1292, 1050, 1048, 1046, 1036, 1034, 1032, 1020, 1018, and 744.The MATmr proxy was calculated to identify the changes that occurred in the mean annual temperature in the Weinan section over the last 430 ka. The calculation was as follows24,51.$${{MAT}}_{{mr}} =7.17+17.1*[{Ia}]+25.9*[{Ib}]+34.4*[{Ic}]-28.6*[{IIa}],(n=222,,{R}^{2} \ =0.68,; {RM}{SE}=4.6 {deg} {{{rm{C}}}},,P ; < ; 0.01)$$ (2) $${{MAT}}_{{mr}}=5.58+17.91*[{Ia}]-18.77*[{IIa}]$$ (3) $${MAT}({SSM})= 20.9-13.4*[{IIa}+{IIa}^{{prime}}]-17.2*[{IIIa}+{IIIa}^{{prime}}]\ -17.5*[{IIb}+{IIb}^{{prime}}]+11.2*[{Ib}]$$ (4) $${MAAT}=0.81-5.67*{CBT}+31.0*{MBT}^{{prime}}$$ (5) The soil pH was calculated using the following formulas24.$${pH}=7.15+1.59*{CBT}^{{prime}}(n=221,,{R}^{2}=0.85,,{RMSE}=0.52,, P , < ,0.0001)$$ (6) $${{CBT}}^{{prime} }=-{{log }}frac{{Ic}+{II}{a}^{{prime}}+{II}{b}^{{prime}}+{{IIc}}^{{prime} }+{{IIIa}}^{{prime} }+{III}{b}^{{prime} }+{{IIIc}}^{{prime} }}{{Ia}+{Ib}+{Ic}}$$ (7) SWC is well correlated with MBT’ when IR6ME  > 0.5, and these proxies were calculated using the following expressions:$${{MBT}}^{{prime} }=frac{({Ia}+{Ib}+{Ic})}{({Ia}+{Ib}+{Ic}+{IIa}+{{IIa}}^{{prime} }+{IIb}+{{IIb}}^{{prime} }+{IIc}+{{IIc}}^{{prime} }+{IIIa}+{IIIa}^{prime} )}$$
(8)
$${{IR}}_{6{ME}}=frac{sum (C6-{methylated; brGDGTs})}{sum {brGDGTs}}$$
(9)
$${{MBT}^{prime} }_{6{ME}}=frac{({Ia}+{Ib}+{Ic})}{({Ia}+{Ib}+{Ic}+{{IIa}}^{{prime} }+{{IIb}}^{{prime} }+{{IIc}}^{{prime} }+{IIIa}^{prime} )}$$
(10)
where the Roman numerals indicate different brGDGTs structures (Supplementary Fig. 1).Principal component analysis (PCA)CANOCO version 5 software was utilized to reveal the relationships among various environmental factors. The first PCA figure (Fig. 3a) was generated for the environmental factors MAT, MAPc, SWC, and pH. These variables are based on the same dataset (238 LSPs samples from Weinan profile) without any data transformation. The second PCA figure (Fig. 3b) was generated for the environmental factors MAT, MAP (based on 10Be), SWC and pH, which were all in the transition of the glacial–interglacial after 430 ka BP on the CLP. As the two LSPs profile had similar sedimentation rates, we obtained the same chronological control through linear interpolation in those transition periods. All datasets passed the Gaussian distribution test in this study.Cross wavelet analysisCompared with wavelet special analysis, cross wavelet analysis is even more complicated. The wavelet cross-spectrum can be defined as follows:$${CS}left(b,, a right)={m}_{1c}(b,, a){m}_{2c}(b,, a)$$
(11)
where ({m}_{1c}) and ({m}_{2c}) describe the covarying fractions of the overall spectra given by:$${m}_{1}left(b,, a right)={m}_{1c}left(b,, a right)+{m}_{1i}(b,, a)$$
(12)
$${m}_{2}left(b,, a right)={m}_{2c}left(b,, a right)+{m}_{2i}left(b,, a right)$$
(13)
where ({m}_{1i}) and ({m}_{2i}) are independent contributions to the variance.Overall, this is a multipart function that may be decomposed into amplitude and phase:$${CS}left(b,, a right)={{{{{rm{|}}}}}}{CS}left(b,, a right){{{{{rm{|}}}}}}{{exp }}(i;{{arg }}({CS}(b,, a)))$$
(14)
In this study, a and b represent the array of reconstructed MAPc and the Sanbao speleothem δ18O, respectively.Multiple regression linear modelTo compare the precision of the DLNN-MAP model we established, we set up a multiple regression linear model based on all 6-methyl brGDGTs except Ib. The basis of the model is defined as:$$y=a+{b}_{1}{x}_{1}+{b}_{2}{x}_{2}ldots+{b}_{n}{x}_{n}$$
(15)
where y represents MAP, x represents all 6-methyl brGDGTs and Ia and Ic, and a, b1, b2…bn represent the partial regression coefficients. n represents the number of 6-methyl we entered into the model (in this study, n = 8).The multiple correlation coefficient (R) was defined as follows:$$R=sqrt{frac{{sum }_{i=1}^{n}{({hat{y}}_{i}-bar{y})}^{2}}{{sum }_{i=1}^{n}{({y}_{i}-bar{y})}^{2}}}$$
(16)
where ({y}_{i}) represents the actual observed value, ({hat{y}}_{i}) represents the calculation value and (bar{y}) represents the mean value of all observed data.The t statistic is used to test the validity of regression coefficients, and it can be defined as follows:$${t}_{{b}_{j}}=frac{{b}_{j}}{{s}_{{b}_{j}}}$$
(17)
$${s}_{{b}_{j}}=sqrt{{p}_{{jj}}}*s$$
(18)
$$s=sqrt{frac{1}{n-m-1}mathop{sum }limits_{i=1}^{n}{({y}_{i}-{hat{y}}_{i})}^{2}}$$
(19)
$$P={({p}_{{jj}})}^{-1}=mathop{sum }limits_{i=1}^{n}({x}_{{ij}}-{bar{x}}_{j})({x}_{{ik}}-{bar{x}}_{k})$$
(20)
where ({b}_{j}) represents the regression coefficient of ({x}_{j}), n represents the number of samples and m represents the number of variables.The F statistic is used to test the linear relationship of variables and can be defined as follows:$$F=frac{1}{m{s}^{2}}mathop{sum }limits_{i=1}^{n}{({hat{y}}_{i}-bar{y})}^{2}$$
(21)
The variance inflation factor (VIF) is used to measure collinearity between variables and can be defined as follows:$${{VIF}}_{j}=frac{1}{1-{R}_{j}^{2}}$$
(22)
As shown in Supplementary Fig. 5, we found no obvious collinearity between different variables. However, there are fewer contributions in IIc’, IIIa’, IIIb’, and IIIc’ in the multiple regression linear model we established, and the value of t cannot attain the 95% confidence level (Table 2). The results of both the training dataset and extrapolated experimental dataset (Supplementary Fig. 6), although they seem good (R2 = 0.44 and 0.46, respectively), still have a considerable gap compared with the DLNN-MAP model. Especially when MAP  > 1500 mm, the multiple linear model is unable to forecast the real MAP. These results all indicate that the MAP influence on the brGDGTs compounds is not a simple linear relationship; instead, we suggest that there are complex pilot processes between them.Table 2 List of the parameters of the multiple linear modelFull size tableDLNN modelsDLNNs usually contain an input layer, a few hidden layers, and an output layer. A conceptual structure of the DLNN model was established for forecasting MAP values. The first layer accepts input signals that are various combinations of brGDGTs. The relationships among different variables are processed and analyzed in the hidden layers. The final class output is presented in the output layer; in this study, the output is the MAP reconstruction at the study site.The rectified linear unit (ReLU) activation function is applied in each neuron of the hidden layer, which is computationally simpler than the traditionally applied sigmoid function. The function of the ReLU activation function is given as follows:$$fleft(xright)=left{begin{array}{c}x,, x , > , 0 \ 0,, x , le , 0end{array}right.={{{{{rm{max }}}}}}(0,, x)$$
(23)
where x represents an input signal to a neuron and f represents the activation function.Furthermore, the bias between the measured and forecasted output values is reflected by the loss function. The loss function applied herein is the MAE (mean absolute error), which is given as follows:$${MAE}=frac{1}{N}mathop{sum }limits_{i=1}^{n}{{{{{rm{|}}}}}}T-Y{{{{{rm{|}}}}}}$$
(24)
where N is the number of training data points, and T and Y represent the measured output value and the forecasted class value, respectively.To realize the backpropagation framework, the derivative of the ReLU activation function needs to be acquired. According to the definition of the ReLU, the derivative is shown as follows.$$f{^prime} left(xright)=left{begin{array}{c}1,; x , > , 0 \ 0,; x , le , 0end{array}right.$$
(25)
Given a minibatch of m training samples obtained from the training set {x(1), x(2)…, x(m)} and their corresponding targets T(i) (i = 1,2…, m), the gradient used in the backpropagation framework is shown as follows:$$f=frac{1}{m}mathop{sum }limits_{i=1}^{n}frac{partial L}{partial w}$$
(26)
where L is the loss function; w represents the network weights; and n = 1 is the number of output values (MAP).In addition, considering that the adaptive moment estimation algorithm (Adam) was proven to be an effective neural network training method with a fast convergence speed and great classification performance, we applied this algorithm to train the DLNN model for MAP forecasting in this study. Adam has two biased equations, which are shown as follows:$$a={rho }_{1}a+(1-{rho }_{1})g$$
(27)
$$b={rho }_{2}b+(1-{rho }_{2})gtimes g$$
(28)
where ({rho }_{1}=0.9) and ({rho }_{2}=0.999) are exponential decay rates; g is the gradient; and (times) represents an elementwise product operator.After this calculation, the correct biases in the above two moments are given as follows:$${a}_{c}=frac{a}{1-{rho }_{1}^{t}}$$
(29)
$${b}_{c}=frac{b}{1-{rho }_{2}^{t}}$$
(30)
where t represents the current time step.Moreover, the update of the network weights is shown as follows:$${triangle }_{w}=-lambda frac{{a}_{c}}{sqrt{{b}_{c}}+epsilon }$$
(31)
where (lambda=0.001) represents the learning rates and (epsilon={10}^{-8}) is a constant used to ensure numerical stability.Eventually, the DLNN parameters can be updated according to the following formula.$$w ,=, w ,+, {triangle }_{w}$$
(32)
brGDGT-MAP modelsWe entered 9 brGDGTs compounds (all 6-methyl brGDGTs; each compound entered in the model is the percentage of all brGDGTs in the surface soil) into the input layer of the DLNN; these compounds are closely related to soil moisture. Then, we selected 533 surface soil samples as the training dataset and 179 surface soil samples as the validation dataset, both of which satisfied the principle of randomness. We assessed the precision of the model using forecast data R2 and root mean square error (RMSE) values.Through several parameters applied in the DLNN model, we found that the frequency of training and the number of neurons play the most significant roles in the brGDGT-MAP models. In addition, four hidden layers containing the other DLNN parameters allow the model to become more stable (detailed parameters are shown in Supplementary Fig. 7). To test the best frequency of training and the number of neurons in each hidden layer, we set a series of gradients to test the model to find the most suitable combination. As shown in Supplementary Fig. 8, for the frequency of training, we set the minimum and maximum training times to 1000 and 1500, respectively, with 100 times as the interval. We also set the numbers of neurons from 160 to 260 with a 20-neuron interval.Testing the weights of different compounds in the DLNN model and determining whether it was essential to eliminate some compounds that may make the dataset redundant were also required. Based on the model in which the Ib parameter was removed, we also set a series of experiments to test the effects of the different 6-methyl isomers on the predicted MAP values. Then, we made seven attempts to test the forecast effect of the brGDGT-MAP models by removing the Ic, IIa’, IIb’, IIc’, IIIa’, IIIb’, and IIIc’ parameters (Supplementary Fig. 9). Then, we obtained the best brGDGT-MAP model (Supplementary Fig. 10).Comparison of various ANN structuresTo improve the accuracy of our brGDGT-MAP models and the models’ universality, we also tested more complex ANN structures and then compared them with our DLNN models.RNNA recurrent neuron network (RNN) is an artificial neural network in which nodes are directionally connected into loops. The essential feature of RNN is that there are both internal feedback connections and feedforward connections between processing units. The inner structure of RNN is similar to that of the human brain, which can learn to transform a lifetime of sensory input streams into an efficient sequence of motor outputs (Supplementary Fig. 11a). Therefore, the basis of the RNN is defined as follows:$${h}_{t}=fleft(U ,*, {X}_{t}+W ,*, {h}_{t-1}right)$$
(33)
$${o}_{t}={softmax}(V ,{h}_{t})$$
(34)
where Xt represents the input value at time t; ot represents the output value at time t; ht represents the memory value at time t; and U, V, and W are the parameters of this network. For the motivative function, we chose softmax.LSTMLong short-term memory networks (LSTM) are a special type of RNN that can learn long-term dependence and contain three gates (forget gate, input gate and output gate) and one memory cell. The horizontal line above the box is called the cell state, and it acts as a conveyor belt to control the flow of information to the next moment (Supplementary Fig. 11b). Therefore, the basis of LSTM is defined as follows:$${C}_{t}={f}_{t}*{C}_{t-1}+{i}_{t}*{widetilde{C}}_{t}$$
(35)
where ({C}_{t-1}) represents the knowledge state of the model at time t − 1 and ({widetilde{C}}_{t}) represents the newly acquired information after entering new observations. ({f}_{t}) and ({i}_{t}) represent the weight parameters of ({C}_{t-1}) and ({widetilde{C}}_{t}), respectively.$${f}_{t}=sigma ({W}_{f}cdot left[{h}_{t-1},, {x}_{t}right]+{b}_{f})$$
(36)
$${i}_{t}=sigma ({W}_{f}cdot left[{h}_{t-1},, {x}_{t}right]+{b}_{i})$$
(37)
$$kern0.9pc {widetilde{C}}_{t}={{tanh }}({W}_{c}cdot left[{h}_{t-1},, {x}_{t}right]+{b}_{c})$$
(38)
where ({h}_{t-1}) represents the output value at time t − 1 and ({x}_{t}) represents the new input value at time t. ({W}_{f}) represents the motivative function in this study. We used tanh as the motivative function when our model was learning. Each new input may not have a positive impact on the machine, but it may also have a negative impact., ({b}_{f}), ({b}_{i}) and ({b}_{c}) represent the random disturbances (white noise).GRUAs mentioned above, the LSTM is proposed to overcome RNN’s inability to address remote dependence and the gate recurrent unit (GRU), a variant of the LSTM, keeps the effect of the LSTM while making the structure simpler.Compared with the LSTM, the GRU only has two gates (update (zt) and reset (rt) gates). The update gate is used to control the degree to which the state information at the previous moment is brought into the current state. The larger the value of the update gate is, the more state information at the previous moment is brought in. The reset gate is used to control the degree to which the state information at the previous moment is ignored (Supplementary Fig. 11c). Therefore, the basis of the LSTM is defined as follows:$${r}_{t}=sigma ({W}_{r}cdot [{h}_{t-1},, {x}_{t}])$$
(39)
$${z}_{t}=sigma ({W}_{z}cdot [{h}_{t-1},, {x}_{t}])$$
(40)
$${widetilde{h}}_{t}={tanh }({W}_{widetilde{h}}cdot [{{r}_{t}*h}_{t-1},, {x}_{t}])$$
(41)
$${h}_{t}=left(1-{z}_{t}right)*{{r}_{t}*h}_{t-1}+{z}_{t}*{widetilde{h}}_{t}$$
(42)
$${y}_{t}=sigma ({W}_{o}cdot {h}_{t})$$
(43)
where [] represents the connection of two vectors and * represents the multiplication of matrix elements. The xt and yt represent the input and output values at time t, respectively.It can be seen from the above formula that the parameters to be learned are the weight parameters of Wr, Wz, Wh, and Wo. The first three weights are spliced; therefore, they need to be segmented during learning. These can be defined as follows:$${W}_{r}={W}_{{rx}}+{W}_{{rh}}$$
(44)
$${W}_{z}={W}_{{zx}}+{W}_{{zh}}$$
(45)
$${W}_{widetilde{h}}={W}_{widetilde{h}x}+{W}_{widetilde{h}h}$$
(46)
As we can find in the RNN, LSTM, and GRU models we reconstructed (Supplementary Fig. 12), the training datasets all show extraordinarily high R2 values (0.99, 0.99, and 0.99, respectively) and low RMSE values (0.36, 0.23, and 0.16, respectively). However, the validation datasets do not show good prediction ability compared with the DLNN. These results indicate that the two ANN structures are not suitable for MAP prediction based on brGDGTs, although their inner structures are more complex than those of the DLNN. The reason we suggested is that the RNN, LSTM and GRU are more appropriate to the massive amounts of data and the data that have obvious spatiotemporal characteristics. The great prediction precision in the training dataset and the poor performance in the extrapolated datasets indicate that the models based on the RNN, LSTM and GRU have significant overfitting. As a result, compared with other ANN structures, we concluded that our DLNN model is the most suitable one to forecast MAP based on brGDGTs.Environmental indicators of n-alkanes proxiesLong-chain n-alkanes in plant leaf waxes are universal in terrestrial environments and can deliver signals of variations in plant sources and past climate. They are widely distributed in surface soils and Quaternary sediments, especially in LPSs. In this study, due to the insufficient samples in Weinan profile, we only analyzed n-alkanes components for 40 LPS samples, which contain ages between 340 and 430 ka BP.Instrumental measurementsFor the apolar fractions, a total of 40 samples in this study, mainly containing n-alkanes, were all investigated utilizing a Shimadzu 2010 gas chromatograph (GC) equipped with a flame ionization detector (FID) and a DB-5 fused silica capillary column (60 m (times) 0.32 mm (times) 0.25 μm film thickness) with helium as the carrier gas. The temperature of the GC oven was enhanced from 70 to 300 °C at a rate of 3 °C/min. Then, this temperature (300 °C) was maintained for 30 min. Finally, the concentrations of the n-alkane homologs were evaluated by assessing the peak area of the n-alkanes to that of the internal standard (cholane).Long-term paleoclimatic changeThe carbon preference index (CPI) evaluates the relative abundances of odd vs. even-numbered n-alkanes. The CPI increases as the environmental aridity increases. The CPI indicated warm–wet periods and cold-dry periods in paleoclimate and corresponded well with the loess–paleosol cycle52. The average chain length (ACL) value is the weighted average of the different carbon chain lengths. The lower ACL value corresponds to the lower temperature in the research of LPSs. The variations in the ACL value have good coordination with the magnetic susceptibility and particle size. The n-alkane CPI53 and ACL54 are calculated as follows:$${CPI}(1)=frac{({C}_{23}+{C}_{25}+{C}_{27}+{C}_{29}+{C}_{31})+({C}_{25}+{C}_{27}+{C}_{29}+{C}_{31}+{C}_{33})}{2({C}_{24}+{C}_{26}+{C}_{28}+{C}_{30}+{C}_{32})}$$
(47)
$${CPI}left(2right)=frac{1}{2}left(frac{{C}_{25}+{C}_{27}+{C}_{29}+{C}_{31}+{C}_{33}}{{C}_{24}+{C}_{26}+{C}_{28}+{C}_{30}+{C}_{32}}+frac{{C}_{25}+{C}_{27}+{C}_{29}+{C}_{31}+{C}_{33}}{{C}_{26}+{C}_{28}+{C}_{30}+{C}_{32}+{C}_{34}}right)$$
(48)
$${ACL}=frac{{23C}_{23}+{25C}_{25}+{27C}_{27}+{29C}_{29}+31{C}_{31}+{33C}_{33}}{{C}_{23}+{C}_{25}+{C}_{27}+{C}_{29}+{C}_{31}{+C}_{33}}$$
(49)
Figure 13 shows the variations in CPI (Supplementary Fig. 13a) and ACL (Supplementary Fig. 13b) values in the Weinan profile from 340 to 430 ka BP. Compared with the MAP (Supplementary Fig. 13c) and SWC (Fig. 2e) reconstructions based on brGDGTs, we found that they all had a peak at ∼350 ka BP, which indicates relatively high soil moisture at approximately 350 ka BP.MAP reconstruction in the XRD sectionIn this section, we test the brGDGT-MAP model in the Xiangride (XRD) profile, which is located in the margin region of the East Asian monsoon (Fig. 1). With robust chronological control, we reconstructed the rainfall changes in 7000 years BP (Supplementary Fig. 14b). We found that MAP was ∼200 mm in the late Holocene, which approaches multiple modern observations in this region (180 mm). Moreover, we suggest that this region experienced the most humid period in the mid-Holocene, when the rainfall reached 600 mm. Afterward, the precipitation declined from 6000 to 4000 years BP and then increased and reached a peak value at ∼3000 years BP. Then, it had a drought trend until modern times.We discovered that our brGDGT-MAP model could precisely capture rainfall dynamics based on the Weinan profile (Supplementary Fig. 14a) and XRD profile (Supplementary Fig. 14b). Combined with the most acceptable rainfall records in the Holocene (i.e., 10Be (Supplementary Fig. 14c), pollen in Gonghai (Fig. 1 and Supplementary Fig. 14d), and Dongge cave δ18O (Fig. 1 and Supplementary Fig. 14e)), we found the same precipitation peak values in the early Holocene and mid-Holocene. In addition, they all revealed a drought trend throughout the whole Holocene. We suggest that brGDGTs can become a robust proxy to reconstruct precipitation in the Holocene. More

.readcube-buybox { display: none !important;}
Chocolate, a serious contender for the world’s most beloved food, is made from the seed kernels of the cacao tree (Theobroma cacao). But despite its popularity, Justine Vansynghel at the University of Würzburg in Germany and her colleagues found that nobody had quantified how species living on small-scale cacao farms collectively affect production1.

Access options

/* style specs start */
style{display:none!important}.LiveAreaSection-193358632 *{align-content:stretch;align-items:stretch;align-self:auto;animation-delay:0s;animation-direction:normal;animation-duration:0s;animation-fill-mode:none;animation-iteration-count:1;animation-name:none;animation-play-state:running;animation-timing-function:ease;azimuth:center;backface-visibility:visible;background-attachment:scroll;background-blend-mode:normal;background-clip:borderBox;background-color:transparent;background-image:none;background-origin:paddingBox;background-position:0 0;background-repeat:repeat;background-size:auto auto;block-size:auto;border-block-end-color:currentcolor;border-block-end-style:none;border-block-end-width:medium;border-block-start-color:currentcolor;border-block-start-style:none;border-block-start-width:medium;border-bottom-color:currentcolor;border-bottom-left-radius:0;border-bottom-right-radius:0;border-bottom-style:none;border-bottom-width:medium;border-collapse:separate;border-image-outset:0s;border-image-repeat:stretch;border-image-slice:100%;border-image-source:none;border-image-width:1;border-inline-end-color:currentcolor;border-inline-end-style:none;border-inline-end-width:medium;border-inline-start-color:currentcolor;border-inline-start-style:none;border-inline-start-width:medium;border-left-color:currentcolor;border-left-style:none;border-left-width:medium;border-right-color:currentcolor;border-right-style:none;border-right-width:medium;border-spacing:0;border-top-color:currentcolor;border-top-left-radius:0;border-top-right-radius:0;border-top-style:none;border-top-width:medium;bottom:auto;box-decoration-break:slice;box-shadow:none;box-sizing:border-box;break-after:auto;break-before:auto;break-inside:auto;caption-side:top;caret-color:auto;clear:none;clip:auto;clip-path:none;color:initial;column-count:auto;column-fill:balance;column-gap:normal;column-rule-color:currentcolor;column-rule-style:none;column-rule-width:medium;column-span:none;column-width:auto;content:normal;counter-increment:none;counter-reset:none;cursor:auto;display:inline;empty-cells:show;filter:none;flex-basis:auto;flex-direction:row;flex-grow:0;flex-shrink:1;flex-wrap:nowrap;float:none;font-family:initial;font-feature-settings:normal;font-kerning:auto;font-language-override:normal;font-size:medium;font-size-adjust:none;font-stretch:normal;font-style:normal;font-synthesis:weight style;font-variant:normal;font-variant-alternates:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-position:normal;font-weight:400;grid-auto-columns:auto;grid-auto-flow:row;grid-auto-rows:auto;grid-column-end:auto;grid-column-gap:0;grid-column-start:auto;grid-row-end:auto;grid-row-gap:0;grid-row-start:auto;grid-template-areas:none;grid-template-columns:none;grid-template-rows:none;height:auto;hyphens:manual;image-orientation:0deg;image-rendering:auto;image-resolution:1dppx;ime-mode:auto;inline-size:auto;isolation:auto;justify-content:flexStart;left:auto;letter-spacing:normal;line-break:auto;line-height:normal;list-style-image:none;list-style-position:outside;list-style-type:disc;margin-block-end:0;margin-block-start:0;margin-bottom:0;margin-inline-end:0;margin-inline-start:0;margin-left:0;margin-right:0;margin-top:0;mask-clip:borderBox;mask-composite:add;mask-image:none;mask-mode:matchSource;mask-origin:borderBox;mask-position:0 0;mask-repeat:repeat;mask-size:auto;mask-type:luminance;max-height:none;max-width:none;min-block-size:0;min-height:0;min-inline-size:0;min-width:0;mix-blend-mode:normal;object-fit:fill;object-position:50% 50%;offset-block-end:auto;offset-block-start:auto;offset-inline-end:auto;offset-inline-start:auto;opacity:1;order:0;orphans:2;outline-color:initial;outline-offset:0;outline-style:none;outline-width:medium;overflow:visible;overflow-wrap:normal;overflow-x:visible;overflow-y:visible;padding-block-end:0;padding-block-start:0;padding-bottom:0;padding-inline-end:0;padding-inline-start:0;padding-left:0;padding-right:0;padding-top:0;page-break-after:auto;page-break-before:auto;page-break-inside:auto;perspective:none;perspective-origin:50% 50%;pointer-events:auto;position:static;quotes:initial;resize:none;right:auto;ruby-align:spaceAround;ruby-merge:separate;ruby-position:over;scroll-behavior:auto;scroll-snap-coordinate:none;scroll-snap-destination:0 0;scroll-snap-points-x:none;scroll-snap-points-y:none;scroll-snap-type:none;shape-image-threshold:0;shape-margin:0;shape-outside:none;tab-size:8;table-layout:auto;text-align:initial;text-align-last:auto;text-combine-upright:none;text-decoration-color:currentcolor;text-decoration-line:none;text-decoration-style:solid;text-emphasis-color:currentcolor;text-emphasis-position:over right;text-emphasis-style:none;text-indent:0;text-justify:auto;text-orientation:mixed;text-overflow:clip;text-rendering:auto;text-shadow:none;text-transform:none;text-underline-position:auto;top:auto;touch-action:auto;transform:none;transform-box:borderBox;transform-origin:50% 50%0;transform-style:flat;transition-delay:0s;transition-duration:0s;transition-property:all;transition-timing-function:ease;vertical-align:baseline;visibility:visible;white-space:normal;widows:2;width:auto;will-change:auto;word-break:normal;word-spacing:normal;word-wrap:normal;writing-mode:horizontalTb;z-index:auto;-webkit-appearance:none;-moz-appearance:none;-ms-appearance:none;appearance:none;margin:0}.LiveAreaSection-193358632{width:100%}.LiveAreaSection-193358632 .login-option-buybox{display:block;width:100%;font-size:17px;line-height:30px;color:#222;padding-top:30px;font-family:Harding,Palatino,serif}.LiveAreaSection-193358632 .additional-access-options{display:block;font-weight:700;font-size:17px;line-height:30px;color:#222;font-family:Harding,Palatino,serif}.LiveAreaSection-193358632 .additional-login >li:not(:first-child)::before{transform:translateY(-50%);content:””;height:1rem;position:absolute;top:50%;left:0;border-left:2px solid #999}.LiveAreaSection-193358632 .additional-login >li:not(:first-child){padding-left:10px}.LiveAreaSection-193358632 .additional-login >li{display:inline-block;position:relative;vertical-align:middle;padding-right:10px}.BuyBoxSection-683559780{display:flex;flex-wrap:wrap;flex:1;flex-direction:row-reverse;margin:-30px -15px 0}.BuyBoxSection-683559780 .box-inner{width:100%;height:100%}.BuyBoxSection-683559780 .readcube-buybox{background-color:#f3f3f3;flex-shrink:1;flex-grow:1;flex-basis:255px;background-clip:content-box;padding:0 15px;margin-top:30px}.BuyBoxSection-683559780 .subscribe-buybox{background-color:#f3f3f3;flex-shrink:1;flex-grow:4;flex-basis:300px;background-clip:content-box;padding:0 15px;margin-top:30px}.BuyBoxSection-683559780 .subscribe-buybox-nature-plus{background-color:#f3f3f3;flex-shrink:1;flex-grow:4;flex-basis:100%;background-clip:content-box;padding:0 15px;margin-top:30px}.BuyBoxSection-683559780 .title-readcube{display:block;margin:0;margin-right:20%;margin-left:20%;font-size:24px;line-height:32px;color:#222;padding-top:30px;text-align:center;font-family:Harding,Palatino,serif}.BuyBoxSection-683559780 .title-buybox{display:block;margin:0;margin-right:29%;margin-left:29%;font-size:24px;line-height:32px;color:#222;padding-top:30px;text-align:center;font-family:Harding,Palatino,serif}.BuyBoxSection-683559780 .title-asia-buybox{display:block;margin:0;margin-right:5%;margin-left:5%;font-size:24px;line-height:32px;color:#222;padding-top:30px;text-align:center;font-family:Harding,Palatino,serif}.BuyBoxSection-683559780 .asia-link{color:#069;cursor:pointer;text-decoration:none;font-size:1.05em;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:1.05em6}.BuyBoxSection-683559780 .access-readcube{display:block;margin:0;margin-right:10%;margin-left:10%;font-size:14px;color:#222;padding-top:10px;text-align:center;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:20px}.BuyBoxSection-683559780 .access-asia-buybox{display:block;margin:0;margin-right:5%;margin-left:5%;font-size:14px;color:#222;padding-top:10px;text-align:center;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:20px}.BuyBoxSection-683559780 .access-buybox{display:block;margin:0;margin-right:30%;margin-left:30%;font-size:14px;color:#222;opacity:.8px;padding-top:10px;text-align:center;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:20px}.BuyBoxSection-683559780 .usps-buybox{display:block;margin:0;margin-right:30%;margin-left:30%;font-size:14px;color:#222;opacity:.8px;text-align:center;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:20px}.BuyBoxSection-683559780 .price-buybox{display:block;font-size:30px;color:#222;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;padding-top:30px;text-align:center}.BuyBoxSection-683559780 .price-from{font-size:14px;padding-right:10px;color:#222;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:20px}.BuyBoxSection-683559780 .issue-buybox{display:block;font-size:13px;text-align:center;color:#222;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:19px}.BuyBoxSection-683559780 .no-price-buybox{display:block;font-size:13px;line-height:18px;text-align:center;padding-right:10%;padding-left:10%;padding-bottom:20px;padding-top:30px;color:#222;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif}.BuyBoxSection-683559780 .vat-buybox{display:block;margin-top:5px;margin-right:20%;margin-left:20%;font-size:11px;color:#222;padding-top:10px;padding-bottom:15px;text-align:center;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;line-height:17px}.BuyBoxSection-683559780 .button-container{display:flex;padding-right:20px;padding-left:20px;justify-content:center}.BuyBoxSection-683559780 .button-container >*{flex:1px}.BuyBoxSection-683559780 .button-container >a:hover,.Button-505204839:hover,.Button-1078489254:hover,.Button-2808614501:hover{text-decoration:none}.BuyBoxSection-683559780 .readcube-button{background:#fff;margin-top:30px}.BuyBoxSection-683559780 .button-asia{background:#069;border:1px solid #069;border-radius:0;cursor:pointer;display:block;padding:9px;outline:0;text-align:center;text-decoration:none;min-width:80px;margin-top:75px}.BuyBoxSection-683559780 .button-label-asia,.ButtonLabel-3869432492,.ButtonLabel-3296148077,.ButtonLabel-1566022830{display:block;color:#fff;font-size:17px;line-height:20px;font-family:-apple-system,BlinkMacSystemFont,”Segoe UI”,Roboto,Oxygen-Sans,Ubuntu,Cantarell,”Helvetica Neue”,sans-serif;text-align:center;text-decoration:none;cursor:pointer}.Button-505204839,.Button-1078489254,.Button-2808614501{background:#069;border:1px solid #069;border-radius:0;cursor:pointer;display:block;padding:9px;outline:0;text-align:center;text-decoration:none;min-width:80px;max-width:320px;margin-top:10px}.Button-505204839 .readcube-label,.Button-1078489254 .readcube-label,.Button-2808614501 .readcube-label{color:#069}
/* style specs end */Subscribe to Nature+Get immediate online access to the entire Nature family of 50+ journals$29.99monthlySubscribe to JournalGet full journal access for 1 year$199.00only $3.90 per issueAll prices are NET prices.VAT will be added later in the checkout.Tax calculation will be finalised during checkout.Buy articleGet time limited or full article access on ReadCube.$32.00All prices are NET prices.

Additional access options:

doi: https://doi.org/10.1038/d41586-022-02908-0

References

Subjects

Latest on: More

Cliffs present a unique flora that has been little studied until now mainly because of the inherent difficulties to access this unique environment, as shown in Fig. 2. The techniques currently used to access plants on steep slopes and cliffs (e.g., abseiling, helicopter) are generally dangerous, costly and time consuming. Using a small aerial manipulator to sample plants on the cliffs can represent many advantages, including safety and portability, as well as the capability of reaching otherwise inaccessible locations easily, quickly and at low cost.Figure 2Examples of the cliff habitats of some critically endangered species on the Kauaʻi Island along with the count of known individuals as of February 2022.Full size imageHowever, several technical challenges make it difficult to develop suitable aerial manipulators for this task. Indeed, the sampling of plants on cliffs necessarily leads to significant collision risks, as well as contact forces and moments during sampling that can destabilize the drone. The samples collected would also need to be accessed from the side of the aerial platform22. Any weight (e.g., sampling tool, collected samples) located horizontally away from the center of mass of the drone creates large additional demands on the propulsion system of most drones. To collect specific plant parts in windy conditions (e.g., scion, flowers, seeds, etc.), precise and fast motion is required even in degraded Global Navigation Satellite System (GNSS) coverage near the cliffs. The great diversity of plant species and morphology found on cliffs, as well as the variety of targeted sections of plant, also represent a major design challenge. Finally, to maximize the adoption of this tool, it is also desirable that scientists with minimal training could use this platform. The next sections describe how these challenges were addressed through the development of the Mamba.Suspended sampling platformThere are a multitude of configurations that could have been explored to sample plants on cliffs. Some drones have manipulators rigidly attached to their structure20,23. However, these manipulators tend to have a limited reach to keep the center of mass within the propeller footprint and to minimize the inertia of the system. This could result in a high collision risk with the propellers in the uneven terrain found on cliffs. The contact forces created during the sampling operation also generate destabilizing moments through manipulators rigidly attached to the drone. To address these challenges, concepts involving a compliant manipulator operated from specialized drones were also explored10. Alternatively, some aerial manipulators were also passively suspended under the drone through a long rod21,24. This keeps the drone above potential obstacles within the environment, significantly reducing the operator’s mental demand and stress while also reducing the disturbances transmitted to the drone to a downward force aligned with the rod and yaw torque. To maintain these advantages while providing better precision, some projects have developed cable suspended platforms equipped with thrusters25,26. As these platforms do not have to counter gravity, the thrusters can be positioned to fight external disturbances more efficiently (e.g., wind, contact forces, drone movements). Existing systems however only stabilize the suspended platform close to its equilibrium point.The chosen concept for the Mamba, illustrated at Fig. 3, consists of a suspended platform that can stabilize itself far from its natural equilibrium to provide a large workspace. The lifting drone in this system stays safely away and above from steep cliff faces, while supporting the platform and providing rough positioning in space through better GNSS coverage. The platform is suspended 10 m below the lifting drone using four attachment points to prevent pitch and roll motions. The cable also acts as a low pass filter, isolating the platform from the fast drone movements required to fight wind disturbances. The suspended platform design can then focus on fast and precise positioning, while also being tolerant to contacts during sampling. To do so, four pairs of bidirectional actuators are used to control the motion in the plane of the pendulum (i.e., x and y translation, as well as yaw). Two pairs of actuators are installed in the x-direction to provide sufficient force to reach plants as far as 4 m from the equilibrium position. This corresponds to roughly 3.3 m from the tip of the lifting drone’s propellers.Figure 3(a) General concept of the Mamba and lifting drone during transit and sampling on cliffs. (b) Side view of the Mamba showing the components and cable installations. (c) Top view showing the antagonist thrusters configuration. (d) Close-up of the sampling tool and 2 degrees of freedom (DOF) wrist specifically designed to sample small fragile plants.Full size imageSince the Mamba is self-powered and has its own communication system, the lifting drone function is simply to lift the platform and hold it in place. This made it possible to select amongst the many commercially available products to accelerate the development of the Mamba. The DJI M300 was chosen as it comes equipped with a 360° optical obstacle avoidance vision system, an IP45 rating, and a flight time of 20 min with the Mamba attached (3.3 kg). It also advertised a four constellation GNSS receiver for better coverage around buildings, structures, and cliffs.Precise control in windsWinds under 20 km/h represent a gentle breeze on the Beaufort scale. At this level, the wind only moves the leaves, and not the branches, which allows for ideal sampling conditions. According to historical weather data from 2020, daily maximum winds are less than 20 km/h for 40 to 70% of the year, depending on the exact location on Kauaʻi Island (i.e., Lihuʻe International airport, as reported by the National Oceanic and Atmospheric Administration, and the Makaha Ridge Weather Station, as reported in the MesoWest database). This also implies that Kauaʻi experiences stronger winds on certain days which would make precise sampling difficult. Wind conditions are also more challenging near cliff faces, with increased turbulence and vertical airflow along the cliff.To allow operations on most days, while providing precise positioning and fast rejection of wind disturbances, the actuators of the Mamba are oriented in the horizontal plane. This allows the actuator forces to directly affect the motion of the suspended platform. Each actuator of the Mamba consists of a pair of brushless DC motors and 23 cm propellers capable of producing 7 N of force. The motors are installed in opposite directions, are always idling at their minimum rotation speed, and are commanded to only create force in their preferred direction. This antagonistic configuration avoids the low-velocity dead zone of a brushless motor during thrust reversal. This makes it possible to quickly revert the direction of the thrust and nearly triples the bandwidth of the actuators to approximately 2.5 Hz27. This configuration, however, comes at the expense of added mass and components.The Mamba is equipped with a flight controller that includes a control system, and a state estimator. To avoid degraded GNSS coverage issues, the state estimator only uses data from a high accuracy inertial measurement unit (IMU) to estimate the attitude of the platform. This provides the relative position of the platform with respect to the drone and is sufficient for teleoperation. Three separated proportional-derivative controllers are used for each of the DOF controlled by the actuators. This control system also provides attitude-hold assistance (i.e., pitch and roll, which correspond to x and y displacements, as well as yaw). This implies that if the user does not send any commands, the suspended platform maintains its current state.Figure 4 illustrates the stabilization accuracy of the Mamba when moving along a representative trajectory when suspended indoors from a 5.7 m cable (limited by ceiling height). This experiment confirmed that the sampling tool can maintain a position at a horizontal reach of 2.25 m with a precision of about 5 cm for 30 s. As the horizontal reach and precision are limited by the cable angular displacements (e.g., component of weight acting on the pendulum, IMU angular resolution), the resulting workspace when operating with a 10 m long cable would reach a radius of 4 m with a positioning accuracy of about 9 cm. To account for potential external disturbances like wind, the sampling tool was designed with an opening of 15 cm. This creates some margin for the pilot to align the target with the sampling mechanism. Field trials detailed below demonstrated that the Mamba actuators and controller could maintain a sufficiently stable position to sample plants in winds During the sampling phase, wind speed averaged 15.7 km/h with a standard deviation of 6.8 km/h, while wind gusts reached an average of 20.1 km/h with a standard deviation of 6.5 km/h. The maximum average wind speed recorded during sampling was 28 km/h with gusts up to 37 km/h. This represents a lower bound of the system performance, as no failure resulted from the wind conditions experienced during the trials. The a ttached Supplementary Video also demonstrates the stability of the system.Figure 4Representative motion of the sampling tool within its workspace based only on feedback from a high accuracy IMU and recorded using a motion capture system. The natural equilibrium point is at (0,0). The experiment starts with a 90° rotation around the z axis, followed by a forward movement along the x-axis of the Mamba and a lateral movement along its y-axis. The system then maintains this position for 30 s without any user inputs. Produced in MATLAB R2021a.Full size imageTeleoperated sampling of cliffs habitatsPlants growing on Kauaʻi cliffs exhibit a wide morphological variety. For this project, targets ranged from small herbaceous plants such as Euphorbia eleanoriae (plants  More

Incentivizing farmers to shift from conventional to regenerative practices could help fulfil the United Nations Food Systems commitments to transform food supply chains — as well as reducing carbon emissions (see L. A. Schulte et al. Nature Sustain. 5, 384–388; 2022).
Competing Interests
The authors declare no competing interests. More

BirdLife International. Red List Update: Parrots of the Americas in Peril. https://www.birdlife.org/news/2021/02/08/red-list-update-parrots-of-the-americas-in-peril/ (2020).Berkunsky, I. et al. Current threats faced by Neotropical parrot populations. Biol. Cons. 214, 278–287. https://doi.org/10.1016/j.biocon.2017.08.016 (2017).Article 

Google Scholar 
ICMBIO—Instituto Chico Mendes de Conservação da Biodiversidade (Org.). Livro Vermelho da Fauna Brasileira Ameaçada de Extinção: Volume III-Aves 709. https://www.icmbio.gov.br/portal/images/stories/comunicacao/publicacoes/publicacoes-diversas/livro_vermelho_2018_vol3.pdf (Ministério do Meio Ambiente, 2018).CBRO—Comitê Brasileiro de Registros Ornitológicos. Listas das Aves do Brasil. 11th ed. http://www.cbro.org.br/wp-content/uploads/2020/06/avesbrasil_2014jan1.pdf (CBRO, 2014).Pacheco, J. F. et al. Annotated checklist of the birds of Brazil by the Brazilian Ornithological Records Committee—second edition. Ornithol. Res. 29(2), 94–105. https://doi.org/10.1007/s43388-021-00058-x (2021).Article 

Google Scholar 
IUCN—International Union for Conservation of Nature. The IUCN Red List of Threatened Species www.iucnredlist.org (2018).Guedes, N. M. R. Biologia reprodutiva da arara azul (Anodorhynchus hyacinthinus) no Pantanal—MS, Brasil. (Dissertação de Mestrado Universidade de São Paulo, São Paulo (1993).Guedes, N. M. R. et al. Technical Report Assessing the Impact of Fire on Blue Macaws, Pantanal, Mato Grosso do Sul, Brazil, p 13, Campo Grande, Instituto Arara Azul (2019).Guedes, N. M. R. Araras azuis: 15 anos de estudos no Pantanal. In Paper presented at IV Simpósio Sobre Recursos Naturais e Sócio-Econômicos do Pantanal, Corumbá: Embrapa Pantanal (2004).Guedes, N. M. R. Sucesso reprodutivo, mortalidade e crescimento de filhotes de araras azuis Anodorhynchus hyacinthinus (Aves, Psittacidae), no Pantanal, Brasil (Tese de doutorado Universidade Estadual Paulista, Botucatu, 2009)Guedes, N. M. R. & Harper, L. H. Hyacinth macaws in the Pantanal. In The Large Macaws (eds Abramson, J. et al.) 394–421 (Raintree Publications, 1995).
Google Scholar 
Vicente, E. C. & Guedes, N. M. Organophosphate poisoning of Hyacinth Macaws in the Southern Pantanal, Brazil. Sci. Rep. 11, 1–6. https://doi.org/10.1038/s41598-021-84228-3 (2021).CAS 
Article 

Google Scholar 
Guedes, N. M. R. et al. Assessment of fire impact on Hyacinth Macaws in Perigara, Pantanal—MT, Brazil, p 35, Campo Grande, Instituto Arara Azul (2020).Guedes, N. M. R. et al. Macaws survive fires and provide hope for resilience—Stubborn survivors. Pantanal Sci. Mag. 6, 36–41 (2021).
Google Scholar 
Oliveira, M. D. R. et al. Lack of protected areas and future habitat loss threaten the Hyacinth Macaw Anodorhynchus hyacinthinus and its main food and nesting resources. Ibis 163, 1217–1234 (2021).Article 

Google Scholar 
Ricklefs, R. E. Patterns of growth in birds. Ibis 110, 419–451. https://doi.org/10.1111/j.1474-919X.1968.tb00058.x (1968).Article 

Google Scholar 
Gebhardt-Henrich, S. & Richner, H. Causes of growth variation and its consequences for fitness. Oxford Ornithol. Ser. 8, 324–339 (1998).
Google Scholar 
Masello, J. F. & Quillfeldt, P. Body size, body condition and ornamental feathers of Burrowing Parrots: Variation between years and sexes, assortative mating and influences on breeding success. Emu Austral Ornithol. 103, 149–161. https://doi.org/10.1071/MU02036 (2003).Article 

Google Scholar 
Renton, K. Influence of environmental variability on the growth of Lilac-crowned Parrot nestlings. Ibis 144, 331–339. https://doi.org/10.1046/j.1474-919X.2002.00015.x (2002).Article 

Google Scholar 
Masello, J. F. & Quillfeldt, P. Chick growth and breeding success of the Burrowing Parrot. Condor 104, 574–586. https://doi.org/10.1650/0010-5422 (2002).Article 

Google Scholar 
Pacheco, M. A., Beissinger, S. R. & Bosque, C. Why grow slowly in a dangerous place? Postnatal growth, thermoregulation, and energetics of nestling green-rumped parrotlets (Forpus passerinus). Auk 127, 558–570. https://doi.org/10.1525/auk.2009.09190 (2010).Article 

Google Scholar 
Vigo, G., Williams, M. & Brightsmith, D. J. Growth of Scarlet Macaw (Ara macao) chicks in southeastern Peru. Neotrop. Ornithol. 22, 143–153 (2011).
Google Scholar 
Lyon, J. P. et al. Reintroduction success of threatened Australian trout cod (Maccullochella macquariensis) based on growth and reproduction. Mar. Freshw. Res. 63, 598–605. https://doi.org/10.1071/MF12034 (2012).Article 

Google Scholar 
Vigo-Trauco, G., Garcia-Anleu, R. & Brightsmith, D. J. Increasing survival of wild macaw chicks using foster parents and supplemental feeding. Diversity 13, 121. https://doi.org/10.3390/d13030121 (2021).Article 

Google Scholar 
Tellería, J. L., De La Hera, I. & Perez-Tris, J. Morphological variation as a tool for monitoring bird populations: A review. Ardeola 60, 191–224. https://doi.org/10.13157/arla.60.2.2013.191 (2013).Article 

Google Scholar 
Silva, J. S. V. Elementos fisiográficos para delimitação do ecossistema Pantanal: Discussão e proposta. Oecol. Brasil. 1, 349–458. https://doi.org/10.4257/OECO.1995.0101.22 (1995).Article 

Google Scholar 
Silva, J. S. V. & Abdon, M. M. Delimitação do Pantanal Brasileiro e suas Sub-Regiões. Pesq. Agropec. Bras. 33, 1703–1711 (1998).
Google Scholar 
Keuroghlian, A., Eaton, D. & Desbiez, A. L. J. The response of a landscape species, white-lipped peccaries, to seasonal resource fluctuations in a tropical wetland, the Brazilian Pantanal. Int. J. Biodivers. Conserv. 1, 87–97 (2009).
Google Scholar 
Donatelli, R. J., Posso, S. R. & Toledo, M. C. B. D. Distribution, composition and seasonality of aquatic birds in the Nhecolândia sub-region of South Pantanal, Brazil. Braz. J. Biol. 74, 844–853 (2014).CAS 
Article 

Google Scholar 
Donatelli, R. J. et al. Temporal and spatial variation of richness and abundance of the community of birds in the Pantanal wetlands of Nhecolândia (Mato Grosso do Sul, Brazil). Rev. Biol. Trop. 65, 1358–1380 (2017).Article 

Google Scholar 
Tomas, W. M. et al. Sustainability agenda for the Pantanal Wetland: Perspectives on a collaborative interface for science, policy, and decision-making. Trop. Conserv. Sci. 12, 1–30. https://doi.org/10.1177/1940082919872634 (2019).ADS 
Article 

Google Scholar 
Harris, M. B. et al. Safeguarding the Pantanal wetlands: Threats and conservation initiatives. Conserv. Biol. 19, 714–720. https://doi.org/10.1111/j.1523-1739.2005.00708.x (2005).Article 

Google Scholar 
Santos Júnior, A. D., Aspectos populacionais de Sterculia apetala (Jacq.) Karst (Sterculiaceae) como subsídios ao plano de conservação da arara-azul no Sul do Pantanal, Mato Grosso do Sul, Brasil. (2006). https://repositorio.ufms.br/handle/123456789/521.Ricklefs, R. E. The optimization of growth rate in altricial birds. Ecology 65, 1602–1616 (1984).Article 

Google Scholar 
Bruford, M. W., Hanotte, O., Brookfield, J. F. Y. & Burke, T. Single-locus and multilocus DNA fingerprinting. In Molecular Genetic Analysis of Populations: A Practical Approach (ed. Hoelzel, A. R.) 225–269 (Oxford University Press, 1992).
Google Scholar 
Miyaki, C. Y. et al. Sex identification of parrots, toucans, and curassows by PCR: Perspectives for wild and captive population studies. Zoo Biol. 17(5), 415–423 (1998).Article 

Google Scholar 
Cavanaugh, J. E. & Neath, A. A. The Akaike information criterion: Background, derivation, properties, application, interpretation, and refinements. Wiley Interdiscip. Rev. Comput. Stat. 11, 1460. https://doi.org/10.1002/wics.1460 (2019).MathSciNet 
Article 

Google Scholar 
Motulsky H. J. GraphPad curve fitting guide. 2021. http://www.graphpad.com/guides/prism/7/curve-fitting/index.htm. Accessed 18 September.Saunders, D. A., Smith, G. T. & Rowley, I. The availability and dimensions of tree hollows that provide nest sites for cockatoos (Psittaciformes) in Western Australia. Wildl. Res. 9, 541–556. https://doi.org/10.1071/WR9820541 (1982).Article 

Google Scholar 
Navarro, J. L. & Bucher, E. H. Growth of monk parakeets. Wilson Bull. 102, 520–525 (1990).
Google Scholar 
Murtaugh, P. A. Performance of several variable-selection methods applied to real ecological data. Ecol. Lett. 12, 1061–1068 (2009).Article 

Google Scholar 
Waltman, J. R. & Beissinger, S. R. Breeding behavior of the Green-rumped Parrotlet. Wilson Bull. 104, 65–84 (1992).
Google Scholar 
Enkerlin-Hoeflich, E. C., Packard, J. M. & González-Elizondo, J. J. Safe field techniques for nest inspections and nestling crop sampling of parrots. J. Field Ornithol. 70, 8–17 (1999).
Google Scholar 
Barros, Y. de M. Biologia comportamental de Propyrrhura maracana (Aves, Psittacidae): Fundamentos para conservação in situ de Cyanopsitta spixii (Aves, Psittacidae) na Caatinga. (Tese de Doutorado Universidade Estadual de São Paulo, Rio Claro, 2001).Seixas, G. H. F. & Mourão, G. M. Growth of nestlings of the BlueFronted Amazon (Amazona aestiva) raised in the wild or in captivity. Ornitol. Neotrop. 14, 295–305 (2003).
Google Scholar 
Vigo-Trauco, G. Crecimiento de pichones de Guacamayo Escarlata, Ara macao (Linneus: 1758) en la Reserva Nacional Tambopata-Madre de Dios-Peru (Tese Universidad Nacional Agraria La Molina, 2007).
Google Scholar 
Tjørve, K. M. & Tjørve, E. The use of Gompertz models in growth analyses, and new Gompertz-model approach: An addition to the Unified-Richards family. PLoS One https://doi.org/10.1371/journal.pone.0178691 (2017).Article 
PubMed 
PubMed Central 
MATH 

Google Scholar 
Reed, J. M. The role of behavior in recent avian extinctions and endangerments. Conserv. Biol. 13, 232–241. https://doi.org/10.1046/j.1523-1739.1999.013002232.x (1999).Article 

Google Scholar 
Tjørve, K. M., Underhill, L. G. & Visser, G. H. Energetics of growth in semi-precocial shorebird chicks in a warm environment: The African black oystercatcher, Haematopus moquini. Zoology 110, 176–188. https://doi.org/10.1016/j.zool.2007.01.002 (2007).Article 
PubMed 

Google Scholar 
Tjørve, K. M., Underhill, L. G. & Visser, G. H. The energetic implications of precocial development for three shorebird species breeding in a warm environment. Ibis 150, 125–138 (2008).Article 

Google Scholar 
Ricklefs, R. E. Weight recession in nestling birds. Auk 85, 30–35. https://doi.org/10.2307/4083621 (1968).Article 

Google Scholar 
Huin, N. & Prince, P. A. Chick growth in albatrosses: Curve fitting with a twist. J. Avian Biol. 31, 418–425. https://doi.org/10.1034/j.1600-048X.2000.310318.x (2000).Article 

Google Scholar 
Corsini, M. et al. Growing in the city: Urban evolutionary ecology of avian growth rates. Evol. Appl. 14, 69–84. https://doi.org/10.1111/eva.13081 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Barbosa, L. T. Avaliação do sucesso reprodutivo da arara-canindé (Ara ararauna—Psittacidae) e o desenvolvimento urbano de Campo Grande, Mato Grosso do Sul (Dissertação de mestrado Universidade Anhanguera Uniderp, Campo Grande, 2015).Giraldo-Deck, L. M. et al. Development of intraspecific size variation in black coucals, white-browed coucals and ruffs from hatching to fledging. J. Avian Biol. 51, e02440. https://doi.org/10.1111/jav.02440 (2020).Article 

Google Scholar 
Guedes et al. Annual Technical Report from the Instituto Arara Azul., Pantanal-MS, Brazil. 35p, Campo Grande, Instituto Arara Azul (2022). More

Cole, J., Findlay, S. & Pace, M. Bacterial production in fresh and saltwater ecosystems: a cross-system overview. Mar. Ecol. Prog. Ser. 43, 1–10 (1988).ADS 
Article 

Google Scholar 
Azam, F. et al. The Ecological Role of Water-Column Microbes in the Sea. Mar. Ecol. Prog. Ser. 10, 257–263 (1983).ADS 
Article 

Google Scholar 
Cotner, J. B. & Biddanda, B. A. Small players, large role: Microbial influence on biogeochemical processes in pelagic aquatic ecosystems. Ecosystems. 5, 105–121 (2002).CAS 
Article 

Google Scholar 
Falkowski, P. G., Fenchel, T. & Delong, E. F. The microbial engines that drive earth’s biogeochemical cycles. Science. 320, 1034–1039 (2008).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Newton, R. J., Jones, S. E., Eiler, A., McMahon, K. D. & Bertilsson, S. A Guide to the Natural History of Freshwater Lake Bacteria. Microbiol. Mol. Biol. Rev. 75, 14–49 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Coleman, M. L. et al. Genomic islands and the ecology and evolution of Prochlorococcus. Science. 311, 1768–1770 (2006).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Franzosa, E. A. et al. Sequencing and beyond: Integrating molecular ‘omics’ for microbial community profiling. Nat. Rev. Microbiol. 13, 360–372 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hanson, C., Fuhrman, J., Horner-Devine, M. & Martiny, J. Beyond biogeographic patterns: processes shaping the microbial landscape. Nat. Rev. Microbiol. 10, 497–506 (2012).CAS 
PubMed 
Article 

Google Scholar 
Dai, A. & Trenberth, K. E. Estimates of freshwater discharge from continents: Latitudinal and seasonal variations. J. Hydrometeorol. 3, 660–687 (2002).ADS 
Article 

Google Scholar 
White, W. R. World water: resources, usage and the role of man-made reservoirs. Report No. FR/R0012. Fundation for Water Research, (2010).Clark, E. A., Sheffield, J., van Vliet, M. T. H., Nijssen, B. & Lettenmaier, D. P. Continental runoff into the oceans (1950–2008). J. Hydrometeorol. 16, 1502–1520 (2015).ADS 
Article 

Google Scholar 
Stevaux, J. C., Paes, R. J., Franco, A. A., Mário, M. L. & Fujita, R. H. Morphodynamics in the confluence of large regulated rivers: The case of Paraná and Paranapanema Rivers. Lat. Am. J. Sedimentol. Basin Anal. 16, 101–109 (2009).
Google Scholar 
Brêda, J. P. L. F. et al. Climate change impacts on South American water balance from a continental-scale hydrological model driven by CMIP5 projections. Clim. Change 159, 503–522 (2020).ADS 
Article 

Google Scholar 
Llames, M. E. & Zagarese, H. E. Lakes and Reservoirs of South America. In Encyclopedia of Inland Waters vol.2 (ed. Linkens, G. E.). (Oxford: Elsevier, 2009).Cabrera, A. L. & Willink, A. Biogeografia De America Latina 2da edn (Organización de los Estados Americanos, 1980).Morrone, J. J. Biogeografía de América Latina y el Caribe 1st edn. (Nature, 2001).Morrone, J. J. Biogeographical regionalisation of the neotropical region. Zootaxa 3782, 1–110 (2014).PubMed 
Article 

Google Scholar 
Antonelli, A. et al. Amazonia is the primary source of Neotropical biodiversity. Proc. Natl. Acad. Sci. USA 115, 6034–6039 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Sarmento, H. New paradigms in tropical limnology: The importance of the microbial food web. Hydrobiologia 686, 1–14 (2012).Article 

Google Scholar 
Meerhoff, M. et al. Environmental Warming in Shallow Lakes. A Review of Potential Changes in Community Structure as Evidenced from Space-for-Time Substitution Approaches. Adv. Ecol. Res. 46, 259–349 (2012).Article 

Google Scholar 
Herlemann, D. P. R. et al. Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. ISME J. 5, 1571–1579 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Metz, S. & Huber, P. et al. A georeferenced rRNA amplicon database of aquatic microbiomes from South America (Dataset), Zenodo, https://doi.org/10.5281/zenodo.6802178 (2022).Callahan, B. J., Sankaran, K., Fukuyama, J. A., McMurdie, P. J. & Holmes, S. P. Bioconductor Workflow for Microbiome Data Analysis: from raw reads to community analyses. F1000 Research 5, 1492 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal 17, 10 (2011).Article 

Google Scholar 
Edgar, R. C. UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing. Preprint at https://www.biorxiv.org/content/10.1101/081257v1 (2016).Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403–410 (1990).CAS 
PubMed 
Article 

Google Scholar 
Quast, C. et al. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2013).CAS 
PubMed 
Article 

Google Scholar 
Griffith, G. E., Omernik, J. M. & Azevedo, S. H. Ecological classification of the Western Hemisphere http://ecologicalregions.info/htm/ecoregions.htm (1998).Salcedo, J. C. R. South America: Argentina, Bolivia, and Peru https://www.worldwildlife.org/ecoregions/nt1002 Accessed (2018).Vidal, J. Geografía del Perú: las ocho regiones naturales, la regionalización transversal, la microregionalización 9th edn (PEISA, 1987).Paruelo, J. M., Beltran, A., Jobbagy, E., Sala, O. E. & Golluscio, R. A. The climate of Patagonia: General patterns and controls on biotic processes. Ecol. Austral 8, 85–101 (1998).
Google Scholar 
Iriondo, M. Quaternary lakes of Argentina. Palaeogeogr. Palaeoclimatol. Palaeoecol. 70, 81–88 (1989).Article 

Google Scholar 
Soto, D. & Campos, H. in Ecología de los bosques templados de Chile vol. 1 (eds. Khalin, J. M. & Villagrán, C.) (Editorial Universitaria, 1995).Modenutti, B. et al. Structure and dynamic of food webs in Andean North Patagonian freshwater systems: Organic matter, light and nutrient relationships. Ecol. Austral 20, 95–114 (2010).
Google Scholar 
Modenutti, B. E. et al. Structure and dynamics of food webs in Andean lakes. Lakes Reserv. Res. Manag. 3, 179–186 (1998).Article 

Google Scholar 
Quirós, R. & Drago, E. The environmental state of Argentinean lakes: An overview. Lakes Reserv. Res. Manag. 4, 55–64 (1999).Article 

Google Scholar 
Morris, D. P. et al. The attenuation of solar UV radiation in lakes and the role of dissolved organic carbon. Limnol. Oceanogr. 40, 1381–1391 (1995).ADS 
CAS 
Article 

Google Scholar 
Bastidas Navarro, M., Balseiro, E. & Modenutti, B. Bacterial Community Structure in Patagonian Andean Lakes Above and Below Timberline: From Community Composition to Community Function. Microb. Ecol. 68, 528–541 (2014).PubMed 
Article 

Google Scholar 
Modenutti, B. et al. Environmental changes affecting light climate in oligotrophic mountain lakes: The deep chlorophyll maxima as a sensitive variable. Aquat. Sci. 75, 361–371 (2013).CAS 
Article 

Google Scholar 
Bastidas Navarro, M., Martyniuk, N., Balseiro, E. & Modenutti, B. Effect of glacial lake outburst floods on the light climate in an Andean Patagonian lake: implications for planktonic phototrophs. Hydrobiologia 816, 39–48 (2018).CAS 
Article 

Google Scholar 
Sioli, H. Hydrochemistry and Geology in the Brazilian Amazon Region. Amazoniana 1, 267–277 (1968).
Google Scholar 
Fick, S. E. & Hijmans, R. J. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315 (2017).Article 

Google Scholar 
Salati, E. & Vose, P. B. Amazon Basin: A system in equilibrium. Science. 225, 129–138 (1984).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Melack, J. M. & Forsberg, B. R. In The Biogeochemistry of the Amazon Basin Vol. 1 (eds. MacCLain, M. E., Victoria, R. & Richey, J. E.). (Oxford Scholarship Online, 2001).Junk, W. J., Bayley, P. B. & Sparks, R. E. The flood pulse concept in river-floodplain systems. Can. J. Fish. Aquat. Sci. 106, 110–127 (1989).
Google Scholar 
Ratter, J. A., Ribeiro, J. F. & Bridgewater, S. The Brazilian cerrado vegetation and threats to its biodiversity. Ann. Bot. 80, 223–230 (1997).Article 

Google Scholar 
Haridasan, M. Nutritional adaptations of native plants of the cerrado biome in acid soils. Braz. J. Plant Physiol. 20, 183–195 (2008).Article 

Google Scholar 
Vasconcelos, V., de Carvalho Júnior, O. A., de Souza Martins, É. & Couto Júnior, A. F. in World Geomorphological Landscapes. Vol. 1 (eds. Vieira, B., Salgado, A. & Santos, L.) (Springer, 2015).Bichsel, D. et al. Water quality of rural ponds in the extensive agricultural landscape of the Cerrado (Brazil). Limnology 17, 239–246 (2016).CAS 
Article 

Google Scholar 
Cunha, D. G. F., Calijuri, M., do, C. & Lamparelli, M. C. A trophic state index for tropical/subtropical reservoirs (TSItsr). Ecol. Eng. 60, 126–134 (2013).Article 

Google Scholar 
Morellato, L. P. C. & Haddad, C. F. B. Introduction: The Brazilian atlantic forest. Biotropica 32, 786–792 (2000).Article 

Google Scholar 
Galindo-Leal, C. & Câmara, I. de G. The Atlantic Forest of South America: Biodiversity status, threats, and outlook 1st edn (Island Press, 2003).Joly, C. A., Metzger, J. P. & Tabarelli, M. Experiences from the Brazilian Atlantic Forest: Ecological findings and conservation initiatives. New Phytologist 204, 459–473 (2014).PubMed 
Article 

Google Scholar 
Caliman, A. et al. Temporal coherence among tropical coastal lagoons: A search for patterns and mechanisms. Brazilian J. Biol. 70, 803–814 (2010).CAS 
Article 

Google Scholar 
Junger, P. C. et al. Salinity Drives the Virioplankton Abundance but Not Production in Tropical Coastal Lagoons. Microb. Ecol. 75, 52–63 (2018).CAS 
PubMed 
Article 

Google Scholar 
Depetris, P. J., Probst, J. L., Pasquini, A. I. & Gaiero, D. M. The geochemical characteristics of the Paraná River suspended sediment load: An initial assessment. Hydrol. Process. 17, 1267–1277 (2003).ADS 
Article 

Google Scholar 
Orfeo, O. & Stevaux, J. Hydraulic and morphological characteristics of middle and upper reaches of the Paraná River (Argentina and Brazil). Geomorphology 44, 309–322 (2002).ADS 
Article 

Google Scholar 
Neiff, J. J. Large rivers of South America: toward the new approach. Verh. Internat. Verein. Limnol 26, 167–180 (1996).
Google Scholar 
Unrein, F. Changes in phytoplankton community along a transversal section of the Lower Paraná floodplain, Argentina. Hydrobiologia 468, 123–134 (2002).Article 

Google Scholar 
Devercelli, M. Changes in phytoplankton morpho-functional groups induced by extreme hydroclimatic events in the Middle Paraná river (Argentina). Hydrobiologia 639, 5–19 (2010).CAS 
Article 

Google Scholar 
Huber, P. et al. Environmental heterogeneity determines the ecological processes that govern bacterial metacommunity assembly in a floodplain river system. ISME J. 14, 2951–2966 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Olson, D. M. et al. Terrestrial ecoregions of the world: A new map of life on Earth. Bioscience 51, 933–938 (2001).Article 

Google Scholar 
Peel, M. C., Finlayson, B. L. & McMahon, T. A. Updated world map of the Köppen-Geiger climate classification. Hydrol. Earth Syst. Sci. 11, 1633–1644 (2007).ADS 
Article 

Google Scholar 
Conde, D., Arocena, R. & Recursos, R.-G. L. acuáticos superficiales de Uruguay: ambientes, algunas problemáticas y desafios para la gestión. Ambios 10, 1–7 (2003).
Google Scholar 
Martin, L. & Suguio, K. Variation of coastal dynamics during the last 7000 years recorded in beach-ridge plains associated with river mouths: example from the central Brazilian coast. Palaeogeogr. Palaeoclimatol. Palaeoecol. 99, 119–140 (1992).Article 

Google Scholar 
Alonso, C. et al. Environmental dynamics as a structuring factor for microbial carbon utilization in a subtropical coastal lagoon. Front. Microbiol. 4, 1664–302X (2013).Article 
CAS 

Google Scholar 
Amaral, V., Graeber, D., Calliari, D. & Alonso, C. Strong linkages between DOM optical properties and main clades of aquatic bacteria. Limnol. Oceanogr. 61, 906–918 (2016).ADS 
Article 

Google Scholar 
Rennella, A. M. M., Quiro, R. & Quirós, R. The effects of hydrology on plankton biomass in shallow lakes of the Pampa Plain. Hydrobiologia 556, 181–191 (2006).Article 

Google Scholar 
Diaz, M., Pedrozo, F. & Baccala, N. Summer classification of Southern Hemisphere temperate lakes (Patagonia, Argentina). Lakes Reserv. Res. Manag. 5, 213–229 (2000).Article 

Google Scholar 
Izaguirre, I. et al. Influence of fish introduction and water level decrease on lakes of the arid Patagonian plateaus with importance for biodiversity conservation. Glob. Ecol. Conserv. 14, e00391 (2018).Article 

Google Scholar 
Porcel, S., Saad, J. F., Sabio y García, C. A. & Izaguirre, I. Microbial planktonic communities in lakes from a Patagonian basaltic plateau: influence of the water level decrease. Aquat. Sci. 81, 51 (2019).Article 
CAS 

Google Scholar 
Bernal, M. C. et al. Spatial variation of picoplankton communities along a cascade reservoir system in Patagonia, Argentina. J. Limnol. 80, 84–99 (2021).
Google Scholar 
Leinonen, R. et al. The European nucleotide archive. Nucleic Acids Res. 39, 44–47 (2011).Article 
CAS 

Google Scholar 
ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA217932 (2013).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA302313 (2017).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA294718 (2022).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA309832 (2016).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA326475 (2016).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJEB48609 (2022).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA289691 (2015).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA414894 (2018).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA323673 (2016).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA356055 (2017).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA310230 (2017).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA390178 (2019).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA411849 (2017).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA725228 (2021).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA292014 (2015).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA310230 (2017).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA411849 (2017).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA316315 (2017).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA406945 (2017).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA515842 (2019).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA310230 (2017).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA321235 (2016).ENA European Nucleotide Archive https://identifiers.org/ena.embl:SAMN07998328 (2015).ENA European Nucleotide Archive https://identifiers.org/ena.embl:SAMN07998330 (2015).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJEB36116 (2020).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJEB29989 (2019).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA788397 (2021).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJEB48353 (2022).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJEB37379 (2020).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJEB46122 (2021).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJEB40710 (2020).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJEB40864 (2020).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJEB40854 (2020).ENA European Nucleotide Archive https://identifiers.org/ena.embl:PRJNA268541 (2015). More