Ecology

Taylor, M. W., Radax, R., Steger, D. & Wagner, M. Sponge-associated microorganisms: Evolution, ecology, and biotechnological potential. Microbiol. Mol. Biol. Rev. 71, 295–347 (2007).CAS 
PubMed 
PubMed Central 

Google Scholar 
Thomas, T. et al. Diversity, structure and convergent evolution of the global sponge microbiome. Nat. Commun. 7, 11870 (2016).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Webster, N. S. et al. Deep sequencing reveals exceptional diversity and modes of transmission for bacterial sponge symbionts. Environ. Microbiol. 12, 2070–2082 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
Sipkema, D. et al. Similar sponge-associated bacteria can be acquired via both vertical and horizontal transmission: Microbial transmission in Petrosia ficiformis. Environ. Microbiol. 17, 3807–3821 (2015).CAS 
PubMed 

Google Scholar 
Cleary, D. F. R. et al. The sponge microbiome within the greater coral reef microbial metacommunity. Nat. Commun. 10, 1644 (2019).Björk, J. R., Díez-Vives, C., Astudillo-García, C., Archie, E. A. & Montoya, J. M. Vertical transmission of sponge microbiota is inconsistent and unfaithful. Nat. Ecol. Evol. 3, 1172–1183 (2019).PubMed 
PubMed Central 

Google Scholar 
Webster, N. S. & Taylor, M. W. Marine sponges and their microbial symbionts: Love and other relationships. Environ. Microbiol. 14, 335–346 (2012).CAS 
PubMed 

Google Scholar 
Kennedy, J. et al. Evidence of a putative deep sea specific microbiome in marine sponges. PLoS ONE 9, e91092 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
Steinert, G. et al. Compositional and quantitative insights into bacterial and archaeal communities of south pacific deep-sea sponges (Demospongiae and Hexactinellida). Front. Microbiol. 11, 716 (2020).Busch, K. et al. On giant shoulders: How a seamount affects the microbial community composition of seawater and sponges. Biogeosciences 17, 3471–3486 (2020).ADS 
CAS 

Google Scholar 
Olson, J. B. & Gao, X. Characterizing the bacterial associates of three Caribbean sponges along a gradient from shallow to mesophotic depths. FEMS Microbiol. Ecol. 85, 74–84 (2013).PubMed 

Google Scholar 
Steinert, G. et al. In four shallow and mesophotic tropical reef sponges from Guam the microbial community largely depends on host identity. PeerJ 4, e1936 (2016).PubMed 
PubMed Central 

Google Scholar 
Morrow, K. M., Fiore, C. L. & Lesser, M. P. Environmental drivers of microbial community shifts in the giant barrel sponge, Xestospongia muta, over a shallow to mesophotic depth gradient. Environ. Microbiol. 18, 2025–2038 (2016).CAS 
PubMed 

Google Scholar 
Ebada, S. S. & Proksch, P. The chemistry of marine sponges. In Handbook of Marine Natural Products (eds Fattorusso, E. et al.) 191–293 (Springer, 2012). https://doi.org/10.1007/978-90-481-3834-0_4.Chapter 

Google Scholar 
Kornprobst, J.-M. Porifera (Sponges). Encyclopedia of Marine Natural Products (Wiley, 2014).
Google Scholar 
Leal, M. C., Puga, J., Serôdio, J., Gomes, N. C. M. & Calado, R. Trends in the discovery of new marine natural products from invertebrates over the last two decades—Where and what are we bioprospecting?. PLoS ONE 7, e30580 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Blunt, J. W., Copp, B. R., Keyzers, R. A., Munro, M. H. G. & Prinsep, M. R. Marine natural products. Nat. Prod. Rep. 34, 235–294 (2017).CAS 
PubMed 

Google Scholar 
Unson, M. D., Holland, N. D. & Faulkner, D. J. A brominated secondary metabolite synthesized by the cyanobacterial symbiont of a marine sponge and accumulation of the crystalline metabolite in the sponge tissue. Mar. Biol. 119, 1–11 (1994).CAS 

Google Scholar 
Bewley, C. A., Holland, N. D. & Faulkner, D. J. Two classes of metabolites from Theonella swinhoei are localized in distinct populations of bacterial symbionts. Experientia 52, 716–722 (1996).CAS 
PubMed 

Google Scholar 
Wilson, M. C. et al. An environmental bacterial taxon with a large and distinct metabolic repertoire. Nature 506, 58–62 (2014).ADS 
CAS 
PubMed 

Google Scholar 
Tianero, M. D., Balaich, J. N. & Donia, M. S. Localized production of defence chemicals by intracellular symbionts of Haliclona sponges. Nat. Microbiol. 4, 1149–1159 (2019).CAS 
PubMed 

Google Scholar 
Ivanišević, J., Thomas, O. P., Lejeusne, C., Chevaldonné, P. & Pérez, T. Metabolic fingerprinting as an indicator of biodiversity: Towards understanding inter-specific relationships among Homoscleromorpha sponges. Metabolomics 7, 289–304 (2011).
Google Scholar 
Pérez, T. et al. Oscarella balibaloi, a new sponge species (Homoscleromorpha: Plakinidae) from the Western Mediterranean Sea: Cytological description, reproductive cycle and ecology: O. balibaloi: Description, reproductive cycle and ecology. Mar. Ecol. (Berl.) 32, 174–187 (2011).ADS 

Google Scholar 
Reveillaud, J. et al. Relevance of an integrative approach for taxonomic revision in sponge taxa: Case study of the shallow-water Atlanto-Mediterranean Hexadella species (Porifera: Ianthellidae: Verongida). Invertebr. Syst. 26, 230–248 (2012).
Google Scholar 
Olsen, E. K. et al. Marine AChE inhibitors isolated from Geodia barretti: Natural compounds and their synthetic analogs. Org. Biomol. Chem. 14, 1629–1640 (2016).CAS 
PubMed 

Google Scholar 
Reverter, M., Perez, T., Ereskovsky, A. V. & Banaigs, B. Secondary metabolome variability and inducible chemical defenses in the Mediterranean Sponge Aplysina cavernicola. J. Chem. Ecol. 42, 60–70 (2016).CAS 
PubMed 

Google Scholar 
Reverter, M., Tribalat, M.-A., Pérez, T. & Thomas, O. P. Metabolome variability for two Mediterranean sponge species of the genus Haliclona: Specificity, time, and space. Metabolomics 14, 114 (2018).Villegas-Plazas, M. et al. Variations in microbial diversity and metabolite profiles of the tropical marine sponge Xestospongia muta with season and depth. Microb. Ecol. 78, 243–256 (2019).CAS 
PubMed 

Google Scholar 
Mohanty, I. et al. Multi-omic profiling of Melophlus sponges reveals diverse metabolomic and microbiome architectures that are non-overlapping with ecological neighbors. Mar. Drugs 18, 124 (2020).CAS 
PubMed Central 

Google Scholar 
Bowerbank, J. S. On the anatomy and physiology of the Spongiadae. Part I. On the spicula. Philos. Trans. R. Soc. Lond. 148, 279–332 (1858).ADS 

Google Scholar 
Vosmaer, G. C. J. The sponges of the ‘Willem Barents’ expedition 1880 and 1881. Bijdragen tot de Dierkunde 12, 1–47 (1885).
Google Scholar 
Radax, R. et al. Metatranscriptomics of the marine sponge Geodia barretti: Tackling phylogeny and function of its microbial community. Environ. Microbiol. 14, 1308–1324 (2012).CAS 
PubMed 

Google Scholar 
Topsent, E. Spongiaires provenant des campagnes scientifiques de la ‘Princesse Alice’ dans les Mers du Nord (1898–1899—1906–1907). Résultats des campagnes scientifiques accomplies par le Prince Albert I. Monaco 45, 1–67 (1913).
Google Scholar 
Yashayaev, I. & Loder, J. W. Further intensification of deep convection in the Labrador Sea in 2016. Geophys. Res. Lett. 44, 1429–1438 (2017).ADS 

Google Scholar 
Gutleben, J. et al. Diversity of tryptophan halogenases in sponges of the genus Aplysina. FEMS Microbiol. Ecol. 95, fiz108 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Indraningrat, A. et al. Cultivation of sponge-associated bacteria from Agelas sventres and Xestospongia muta collected from different depths. Mar. Drugs 17, 578 (2019).CAS 
PubMed Central 

Google Scholar 
Ramiro-Garcia, J. et al. NG-Tax, a highly accurate and validated pipeline for analysis of 16S rRNA amplicons from complex biomes. F1000 Res. 5, 1791 (2018).
Google Scholar 
Yilmaz, P. et al. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucl. Acids Res. 42, D643–D648 (2014).CAS 
PubMed 

Google Scholar 
Erngren, I., Smit, E., Pettersson, C., Cárdenas, P. & Hedeland, M. The effects of sampling and storage conditions on the metabolite profile of the marine sponge Geodia barretti. Front. Chem. 9:662659 (2021)Smith, C. A., Want, E. J., O’Maille, G., Abagyan, R. & Siuzdak, G. XCMS: Processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal. Chem. 78, 779–787 (2006).CAS 
PubMed 

Google Scholar 
Kuhl, C., Tautenhahn, R., Böttcher, C., Larson, T. R. & Neumann, S. CAMERA: An integrated strategy for compound spectra extraction and annotation of liquid chromatography/mass spectrometry data sets. Anal. Chem. 84, 283–289 (2012).CAS 
PubMed 

Google Scholar 
Oksanen, J. et al. vegan: Community Ecology Package (2017).Dat, T. T. H., Steinert, G., Thi Kim Cuc, N., Smidt, H. & Sipkema, D. Archaeal and bacterial diversity and community composition from 18 phylogenetically divergent sponge species in Vietnam. PeerJ 6, e4970 (2018).PubMed 
PubMed Central 

Google Scholar 
Miller, M. A., Pfeiffer, W. & Schwartz, T. Creating the CIPRES science gateway for inference of large phylogenetic trees. In 2010 Gateway Computing Environments Workshop (GCE) 1–8 (IEEE, 2010). https://doi.org/10.1109/GCE.2010.5676129.Letunic, I. & Bork, P. Interactive Tree Of Life (iTOL) v4: Recent updates and new developments. Nucl. Acids Res. 47, W256–W259 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Thévenot, E. A., Roux, A., Xu, Y., Ezan, E. & Junot, C. Analysis of the human adult urinary metabolome variations with age, body mass index, and gender by implementing a comprehensive workflow for univariate and OPLS statistical analyses. J. Proteome Res. 14, 3322–3335 (2015).PubMed 

Google Scholar 
Weiss, S. et al. Correlation detection strategies in microbial data sets vary widely in sensitivity and precision. ISME J. 10, 1669–1681 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Deng, Y. et al. Molecular ecological network analyses. BMC Bioinform. 13, 113 (2012).
Google Scholar 
Friedman, J. & Alm, E. J. Inferring correlation networks from genomic survey data. PLoS Comput. Biol. 8, e1002687 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Durno, W. E., Hanson, N. W., Konwar, K. M. & Hallam, S. J. Expanding the boundaries of local similarity analysis. BMC Genom. 14, S3 (2013).
Google Scholar 
Reshef, D. N. et al. Detecting novel associations in large data sets. Science 334, 1518–1524 (2011).ADS 
CAS 
PubMed 
PubMed Central 
MATH 

Google Scholar 
Hall, M. M., Torres, D. J. & Yashayaev, I. Absolute velocity along the AR7W section in the Labrador Sea. Deep Sea Res. Part 1 Oceanogr. Res. Pap. 72, 72–87 (2013).
Google Scholar 
Reveillaud, J. et al. Host-specificity among abundant and rare taxa in the sponge microbiome. ISME J. 8, 1198–1209 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
Moitinho-Silva, L. et al. Predicting the HMA-LMA status in marine sponges by machine learning. Front. Microbiol. 8, 752 (2017).Lidgren, G., Bohlin, L. & Bergman, J. Studies of Swedish marine organisms VII. A novel biologically active indole alkaloid from the sponge Geodia barretti. Tetrahedron Lett. 27, 3283–3284 (1986).CAS 

Google Scholar 
Sjögren, M. et al. Antifouling activity of brominated cyclopeptides from the marine sponge Geodia barretti. J. Nat. Prod. 67, 368–372 (2004).PubMed 

Google Scholar 
Sölter, S. Identifizierung und Synthese von Naturstoffen aus Borealen Schwämmen (Universität Hamburg, 2004).
Google Scholar 
Di, X. et al. 6-Bromoindole derivatives from the Icelandic marine sponge Geodia barretti: Isolation and anti-inflammatory activity. Mar. Drugs 16, 437 (2018).CAS 
PubMed Central 

Google Scholar 
Carstens, B. B. et al. Isolation, characterization, and synthesis of the barrettides: Disulfide-containing peptides from the marine sponge Geodia barretti. J. Nat. Prod. 78, 1886–1893 (2015).CAS 
PubMed 

Google Scholar 
Hedner, E. et al. Brominated cyclodipeptides from the marine sponge Geodia barretti as selective 5-HT ligands. J. Nat. Prod. 69, 1421–1424 (2006).CAS 
PubMed 

Google Scholar 
Hedner, E. et al. Antifouling activity of a dibrominated cyclopeptide from the marine sponge Geodia barretti. J. Nat. Prod. 71, 330–333 (2008).CAS 
PubMed 

Google Scholar 
Erwin, P. M., Pita, L., López-Legentil, S. & Turon, X. Stability of sponge-associated bacteria over large seasonal shifts in temperature and irradiance. Appl. Environ. Microbiol. 78, 7358–7368 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Cárdenas, C. A., Bell, J. J., Davy, S. K., Hoggard, M. & Taylor, M. W. Influence of environmental variation on symbiotic bacterial communities of two temperate sponges. FEMS Microbiol. Ecol. 88, 516–527 (2014).PubMed 

Google Scholar 
Glasl, B., Smith, C. E., Bourne, D. G. & Webster, N. S. Exploring the diversity-stability paradigm using sponge microbial communities. Sci. Rep. 8, 8425 (2018).Schöttner, S. et al. Relationships between host phylogeny, host type and bacterial community diversity in cold-water coral reef sponges. PLoS ONE 8, e55505 (2013).ADS 
PubMed 
PubMed Central 

Google Scholar 
Lurgi, M., Thomas, T., Wemheuer, B., Webster, N. S. & Montoya, J. M. Modularity and predicted functions of the global sponge-microbiome network. Nat. Commun. 10, 992 (2019).ADS 
PubMed 
PubMed Central 

Google Scholar 
Luter, H. M. et al. Microbiome analysis of a disease affecting the deep-sea sponge Geodia barretti. FEMS Microbiol. Ecol. 93, fix074 (2017).Thistle, D. Ecosystems of the Deep Oceans (Elsevier, 2003).
Google Scholar 
Pita, L., Erwin, P. M., Turon, X. & López-Legentil, S. Till death do us part: Stable sponge-bacteria associations under thermal and food shortage stresses. PLoS ONE 8, e80307 (2013).ADS 
PubMed 
PubMed Central 

Google Scholar 
Webster, N. S., Cobb, R. E. & Negri, A. P. Temperature thresholds for bacterial symbiosis with a sponge. ISME J. 2, 830–842 (2008).CAS 
PubMed 

Google Scholar 
Gerringer, M. E., Drazen, J. C. & Yancey, P. H. Metabolic enzyme activities of abyssal and hadal fishes: Pressure effects and a re-evaluation of depth-related changes. Deep Sea Res. Part 1 Oceanogr. Res. Pap. 125, 135–146 (2017).CAS 

Google Scholar 
Yashayaev, I. Hydrographic changes in the Labrador Sea, 1960–2005. Prog. Oceanogr. 73, 242–276 (2007).ADS 

Google Scholar 
Rhein, M., Steinfeldt, R., Kieke, D., Stendardo, I. & Yashayaev, I. Ventilation variability of Labrador Sea Water and its impact on oxygen and anthropogenic carbon: A review. Philos. Trans. A Math. Phys. Eng. Sci. 375, 20160321 (2017).ADS 
PubMed 
PubMed Central 

Google Scholar 
Galand, P. E., Potvin, M., Casamayor, E. O. & Lovejoy, C. Hydrography shapes bacterial biogeography of the deep Arctic Ocean. ISME J. 4, 564–576 (2010).PubMed 

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

Google Scholar 
Agogué, H., Lamy, D., Neal, P. R., Sogin, M. L. & Herndl, G. J. Water mass-specificity of bacterial communities in the North Atlantic revealed by massively parallel sequencing. Mol. Ecol. 20, 258–274 (2011).PubMed 

Google Scholar 
Djurhuus, A., Boersch-Supan, P. H., Mikalsen, S.-O. & Rogers, A. D. Microbe biogeography tracks water masses in a dynamic oceanic frontal system. R. Soc. Open Sci. 4, 170033 (2017).ADS 
PubMed 
PubMed Central 

Google Scholar 
Müller, O. et al. Spatiotemporal dynamics of ammonia-oxidizing Thaumarchaeota in distinct Arctic water masses. Front. Microbiol. 9, 1–13 (2018).ADS 

Google Scholar 
Kraemer, S., Ramachandran, A., Colatriano, D., Lovejoy, C. & Walsh, D. A. Diversity and biogeography of SAR11 bacteria from the Arctic Ocean. ISME J. https://doi.org/10.1038/s41396-019-0499-4 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Monier, A. et al. Upper Arctic Ocean water masses harbor distinct communities of heterotrophic flagellates. Biogeosciences 10, 4273–4286 (2013).ADS 

Google Scholar 
Monier, A. et al. Oceanographic structure drives the assembly processes of microbial eukaryotic communities. ISME J. 9, 990–1002 (2015).CAS 
PubMed 

Google Scholar 
Corrège, T. The relationship between water masses and benthic ostracod assemblages in the western Coral Sea, Southwest Pacific. Palaeogeogr. Palaeoclimatol. Palaeoecol. 105, 245–266 (1993).
Google Scholar 
Muhling, B. A., Beckley, L. E., Koslow, J. A. & Pearce, A. F. Larval fish assemblages and water mass structure off the oligotrophic south-western Australian coast: SW Australian larval fish assemblages. Fish. Oceanogr. 17, 16–31 (2007).
Google Scholar 
Eerkes-Medrano, D. et al. A community assessment of the demersal fish and benthic invertebrates of the Rosemary Bank Seamount Marine Protected Area (NE Atlantic). Deep Sea Res. Part 1 Oceanogr. Res. Pap. https://doi.org/10.1016/j.dsr.2019.103180 (2019).Article 

Google Scholar 
Puerta, P. et al. Influence of water masses on the biodiversity and biogeography of deep-sea benthic ecosystems in the North Atlantic. Front. Mar. Sci. 7, 239 (2020).Roberts, E. et al. Water masses constrain the distribution of deep-sea sponges in the North Atlantic Ocean and Nordic Seas. Mar. Ecol. Prog. Ser. 659, 75–96 (2021).ADS 

Google Scholar 
Kenchington, E. et al. Connectivity modelling of areas closed to protect vulnerable marine ecosystems in the northwest Atlantic. Deep Sea Res. Part 1 Oceanogr. Res. Pap. 143, 85–103 (2019).
Google Scholar 
Louca, S. et al. Function and functional redundancy in microbial systems. Nat. Ecol. Evol. 2, 936–943 (2018).PubMed 

Google Scholar 
McCauley, M., Chiarello, M., Atkinson, C. L. & Jackson, C. R. Gut microbiomes of freshwater mussels (Unionidae) are taxonomically and phylogenetically variable across years but remain functionally stable. Microorganisms 9, 411 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Page, M., West, L., Northcote, P., Battershill, C. & Kelly, M. Spatial and temporal variability of cytotoxic metabolites in populations of the New Zealand Sponge Mycale hentscheli. J. Chem. Ecol. 31, 1161–1174 (2005).CAS 
PubMed 

Google Scholar 
Ternon, E., Perino, E., Manconi, R., Pronzato, R. & Thomas, O. P. How environmental factors affect the production of guanidine alkaloids by the Mediterranean sponge Crambe crambe. Mar. Drugs 15, 181 (2017).PubMed Central 

Google Scholar 
Sacristán-Soriano, O., Banaigs, B. & Becerro, M. A. Temporal trends in the secondary metabolite production of the sponge Aplysina aerophoba. Mar. Drugs 10, 677–693 (2012).PubMed 
PubMed Central 

Google Scholar 
Ivanisevic, J. et al. Biochemical trade-offs: Evidence for ecologically linked secondary metabolism of the sponge Oscarella balibaloi. PLoS ONE 6, e28059 (2011).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Burg, M. B. & Ferraris, J. D. Intracellular organic osmolytes: Function and regulation. J. Biol. Chem. 283, 7309–7313 (2008).CAS 
PubMed 
PubMed Central 

Google Scholar 
Nau-Wagner, G., Boch, J., Le Good, J. A. & Bremer, E. High-affinity transport of choline-O-sulfate and its use as a compatible solute in Bacillus subtilis. Appl. Environ. Microbiol. 65, 560–568 (1999).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Popowich, A., Zhang, Q. & Le, X. C. Arsenobetaine: The ongoing mystery. Natl. Sci. Rev. 3, 451–458 (2016).CAS 

Google Scholar 
Connor, K. M. & Gracey, A. Y. High-resolution analysis of metabolic cycles in the intertidal mussel Mytilus californianus. Am. J. Physiol. Regul. Integr. Comp. Physiol. 302, R103–R111 (2012).CAS 
PubMed 

Google Scholar 
Cárdenas, P. Who produces Ianthelline? The Arctic sponge Stryphnus fortis or its sponge Epibiont Hexadella dedritifera: A probable case of sponge–sponge contamination. J. Chem. Ecol. 42, 339–347 (2016).PubMed 

Google Scholar 
Steffen, K. et al. Barrettides: A peptide family specifically produced by the deep-sea sponge Geodia barretti. J. Nat. Prod. 84, 3138–3146 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Abbamondi, G. R., De Rosa, S., Iodice, C. & Tommonaro, G. Cyclic dipeptides produced by marine sponge-associated bacteria as quorum sensing signals. Nat. Prod. Commun. 9, 229–232 (2014).CAS 
PubMed 

Google Scholar 
Kasheverov, I. et al. 6-Bromohypaphorine from Marine Nudibranch Mollusk Hermissenda crassicornis is an agonist of human α7 nicotinic acetylcholine receptor. Mar. Drugs 13, 1255–1266 (2015).CAS 
PubMed 
PubMed Central 

Google Scholar 
Moitinho-Silva, L. et al. The sponge microbiome project. Gigascience 6, 1–7 (2017).CAS 
PubMed 

Google Scholar 
Kielak, A. M., Barreto, C. C., Kowalchuk, G. A., van Veen, J. A. & Kuramae, E. E. The ecology of acidobacteria: Moving beyond genes and genomes. Front. Microbiol. 7, 744 (2016).Crits-Christoph, A., Diamond, S., Butterfield, C. N., Thomas, B. C. & Banfield, J. F. Novel soil bacteria possess diverse genes for secondary metabolite biosynthesis. Nature 558, 440–444 (2018).ADS 
CAS 
PubMed 

Google Scholar  More

To our knowledge, this study provides the first recorded information on the active swimming of Japanese eels and on their transport by currents in the open ocean. Specifically, the strong flow of the KC largely dominated the movements of the eels and transported them northeastward while they swam mainly southward, and active swimming contributed a little to their travel trajectories. In contrast, the swimming of eels made a relatively higher contribution to their travel trajectories in the TS area.Our in situ estimates of the mean swimming speeds of Japanese eels (26–41 cm/s) were similar or slightly lower than those of European eels. In the acoustic tracking experiment of European eels considering environmental current vectors, their swimming speeds were 35–58 cm/s in the coastal midwater26. In a laboratory experiment using stamina tunnels with stable temperatures, the optimal swimming speeds of European eels were estimated to be 61–68 cm/s (0.74–1.02 BL/s)56, which were higher than the in situ estimates. The minimum swimming speed of European eels is considered to be 40 cm/s if they will arrive at their spawning area in the Sargasso Sea (distance of 5500 km) in 6 months, and their optimal swimming speeds were sufficient to migrate over the long distance in time for the near-spawning period after escape from their growth habitats56. However, field studies using PSAT tagging also reported that in situ migration speeds (including transport by currents) were less than the optimal swimming speeds and suggested that some European eels could reach their spawning area within the near-spawning periods and that others only arrive in time for the following spawning season19.Our estimated effective swimming speed of Japanese eels, all day and all night over the tracking periods, ranged from 3 to 30 cm/s with individual variations. These estimates were consistent with the swimming speeds (excluding transport by currents) of 2.2–15.1 km/day (2–18 cm/s) estimated in the PSAT study of Japanese eels14. Silver-phase Japanese eels start migrating from their coastal growth habitats in Japan primarily in October to December57, 58, and spawning near the West Mariana Ridge occurs in April to August33, 35. Numerical models assuming that migrating eels use true navigation (readjusted compass) or a constant compass heading (fixed compass from the departure place to the spawning site) indicate that the minimal swimming speed required to arrive at the spawning area within 8 months is 10–12 cm/s37. Our estimated effective swimming speeds of five out of ten eels during the day and eight out of ten eels during the night were similar or higher than these minimal speeds. The low effective swimming speeds frequently observed during the day might be due to the relatively low values observed in the swimming speed at 10 min intervals and the swimming directions often varying during the day. When eels swim with stable orientation, as observed in three of the eels (WE2999_TS, WE3001_TS, and WE3002_TS) during the night, the effective swimming speeds exceeded 25 cm/s. If such a stable orientation is maintained and compensate the low speeds during the day, the eels that leave during autumn and winter will be able to arrive at the spawning area during the next spring to summer.It should also be noted that the swimming speeds in body length per second were significantly higher in shallow water during the night than in deep water during the day. In the open ocean, anguillid species exhibit DVMs during oceanic migration, swimming at depth during the day and in the shallows during the night9,10,11,12,13,14,15,16,17,18,19,20,21,22. These DVMs are likely related to the possible avoidance from visual predators under light conditions19 or maturation control59. Essentially, through the DVMs, the eels encounter low temperatures ( 20 °C) during the same day. Generally, the swimming speeds of fishes are restricted by the ambient water temperature60, and the water temperature encountered through DVMs might influence the horizontal-swimming speeds of Japanese eels.Other factors besides swimming speed are important for the success of eel migrations, such as adapting to mesopelagic zones that silver eels undergo during their spawning migrations. The most important and obvious morphological adaptation in mesopelagic fish is their well-developed eyes, and migrating eels also seem to use this strategy. These fish often have relatively large pupils61, high photosensitive structures, such as tubular eyes62, a pure rod multibank retina63, and maximum rhodopsin absorption to adapt to the blue-green light in the deep sea64. The eyes of catadromous eels displayed enlargement during their transformation into migrating silver-phase eels65, 66 and potentially increase their retinal surface area, which results in the possibility of increased photon capture. In addition, the rhodopsins in the eyes change from a freshwater type with a maximum absorption of ~ 500 nm to a deep-sea type with a maximum absorption of ~ 480 nm67,68,69. Their extreme sensitivity to light is evident through their DVM in mesopelagic water, where the timing of a large descent and ascent in the DVM demonstrated by migrating catadromous eels is precisely synchronized with sunrise and sunset. Furthermore, eels alter their swimming depth in response to the phase of the Moon9, 15, 20, 21, appearing to be capable of perceiving extremely low-intensity moonlight.This study showed that three eels released in the TS area (mainly 300–400-m depth) and one eel in the KC area (near surface) were found to change their swimming direction around the time of the solar culmination when the Sun’s bearing changed. The clockwise and counterclockwise trajectories of these eels corresponded to whether the Sun moved from the east to west in the southern and northern sky, suggesting that they demonstrated horizontal negative phototaxis swimming to avoid sunlight. They might move to avoid high-intensity sunlight horizontally, not vertically, as they gradually increase the swimming depths possibly due to acclimation to cold deep water after release. The daytime swimming depths of the eels became deeper day-by-day after their release (Fig. 4); a similar phenomenon was observed in European eels12, American eels17, and long fin eels13. Recently, Higuchi et al.20 observed that the daytime swimming depths of Japanese eels released in the TS area gradually became deeper until 13 days after their release. These facts indicate that they gradually acclimate to the cold water at the deep depths after release. Since this tracking study was conducted 2–8 days after their release, the daytime swimming depth of eels would not have reached a steady state yet. The relatively high intensity from sunlight at the shallow depths where eels swam immediately after release in the TS area might cause horizontal avoidance behavior from the light.In other cases, many eels, especially those released in the KC area, did not demonstrate the rotational behavior. The eels in the KC area mostly stayed deeper (500–800 m) during the day than the eels in the TS area (stayed at depths of 300–600 m) even during the periods shortly after their release. This is possibly due to higher water temperatures even at the deeper depths in the KC area (Fig. 4). The eels in the TS area did not demonstrate clear rotational behavior at depths of more than 400 m. The PSAT studies have reported that the steady swimming depths during the day were 500–800 m14, 20. Therefore, it was assumed that the rotational behavior observed in some eels was not a regular behavior during their migration. However, the rotational behavior observed in this study suggests that they surely perceive the horizontal direction of Sun’s bearing at 400 m depths at least. Generally, they exhibit DVM precisely synchronizing with sunrise and sunset and surely perceive the change in sunlight intensity at deeper depths9,10,11,12,13,14,15,16,17,18,19,20,21,22. Even though the rotational behavior were not observed below 400 m, it remains unknown whether the eels could not perceive the Sun’s bearing from the light penetrated at depth; thus, further investigation of response to underwater light is required in future.While possible negative phototaxis behaviors were observed in some eels after release around the time of solar culmination, the trajectories of ten eels during the entire period of tracking experiments implied that each eel tended to swim meridionally toward the bearing of the Sun at culmination. We observed that eels released at middle (20°–34° N) and low (12°–13°N) latitudes tended to swim southward and northward in the meridional direction, respectively (Fig. 6A, B). The tendency to move in a north–south swimming direction corresponded to whether the Sun culminated to the north or south: eels swam southward if the culmination occurred in the southern sky, but they swam northward if it occurred in the northern sky (Fig. 6). In the KC area (33°–35° N), the Sun rose in the southeast, passed celestial meridian in the southern sky, and set in the southwest (Fig. 6C). At 20° N in the summer time when the tracking study was conducted, the Sun also passed a celestial meridian in the southern sky, but rose in the northeast and set in the northwest (Fig. 6C). When Sun culmination occurred in the southern sky, the meridional swimming directions tended to be southward (Fig. 6A). However, at 12° to 13° N in the summer time, the Sun rose in the northeast, passed the celestial meridian in the northern sky, and set in the northwest (Fig. 6D). When the Sun at culmination appeared in the northern sky, the meridional swimming directions tended to be northward (Fig. 6B). Furthermore, the swimming behavior by one eel (WE4264_TS) that was released slightly south (14° 15′ N) from the latitude with the Sun passing through the zenith was also indicative of the meridional swimming traits. This eel moved in a northerly direction on the first day, but then it lost its north–south bias in swimming around 14° 30′ N, where the Sun nearly passed through the zenith (Figs. 1 and 6D). These observations imply that the eels might move toward the latitude with the Sun passing through the zenith.Figure 6Swimming trajectories of eels and solar paths in the celestial sphere viewed from east during each tracking period. Swimming trajectories of eels released at (A) 20°N in the tropical–subtropical area and the Kuroshio Current area, and (B) 12°–14°15′N in the tropical–subtropical area. Solar paths through the north (N)–south (S) axis and the zenith at the time of tracking in (C) 20°N in the tropical–subtropical area and the Kuroshio Current area, and (D) 12°–14°15′N in the tropical–subtropical area.Full size imageTheoretically, it is possible for mesopelagic animals to use solar cues for navigation at depths shallower than the asymptotic depth, below which penetrating light rays are symmetrical around the vertical axis and the polarization plane becomes horizontal. For example, the Sargasso Sea, where the two Atlantic catadromous eels spawn1, 3, has extremely transparent water70, and the major axis of radiance distribution still remains tilted in the mesopelagic zone. The angle of maximum radiance of sunlight at 475 nm was 13° at depths of 400 m when the Sun’s elevation was 60° (Fig. 7)52, 53. In highly transparent water, the asymptotic depth could be as high as 1000–1200 m, and the depths below this cannot be utilized for compass use53. Currently, it is not possible to verify whether the Sun culminating to north or south caused the meridional swimming tendencies of eels in this study. Potentially, these meridional swimming tendencies could be due to other orientation clues, such as the geomagnetic field, as discussed for temperate anguillid eels17, 45. Nevertheless, in future studies, it would be worthwhile considering solar cues as a possible candidate factor in the orientation of eels, even when under faint underwater light conditions.Figure 7Optical features of underwater sunlight. (A) Schematic diagram of sunlight penetrating the deep ocean at 90° to the solar bearing. The line of arrows indicates the major axis of the incident beam in a vertical plane perpendicular to the Sun’s bearing. Blue light (around 475 nm) reaches the lowest depths. With increasing depth, the light field alters its character into a less directed distribution and a lower energetic level through scattering and absorption processes. Penetrating light rays are symmetrical around the region below the asymptotic depth. (B) An example of spectral radiance distribution (e. g. 475 nm) at a certain depth. The radiance distribution is shown by an ellipsoid and the major axis is drawn by a line with arrow. The refracted angular deviation (a) of the major axis of underwater radiance distribution from the vertical axis equals the tilt of the electric vector (ee bar) from the horizontal axis53. When the Sun’ s elevation was 60° in the Sargasso Sea, the radiance distributions were measured at three different depths and the tilt of the electric vector were estimated to be 24° at depths of 100 m and 200 m and 13° at depth of 400 m52, 53.Full size imageGiven that eels might be able to use the Sun’s bearing at culmination to orient their meridional swimming direction, this orientation scheme could support a clockwise eel migration route following a partial subtropical gyre2, 37. Japanese eels that departed from the nursery area first transported northeastward via the strong KC. Maintaining southward swimming in the current, they eventually crossed the current and shifted to the southward migration course. When they enter the KC, movement to the left of the bearing of the Sun at culmination (i.e., south) is the typical pattern for the early migration of eels from Japan. The movements of eels observed in the KC were consistent with the expected route; however, eels released at low latitudes of the TS area often swam northward but also westward, which resulted in their traveling an unreasonable distance from the spawning area. This might be due to their behavior during early migration. In this study, eels were transported from Japan and released into the open ocean at low latitudes. They might have swum toward the expected bearing of the Sun at culmination as if they were in the north and moved to the left of the Sun’s bearing along with the North Equatorial Current, which would mimic the early migration of eels leaving Japan and moving along the KC.Among the eels tracked in this study were individuals with impaired swim bladders, yellow-phase eels in the process of hormone-treatment maturation, and silver-phase eels collected from different rivers in different years. Despite these variations, the swimming characteristics of the eels did not differ in terms of their DVM behavior16 and swimming speed. Nevertheless, confirmation of our results using samples with a uniform status in future research would be highly desirable. In this study, the tracked eel position was assumed to be identical to that of the tracking ship and the errors between these two positions could not be evaluated; thus, the positioning of tracked fish also may need to be improved in future studies. Experimental studies, such as tracking of blind, magnetically disturbed, or olfactory-blocked eels, could help obtain or eliminate alternative candidate clues and enhance our understanding of the navigational system of anguillid eels. Controlled laboratory experiments are required to directly quantify the ability of eels to perceive radiance distribution or polarization, along with any associated behaviors. In addition, the internal clock of eels required to perform celestial navigation should be investigated. Meanwhile, the results obtained from this study can enhance our knowledge of the mechanisms underlying the migratory behaviors of eels in the open ocean. More

Wilson, A. D. M. et al. Social networks in changing environments. Behav. Ecol. Sociobiol. 69, 1617–1629 (2015).
Google Scholar 
Ward, A. J. W. et al. Association patterns and shoal fidelity in the three–spined stickleback. Proc. R. Soc. Lond. Ser. B Biol. Sci. 269, 2451–2455 (2002).
Google Scholar 
Croft, D. P. et al. Assortative interactions and social networks in fish. Oecologia 143, 211–219 (2005).CAS 
PubMed 

Google Scholar 
Helfman, G. S. & Schultz, E. T. Social transmission of behavioural traditions in a coral reef fish. Anim. Behav. 32, 379–384 (1984).
Google Scholar 
Wong, M. Y. L., Buston, P. M., Munday, P. L. & Jones, G. P. The threat of punishment enforces peaceful cooperation and stabilizes queues in a coral-reef fish. Proc. R. Soc. B Biol. Sci. 274, 1093–1099 (2007).
Google Scholar 
King, A. J., Fehlmann, G., Biro, D., Ward, A. J. & Fürtbauer, I. Re-wilding collective behaviour: an ecological perspective. Trends Ecol. Evol. 33, 347–357 (2018).PubMed 

Google Scholar 
Bro-Jørgensen, J., Franks, D. W. & Meise, K. Linking behaviour to dynamics of populations and communities: application of novel approaches in behavioural ecology to conservation. Philos. Trans. R. Soc. B Biol. Sci. 374, 20190008 (2019).
Google Scholar 
Rose, G. A. Cod spawning on a migration highway in the north-west Atlantic. Nature 366, 458 (1993).
Google Scholar 
Wilson, A. D. M., Croft, D. P. & Krause, J. Social networks in elasmobranchs and teleost fishes. Fish Fish. 15, 676–689 (2014). This study reviewed the state of knowledge of the mechanisms and functions underpinning social network structure in fishes, including a discussion on methodological issues and developments in this area of research.Taborsky, M. & Wong, M. In Comparative Social Evolution (eds. Rubenstein, D. R., Abbot, P.) 354–389 (Cambridge University Press, 2017).Lusseau, D. Evidence for social role in a dolphin social network. Evol. Ecol. 21, 357–366 (2007).
Google Scholar 
Krause, J., James, R., Franks, D. W. & Croft, D. P. Animal social networks. (Oxford University Press, 2015).Smith, J. E. & Pinter‐Wollman, N. Observing the unwatchable: Integrating automated sensing, naturalistic observations and animal social network analysis in the age of big data. J. Anim. Ecol. 90, 62–75 (2021).PubMed 

Google Scholar 
Webber, Q. M. R. & Vander Wal, E. Trends and perspectives on the use of animal social network analysis in behavioural ecology: a bibliometric approach. Anim. Behav. 149, 77–87 (2019).
Google Scholar 
Aspillaga, E., Arlinghaus, R., Martorell-Barceló, M., Barcelo-Serra, M. & Alós, J. High-throughput tracking of social networks in marine fish populations. Front. Mar. Sci. 8, 794 (2021). This original and pioneering study demonstrated the use of high-resolution tracking to infer social behaviour and social structure in the marine environment.Silk, M. J., Jackson, A. L., Croft, D. P., Colhoun, K. & Bearhop, S. The consequences of unidentifiable individuals for the analysis of an animal social network. Anim. Behav. 104, 1–11 (2015).
Google Scholar 
Hughey, L. F., Hein, A. M., Strandburg-Peshkin, A. & Jensen, F. H. Challenges and solutions for studying collective animal behaviour in the wild. Philos. Trans. R. Soc. B Biol. Sci. 373, 20170005 (2018).
Google Scholar 
Hussey, N. E. et al. Aquatic animal telemetry: a panoramic window into the underwater world. Science 348, 1255642 (2015).PubMed 

Google Scholar 
Barkley, A. N. et al. A framework to estimate the likelihood of species interactions and behavioural responses using animal-borne acoustic telemetry transceivers and accelerometers. J. Anim. Ecol. 89, 146–160 (2020).PubMed 

Google Scholar 
Baktoft, H., Gjelland, K. Ø., Økland, F. & Thygesen, U. H. Positioning of aquatic animals based on time-of-arrival and random walk models using YAPS (Yet Another Positioning Solver). Sci. Rep. 7, 14294 (2017).PubMed 
PubMed Central 

Google Scholar 
Aspillaga, E. et al. Performance of a novel system for high-resolution tracking of marine fish societies. Anim. Biotelemetry 9, 1 (2021).
Google Scholar 
Jacoby, D. M. P., Papastamatiou, Y. P. & Freeman, R. Inferring animal social networks and leadership: applications for passive monitoring arrays. J. R. Soc. Interface 13, 20160676 (2016).PubMed 
PubMed Central 

Google Scholar 
Papastamatiou, Y. P., Meyer, C. G., Watanabe, Y. & Heithaus, M. in Shark Research: Emerging Technologies and Applications for the Field and Laboratory, (eds. Carrier, J. C., Heithaus, M. R., Simpfendorfer, C. A.) 83–119 (C. R. C. Press, 2018).Butcher, P. A. et al. The drone revolution of shark. Sci. A Rev. Drones 5, 8 (2021).
Google Scholar 
Hamede, R. K., Bashford, J., McCallum, H. & Jones, M. Contact networks in a wild Tasmanian devil (Sarcophilus harrisii) population: using social network analysis to reveal seasonal variability in social behaviour and its implications for transmission of devil facial tumour disease. Ecol. Lett. 12, 1147–1157 (2009).PubMed 

Google Scholar 
Sih, A., Spiegel, O., Godfrey, S., Leu, S. & Bull, C. M. Integrating social networks, animal personalities, movement ecology and parasites: a framework with examples from a lizard. Anim. Behav. 136, 195–205 (2018).
Google Scholar 
Carne, C., Semple, S., Morrogh-Bernard, H., Zuberbühler, K. & Lehmann, J. Predicting the vulnerability of great apes to disease: the role of superspreaders and their potential vaccination. PLoS ONE 8, e84642 (2013).PubMed 
PubMed Central 

Google Scholar 
Fielding, H. R. et al. Spatial and temporal variation in proximity networks of commercial dairy cattle in Great Britain. Prev. Vet. Med. 194, 105443 (2021).PubMed 
PubMed Central 

Google Scholar 
Haulsee, D. E. et al. Social network analysis reveals potential fission-fusion behavior in a shark. Sci. Rep. 6, 34087 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Merrick, M. J. & Koprowski, J. L. Should we consider individual behavior differences in applied wildlife conservation studies? Biol. Conserv. 209, 34–44 (2017).
Google Scholar 
Kressler, M. M., Gerlam, A., Spence-Jones, H. & Webster, M. M. Passive traps and sampling bias: Social effects and personality affect trap entry by sticklebacks. Ethology 127, 446–452 (2021).
Google Scholar 
Blumstein, D. T. In Social Behaviour (eds. Szekely, T., Moore, A. J., Komdeur, J.) 520–534 (Cambridge University Press, 2010).Berger-Tal, O. et al. A systematic survey of the integration of animal behavior into conservation. Conserv. Biol. 30, 744–753 (2016).PubMed 

Google Scholar 
Mucientes, G. R., Queiroz, N., Sousa, L. L., Tarroso, P. & Sims, D. W. Sexual segregation of pelagic sharks and the potential threat from fisheries. Biol. Lett. 5, 156–159 (2009).PubMed 
PubMed Central 

Google Scholar 
Mourier, J., Vercelloni, J. & Planes, S. Evidence of social communities in a spatially structured network of a free-ranging shark species. Anim. Behav. 83, 389–401 (2012).
Google Scholar 
Perryman, R. J. Y. et al. Social preferences and network structure in a population of reef manta rays. Behav. Ecol. Sociobiol. 73, 114 (2019).
Google Scholar 
He, P., Maldonado-Chaparro, A. A. & Farine, D. R. The role of habitat configuration in shaping social structure: a gap in studies of animal social complexity. Behav. Ecol. Sociobiol. 73, 9 (2019).
Google Scholar 
Mourier, J., Lédée, E. J. I. & Jacoby, D. M. P. A multilayer perspective for inferring spatial and social functioning in animal movement networks. bioRxiv https://www.biorxiv.org/content/10.1101/749085v1.full (2019).Snijders, L., Blumstein, D. T., Stanley, C. R. & Franks, D. W. Animal social network theory can help wildlife conservation. Trends Ecol. Evol. 32, 567–577 (2017). This review paper outlines how understanding of direct and indirect relationships between animals can be profitably applied by wildlife managers and conservationists.Beyer, K., Gozlan, R. E. & Copp, G. H. Social network properties within a fish assemblage invaded by non-native sunbleak Leucaspius delineatus. Ecol. Modell. 221, 2118–2122 (2010).
Google Scholar 
Hasenjager, M. J., Leadbeater, E. & Hoppitt, W. Detecting and quantifying social transmission using network-based diffusion analysis. J. Anim. Ecol. 90, 8–26 (2021).PubMed 

Google Scholar 
Fritzsche McKay, A. & Hoye, B. J. Are migratory animals superspreaders of infection? Integr. Comp. Biol. 56, 260–267 (2016).PubMed 

Google Scholar 
Albery, G. F., Kirkpatrick, L., Firth, J. A. & Bansal, S. Unifying spatial and social network analysis in disease ecology. J. Anim. Ecol. 90, 45–61 (2021).PubMed 

Google Scholar 
Salvanes, A. & Braithwaite, V. The need to understand the behaviour of fish reared for mariculture or restocking. ICES J. Mar. Sci. 63, 345–354 (2006).
Google Scholar 
Andrew, J. E., Holm, J., Kadri, S. & Huntingford, F. A. The effect of competition on the feeding efficiency and feed handling behaviour in gilthead sea bream (Sparus aurata L.) held in tanks. Aquaculture 232, 317–331 (2004).
Google Scholar 
Muñoz, L., Aspillaga, E., Palmer, M., Saraiva, J. L. & Arechavala-Lopez, P. Acoustic telemetry: a tool to monitor fish swimming behavior in sea-cage aquaculture. Front. Mar. Sci. 7, 645 (2020).
Google Scholar 
Macaulay, G., Bui, S., Oppedal, F. & Dempster, T. Challenges and benefits of applying fish behaviour to improve production and welfare in industrial aquaculture. Rev. Aquac. 13, 934–948 (2021).
Google Scholar 
Jacoby, D. M. P. et al. Social network analysis reveals the subtle impacts of tourist provisioning on the social behavior of a generalist marine apex predator. Front. Mar. Sci. 8, 1202 (2021).
Google Scholar 
Shizuka, D. & Johnson, A. E. How demographic processes shape animal social networks. Behav. Ecol. 31, 1–11 (2020).
Google Scholar 
Guerra, A. S., Kao, A. B., McCauley, D. J. & Berdahl, A. M. Fisheries-induced selection against schooling behaviour in marine fishes. Proc. R. Soc. B Biol. Sci. 287, 20201752 (2020).
Google Scholar 
Frisch, A. Sex-change and gonadal steroids in sequentially-hermaphroditic teleost fish. Rev. Fish. Biol. Fish. 14, 481–499 (2004).
Google Scholar 
Webber, Q. M. R. & Vander Wal, E. An evolutionary framework outlining the integration of individual social and spatial ecology. J. Anim. Ecol. 87, 113–127 (2018).PubMed 

Google Scholar 
Staveley, T. A. B. et al. Sea surface temperature dictates movement and habitat connectivity of Atlantic cod in a coastal fjord system. Ecol. Evol. 9, 9076–9086 (2019).PubMed 
PubMed Central 

Google Scholar 
Sosa, S., Jacoby, D. M. P., Lihoreau, M. & Sueur, C. Animal social networks: towards an integrative framework embedding social interactions, space and time. Methods Ecol. Evol. 12, 4–9 (2021).
Google Scholar 
Albery, G. F. et al. Multiple spatial behaviours govern social network positions in a wild ungulate. Ecol. Lett. 24, 676–686 (2021).PubMed 

Google Scholar 
Ellis, S. et al. Mortality risk and social network position in resident killer whales: sex differences and the importance of resource abundance. Proc. R. Soc. B Biol. Sci. 284, 20171313 (2017).
Google Scholar 
Ellis, S., Snyder-Mackler, N., Ruiz-Lambides, A., Platt, M. L. & Brent, L. J. N. Deconstructing sociality: the types of social connections that predict longevity in a group-living primate. Proc. R. Soc. B Biol. Sci. 286, 20191991 (2019).
Google Scholar 
Kohn, G. M. Friends give benefits: autumn social familiarity preferences predict reproductive output. Anim. Behav. 132, 201–208 (2017).
Google Scholar 
Villegas-Ríos, D., Freitas, C., Moland, E., Thorbjørnsen, S. H. & Olsen, E. M. Inferring individual fate from aquatic acoustic telemetry data. Methods Ecol. Evol. 11, 1186–1198 (2020).
Google Scholar 
Mourier, J., Bass, N. C., Guttridge, T. L., Day, J. & Brown, C. Does detection range matter for inferring social networks in a benthic shark using acoustic telemetry? R. Soc. open Sci. 4, 170485 (2017).PubMed 
PubMed Central 

Google Scholar 
Vanovac, S., Howard, D., Monk, C. T., Arlinghaus, R. & Giabbanelli, P. J. Network analysis of intra- and interspecific freshwater fish interactions using year-around tracking. J. R. Soc. Interface 18, 20210445 (2021).PubMed 

Google Scholar 
Dahl, K. A., Patterson, W. F. & Snyder, R. A. Experimental assessment of lionfish removals to mitigate reef fish community shifts on northern Gulf of Mexico artificial reefs. Mar. Ecol. Prog. Ser. 558, 207–221 (2016).
Google Scholar 
Fitzpatrick, J. L. et al. Female-mediated causes and consequences of status change in a social fish. Proc. R. Soc. B Biol. Sci. 275, 929–936 (2008).CAS 

Google Scholar 
Mourier, J., Brown, C. & Planes, S. Learning and robustness to catch-and-release fishing in a shark social network. Biol. Lett. 13, 20160824 (2017).PubMed 
PubMed Central 

Google Scholar 
Rutledge, L. Y. et al. Protection from harvesting restores the natural social structure of eastern wolf packs. Biol. Conserv. 143, 332–339 (2010).
Google Scholar 
Jacoby, D. M. P. et al. Synergistic patterns of threat and the challenges facing global anguillid eel conservation. Glob. Ecol. Conserv 4, 321–333 (2015).
Google Scholar 
Geffroy, B., Bru, N., Dossou-Gbété, S., Tentelier, C. & Bardonnet, A. The link between social network density and rank-order consistency of aggressiveness in juvenile eels. Behav. Ecol. Sociobiol. 68, 1073–1083 (2014).
Google Scholar  More

Lombard, M. et al. South African and Lesotho Stone Age sequence updated. S. Afr. Archaeol. Bull. 67, 120–144 (2012).
Google Scholar 
Kandel, A. W. et al. Increasing behavioral flexibility? An integrative macro-scale approach to understanding the Middle Stone Age of southern Africa. J. Archaeol. Method Theory 23, 623–628 (2015).
Google Scholar 
Porraz, G. et al. Experimentation preceding innovation in a MIS5 Pre-Still Bay layer from Diepkloof Rock Shelter (South Africa): emerging technologies and symbols. Preprint at EcoEvoRxiv https://ecoevorxiv.org/ch53r/ (2020).Texier, P. J. et al. A Howiesons Poort tradition of engraving ostrich eggshell containers dated to 60,000 years ago at Diepkloof Rock Shelter, South Africa. Proc. Natl Acad. Sci. USA 107, 6180–6185 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
Henshilwood, C. S. et al. Klipdrift Shelter, southern Cape, South Africa: preliminary report on the Howiesons Poort layers. J. Archaeol. Sci. 45, 284–303 (2014).
Google Scholar 
Powell, A., Shennan, S. & Thomas, M. G. Late Pleistocene demography and the appearance of modern human behavior. Science 324, 1298–1301 (2009).CAS 
PubMed 

Google Scholar 
Marean, C. W. The transition to foraging for dense and predictable resources and its impact on the evolution of modern humans. Phil. Trans. R. Soc. B 371, 20150239 (2016).PubMed 
PubMed Central 

Google Scholar 
Mackay, A., Stewart, B. A. & Chase, B. M. Coalescence and fragmentation in the late Pleistocene archaeology of southernmost Africa. J. Hum. Evol. 72, 26–51 (2014).PubMed 

Google Scholar 
Wilkins, J. et al. Innovative Homo sapiens behaviours 105,000 years ago in a wetter Kalahari. Nature 592, 248–252 (2021).CAS 
PubMed 

Google Scholar 
Dewar, G. & Stewart, B. A. Preliminary results of excavations at Spitzkloof Rockshelter, Richtersveld, South Africa. Quat. Int. 270, 30–39 (2012).
Google Scholar 
Cowling, R. M. & Pierce, S. Namaqualand: A Succulent Desert (Fernwood Press, 1999).Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G. & Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978 (2005).
Google Scholar 
Mucina, L. et al. in The Vegetation of South Africa, Lesotho and Swaziland (eds Mucina, L. & Rutherford, M. C.) 221–299 (SANBI, 2006).Rebelo, A. G., Boucher, C., Helme, N., Mucina, L. & Rutherford, M. C. in The Vegetation of South Africa, Lesotho and Swaziland (eds Mucina, L. & Rutherford, M. C.) 53–219 (SANBI, 2006).Marean, C. W. et al. in Fynbos: Ecology, Evolution, and Conservation of a Megadiverse Region (eds Allsopp, N. et al.) 164–199 (Oxford Univ. Press, 2014).Carr, A. S., Chase, B. M. & Mackay, A. in Africa from MIS 6-2: Population Dynamics and Paleoenvironments (eds Jones, S. & Stewart, B. A.) 23–47 (Springer, 2016).Chase, B. M. & Meadows, M. E. Late Quaternary dynamics of southern Africa’s winter rainfall zone. Earth Sci. Rev. 84, 103–138 (2007).
Google Scholar 
Steele, T. E. et al. Varsche Rivier 003: a Middle and Later Stone Age site with Still Bay and Howiesons Poort assemblages in southern Namaqualand, South Africa. Paleoanthropology 2016, 100–163 (2016).
Google Scholar 
Sharp, W. D. et al. 230Th/U burial dating of ostrich eggshell. Quat. Sci. Rev. 219, 263–276 (2019).
Google Scholar 
Chase, B. M. et al. South African speleothems reveal influence of high- and low-latitude forcing over the last 113.5 kyr. Geology 49, 1353–1357 (2021).CAS 

Google Scholar 
Chase, B. M. et al. Influence of tropical easterlies in southern Africa’s winter rainfall zone during the Holocene. Quat. Sci. Rev. 107, 138–148 (2015).
Google Scholar 
Manning, J. Namaqualand (Briza, 2008).Skinner, J. D. & Chimimba, C. T. The Mammals of the Southern African Subregion 3rd edn (Cambridge Univ. Press, 2005).Skead, C. J. Historical Mammal Incidence in the Cape Province Vol 1: The Western and Northern Cape (Cape Town Department of Nature and Environmental Conservation, 1980).Churcher, C. S. Distribution and history of the Cape zebra (Equus capensis) in the Quaternary of Africa. Trans. R. Soc. S. Afr. 61, 89–95 (2006).
Google Scholar 
Spratt, R. M. & Lisiecki, L. E. A Late Pleistocene sea level stack. Climate 12, 1079–1092 (2016).
Google Scholar 
De Wet, W. Bathymetry of the South African Continental Shelf. MSc thesis, Univ. Cape Town (2013).Jerardino, A. & Marean, C. W. Shellfish gathering, marine paleoecology and modern human behavior: perspectives from Cave PP13B, Pinnacle Point, South Africa. J. Hum. Evol. 59, 412–424 (2010).PubMed 

Google Scholar 
Marean, C. W. Pinnacle Point Cave 13B (Western Cape Province, South Africa) in context: the Cape Floral kingdom, shellfish, and modern human origins. J. Hum. Evol. 59, 425–443 (2010).PubMed 

Google Scholar 
Kandel, A. W. Modification of ostrich eggs by carnivores and its bearing on the interpretation of archaeological and paleontological find. J. Archaeol. Sci. 31, 377–391 (2004).
Google Scholar 
Steele, T. E. & Klein, R. G. The Middle and Later Stone Age faunal remains from Diepkloof Rock Shelter, Western Cape, South Africa. J. Archaeol. Sci. 40, 3453–3462 (2013).
Google Scholar 
Klein, R. G. et al. The Ysterfontein 1 Middle Stone Age site, South Africa, and early human exploitation of coastal resources. Proc. Natl Acad. Sci. USA 101, 5708–5715 (2004).CAS 
PubMed 
PubMed Central 

Google Scholar 
Vogelsang, R. et al. New excavations of Middle Stone Age deposits at Apollo 11 Rockshelter, Namibia: stratigraphy, archaeology, chronology and past environments. J. Afr. Archaeol. 8, 185–218 (2010).Marean, C. W. et al. Early human use of marine resources and pigment in South Africa during the Middle Pleistocene. Nature 449, 905–908 (2007).CAS 
PubMed 

Google Scholar 
Schmidt, I. et al. New investigations at the Middle Stone Age site of Pockenbank Rockshelter, Namibia. Antiquity 90, e2 (2016).
Google Scholar 
Vogelsang, R. Middle Stone Age Fundstellen in Südwest-Namibia, Africa (Heinrich-Barth-Institut, 1998).Plug, I. Aquatic animals and their associates from the Middle Stone Age levels at Sibudu. South. Afr. Humanit. 18, 289–299 (2006).
Google Scholar 
Wurz, S. Technological trends in the Middle Stone Age of South Africa between MIS 7 and MIS 3. Curr. Anthropol. 54, S305–S319 (2013).
Google Scholar 
Volman, T. P. The Middle Stone Age in the Southern Cape. PhD thesis, Univ. Chicago (1981).Schmidt, P. & Mackay, A. Why was silcrete heat-treated in the Middle Stone Age? An early transformative technology in the context of raw material use at Mertenhof Rock Shelter, South Africa. PloS ONE 11, e0149243 (2016).PubMed 
PubMed Central 

Google Scholar 
Porraz, G. et al. Technological successions in the Middle Stone Age sequence of Diepkloof Rock Shelter, Western Cape, South Africa. J. Archaeol. Sci. 40, 3376–3400 (2013).
Google Scholar 
Schmid, V., Conard, N. J., Parkington, J., Texier, P. J. & Porraz, G. The ‘MSA 1’ of Elands Bay Cave (South Africa) in the context of the southern African early MSA technologies. South. Afr. Humanit. 29, 153–201 (2016).
Google Scholar 
Evans, U. Hollow Rock Shelter, a Middle Stone Age site in the Cederberg. South. Afr. Field Archaeol. 3, 63–73 (1994).
Google Scholar 
Mackay, A., Jacobs, Z. & Steele, T. E. Pleistocene archaeology and chronology of Putslaagte 8 (PL8) rockshelter, Western Cape, South Africa. J. Afr. Archaeol. 13, 71–98 (2015).
Google Scholar 
Thompson, J. C. et al. Ecological risk, demography and technological complexity in the Late Pleistocene of northern Malawi: implications for geographical patterning in the Middle Stone Age. J. Quat. Sci. 33, 261–284 (2018).
Google Scholar 
Vaesen, K. & Houkes, W. Is human culture cumulative? Curr. Anthropol. 62, 218–238 (2021).
Google Scholar 
Boyd, R., Richerson, P. J. & Henrich, J. The cultural niche: why social learning is essential for human adaptation. Proc. Natl Acad. Sci. USA 108, 10918–10925 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
Sterelny, K. From hominins to humans: how sapiens became behaviourally modern. Phil. Trans. R. Soc. B 366, 809–822 (2011).PubMed 
PubMed Central 

Google Scholar 
Gärdenfors, P. & Högberg, A. The archaeology of teaching and the evolution of Homo docens. Curr. Anthropol. 58, 188–208 (2017).
Google Scholar 
Marwick, B. Pleistocene exchange networks as evidence for the evolution of language. Camb. Archaeol. J. 13, 67–81 (2003).
Google Scholar 
Blegen, N. The earliest long-distance obsidian transport: evidence from the ∼200 ka Middle Stone Age Sibilo School Road Site, Baringo, Kenya. J. Hum. Evol. 103, 1–19 (2017).PubMed 

Google Scholar 
McBrearty, S. & Brooks, A. S. The revolution that wasn’t: a new interpretation of the origin of modern human behavior. J. Hum. Evol. 39, 453–563 (2000).CAS 
PubMed 

Google Scholar 
Klein, R. G. Archeology and the evolution of human behavior. Evol. Anthropol. 9, 17–36 (2000).
Google Scholar 
Wynn, T. & Coolidge, F. L. Archeological insights into hominin cognitive evolution. Evol. Anthropol. 25, 200–213 (2016).PubMed 

Google Scholar 
Derex, M. & Mesoudi, A. Cumulative cultural evolution within evolving population structures. Trends Cogn. Sci. 24, 654–667 (2020).PubMed 

Google Scholar 
Sterelny, K. Adaptable individuals and innovative lineages. Phil. Trans. R. Soc. B 371, 20150196 (2016).PubMed 
PubMed Central 

Google Scholar 
Henshilwood, C. S. & Marean, C. W. The origin of modern human behavior: critique of the models and their test implications. Curr. Anthropol. 44, 627–651 (2003).PubMed 

Google Scholar 
Stoops, G., Marcelino, V. & Mees, F. (eds) Interpretation of Micromorphological Features of Soils and Regoliths (Elsevier, 2010).Stoops, G. Guidelines for Analysis and Description of Soil and Regolith Thin Sections (Soil Science Society of America, 2003).McPherron, S. P. Additional statistical and graphical methods for analyzing site formation processes using artifact orientations. PLoS ONE 13, e0190195 (2018).PubMed 
PubMed Central 

Google Scholar 
Thomsen, K. J., Murray, A. S., Jain, M. & Bøtter-Jensen, L. Laboratory fading rates of various luminescence signals from feldspar-rich sediment extracts. Radiat. Meas. 43, 1474–1486 (2008).CAS 

Google Scholar 
Armitage, S. J. & Bailey, R. M. The measured dependence of laboratory beta dose rates on sample grain size. Radiat. Meas. 39, 123–127 (2005).CAS 

Google Scholar 
Huntley, D. J. & Lamothe, M. Ubiquity of anomalous fading in K-feldspars and the measurement and correction for it in optical dating. Can. J. Earth Sci. 38, 1093–1106 (2001).CAS 

Google Scholar 
Jaffey, A. H., Flynn, K. F., Glendenin, L. E. & Essling, A. M. Precision measurement of half-lives and specific activities of 235U and 238U. Phys. Rev. C 4, 1889–1906 (1971).
Google Scholar 
Cheng, H. et al. Improvements in 230Th dating, 230Th and 234U half-life values, and U–Th isotopic measurements by multi-collector inductively coupled plasma mass spectrometry. Earth Planet. Sci. Lett. 371–372, 82–91 (2013).
Google Scholar 
Holden, N. E. Total half-lives for selected nuclides. Pure Appl. Chem. 62, 941–958 (1990).CAS 

Google Scholar 
Ludwig, K. R. Isoplot/Ex Version 3.75: A Geochronological Toolkit for Microsoft Excel (Berkeley Geochronology Center Special Publication, 2010).Collins, B. & Steele, T. E. An often overlooked resource: ostrich (Struthio spp.) eggshell in the archaeological record. J. Archaeol. Sci. Rep. 13, 121–131 (2017).
Google Scholar 
Schmidt, P. How reliable is the visual identification of heat treatment on silcrete? A quantitative verification with a new method. Archaeol. Anthropol. Sci. 11, 713–726 (2017).
Google Scholar 
Roberts, D. L. Age, Genesis and Significance of South African Coastal Belt Silcretes (Council for Geoscience, South Africa, 2003).Schmidt, P. et al. A previously undescribed organic residue sheds light on heat treatment in the Middle Stone Age. J. Hum. Evol. 85, 22–34 (2015).PubMed 

Google Scholar 
Trabucco, A. & Zomer, R. Global Aridity Index and Potential Evapotranspiration (ET0) Climate Database v2 (figshare, 2019); https://doi.org/10.6084/m9.figshare.7504448.v3World Atlas of Desertification 2nd edn (UNEP, 1997).Mucina, L. & Rutherford, M. C. The Vegetation of South Africa, Lesotho and Swaziland (South African National Biodiversity Institute, 2006).Cordova, C. E. C3 Poaceae and Restionaceae phytoliths as potential proxies for reconstructing winter rainfall in South Africa. Quat. Int. 287, 121–140 (2013).
Google Scholar 
Esteban, I. et al. Modern soil phytolith assemblages used as proxies for paleoscape reconstruction on the south coast of South Africa. Quat. Int. 434, 160–179 (2017).
Google Scholar 
Laskar, J. et al. A long-term numerical solution for the insolation quantities of the Earth. Astron. Astrophys. 428, 261–285 (2004).
Google Scholar 
Chase, B. M. et al. Orbital controls on Namib Desert hydroclimate over the past 50,000 years. Geology 47, 867–871 (2019).
Google Scholar 
Farmer, E. C., deMenocal, P. B. & Marchitto, T. M. Holocene and deglacial ocean temperature variability in the Benguela upwelling region: implications for low‐latitude atmospheric circulation. Paleoceanography 20, PA2018 (2005).
Google Scholar 
Pichevin, L., Cremer, M., Giraudeau, J. & Bertrand, P. A 190 kyr record of lithogenic grain size on the Namibian slope: forging a tight link between past wind‐strength and coastal upwelling dynamics. Mar. Geol. 218, 81–96 (2005).
Google Scholar 
Little, M. G. et al. Trade wind forcing of upwelling, seasonality, and Heinrich events as a response to sub‐Milankovitch climate variability. Paleoceanography 12, 568–576 (2005).
Google Scholar 
Stuut, J.-B. et al. A 300‐kyr record of aridity and wind strength in southwestern Africa: inferences from grain‐size distributions of sediments on Walvis Ridge, SE Atlantic. Mar. Geol. 180, 221–233 (2002).
Google Scholar 
Kandel, A. W. & Conard, N. J. Production sequences of ostrich eggshell beads and settlement dynamics in the Geelbek Dunes of the Western Cape, South Africa. J. Archaeol. Sci. 32, 1711–1721 (2005).
Google Scholar  More

Trenberth, K. E. & Jones, P. D. in Climate Change 2007: The Physical Science Basis (eds Solomon, S. et al.) 235–335 (Cambridge Univ. Press, 2007).Linderholm, H. W. Growing season changes in the last century. Agr. For. Meteorol. 137, 1–14 (2006).
Google Scholar 
Yang, B. et al. New perspective on spring vegetation phenology and global climate change based on Tibetan Plateau tree-ring data. Proc. Natl Acad. Sci. USA 114, 6966–6971 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
Shen, M., Tang, Y., Chen, J. & Yang, W. Specification of thermal growing season in temperate China from 1960 to 2009. Clim. Change 114, 783–798 (2012).
Google Scholar 
Zhou, B., Zhai, P., Chen, Y. & Yu, R. Projected changes of thermal growing season over Northern Eurasia in a 1.5 °C and 2 °C warming world. Environ. Res. Lett. 13, 35004 (2018).
Google Scholar 
Barichivich, J., Briffa, K. R., Osborn, T. J., Melvin, T. M. & Caesar, J. Thermal growing season and timing of biospheric carbon uptake across the Northern Hemisphere. Glob. Biogeochem. Cycles 26, B4015 (2012).
Google Scholar 
Buitenwerf, R., Rose, L. & Higgins, S. I. Three decades of multi-dimensional change in global leaf phenology. Nat. Clim. Change 5, 364–368 (2015).
Google Scholar 
Gonsamo, A., Chen, J. M. & Ooi, Y. W. Peak season plant activity shift towards spring is reflected by increasing carbon uptake by extratropical ecosystems. Glob. Change Biol. 24, 2117–2128 (2018).
Google Scholar 
Menzel, A. et al. European phenological response to climate change matches the warming pattern. Glob. Change Biol. 12, 1969–1976 (2006).
Google Scholar 
Montgomery, R. A., Rice, K. E., Stefanski, A., Rich, R. L. & Reich, P. B. Phenological responses of temperate and boreal trees to warming depend on ambient spring temperatures, leaf habit, and geographic range. Proc. Natl Acad. Sci. USA 117, 10397–10405 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Piao, S. et al. Leaf onset in the northern hemisphere triggered by daytime temperature. Nat. Commun. 6, 6911 (2015).CAS 
PubMed 

Google Scholar 
Barichivich, J. et al. Large-scale variations in the vegetation growing season and annual cycle of atmospheric CO2 at high northern latitudes from 1950 to 2011. Glob. Change Biol. 19, 3167–3183 (2013).
Google Scholar 
Peñuelas, J., Rutishauser, T. & Filella, I. Phenology feedbacks on climate change. Science 324, 887–888 (2009).PubMed 

Google Scholar 
Piao, S. et al. Weakening temperature control on the interannual variations of spring carbon uptake across northern lands. Nat. Clim. Change 7, 359–363 (2017).CAS 

Google Scholar 
Richardson, A. D. et al. Influence of spring and autumn phenological transitions on forest ecosystem productivity. Philos. Trans. R. Soc. B 365, 3227–3246 (2010).
Google Scholar 
Bonan, G. B. Forests and climate change: forcings, feedbacks, and the climate benefits of forests. Science 320, 1444–1449 (2008).CAS 
PubMed 

Google Scholar 
Piao, S. et al. Plant phenology and global climate change: current progresses and challenges. Glob. Change Biol. 25, 1922–1940 (2019).
Google Scholar 
Fu, Y. H. et al. Declining global warming effects on the phenology of spring leaf unfolding. Nature 526, 104–107 (2015).CAS 
PubMed 

Google Scholar 
Park, T. et al. Changes in timing of seasonal peak photosynthetic activity in northern ecosystems. Glob. Change Biol. 25, 2382–2395 (2019).
Google Scholar 
Xu, C., Liu, H., Williams, A. P., Yin, Y. & Wu, X. Trends toward an earlier peak of the growing season in Northern Hemisphere mid-latitudes. Glob. Change Biol. 22, 2852–2860 (2016).
Google Scholar 
Wang, X. et al. No trends in spring and autumn phenology during the global warming hiatus. Nat. Commun. 10, 2389 (2019).PubMed 
PubMed Central 

Google Scholar 
Piao, S., Friedlingstein, P., Ciais, P., Viovy, N. & Demarty, J. Growing season extension and its impact on terrestrial carbon cycle in the Northern Hemisphere over the past 2 decades. Glob. Biogeochem. Cycles 21, B3018 (2007).
Google Scholar 
Buermann, W., Bikash, P. R., Jung, M., Burn, D. H. & Reichstein, M. Earlier springs decrease peak summer productivity in North American boreal forests. Environ. Res. Lett. 8, 24027 (2013).
Google Scholar 
Buermann, W. et al. Widespread seasonal compensation effects of spring warming on northern plant productivity. Nature 562, 110–114 (2018).CAS 
PubMed 

Google Scholar 
Lian, X. et al. Summer soil drying exacerbated by earlier spring greening of northern vegetation. Sci. Adv. 6, eaax0255 (2020).PubMed 
PubMed Central 

Google Scholar 
Piao, S. et al. Net carbon dioxide losses of northern ecosystems in response to autumn warming. Nature 451, 49–52 (2008).CAS 
PubMed 

Google Scholar 
Wang, H. et al. Alpine grassland plants grow earlier and faster but biomass remains unchanged over 35 years of climate change. Ecol. Lett. 23, 701–710 (2020).PubMed 
PubMed Central 

Google Scholar 
Pan, Y. et al. A large and persistent carbon sink in the world’s forests. Science 333, 988–993 (2011).CAS 
PubMed 

Google Scholar 
Delpierre, N. et al. Temperate and boreal forest tree phenology: from organ-scale processes to terrestrial ecosystem models. Ann. For. Sci. 73, 5–25 (2016).
Google Scholar 
Huang, J. et al. Photoperiod and temperature as dominant environmental drivers triggering secondary growth resumption in Northern Hemisphere conifers. Proc. Natl Acad. Sci. USA 117, 20645–20652 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Li, X. et al. Critical minimum temperature limits xylogenesis and maintains treelines on the southeastern Tibetan Plateau. Sci. Bull. 62, 804–812 (2017).
Google Scholar 
Rossi, S. et al. Critical temperatures for xylogenesis in conifers of cold climates. Glob. Ecol. Biogeogr. 17, 696–707 (2008).
Google Scholar 
Lenz, A., Vitasse, Y., Hoch, G. & Körner, C. Growth and carbon relations of temperate deciduous tree species at their upper elevation range limit. J. Ecol. 102, 1537–1548 (2014).
Google Scholar 
Zeng, Q., Rossi, S., Yang, B., Qin, C. & Li, G. Environmental drivers for cambial reactivation of Qilian junipers (Juniperus przewalskii) in a semi-arid region of northwestern China. Atmosphere 11, 232 (2020).
Google Scholar 
Ren, P. et al. Growth rate rather than growing season length determines wood biomass in dry environments. Agr. For. Meteorol. 271, 46–53 (2019).
Google Scholar 
Sanginés De Cárcer, P. et al. Vapor-pressure deficit and extreme climatic variables limit tree growth. Glob. Change Biol. 24, 1108–1122 (2017).
Google Scholar 
Zhang, J. et al. Drought limits wood production of Juniperus przewalskii even as growing seasons lengthens in a cold and arid environment. Catena 196, 104936 (2021).
Google Scholar 
Huang, J., Deslauriers, A. & Rossi, S. Xylem formation can be modeled statistically as a function of primary growth and cambium activity. New Phytol. 203, 831–841 (2014).CAS 
PubMed 

Google Scholar 
Rossi, S., Morin, H. & Deslauriers, A. Causes and correlations in cambium phenology: towards an integrated framework of xylogenesis. J. Exp. Bot. 63, 2117–2126 (2012).CAS 
PubMed 

Google Scholar 
Rossi, S., Girard, M. J. & Morin, H. Lengthening of the duration of xylogenesis engenders disproportionate increases in xylem production. Glob. Change Biol. 20, 2261–2271 (2014).
Google Scholar 
Cuny, H. E. et al. Woody biomass production lags stem-girth increase by over one month in coniferous forests. Nat. Plants 1, 15160 (2015).CAS 
PubMed 

Google Scholar 
Pasho, E., Camarero, J. J. & Vicente-Serrano, S. M. Climatic impacts and drought control of radial growth and seasonal wood formation in Pinus halepensis. Trees 26, 1875–1886 (2012).
Google Scholar 
Keenan, T. F. et al. Net carbon uptake has increased through warming-induced changes in temperate forest phenology. Nat. Clim. Change 4, 598–604 (2014).CAS 

Google Scholar 
Chen, L. et al. Leaf senescence exhibits stronger climatic responses during warm than during cold autumns. Nat. Clim. Change 10, 777–780 (2020).CAS 

Google Scholar 
Körner, C. Paradigm shift in plant growth control. Curr. Opin. Plant Biol. 25, 107–114 (2015).PubMed 

Google Scholar 
Muller, B. et al. Water deficits uncouple growth from photosynthesis, increase C content, and modify the relationships between C and growth in sink organs. J. Exp. Bot. 62, 1715–1729 (2011).CAS 
PubMed 

Google Scholar 
Charney, N. D. et al. Observed forest sensitivity to climate implies large changes in 21st century North American forest growth. Ecol. Lett. 19, 1119–1128 (2016).PubMed 

Google Scholar 
Liu, Q. et al. Extension of the growing season increases vegetation exposure to frost. Nat. Commun. 9, 426 (2018).PubMed 
PubMed Central 

Google Scholar 
Deslauriers, A. & Morin, H. Intra-annual tracheid production in balsam fir stems and the effect of meteorological variables. Trees 19, 402–408 (2005).
Google Scholar 
Piao, S. et al. Characteristics, drivers and feedbacks of global greening. Nat. Rev. Earth Environ. 1, 14–27 (2020).
Google Scholar 
Huang, M. et al. Air temperature optima of vegetation productivity across global biomes. Nat. Ecol. Evol. 3, 772–779 (2019).PubMed 
PubMed Central 

Google Scholar 
Keenan, T. F. & Riley, W. J. Greening of the land surface in the world’s cold regions consistent with recent warming. Nat. Clim. Change 8, 825–828 (2018).CAS 

Google Scholar 
Camarero, J. J., Olano, J. M. & Parras, A. Plastic bimodal xylogenesis in conifers from continental Mediterranean climates. New Phytol. 185, 471–480 (2010).PubMed 

Google Scholar 
Fu, Y. H. et al. Unexpected role of winter precipitation in determining heat requirement for spring vegetation green-up at northern middle and high latitudes. Glob. Change Biol. 20, 3743–3755 (2014).
Google Scholar 
Wu, X. et al. Uneven winter snow influence on tree growth across temperate China. Glob. Change Biol. 25, 144–154 (2018).
Google Scholar 
Wang, X. et al. Disentangling the mechanisms behind winter snow impact on vegetation activity in northern ecosystems. Glob. Change Biol. 24, 1651–1662 (2018).
Google Scholar 
Adams, H. D. et al. Experimental drought and heat can delay phenological development and reduce foliar and shoot growth in semiarid trees. Glob. Change Biol. 21, 4210–4220 (2015).
Google Scholar 
He, W., Liu, H., Qi, Y., Liu, F. & Zhu, X. Patterns in nonstructural carbohydrate contents at the tree organ level in response to drought duration. Glob. Change Biol. 26, 3627–3638 (2020).
Google Scholar 
Williams, A. P. et al. Temperature as a potent driver of regional forest drought stress and tree mortality. Nat. Clim. Change 3, 292–297 (2012).
Google Scholar 
Vitasse, Y. et al. Contrasting resistance and resilience to extreme drought and late spring frost in five major European tree species. Glob. Change Biol. 25, 3781–3792 (2019).
Google Scholar 
Zhao, S. et al. The International Tree-Ring Data Bank (ITRDB) revisited: data availability and global ecological representativity. J. Biogeogr. 46, 355–368 (2019).
Google Scholar 
Babst, F., Poulter, B., Bodesheim, P., Mahecha, M. D. & Frank, D. C. Improved tree-ring archives will support earth-system science. Nat. Ecol. Evol. 1, 8 (2017).PubMed 

Google Scholar 
Elmore, A. J., Guinn, S. M., Minsley, B. J. & Richardson, A. D. Landscape controls on the timing of spring, autumn, and growing season length in mid-Atlantic forests. Glob. Change Biol. 18, 656–674 (2012).
Google Scholar 
Kannenberg, S. A. et al. Drought legacies are dependent on water table depth, wood anatomy and drought timing across the eastern US. Ecol. Lett. 22, 119–127 (2018).PubMed 

Google Scholar 
Rossi, S., Deslauriers, A., Anfodillo, T. & Carraro, V. Evidence of threshold temperatures for xylogenesis in conifers at high altitudes. Oecologia 152, 1–12 (2007).PubMed 

Google Scholar 
Gao, S. et al. Dynamic responses of tree-ring growth to multiple dimensions of drought. Glob. Change Biol. 24, 5380–5390 (2018).
Google Scholar 
Peltier, D. M. P. & Ogle, K. Tree growth sensitivity to climate is temporally variable. Ecol. Lett. 23, 1561–1572 (2020).PubMed 

Google Scholar 
Wilmking, M. et al. Global assessment of relationships between climate and tree growth. Glob. Change Biol. 26, 3212–3220 (2020).
Google Scholar 
Seftigen, K., Frank, D. C., Björklund, J., Babst, F. & Poulter, B. The climatic drivers of normalized difference vegetation index and tree-ring-based estimates of forest productivity are spatially coherent but temporally decoupled in Northern Hemispheric forests. Glob. Ecol. Biogeogr. 27, 1352–1365 (2018).
Google Scholar 
Bunn, A. G. A dendrochronology program library in R (dplR). Dendrochronologia 26, 115–124 (2008).
Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2019); https://www.R-project.org/Sheffield, J., Goteti, G. & Wood, E. F. Development of a 50-year high-resolution global dataset of meteorological forcings for land surface modeling. J. Clim. 19, 3088–3111 (2006).
Google Scholar 
Frich, P. L. et al. Observed coherent changes in climatic extremes during the second half of the twentieth century. Clim. Res. 19, 193–212 (2002).
Google Scholar 
Selyaninov, G. T. About climate agricultural estimation (in Russian). Proc. Agric. Meteorol. 20, 165–177 (1928).
Google Scholar 
Streiner, D. L. Finding our way: an introduction to path analysis. Can. J. Psychiatry 50, 115–122 (2005).PubMed 

Google Scholar 
Fox, J., Nie, Z. & Byrnes, J. sem: Structural equation models. R package version 3.1-9 https://CRAN.R-project.org/package=sem (2017).Iturbide, M. et al. An update of IPCC climate reference regions for subcontinental analysis of climate model data: definition and aggregated datasets. Earth Syst. Sci. Data 12, 2959–2970 (2020).
Google Scholar 
Bagozzi, R. P. & Yi, Y. Specification, evaluation, and interpretation of structural equation models. J. Acad. Mark. Sci. 40, 8–34 (2012).
Google Scholar  More

Global map of forest cover change and its validationWe use a high-resolution map of global forest cover change15 (annual intervals at 30 m spatial resolution; version 1.7) to quantify forest loss across the tropics. The GFC dataset maps where and when forests were converted (naturally and anthropogenically) from 2001 to 2019. Trees are defined as all vegetation taller than 5 m in height, and forests are defined with a tree canopy threshold of at least 30%. The definition of forests includes plantations and tree crops such as oil palm. Forest loss is the mortality or removal of all tree cover within a pixel.Our study uses the v.1.7 product that spans the period 2001–2019 for the analysis. Different methods are used for detecting forest cover loss in two periods (2001–2010 and 2011–2019). This change in detection method as well as in satellite data (Landsat 7 and Landsat 8) might result in inconsistencies of data during the two periods. Therefore, we perform an independent assessment of the v.1.7 product throughout the study period (2001–2019) using stratified random-sample reference data. We randomly sample 18,000 pixels, a much larger sample population than the assessment of the original product (v.1.0, 628 pixels across the tropics; ref. 15), and visually interpret forest loss using Landsat imagery. Specifically, we randomly select 50 path/row locations (World Reference System II) of Landsat imagery in each tropical continent (150 path/row locations in total in the three tropical continents; Supplementary Fig. 6). For each path/row location, we randomly select 20 loss pixels and 10 non-loss pixels in each period of 2001–2005, 2006–2010, 2011–2014 and 2015–2019, with total sampling pixels of ~18,000 ((20 loss + 10 non-loss) × 50 path/row locations × 4 periods × 3 continents). Our study compares the increase in forest carbon loss from the start (~2001) to the end (~2019) of the study period. Thus, we divided the whole 19 yr study period into four subperiods (2001–2005, 2006–2010, 2011–2014 and 2015–2019), with the first five years considered as the start period and the last five years considered as the end period for the comparison. Some path/row locations do not have 20 loss pixels in a specific period; for example, there is no loss detected in some locations of the Sahara. Therefore, we sample 11,198 loss pixels and 6,000 non-loss pixels (Supplementary Data). Finally, we download time-series Landsat imagery covering 1999–2020 to visually interpret these pixels as reference data.Following the suggestion of Global Forest Watch19 and per best practice guidance of ref. 18, we use a stratified random-sample approach for area estimation, which is independent of the method and satellite changes in the GFC data. The sampling reference data (Supplementary Data) are used to estimate loss area:$${it{p}}_{hi} = w_hfrac{{mathop {sum}limits_{j = 1}^{150} {n_{hij} times A_{hi}} }}{{mathop {sum}limits_{j = 1}^{150} {n_{hj} times A_i} }}$$
(1)
where phi is the stratified random-sample estimated area for GFC map class h that is classified as reference class i; wh is the proportion of the total area GFC map class h; nhij is the number of pixels in GFC map class h that is classified as reference class i in jth Landsat path/row location; nhj is the total number of pixels in GFC map class h in jth Landsat path/row location; Ahj is the total area in GFC map class h in jth Landsat path/row location; and Aj is the total land area of jth Landsat path/row location.We then calculate the error matrix, which includes overall accuracy (OA), user’s accuracy (UA) and producer’s accuracy (PA), as follows:$${mathrm{OA}} = mathop {sum}limits_{h = 1}^H {p_{hh}}$$
(2)
$${mathrm{UA}}_h = frac{{p_{hh}}}{{p_h}}$$
(3)
$${mathrm{PA}}_j = frac{{p_{jj}}}{{p_j}}$$
(4)
where UAh is the UA for stratum h; ph and pj are the total area in stratum h and j, respectively; and PAj is the PA in stratum j.Using the preceding equations, we estimate OA, UA and PA in the four periods (Supplementary Table 2). In general, OAs are >99%, UAs are >88% and PAs are >72% in each period. The stratified random-sample approach for area estimation is considered the most robust method to investigate loss trends in GFC product and can avoid inconsistencies of the dataset due to changes in detection model and satellite sensors19. Our results show that the forest loss from the stratified random-sample approach is similar to GFC mapped loss, both of which show a consistent increase during the four periods of 2001–2019 (Fig. 1a), confirming the increasing forest carbon loss across the tropics during the early twenty-first century.Forest carbon stocksWe estimate forest (aboveground and belowground) biomass carbon losses by co-locating GFC loss data with corresponding biomass data. Forest biomass maps are not universally reliable, owing to uncertainties and some degree of bias. We use four biomass maps to quantify forest carbon stocks, which helps reduce the uncertainties and bias44. The four maps were developed by refs. 9,48,49,20 and are hereafter referred to as ‘Baccini’, ‘Saatchi’, ‘Avitabile’ and ‘Zarin’ maps, respectively. The Baccini map, derived from Moderate Resolution Imaging Spectroradiometer (MODIS) data, presents aboveground live woody biomass (AGB) across the tropics at 500 m spatial resolution. The Saatchi map, also derived from MODIS data, presents total forest carbon stocks at 1 km spatial resolution across the tropics. The Avitabile map, integrated from the Baccini and Saatchi maps, shows AGB at 1 km resolution across the tropics. The Zarin map, derived from Landsat data, presents AGB across the globe at 30 m resolution.Belowground root biomass (BGB) data are sparse because measurements of BGB are time consuming, laborious and technically challenging50. Thus, we calculate BGB (in Mg ha−1 biomass) from AGB maps (Baccini, Avitabile and Zarin) using an empirical model at the pixel level51:$${mathrm{BGB}} = 0.489 times {mathrm{AGB}}^{0.89}$$
(5)
Total forest biomass is calculated as the sum of AGB and BGB. Finally, total forest carbon stocks (MgC ha−1) in live woody forest are estimated as 50% of total biomass20,50. The Saatchi map provides total forest carbon stocks rather than AGB. The total forest carbon stocks are calculated from AGB using the same method mentioned in the preceding50. Thus, we estimate AGB and BGB from total forest carbon stocks in the Saatchi map using the preceding method to separate aboveground and belowground parts of forest carbon stocks.There are inconsistencies in the MODIS-derived biomass maps (Baccini, Saatchi and Avitabile) and Landsat-derived GFC data41, which may underestimate forest carbon loss by the three coarse forest biomass maps (Supplementary Fig. 7). The Zarin map is derived from Landsat data and considers tree cover using GFC data, and tree loss can be co-located with the corresponding biomass20, indicating the consistencies of the Zarin map and GFC. Therefore, to correct the three biomass maps with coarse resolution and reduce the inconsistencies, we resample the Zarin map from 30 m to 500 m (the resolution of the Baccini map) and 1 km (the resolutions of the Avitabile and Baccini maps) and calculate forest carbon loss using the two resampled biomass maps. We then estimate the ratios of forest carbon loss derived from the 30 m biomass map to the forest carbon loss derived from resampled biomass maps in each GFC tile (10° × 10°). The ratios are then used as a scale factor to correct the three biomass maps (Baccini, Saatchi and Avitabile). For the resampling, forest carbon density at 500 m or 1 km is averaged from all 30 m pixels in the corresponding locations.Soil organic carbon stocksDeforestation not only causes forest biomass carbon loss, but also results in loss of SOC10. We calculate SOC loss at 0–30 cm soil depth as:$${mathrm{SOC}}_{{{{mathrm{loss}}}}} = {mathrm{OCS}} times theta$$
(6)
where SOCloss is SOC loss resulting from forest loss measured in MgC ha−1, OCS is SOC stocks at 0–30 cm depth measured in MgC ha−1 and θ is the SOC loss rate.We obtain SOC stocks at 0–30 cm depth from SoilGrids (version 2.0), created by the International Soil Reference and Information Centre52. SOC stocks are calculated using a calibrated quantile random forest model at a spatial resolution of 250 m. We further resample the data from 250 m to 30 m using the nearest-neighbour method to match the scale of forest cover loss data.SOC loss rate is affected predominantly by land-use types following forest loss and tree species53. The loss rate data are compiled from a previous meta-analysis, which summarizes the rate of SOC loss resulting from forest loss over the tropics50. SOC losses resulting from primary and secondary forest loss differ (Supplementary Table 3). We use a map of primary humid tropical forests in 2001 developed by ref. 54 to classify primary and secondary forests. The land covers following forest loss are determined according to the driver of forest loss (see Drivers of forest carbon loss). The spatial resolution of the SOC data is coarse compared with GFC data, which may result in inconsistencies in the calculation. However, as SOC loss resulting from forest loss accounts for a small proportion of total forest carbon loss (8%), we ignore the potential inconsistencies.Drivers of forest carbon lossWe determine drivers of tree-cover loss using the dataset generated by ref. 27. This dataset shows the dominant driver of tree-cover loss at each 10 km grid cell for 2001–2019. There are five categories of drivers of tree-cover loss: commodity-driven deforestation, defined as permanent and/or long-term clearing of trees to other land uses (for example, commodity croplands), shifting agriculture, forestry, wildfire and urbanization. Commodity-driven deforestation, shifting agriculture and forestry dominate tropical forest loss27. Thus, wildfire and urbanization are combined and categorized as ‘others’. In tropical Africa, spatial patterns of commodity-driven deforestation are almost similar to that of shifting agriculture, as pointed out by ref. 27, resulting in large uncertainties in separating commodity agriculture from shifting agriculture. Since commodity-driven deforestation is usually for large-scale agricultural plantations, we treat commodity-driven deforestation as loss for large-scale agriculture. Because shifting agriculture is usually smallholder and/or patchy farming systems, we term shifting agriculture as small-scale agriculture, which may include commodity agriculture with similar spatial patterns to shifting agriculture in some regions such as tropical Africa. The driver data are created using decision-tree models trained by ~5,000 high-resolution Google Earth imagery cells, showing overall accuracy of 89 ± 3% from a separate validation of more than 1,500 randomly selected cells. We resample the data from 10 km resolution to 30 m using the nearest-neighbour method to match the scale of forest cover loss data.Interpretation of post-forest-loss land covers in 2020We collect cloudless and very-high-resolution satellite imagery in 2020 from Planet to determine the fate of the agriculture-driven forest loss during 2001–2019. Planet provides two products, RapidEye (at a spatial resolution of 5 m) and Doves (at a spatial resolution of 3 m; 4-band PlanetScope Scene). We randomly sample 500 pixels that show forest loss owing to agriculture expansion during 2001–2005 and 500 pixels that show forest loss owing to agriculture expansion during 2015–2019, then check cloudless satellite imagery in 2020 to visually interpret the land cover of each sampled pixel in 2020. We classify three types of land cover: agricultural land, forest/shrubland and others (Supplementary Fig. 4). Rubber and oil-palm plantations are classified as agricultural lands.Uncertainty and methods for analysisWe use committed emissions of forest carbon, even though some of this carbon will be lost only in later years or transited to other carbon pools or stored as wood products12. Forest carbon loss is defined as gross carbon loss due to forest removal (as indicated by GFC product), including (aboveground and belowground) forest biomass carbon and SOC losses. We calculate only the gross loss of forest carbon stocks while gain of carbon via reforestation and afforestation is not considered.We first estimate reference sample-based forest loss area using sample data:$${mathrm{AS}}_j = {mathrm{AM}}_jmathop {sum}limits_{h = 1}^H {p_{hj}}$$
(7)
where ASj and AMj are reference sample-based forest loss areas and mapped forest loss areas in stratum j, respectively.To estimate forest carbon loss (aboveground, belowground and soil carbon) using sample-based forest loss area, we apply a ‘stratify and multiply’ approach21,22 by assigning mean forest carbon density for each stratum. Aboveground and belowground forest carbon losses are estimated using four biomass maps (Baccini, Saatchi, Avitabile and Zarin), and we report ensemble mean ± s.d. from the four maps as our best estimate.Mountain forest carbon loss at different elevations is calculated by overlaying mountain polygon from the Global Mountain Biodiversity Assessment inventory55 (version 1.2) and the 30 m ASTER Global Digital Elevation Model56 (version 3).Although our validation shows that the GFC (v.1.7) product could accurately map forest loss, we cannot reduce the omission and commission errors. In addition, the accuracy of the disturbance year is 75.2%, with 96.7% of the disturbance occurring within one year before or after the estimated disturbance year in GFC product15. Therefore, we calculate 3 yr moving averages of annual forest and related carbon losses for time-series analysis, following the suggestion of the Global Forest Watch19 and refs. 28,38. We use a non-parametric Theil–Sen estimator regression method53 to detect trends in time-series results and test the significance of the trend by Mann–Kendall test57.Reporting SummaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

Kim, I. R. et al. Genetic diversity and population structure of nutria (Myocastor coypus) in South Korea. Animals 9, 1164. https://doi.org/10.3390/ani9121164 (2019).Article 
PubMed Central 

Google Scholar 
GISD. Of the World’s Worst Invasive Alien Species. Global Invasive Species Database. http://www.iucngisd.org/gisd/100_worst.php. 100, (2021).Hong, S., Do, Y., Kim, J. Y., Kim, D. & Joo, G. Distribution, spread and habitat preferences of nutria (Myocastor coypus) invading the lower Nakdong River, South Korea. Biol. Invas. 17, 1485–1496. https://doi.org/10.1007/s10530-014-0809-8 (2015).Article 

Google Scholar 
Ojeda, R., Bidau, C. & Emmons, L. Myocastor coypus (errata version published in 2017). The IUCN Red List Threat. Species (2016): e.T14085A121734257.Tsiamis, K. et al. Baseline Distribution of Invasive Alien Species of Union Concern (Publications Office of the European Union, 2017).
Google Scholar 
Carter, J. & Leonard, B. P. A review of the literature on the worldwide distribution, spread of, and efforts to eradicate the coypu (Myocastor coypus). Wildl. Soc. Bull. 30, 162–175 (2002).
Google Scholar 
Kim, Y. C. et al. Distribution and management of nutria (Myocastor coypus) populations in South Korea. Sustainability 11, 4169. https://doi.org/10.3390/su11154169 (2019).Article 

Google Scholar 
Park, J. H. et al. The first case of Capillaria hepatica infection in a nutria (Myocastor coypus) in Korea. Korean J. Parasitol. 52, 527–529. https://doi.org/10.3347/kjp.2014.52.5.527 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
Fratini, F., Turchi, B. E., Ebani, V. V. & Bertelloni, F. The presence of Leptospira in coypus (Myocastor coypus) and rats (Rattus norvegicus) living in a protected wetland in Tuscany (Italy). Vet. Arh. 85, 407–414 (2015).
Google Scholar 
Lee, D. H., Kil, J. H. & Kim, D. E. The study on the distribution and inhabiting status of nutria (Myocastor coypus) in Korea. Korean J. Environ. Ecol. 27, 316–326 (2013).CAS 

Google Scholar 
Guichón, M. L., Doncaster, C. P. & Cassini, M. H. Population structure of coypus (Myocastor coypus) in their region of origin and comparison with introduced populations. J. Zool. 261, 265–272. https://doi.org/10.1017/S0952836903004187 (2003).Article 

Google Scholar 
Bertolino, S., Perrone, A. & Gola, L. Effectiveness of coypu control in small Italian Wetland areas. Wildl. Soc. Bull. 33, 714–720. https://doi.org/10.2193/0091-7648(2005)33[714:EOCCIS]2.0.CO;2 (2005).Article 

Google Scholar 
Schertler, A. et al. The potential current distribution of the coypu (Myocastor coypus) in Europe and climate change induced shifts in the near future. NeoBiota 58, 129–160. https://doi.org/10.3897/neobiota.58.33118 (2020).Article 

Google Scholar 
Hilts, D. J., Belitz, M. W., Gehring, T. M., Pangle, K. L. & Uzarski, D. G. Climate change and nutria range expansion in the Eastern United States. J. Wild. Manaag. 83, 591–598. https://doi.org/10.1002/jwmg.21629’ (2019).Article 

Google Scholar 
Jarnevich, C. et al. Evaluating simplistic methods to understand current distributions and forecast distribution changes under climate change scenarios: An example with coypu (Myocastor coypus). NeoBiota 32, 107–125. https://doi.org/10.3897/neobiota.32.8884 (2017).Article 

Google Scholar 
Korean Metrological Administration, (2020). Korean Climate Change Assessment Report 2020.Guillera-Arroita, G. et al. Is my species distribution model fit for purpose? Matching data and models to applications. Glob. Ecol. Biogeogr. 24, 276–292. https://doi.org/10.1111/geb.12268 (2015).Article 

Google Scholar 
Phillips, S. J., Anderson, R. P. & Schapire, R. E. Maximum entropy modeling of species geographic distributions. Ecol. Modell. 190, 231–259. https://doi.org/10.1016/j.ecolmodel.2005.03.026 (2006).Article 

Google Scholar 
Hong, S., Cowan, P., Do, Y. & Gim, J. S. Seasonal feeding habits of coypu (Myocastor coypus) in South Korea. Hystrix 27, 123–128 (2016).
Google Scholar 
Kim, H. S., Kong, J. Y., Kim, J. H., Yeon, S. C. & Hong, I. H. A Case of Fascioliasis in A Wild Nutria, Myocastor coypus Republic of Korea. Korean J. Parasitol. 56, 375–378. https://doi.org/10.3347/kjp.2018.56.4.375 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Do, Y., Kim, J. Y., Im, R. Y. & Kim, S. B. Spatial distribution and social characteristics for wetlands in Gyeongsangnam-do Province. Korean J. Limnol. 45, 252–260 (2012).
Google Scholar 
IPCC. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. (2013).Sheffels, T. R. Status of Nutria (Myocastor coypus) Populations in the Pacific Northwest and Development of Associated Control and Management Strategies, with an Emphasis on Metropolitan Habitats, PhD Thesis (Portland State Univ., 2013).Doncaster, C. P. & MlCOL, T. Annual cycle of a coypu (Myocastor coypus) population: Male and female strategies. J. Zool. 217, 227–240. https://doi.org/10.1111/j.1469-7998.1989.tb02484.x (1989).Article 

Google Scholar 
Reggiani, G., Boitani, L. & Stefano, R. Population dynamics and regulation in the coypu Myocastor coypus in Central Italy. Ecography 18, 138–146. https://doi.org/10.1111/j.1600-0587.1995.tb00334.x (1995).Article 

Google Scholar 
Cha, Y., Cho, K. H., Lee, H., Kang, T. & Kim, J. H. The relative importance of water temperature and residence time in predicting cyanobacteria abundance in regulated rivers. Water Res. 124, 11–19. https://doi.org/10.1016/j.watres.2017.07.040 (2017).CAS 
Article 
PubMed 

Google Scholar 
Hellmann, J. J., Byers, J. E., Bierwagen, B. G. & Dukes, J. S. Five potential consequences of climate change for invasive species. Conserv. Biol. 22, 534–543. https://doi.org/10.1111/j.1523-1739.2008.00951.x (2008).Article 
PubMed 

Google Scholar 
Pereira, A. D. et al. Modeling the geographic distribution of Myocastor coypus (Mammalia, Rodentia) in Brazil: Establishing priority areas for monitoring and an alert about the risk of invasion. Stud. Neotrop. Fauna Environ. 55, 139–148. https://doi.org/10.1080/01650521.2019.1707419 (2020).Article 

Google Scholar 
Pearson, R. G. & Dawson, T. P. Predicting the impacts of climate change on the distribution of species: Are bioclimate envelope models useful?. Glob. Ecol. Biogeogr. 12, 361–371. https://doi.org/10.1046/j.1466-822X.2003.00042.x (2003).Article 

Google Scholar 
Rogers, C. E. & McCarty, J. P. Climate change and ecosystems of the mid-atlantic region. Clim. Res. 14, 235–244. https://doi.org/10.3354/cr014235 (2000).Article 

Google Scholar 
Adhikari, P. et al. Potential impact of climate change on plant invasion in the Republic of Korea. J. Ecol. Environ. 43, 36. https://doi.org/10.1186/s41610-019-0134-3 (2019).Article 

Google Scholar 
Welsch, D. J., Smart, D. L., Boyer, J. N. & Minkin, P. Forested Wetlands: Functions, Benefits and the Use of Best Management Practices (US Dept of the Interior Fish and Wildlife Service, 2021).
Google Scholar 
Borgnia, M., Galante, M. L. & Cassini, M. H. Diet of the coypu (nutria, Myocastor coypus) in agro-systems of Argentinean pampas. J. Wildl. Manag. 64, 354–361. https://doi.org/10.2307/3803233 (2000).Article 

Google Scholar 
Colares, I. G., Oliveira, R. N. V., Liveira, R. M. & Colares, E. P. Feeding habits of coypu (Myocastor coypus Molina 1978) in the wetlands of the Southern region of Brazil. An. Acad. Bras. Cienc. 82, 671–678. https://doi.org/10.1590/s0001-37652010000300015 (2010).Article 
PubMed 

Google Scholar 
Corriale, M. J., Arias, S. M., Bó, R. F. & Porini, G. Habitat-use patterns of the coypu (Myocastor coypus) in an urban wetland of its original distribution. Acta Theriol. 51, 295–302. https://doi.org/10.1007/BF03192681 (2006).Article 

Google Scholar 
Linscombe, G., Kinler, N. & Wright, V. Nutria population density and vegetative changes in brackish marsh in coastal Louisiana. In Worldwide Furbearer Conference Proceedings (eds Chapman, J. A. & Pursley, D.) 129–141 (Worlwide Furbearer Conference Inc, 1981).
Google Scholar 
Aliev, F. Contribution to the study of nutria migrations (Myocastor coypus). Saugetierkd. Mitt. 16, 301–303 (1968).
Google Scholar 
Farashi, A. & Najafabadi, M. S. A model to predict dispersion of the alien nutria, Myocastor coypus Molina, 1782 (Rodentia) Northern Iran. Acta Zool. Bulg. 69, 65–70 (2017).
Google Scholar 
Vilà, M. et al. How well do we understand the impacts of alien species on ecosystem services? A pan-European, cross-taxa assessment. Front. Ecol. Environ. 8, 135–144. https://doi.org/10.1890/080083 (2010).Article 

Google Scholar 
Adhikari, P. et al. Seasonal and altitudinal variation in roe deer (Capreolus pygargus tianschanicus) diet on Jeju Island, South Korea. J. Asia Pac. Biodivers. 9, 422–428. https://doi.org/10.1016/j.japb.2016.09.001 (2016).Article 

Google Scholar 
Koo, K. A., Kong, W. S., Nibbelink, N. P., Hopkinson, C. S. & Lee, J. H. Potential effects of climate change on the distribution of cold-tolerant evergreen broadleaved woody plants in the Korean Peninsula. PLoS ONE 10, e0134043. https://doi.org/10.1371/journal.pone.0134043 (2015).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
National Institute of Biological Research. Korean Red List of Threatened Species 2nd edn. (Ministry of Environement of Korea, 2014).
Google Scholar 
Kil, J. et al. Monitoring of Invasive Alien Species Designated by the Wildlife Protection Act (VII) (Natl Inst. of Environmental Research, 2013).
Google Scholar 
Busby, J. R. In Bioclim, a Bioclimatic Analysis and Prediction System in Nature Conservation: Cost Effective Biological Surveys and Data Analysis (eds Margules, C. R. & Austin, M. P.) 64–68 (CSIRO, 1991).
Google Scholar 
Lee I. H., Park S. H., Kang, H. S. & Cho C. H. Regional climate projections using the HadGEM3-RA in Proceedings of the 3rd International Conference on Earth System Modelling; Hamburg, Germany. 17–21 September 2012. (2012).Robert, J. H., Phillips, S., Leathwick, J. & Elith, J. Package ‘dismo’ version 1.3. , https://cran.rproject.org/web/packages/dismo.pdf (2021).Jeon, J. Y., Adhikari, P. & Seo, C. Impact of climate change on potential dispersal of Paeonia obovata (Paeoniaceae), a critically endangered medicinal plant of South Korea. Ecol. Environ. Conserv. 26, S145–S155 (2020).
Google Scholar 
Dormann, C. F. et al. Collinearity: A review of methods to deal with it and a simulation study evaluating their performance. Ecography 36, 27–46. https://doi.org/10.1111/j.1600-0587.2012.07348.x (2013).Article 

Google Scholar 
Shin, M. S., Seo, C., Lee, M. & Kim, J. Y. Prediction of potential species richness of plants adaptable to climate change in the Korean Peninsula. J. Environ. Impact Assess. 27, 562–581 (2018).
Google Scholar 
Adhikari, P. et al. Northward range expansion of southern butterflies according to climate change in South Korea. KSCCR 11, 643–656. https://doi.org/10.15531/KSCCR.2020.11.6.643 (2020).Article 

Google Scholar 
Song, C. et al. Estimation of future land cover considering shared socioeconomic pathways using scenario generators. KSCCR 9, 223–234. https://doi.org/10.15531/KSCCR.2018.9.3.223 (2018).Article 

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. https://doi.org/10.1002/joc.5086 (2017).Article 

Google Scholar 
Dukes, J. S. & Mooney, H. A. Does global change increase the success of biological invaders?. Trends Ecol. Evol. 14, 135–139. https://doi.org/10.1016/s0169-5347(98)01554-7 (1999).CAS 
Article 
PubMed 

Google Scholar 
Thuiller, W., Georges, D., Gueguen, M., Engler, R. & Breiner, F. Package ‘biomod2’: Ensemble Platform for Species Distribution Modeling, version 3.5.1. https://cran.r-project.org/web/packages/biomod2/biomod2.pdf (2021).Elith, J. et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography 29, 129–151. https://doi.org/10.1111/j.2006.0906-7590.04596.x (2006).Article 

Google Scholar 
Barbet-Massin, M., Jiguet, F., Albert, C. H. & Thuiller, W. Selecting pseudo-absences for species distribution models: how, where and how many?. Methods Ecol. Evol. 3, 327–338. https://doi.org/10.1111/j.2041-210X.2011.00172.x (2012).Article 

Google Scholar 
Brown, J. L. SDM toolbox: A python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. Methods Ecol. Evol. 5, 694–700. https://doi.org/10.1111/2041-210X.12200 (2014).Article 

Google Scholar 
Veloz, S. D. Spatially autocorrelated sampling falsely inflates measures of accuracy for presence–only niche models. J. Biogeogr. 36, 2290–2299. https://doi.org/10.1111/j.1365-2699.2009.02174.x (2009).Article 

Google Scholar 
Adhikari, P., Lee, Y. H., Park, Y.-S. & Hong, S. H. Assessment of the spatial invasion risk of intentionally introduced alien plant species (IIAPS) under environmental change in South Korea. Biology 10, 1169 (2021).Article 

Google Scholar 
Hong, S. H., Lee, Y. H., Lee, G., Lee, D. H. & Adhikari, P. Predicting impacts of climate change on northward range expansion of invasive weeds in South Korea. Plants 10, 1604. https://doi.org/10.3390/plants10081604 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Pearsons, R. G. Species distribution modeling for conservation educators and practitioners. Lessons Conserv. 3, 54–58 (2010).
Google Scholar 
Allouche, O., Tsoar, A. & Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 43, 1223–1232. https://doi.org/10.1111/j.1365-2664.2006.01214.x (2006).Article 

Google Scholar 
Thuiller, W., Lavorel, S. & Araújo, M. B. Niche properties and geographical extent as predictors of species sensitivity to climate change. Glob. Ecol. Biogeogr. 14, 347–357. https://doi.org/10.1111/j.1466-822X.2005.00162.x (2005).Article 

Google Scholar 
Lobo, J. M., Jiménez-Valverde, A. & Real, R. AUC: A misleading measure of the performance of predictive distribution models. Global. Ecol. Biogeography. 17, 145–151. https://doi.org/10.1111/j.1466-8238.2007.00358.x (2008).Article 

Google Scholar 
Swets, J. A. Measuring the accuracy of diagnostic systems. Science 240, 1285–1293 (1988).MathSciNet 
CAS 
Article 
ADS 

Google Scholar 
Baldwin, R. Use of maximum entropy modeling in wildlife research. Entropy 11, 854–866. https://doi.org/10.3390/e11040854 (2009).Article 
ADS 

Google Scholar 
Adhikari, P. et al. Potential impact of climate change on the species richness of subalpine plant species in the mountain national parks of South Korea. J. Ecol. Environ. 42, 36. https://doi.org/10.1186/s41610-018-0095-y (2018).Article 

Google Scholar 
Hijmans, R. J. et al. Package ‘raster’ v 3.5: geographical data analysis and modeling. https://cran.r-project.org/web/packages/raster/raster.pdf, (2021). More

We used publicly available, daily, county-level COVID-19 cases and deaths from the Pennsylvania Department of Health (PA DOH) (https://www.health.pa.gov/topics/disease/coronavirus/pages/cases.aspx)13,14 for Centre County and the six neighboring counties with which it shares borders: Blair, Clearfield, Clinton, Huntingdon, Mifflin, and Union (Table 1, Fig. 1). Official COVID-19 reporting for these counties began on March 1, 2020 and is ongoing.Table 1 Summary statistics. COVID-19 reporting, census data, SafeGraph mobile-device derived data.Full size tableFigure 1(a) The cumulative COVID-19 case trajectory for Centre County minus the student cases (red line) has the same shape as the outbreak for the neighboring counties. When looking at student cases only (blue line), the curve leads other counties. Centre County cumulative cases including the university (purple line) take on the shape of an early increase because of the student cases. (b) When aggregating cases from students and non-students, Centre County (purple dot) reported about the number of cases expected for its population size, relative to the neighboring counties (black dots). When the university-reported student cases are separated from the non-student residents of the county, cases reported in Centre County non-students (red dots show possible range of total cases) fall below the number of cases we would expect for the population size. Student cases only (blue dot) are high for the student population size.Full size imageWithin Centre County, PSU provided COVID-19 testing for UP students from August 7, 2020 onward and reported anonymized weekly (2020) and daily (2021) confirmed cases, negative test results, and total tests completed for each campus in a public dashboard (Figs. 1a, S1) (https://virusinfo.psu.edu/covid-19-dashboard/)8. Two types of testing were conducted: students who were enrolled in on-campus classes were randomly selected for surveillance testing and all students could use on-demand testing. Through March 23, 2021, a total of 45,092 random tests were conducted for surveillance, of which 462, or 1.0%, were infected. Surveillance testing efforts ranged from 2440 to 4020 weekly tests through the Fall 2020 semester and were designed to consistently test approximately 1% of students throughout the school year.During the same time period, 75,436 on-demand tests were conducted, of which 6093, or 8.1%, were infected. Students living in both on-campus dorms and off-campus apartments had equal access to university-provided testing. Both on-campus and off-campus residences are within Centre County so positive and negative tests results were also included in the overall Centre County reports of COVID-19 cases.Pre-arrival testing was required for students returning to campus from transmission hotspots. Students with positive tests from pre-arrival testing were required to isolate for 10–14 days after their positive test before arriving on campus. Results from pre-arrival testing for students returning to campus in the Fall of 2020 are not included in these data.At the county level, PA DOH reports the total positive, probable, and negative tests for each county. Because PSU is within Centre County, we estimated the number of total positive and negative tests for non-student Centre County residents by subtracting the PSU estimates (from the PSU dashboard) from the Centre County estimates provided by PA DOH. However, not all student tests were reported to DOH. A portion of the on-demand tests conducted for PSU UP students were completed by a third-party vendor, which required student registration. At the time of student registration, an estimated 0–25% of students registered with an address for a family home that did not reflect their residence in Centre County. Their test results were reported to the county of their registered address. This impacts a maximum of 1,166 positive student test results and 10,760 negative student tests.We conducted a sensitivity analysis to assess the uncertainty in reporting around the negative and positive students tests that may have been misallocated due to the reported residence of student tests. We have calculated the minimum and maximum number of affected positive and negative student tests. This uncertainty from student tests impacts non-student values, which are calculated by subtracting student values from county level reports. The calculations are based on a range of a possible 0–1166 positive student tests misallocated to other counties and up to 10,760 misallocated negative student tests. We have used the ranges of misallocated student tests to calculate, for non-student Centre County residents, the full possible range of (1) total cases, (2) reported cases per capita, and (3) tests per capita (Table 1, Fig. 1b). As a result, our estimates of cases and per capita testing among non-student residents in Centre County are imprecise (Table 1).We also used publicly available data from PA DOH data and PSU to calculate COVID-19 deaths per 100,000 for Centre County, the six neighboring counties, and PSU UP.We acquired county-level data on median household income, population size, and college enrollment status from the 2019 United States Census Bureau’s American Community Survey (ACS) 5-year data (https://www.census.gov/data/developers/data-sets/acs-5year.html) for all seven previously mentioned counties in central PA15.We divide the census block groups (CBG) of Centre County into two categories. We first designated ‘student-dominated CBGs’ as CBGs where  > 50% of ACS responses report enrollment as undergraduate students. We consider data from the 19 student-dominated CBGs in Centre County to be representative of the student population in Centre County. In addition to off-campus locations, the 19 student-dominated CBGs include all on-campus dorms. These 19 CBGs are either on or adjacent to PSU’s UP campus and occupy exactly 6 census tracts. The remaining 25 county census tracts were designated as non-student dominated areas.SafeGraph16 receives geolocation data from anonymized mobile devices collected from numerous applications. We analyzed SafeGraph’s mobile device-derived daily visit counts to points of interest (POI), which are fixed locations, such as businesses or attractions. SafeGraph data provide daily counts for total numbers of visits by mobile devices while using at least one application that provides geolocation data to SafeGraph. A “visit” indicates that the device entered the building or spatial perimeter designated as a POI. We acquired daily visit counts for POIs in the seven previously mentioned counties in central PA from January 1, 2019 forward (Table 1) and within Centre County grouped counts into student-dominated CBGs and non-student dominated CBGs. From January 1, 2020 forward, we used SafeGraph data on the median daily minutes that devices spent outside of their home in each county and the student- and non-student dominated CBG divisions in Centre County. The “home location” of each device is defined by its location overnight. Finally, we used SafeGraph’s weekly calculated number of devices residing in each county and the CBGs of Centre County for 2019 to measure SafeGraph’s data representation across the seven counties and the CBGs of Centre County.No administrative permissions were required to obtain these data. Academic researchers can register to receive access to SafeGraph data at no charge for non-commercial purposes only. See Data Availability statement below for details. More