Ecology

Martin WF, Bryant DA, Beatty JT. A physiological perspective on the origin and evolution of photosynthesis. FEMS Microbiol Rev. 2018;42:205–31.CAS 
PubMed 
Article 

Google Scholar 
Partensky F, Blanchot J, Vaulot D. Differential distribution and ecology of Prochlorococcus and Synechococcus in oceanic waters: a review. Bull Oceanogr Monaco, no Spec. 1999;19:457–76.
Google Scholar 
Zwirglmaier K, Jardillier L, Ostrowski M, Mazard S, Garczarek L, Vaulot D, et al. Global phylogeography of marine Synechococcus and Prochlorococcus reveals a distinct partitioning of lineages among oceanic biomes. Environ Microbiol. 2008;10:147–61.PubMed 

Google Scholar 
Callieri C. Picophytoplankton in freshwater ecosystems: the importance of small-sized phototrophs. Freshw Rev. 2008;1:1–28.Article 

Google Scholar 
Stal LJ. Physiological ecology of cyanobacteria in microbial mats and other communities. New Phytol. 1995;131:1–32.CAS 
PubMed 
Article 

Google Scholar 
Rikkinen J. Cyanobacteria in terrestrial symbiotic systems. In: Hallenbeck PC editor. Modern topics in the phototrophic prokaryotes. Switzerland: Springer; 2017. p. 243–94.Badger MR, Price GD, Long BM, Woodger FJ. The environmental plasticity and ecological genomics of the cyanobacterial CO2 concentrating mechanism. J Exp Bot. 2006;57:249–65.CAS 
PubMed 
Article 

Google Scholar 
Rae BD, Long BM, Badger MR, Price GD. Functions, compositions, and evolution of the two types of carboxysomes: polyhedral microcompartments that facilitate CO2 fixation in cyanobacteria and some proteobacteria. Microbiol Mol Biol Rev. 2013;77:357–79.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bar-On YM, Phillips R, Milo R. The biomass distribution on Earth. Proc Natl Acad Sci USA. 2018;115:6506–11.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Buitenhuis ET, Li WKW, Vaulot D, Lomas MW, Landry MR, Partensky F, et al. Picophytoplankton biomass distribution in the global ocean. Earth Syst Sci Data. 2012;4:37–46.Article 

Google Scholar 
Flombaum P, Gallegos JL, Gordillo RA, Rincón J, Zabala LL, Jiao N, et al. Present and future global distributions of the marine Cyanobacteria Prochlorococcus and Synechococcus. Proc Natl Acad Sci USA. 2013;110:9824–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Field CB, Behrenfeld MJ, Randerson JT, Falkowski P. Primary production of the biosphere: integrating terrestrial and oceanic components. Science. 1998;281:237–40.CAS 
PubMed 
Article 

Google Scholar 
Garcia-Pichel F, Belnap J, Neuer S, Schanz F. Estimates of global cyanobacterial biomass and its distribution. Arch Hydrobiol Suppl Algol Stud. 2003;109:213.
Google Scholar 
Scanlan DJ, Ostrowski M, Mazard S, Dufresne A, Garczarek L, Hess WR, et al. Ecological genomics of marine picocyanobacteria. Microbiol Mol Biol Rev. 2009;73:249–99.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Doré H, Farrant GK, Guyet U, Haguait J, Humily F, Ratin M, et al. Evolutionary mechanisms of long-term genome diversification associated with niche partitioning in marine picocyanobacteria. Front Microbiol. 2020;11:2129.Article 

Google Scholar 
Dufresne A, Ostrowski M, Scanlan DJ, Garczarek L, Mazard S, Palenik BP, et al. Unraveling the genomic mosaic of a ubiquitous genus of marine cyanobacteria. Genome Biol. 2008;9:R90.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Badger MR, Hanson D, Price GD. Evolution and diversity of CO2 concentrating mechanisms in cyanobacteria. Funct Plant Biol. 2002;29:161–73.CAS 
PubMed 
Article 

Google Scholar 
Whitehead L, Long BM, Price GD, Badger MR. Comparing the in vivo function of α-carboxysomes and β-carboxysomes in two model cyanobacteria. Plant Physiol. 2014;165:398–411.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Castenholz RW, Wilmotte A, Herdman M, Rippka R, Waterbury JB, Iteman I, et al. Phylum BX. cyanobacteria. In: Boone DR, Castenholz RW, Garrity GM editors. Bergey’s manual of systematic bacteriology. New York, NY: Springer; 2001. p. 473–599.Cabello‐Yeves PJ, Picazo A, Camacho A, Callieri C, Rosselli R, Roda-Garcia JJ, et al. Ecological and genomic features of two widespread freshwater picocyanobacteria. Environ Microbiol. 2018;20:3757–71.PubMed 
Article 
CAS 

Google Scholar 
Di Cesare A, Cabello-Yeves PJ, Chrismas NAM, Sánchez-Baracaldo P, Salcher MM, Callieri C, et al. Genome analysis of the freshwater planktonic Vulcanococcus limneticus sp. nov. reveals horizontal transfer of nitrogenase operon and alternative pathways of nitrogen utilization. BMC Genomics. 2018;19:259.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Sánchez-Baracaldo P, Bianchini G, Di Cesare A, Callieri C, Chrismas NAM. Insights into the evolution of picocyanobacteria and phycoerythrin genes (mpeBA and cpeBA). Front Microbiol. 2019;10:a45.Article 

Google Scholar 
Callieri C, Mandolini E, Bertoni R, Lauceri R, Picazo A, Camacho A, et al. Atlas of picocyanobacteria monoclonal strains from the collection of CNR-IRSA, Italy. J Limnol. 2021;80:2002.Article 

Google Scholar 
Herdman M, Castenholz RW, Iteman I, Waterbury JB, Rippka R. Subsection I (Formerly Chroococcales Wettstein 1924, emend. Rippka, Deruelles, Waterbury, Herdman and Stanier 1979). In: Boone DR, Castenholz RW, Garrity GM, editors. Bergey’s manual of systematic bacteriology, Vol 1, 2nd ed., The archaea and the deeply branching and phototrophic bacteria. New York: Springer; 2001. p. 493–514.Cabello-Yeves PJ, Haro-Moreno JM, Martin-Cuadrado A, Ghai R, Picazo A, Camacho A, et al. Novel Synechococcus genomes reconstructed from freshwater reservoirs. Front Microbiol. 2017;8:1151.PubMed 
PubMed Central 
Article 

Google Scholar 
Badger MR, Price GD. CO2 concentrating mechanisms in cyanobacteria: molecular components, their diversity and evolution. J Exp Bot. 2003;54:609–22.CAS 
PubMed 
Article 

Google Scholar 
Wheatley NM, Sundberg CD, Gidaniyan SD, Cascio D, Yeates TO. Structure and identification of a pterin dehydratase-like protein as a ribulose-bisphosphate carboxylase/oxygenase (RuBisCO) assembly factor in the α-carboxysome. J Biol Chem. 2014;289:7973–81.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Huang F, Kong WW, Sun Y, Chen T, Dykes GF, Jiang Y, et al. Rubisco accumulation factor 1 (Raf1) plays essential roles in mediating Rubisco assembly and carboxysome biogenesis. Proc Natl Acad Sci USA. 2020;117:17418–28.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kerfeld CA, Melnicki MR. Assembly, function and evolution of cyanobacterial carboxysomes. Curr Opin Plant Biol. 2016;31:66–75.CAS 
PubMed 
Article 

Google Scholar 
Kupriyanova E, Pronina N, Los D. Carbonic anhydrase—a universal enzyme of the carbon-based life. Photosynthetica. 2017;55:3–19.CAS 
Article 

Google Scholar 
DiMario RJ, Machingura MC, Waldrop GL, Moroney JV. The many types of carbonic anhydrases in photosynthetic organisms. Plant Sci. 2018;268:11–17.Heinhorst S, Cannon GC. A novel evolutionary lineage of carbonic anhydrase (epsilon class) is a component of the carboxysome shell. J Bacteriol. 2004;186:623–30.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Sawaya MR, Cannon GC, Heinhorst S, Tanaka S, Williams EB, Yeates TO, et al. The structure of beta-carbonic anhydrase from the carboxysomal shell reveals a distinct subclass with one active site for the price of two. J Biol Chem. 2006;281:7546–55.CAS 
PubMed 
Article 

Google Scholar 
Heinhorst S, Williams EB, Cai F, Murin CD, Shively JM, Cannon GC. Characterization of the carboxysomal carbonic anhydrase CsoSCA from Halothiobacillus neapolitanus. J Bacteriol. 2006;188:8087–94.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Omata T, Takahashi Y, Yamaguchi O, Nishimura T. Structure, function and regulation of the cyanobacterial high-affinity bicarbonate transporter, BCT1. Funct Plant Biol. 2002;29:151–9.CAS 
PubMed 
Article 

Google Scholar 
Omata T, Price GD, Badger MR, Okamura M, Gohta S, Ogawa T. Identification of an ATP-binding cassette transporter involved in bicarbonate uptake in the cyanobacterium Synechococcus sp. strain PCC 7942. Proc Natl Acad Sci USA. 1999;96:13571–6.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Shelden MC, Howitt SM, Price GD. Membrane topology of the cyanobacterial bicarbonate transporter, BicA, a member of the SulP (SLC26A) family. Mol Membr Biol. 2010;27:12–22.PubMed 
Article 
CAS 

Google Scholar 
Price GD, Howitt SM. The cyanobacterial bicarbonate transporter BicA: its physiological role and the implications of structural similarities with human SLC26 transporters. Biochem Cell Biol. 2011;89:178–88.CAS 
PubMed 
Article 

Google Scholar 
Price GD, Woodger FJ, Badger MR, Howitt SM, Tucker L. Identification of a SulP-type bicarbonate transporter in marine cyanobacteria. Proc Natl Acad Sci USA. 2004;101:18228–33.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bonfil DJ, Ronen-Tarazia M, Sültemeyer D, Lieman-Hurwitza J, Schatz D, Kaplan A. A putative HCO−3 transporter in the cyanobacterium Synechococcus sp. strain PCC 7942. FEBS Lett. 1998;430:236–40.CAS 
PubMed 
Article 

Google Scholar 
Price GD, Shelden MC, Howitt SM. Membrane topology of the cyanobacterial bicarbonate transporter, SbtA, and identification of potential regulatory loops. Mol Membr Biol. 2011;28:265–75.CAS 
PubMed 
Article 

Google Scholar 
Shibata M, Katoh H, Sonoda M, Ohkawa H, Shimoyama M, Fukuzawa H. Genes essential to sodium-dependent bicarbonate transport in cyanobacteria: function and phylogenetic analysis. J Biol Chem. 2002;277:18658–64.CAS 
PubMed 
Article 

Google Scholar 
Zhang P, Battchikova N, Jansen T, Appel J, Ogawa T, Aro EM. Expression and functional roles of the two distinct NDH-1 complexes and the carbon acquisition complex NdhD3/NdhF3/CupA/Sll1735 in Synechocystis sp PCC 6803. Plant Cell. 2004;16:3326–40.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Price GD, Badger MR, Woodger FJ, Long BM. Advances in understanding the cyanobacterial CO2-concentrating-mechanism (CCM): functional components, Ci transporters, diversity, genetic regulation and prospects for engineering into plants. J Exp Bot. 2008;59:1441–61.CAS 
PubMed 
Article 

Google Scholar 
Battchikova N, Eisenhut M, Aro E-M. Cyanobacterial NDH-1 complexes: novel insights and remaining puzzles. Biochim Biophys Acta Bioenerg. 2011;1807:935–44.CAS 
Article 

Google Scholar 
Koester RP, Pignon CP, Kesler DC, Willison RS, Kang M, Shen Y. Transgenic insertion of the cyanobacterial membrane protein ictB increases grain yield in Zea mays through increased photosynthesis and carbohydrate production. PLoS ONE. 2021;16:e0246359.Johnson ZI, Zinser ER, Coe A, McNulty NP, Woodward EMS, Chisholm SW. Niche partitioning among Prochlorococcus ecotypes along ocean-scale environmental gradients. Science. 2006;311:1737–40.CAS 
PubMed 
Article 

Google Scholar 
Callieri C, Coci M, Corno G, Macek M, Modenutti B, Balseiro E. Phylogenetic diversity of nonmarine picocyanobacteria. FEMS Microbiol Ecol. 2013;85:293–301.CAS 
PubMed 
Article 

Google Scholar 
Schallenberg LA, Pearman JK, Burns CW, Wood SA. Spatial abundance and distribution of picocyanobacterial communities in two contrasting lakes revealed using environmental DNA metabarcoding. FEMS Microbiol Ecol. 2021;97:fiab075.CAS 
PubMed 
Article 

Google Scholar 
Mózes A, Présing M, Vörös L. Seasonal dynamics of picocyanobacteria and picoeukaryotes in a large shallow lake (Lake Balaton, Hungary). Int Rev Hydrobiol. 2006;91:38–50.Article 
CAS 

Google Scholar 
Vörös L, Callieri C, Balogh KV, Bertoni R. Freshwater picocyanobacteria along a trophic gradient and light quality range. Hydrobiologia. 1998;369/370:117–25.Watanabe MF, Harada K, Carmichael WW, Fujiki H. Toxic microcystis. Boca Raton, FL: CRC Press; 1995.Stockner J, Callieri C, Cronberg G. Picoplankton and other non-bloom-forming cyanobacteria in lakes. In: Whitton BA, Potts M editors. The ecology of cyanobacteria. The Netherlands: Springer; 2000. p. 195–231.Flamholz AI, Prywes N, Moran U, Davidi D, Bar-On YM, Oltrogge LM. Revisiting trade-offs between Rubisco kinetic parameters. Biochemistry. 2019;58:3365–76.CAS 
PubMed 
Article 

Google Scholar 
Filazzola A, Mahdiyan O, Shuvo A, Ewins C, Moslenko L, Sadid T. A database of chlorophyll and water chemistry in freshwater lakes. Sci Data. 2020;7:1–10.Article 

Google Scholar 
Cabello‐Yeves PJ, Zemskaya TI, Zakharenko AS, Sakirko MV, Ivanov VG, Ghai R, et al. Microbiome of the deep Lake Baikal, a unique oxic bathypelagic habitat. Limnol Oceanogr. 2019.Rodrigo MA, Miracle MR, Vicente E. The meromictic Lake La Cruz (Central Spain). Patterns of stratification. Aquat Sci. 2001;63:406–16.Article 

Google Scholar 
Camacho A, Picazo A, Miracle MR, Vicente E. Spatial distribution and temporal dynamics of picocyanobacteria in a meromictic karstic lake. Arch Hydrobiol Suppl Algol Stud. 2003;109:171–84.
Google Scholar 
Vicente E, Camacho A, Rodrigo MA. Morphometry and physico-chemistry of the crenogenic meromictic Lake El Tobar (Spain). Int Ver für Theor und Angew Limnol Verhandlungen. 1993;25:698–704.CAS 

Google Scholar 
Camacho A, Miracle MR, Vicente E. Which factors determine the abundance and distribution of picocyanobacteria in inland waters? A comparison among different types of lakes and ponds. Arch für Hydrobiol. 2003;157:321–38.Article 

Google Scholar 
Kaźmierczak J, Kempe S, Kremer B, López-García P, Moreira D, Tavera R. Hydrochemistry and microbialites of the alkaline crater lake Alchichica, Mexico. Facies. 2011;57:543–70.Article 

Google Scholar 
Ghai R, Mizuno CM, Picazo A, Camacho A, Rodriguez‐Valera F. Key roles for freshwater Actinobacteria revealed by deep metagenomic sequencing. Mol Ecol. 2014;23:6073–90.CAS 
PubMed 
Article 

Google Scholar 
Cabello-Yeves PJ, Ghai R, Mehrshad M, Picazo A, Camacho A, Rodriguez-Valera F. Reconstruction of diverse Verrucomicrobial genomes from metagenome datasets of freshwater reservoirs. Front Microbiol. 2017;8:2131.PubMed 
PubMed Central 
Article 

Google Scholar 
de Hoyos C, Negro AI, Aldasoro JJ. Cyanobacteria distribution and abundance in the Spanish water reservoirs during thermal stratification. Limnetica. 2004;23:119–32.Article 

Google Scholar 
Raven J, Caldeira K, Eldefield H, Hoegh-Guldberg O, Liss P, Riebesell U, et al. Ocean acidification due to increasing atmospheric carbon dioxide. London: The Royal Society; 2005.Mangan NM, Flamholz A, Hood RD, Milo R, Savage DF. pH determines the energetic efficiency of the cyanobacterial CO2 concentrating mechanism. Proc Natl Acad Sci USA. 2016;113:E5354–62.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tadesse I, Green FB, Puhakka JA. Seasonal and diurnal variations of temperature, pH and dissolved oxygen in advanced integrated wastewater pond system® treating tannery effluent. Water Res. 2004;38:645–54.CAS 
PubMed 
Article 

Google Scholar 
Gao Y, Zhang Z, Liu X, Yi N, Zhang L, Song W, et al. Seasonal and diurnal dynamics of physicochemical parameters and gas production in vertical water column of a eutrophic pond. Ecol Eng. 2016;87:313–23.Article 

Google Scholar 
Schindler DW. Recent advances in the understanding and management of eutrophication. Limnol Oceanogr. 2006;51:356–63.Article 

Google Scholar 
Bosak T, Bush JWM, Flynn MR, Liang B, Ono S, Petroff AP, et al. Formation and stability of oxygen‐rich bubbles that shape photosynthetic mats. Geobiology. 2010;8:45–55.CAS 
PubMed 
Article 

Google Scholar 
Zeebe RE, Wolf-Gladrow D. CO2 in seawater: equilibrium, kinetics, isotopes. Amsterdam: Elsevier Science B.V.; 2001.Rae BD, Förster B, Badger MR, Price GD. The CO2-concentrating mechanism of Synechococcus WH5701 is composed of native and horizontally-acquired components. Photosynth Res. 2011;109:59–72.CAS 
PubMed 
Article 

Google Scholar 
Rippka R, Deruelles J, Waterbury JB, Herdman M, Stanier RY. Generic assignments, strain histories and properties of pure cultures of cyanobacteria. Microbiology. 1979;111:1–61.Article 

Google Scholar 
Martín-Cuadrado A-B, López-García P, Alba J-C, Moreira D, Monticelli L, Strittmatter A, et al. Metagenomics of the deep Mediterranean, a warm bathypelagic habitat. PLoS ONE. 2007;2:e914.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:455–77.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ, et al. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinforma. 2010;11:1.Article 
CAS 

Google Scholar 
Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12:59–60.CAS 
PubMed 
Article 

Google Scholar 
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2013;42:D206–14.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Tatusov RL, Natale DA, Garkavtsev IV, Tatusova TA, Shankavaram UT, Rao BS, et al. The COG database: new developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Res. 2001;29:22–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Haft DH, Loftus BJ, Richardson DL, Yang F, Eisen JA, Paulsen IT, et al. TIGRFAMs: a protein family resource for the functional identification of proteins. Nucleic Acids Res. 2001;29:41–3.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kang D, Li F, Kirton E, Thomas A, Egan R, An H, et al. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ Prepr. 2019;7:e27522v1.
Google Scholar 
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Parks DH, Chuvochina M, Waite DW, Rinke C, Skarshewski A, Chaumeil P-A, et al. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat Biotechnol. 2018;36:996–1004.CAS 
PubMed 
Article 

Google Scholar 
Asnicar F, Thomas AM, Beghini F, Mengoni C, Manara S, Manghi P, et al. Precise phylogenetic analysis of microbial isolates and genomes from metagenomes using PhyloPhlAn 3.0. Nat Commun. 2020;11:1–10.Article 
CAS 

Google Scholar 
Garczarek L, Guyet U, Doré H, Farrant GK, Hoebeke M, Brillet-Guéguen L, et al. Cyanorak v2.1: a scalable information system dedicated to the visualization and expert curation of marine and brackish picocyanobacteria genomes. Nucleic Acids Res. 2021;49:D667–76.CAS 
PubMed 
Article 

Google Scholar 
Erwin PM, Thacker RW. Cryptic diversity of the symbiotic cyanobacterium Synechococcus spongiarum among sponge hosts. Mol Ecol. 2008;17:2937–47.CAS 
PubMed 
Article 

Google Scholar 
Usher KM, Toze S, Fromont J, Ku J, Sutton DC. A new species of cyanobacterial symbiont from the marine sponge Chondrilla nucula. Symbiosis. 2004.Holtman CK, Chen Y, Sandoval P, Gonzales A, Nalty MS, Thomas TL, et al. High-throughput functional analysis of the Synechococcus elongatus PCC 7942 genome. DNA Res. 2005;12:103–15.CAS 
PubMed 
Article 

Google Scholar 
Chen M-Y, Teng W-K, Zhao L, Hu C-X, Zhou Y-K, Han B-P, et al. Comparative genomics reveals insights into cyanobacterial evolution and habitat adaptation. ISME J. 2021;15:211–27.PubMed 
Article 

Google Scholar 
Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30:3059–66.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Price MN, Dehal PS, Arkin AP. FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS ONE. 2010;5:e9490.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods Mol. Biol. 2000;132:365–86.CAS 
PubMed 

Google Scholar  More

Abiotic parametersThe temperature regime experienced by the embryos was purposefully natural and therefore varied between the three air exposure treatments. The subtidal treatment, where embryos were continuously submerged in water, remained around 9.5 °C for the duration of the experiment, while the intertidal treatments experienced dips in temperature during outside air exposure down to 2.5 °C and 0.8 °C, for the low and high intertidal respectively (Fig. 3A). Accumulated thermal units (ATU; days × temperature post collection until hatch) for each air exposure treatment were 79.6, 75.4 and 65.3 for subtidal, low intertidal and high intertidal, respectively. Despite differences in thermal regime, peak hatch was on the same day (March 14, 2021) for all air exposure and CO2 treatments, estimated at 11 dpf.Figure 3Temperature and pH experienced by the herring embryos and larvae throughout the experiment. Hourly measurements of (A) air/water temperature experienced by herring embryos in each of the tidal treatments (subtidal: continuous immersion in 9.5 °C water; low intertidal: 2 × 2 h air exposure; high intertidal: 5 + 9 h air exposure) and (B) pH levels in the tanks for each of the CO2 treatments (greens = control, 400 µatm CO2, yellows = medium, 1600 µatm CO2, reds = high, 3000 µatm CO2); dots are pH levels measured in the individual jars during larval incubation.Full size imageThe pH levels in the tanks were measured hourly and were stable over the course of the embryonic incubation period, with no overlap between treatments, although there was some overlap between individual jars. Control treatment was consistently around a pH of 8, the medium treatment had a pH of 7.4 and the high CO2 treatment had a mean pH of 7.1 (Table 1). After hatch, when the larvae were transferred to the jars, circulation and gas exchange between jars and tank were not as high and CO2 accumulated in the jars over time, leading to pH levels deviating from tank pH levels (Fig. 3B). Although oxygen levels remained high (7–9 mg/L), the pH dropped from a mean 8–7.6 in the control on two occasions, and was brought back up with a partial water exchange from the incubation tank water. The pH in the medium and high CO2 treatments were not as affected (Fig. 3B), however, final water chemistry measurements after completion of the experiment (2 days post water exchange) revealed much higher CO2 levels in all treatments (Table 1: day 15).Table 1 Mean water parameters for each treatment (mean of 3 tanks ± S.D.) at the beginning (day 1, 2021-03-06) and end (day 6, 2021-03-12) of embryonic incubation and mean parameters in the jars (N = 9) at the end of larval incubation (day 15, 2021-03-19); Temperature (T), salinity, pCO2, total CO2 (TCO2) measured at distinct sampling intervals with the BoL; total alkalinity (TA) and pH (on the total scale) calculated with CO2SYS.Full size tableEffect of air exposure and CO2 treatment during embryonic developmentNeither embryonic survival nor growth were significantly affected by treatment in our experiment. Percent daily embryonic mortality was low and not significantly affected by CO2 treatment or air exposure (CO2: p = 0.088, F2 = 2.45; Tide: p = 0.11, F2 = 2.19; CO2*Tide: p = 0.18, F2 = 1.59) . Egg diameter at 6 dpf was also not significantly affected by treatment (CO2: p = 0.38, X2 (2, N = 30) = 1.92; Tide: p = 0.83, X2 (2, N = 30) = 0.33; CO2*Tide: p = 0.08, X2 (2, N = 30) = 8.25). Metabolic rate, as indicated by embryonic heart rate, was significantly affected by air exposure at 6 dpf (p  *; 0.1  >).Full size image More

Pfennig, D. The adaptive significance of an environmentally-cued developmental switch in an anuran tadpole. Oecologia 85, 101–107 (1990).ADS 
PubMed 
Article 

Google Scholar 
Brönmark, C. & Miner, J. G. Predator-induced phenotypical change in body morphology in Crucian carp. Science 258, 1348–1350 (1992).ADS 
PubMed 
Article 

Google Scholar 
Wikelski, M. & Thom, C. Marine iguanas shrink to survive El Niño. Nature 403, 37–38 (2000).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Agrawal, A. A. Phenotypic plasticity in the interactions and evolution of species. Science 294, 321–326 (2001).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Huchard, E., English, S., Bell, M. B. V., Thavarajah, N. & Clutton-Brock, T. Competitive growth in a cooperative mammal. Nature 533, 532–534 (2016).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Travis, J. Evaluating the adaptive role of morphological plasticity. In: Ecological morphology (pp. 99–122) (The University of Chicago Press, Chicago, 1994).Lázaro, J., Dechmann, D. K. N., LaPoint, S., Wikelski, M. & Hertel, M. Profound reversible seasonal changes of individual skull size in a mammal. Curr. Biol. 27, R1106–R1107 (2017).PubMed 
Article 
CAS 

Google Scholar 
Lázaro, J. & Dechmann, D. K. Dehnel’s phenomenon. Ecol. Evol. 31, R463–R465 (2021).
Google Scholar 
Bronstein, J. L. The evolution of facilitation and mutualism. J. Ecol. 97, 1160–1170 (2009).Article 

Google Scholar 
Leigh, J. The evolution of mutualism. J. Environ. Biol. 23, 2507–2528 (2010).
Google Scholar 
Liu, C., Yang, D. R. & Peng, Y. Q. Body size in a pollinating fig wasp and implications for stability in a fig-pollinator mutualism. Entomol. Exper. Appl. 138, 249–255 (2011).Article 

Google Scholar 
Pires, M. M., Guimarães, P. R., Galetti, M. & Jordano, P. Pleistocene megafaunal extinctions and the functional loss of long-distance seed-dispersal services. Ecography 41, 153–163 (2018).Article 

Google Scholar 
Boucher, D., James, S. & Keeler, K. The ecology of mutualism. Annu. Rev. Ecol. Syst. 13, 315–347 (1982).Article 

Google Scholar 
Irwin, R. E. & Brody, A. K. Nectar robbing in Ipomopsis aggregata: effects on pollinator behavior and plant fitness. Oecologia 116, 519–527 (1998).ADS 
PubMed 
Article 

Google Scholar 
Allen, G. The Anemonefishes: their classification and biology (T.F.H. Publications, 1972).
Google Scholar 
Fautin, D.G. & Allen, G.R. Field guide to anemonefishes and their host sea anemones. (Western Australian Museum, Perth, 1992).Ollerton, J., McCollin, D., Fautin, D. G. & Allen, G. R. Finding NEMO: nestedness engendered by mutualistic organization in anemonefish and their hosts. Proc. R. Soc. B Biol. Sci. 274, 591–598 (2006).Article 

Google Scholar 
Fricke, H. & Fricke, S. Monogamy and sex change by aggressive dominance in coral reef fish. Nature 266, 830–832 (1977).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Buston, P. M. Size and growth modification in clownfish. Nature 424, 145–146 (2003).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Mariscal, R. N. The nature of the symbiosis between Indo-Pacific anemone fishes and sea anemones. Mar. Biol. 6, 58–65 (1970).Article 

Google Scholar 
Elliott, J. K., Elliott, J. M. & Mariscal, R. N. Host selection, location, and association behaviors of anemonefishes in field settlement experiments. Mar. Biol. 122, 377–389 (1995).Article 

Google Scholar 
Verde, A. E., Cleveland, A. & Lee, R. W. Nutritional exchange in a tropical tripartite symbiosis II: direct evidence for the transfer of nutrients from host anemone and zooxanthellae to anemonefish. Mar. Biol. 162, 2409–2429 (2015).Article 
CAS 

Google Scholar 
Cleveland, A., Verde, E. A. & Lee, R. W. Nutritional exchange in a tropical tripartite symbiosis: direct evidence for the transfer of nutrients from anemonefish to host anemone and zooxanthellae. Mar. Biol. 158, 589–602 (2011).Article 

Google Scholar 
Sale, P. F. Effect of cover on agonistic behavior of a reef fish: a possible spacing mechanism. Ecology 53, 753–758 (1972).Article 

Google Scholar 
Fricke, H. W. & Holzberg, S. Social units and hermaphroditism in a pomacentrid fish. Naturwissenschaften 61, 367–368 (1974).ADS 
Article 

Google Scholar 
Fricke, H. W. Control of different mating systems in a coral reef fish by one environmental factor. Anim. Behav. 28, 561–569 (1980).Article 

Google Scholar 
Mitchell, J. S. & Dill, L. M. Why is group size correlated with the size of the host sea anemone in the false clown anemonefish?. Canad. J. Zool. 83, 372–376 (2005).Article 

Google Scholar 
Chausson, J., Srinivasan, M. & Jones, G. P. Host anemone size as a determinant of social group size and structure in the orange clownfish (Amphiprion percula). PeerJ 6, e5841 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Reed, C., Branconi, R., Majoris, J., Johnson, C. & Buston, P. Competitive growth in a social fish. Biol. Lett. 15, 20180737 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Buston, P. M. Mortality is associated with social rank in the clown anemonefish (Amphiprion percula). Mar. Biol. 143, 811–815 (2003).Article 

Google Scholar 
Branconi, R. et al. Ecological and social constraints combine to promote evolution of non-breeding strategies in clownfish. Comm. Biol. 3, 1–7 (2020).Article 
CAS 

Google Scholar 
Schmiege, P. F., D’Aloia, C. C. & Buston, P. M. Anemonefish personalities influence the strength of mutualistic interactions with host sea anemones. Mar. Biol. 164, 24 (2017).Article 

Google Scholar 
Barbasch, T. A. & Buston, P. M. Plasticity and personality of parental care in the clown anemonefish. Anim. Behav. 136, 65–73 (2018).Article 

Google Scholar 
Abramoff, M. D., Magalhaes, P. J. & Ram, S. J. Image PROcessing with ImageJ. Biophoto. Int. 11, 36–42 (2004).
Google Scholar 
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/ (2020).Goodrich, B., Gabry, J., Ali I. & Brilleman, S. Rstanarm: Bayesian applied regression modeling via Stan. R package version 2.21.1 https://mc-stan.org/rstanarm (2020).Weatherley, A. H. Approaches to understanding fish growth. Trans. Am. Fish. Soc. 119, 662–672 (1990).Article 

Google Scholar 
Gabry, J. Shinystan: interactive visual and numerical diagnostics and posterior analysis for Bayesian models. R package version 2.5.0. https://CRAN.R-project.org/package=shinystan (2018).Gelman, A., Goodrich, B., Gabry, J. & Vehtari, A. R-squared for Bayesian regression models. Am. Stat. 3, 307–309 (2018).MathSciNet 

Google Scholar 
Lüdecke, D., Ben-Shachar, M. S., Patil, I., Waggoner, P. & Makowski, D. Performance: an R package for assessment, comparison and testing of statistical models. J. Open Sour. Softw. 6, 60 (2021).
Google Scholar 
Gabry, J., Simpson, D., Vehtari, A., Betancourt, M. & Gelman, A. Visualization in Bayesian workflow. J. R. Stat. Soc. Ser. A Stat. Soc. 182, 389–402 (2019).MathSciNet 
Article 

Google Scholar 
Gabry, J. & Mahr, T. Bayesplot: plotting for bayesian models. R package version 1.8.0. https://mc-stan.org/bayesplot/ (2021).Elliott, J. K. & Mariscal, R. N. Coexistence of nine anemonefish species: differential host and habitat utilization, size and recruitment. Mar. Biol. 138, 23–36 (2001).Article 

Google Scholar 
Buston, P. M. Forcible eviction and prevention of recruitment in the clown anemonefish. Behav. Ecol. 14, 576–582 (2003).Article 

Google Scholar 
Fautin, D. G. & Allen, G. R. Anemone fishes and their host sea anemones: a guide for aquarists and divers. Sea Challengers (1997).Beldade, R., Blandin, A., O’Donnell, R. & Mills, S. C. Cascading effects of thermally-induced anemone bleaching on associated anemonefish hormonal stress response and reproduction. Nat. Commun. 8, 1–9 (2017).ADS 
CAS 
Article 

Google Scholar 
Cortese, D. et al. Physiological and behavioural effects of anemone bleaching on symbiont anemonefish in the wild. Funct. Ecol. 35, 663–674 (2021).Article 

Google Scholar 
Scherbatskoy, E. C. et al. Characterization of a novel picornavirus isolated from moribund aquacultured clownfish. J. Gen. Virol. 101, 735–745 (2020).CAS 
PubMed 
Article 

Google Scholar 
Saenz-Agudelo, P., Jones, G. P., Thorrold, S. R. & Planes, S. Mothers matter: contribution to local replenishment is linked to female size, mate replacement and fecundity in a fish metapopulation. Mar. Biol. 162, 3–14 (2014).Article 

Google Scholar 
Barbasch, T. A. et al. Substantial plasticity of reproduction and parental care in response to local resource availability in a wild clownfish population. Oikos 129, 1844–1855 (2020).Article 

Google Scholar 
Sebens, K. P. The ecology of indeterminate growth in animals. A. Rev. Ecol. Syst. 18, 371–407 (1987).Article 

Google Scholar 
Buston, P. M. & García, M. B. An extraordinary life span estimate for the clown anemonefish Amphiprion percula. J. Fish Biol. 70, 1710–1719 (2007).Article 

Google Scholar 
Chamberlain, S. A., Kilpatrick, J. R. & Holland, J. N. Do extrafloral nectar resources, species abundances, and body sizes contribute to the structure of ant–plant mutualistic networks?. Oecologia 164, 741–750 (2010).ADS 
PubMed 
Article 

Google Scholar 
Marting, P. R., Kallman, N. M., Wcislo, W. T. & Pratt, S. C. Ant-plant sociometry in the Azteca-Cecropia mutualism. Sci. Rep. 8, 1–15 (2018).Article 
CAS 

Google Scholar 
Fordyce, J. A. The evolutionary consequences of ecological interactions mediated through phenotypic plasticity. J. Exp. Biol. 209, 2377–2383 (2006).PubMed 
Article 

Google Scholar 
West-Eberhard, M. J. Developmental plasticity and evolution (Oxford University Press, 2003).Book 

Google Scholar 
West-Eberhard, M. J. Phenotypic accommodation: adaptive innovation due to developmental plasticity. J. Exp. Zool. B Mol. Develop. Evol. 304, 610–618 (2005).Article 

Google Scholar 
Moczek, A. P. et al. The role of developmental plasticity in evolutionary innovation. Proc. R. Soc. B Biol. Sci. 278, 2705–2713 (2011).Article 

Google Scholar  More

Angilletta, M. J., Niewiarowski, P. H. & Navas, C. A. The evolution of thermal physiology in ectotherms. J. Therm. Biol 27, 249–268. https://doi.org/10.1016/S0306-4565(01)00094-8 (2002).Article 

Google Scholar 
Ebersole, J. L., Liss, W. J. & Frissell, C. A. Cold water patches in warm streams: physicochemical characteristics and the influence of shading. JAWRA J. Am. Water Resour. Assoc. 39, 355–368. https://doi.org/10.1111/j.1752-1688.2003.tb04390.x (2003).ADS 
Article 

Google Scholar 
Comte, L. & Grenouillet, G. Do stream fish track climate change? Assessing distribution shifts in recent decades. Ecography 36, 1236–1246. https://doi.org/10.1111/j.1600-0587.2013.00282.x (2013).Article 

Google Scholar 
Kurylyk, B. L., MacQuarrie, K. T. B., Linnansaari, T., Cunjak, R. A. & Curry, R. A. Preserving, augmenting, and creating cold-water thermal refugia in rivers: Concepts derived from research on the Miramichi River, New Brunswick (Canada). Ecohydrology 8, 1095–1108. https://doi.org/10.1002/eco.1566 (2015).Article 

Google Scholar 
Ebersole, J. L., Quiñones, R. M., Clements, S. & Letcher, B. H. Managing climate refugia for freshwater fishes under an expanding human footprint. Front. Ecol. Environ. 18, 271–280. https://doi.org/10.1002/fee.2206 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Caissie, D. The thermal regime of rivers: a review. Freshw. Biol. 51, 1389–1406. https://doi.org/10.1111/j.1365-2427.2006.01597.x (2006).Article 

Google Scholar 
Dick, J. J., Tetzlaff, D. & Soulsby, C. Landscape influence on small-scale water temperature variations in a moorland catchment. Hydrol. Process. 29, 3098–3111. https://doi.org/10.1002/hyp.10423 (2015).ADS 
Article 

Google Scholar 
Fullerton, A. H. et al. Rethinking the longitudinal stream temperature paradigm: region-wide comparison of thermal infrared imagery reveals unexpected complexity of river temperatures. Hydrol. Process. 29, 4719–4737. https://doi.org/10.1002/hyp.10506 (2015).ADS 
Article 

Google Scholar 
Fullerton, A. H. et al. Longitudinal thermal heterogeneity in rivers and refugia for coldwater species: Effects of scale and climate change. Aquatic Sci. 80, 3. https://doi.org/10.1007/s00027-017-0557-9 (2018).Article 

Google Scholar 
Segura, C., Caldwell, P., Sun, G., McNulty, S. & Zhang, Y. A model to predict stream water temperature across the conterminous USA. Hydrol. Process. 29, 2178–2195. https://doi.org/10.1002/hyp.10357 (2015).ADS 
Article 

Google Scholar 
Jonkers, A. R. T. & Sharkey, K. J. The differential warming response of Britain’s rivers (1982–2011). PLOS One 11, e0166247. https://doi.org/10.1371/journal.pone.0166247 (2016).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Jackson, F. L., Hannah, D. M., Fryer, R. J., Millar, C. P. & Malcolm, I. A. Development of spatial regression models for predicting summer river temperatures from landscape characteristics: Implications for land and fisheries management. Hydrol. Process. 31, 1225–1238. https://doi.org/10.1002/hyp.11087 (2017).ADS 
Article 

Google Scholar 
Maheu, A., Poff, N. L. & St-Hilaire, A. A classification of stream water temperature regimes in the conterminous USA. River Res. Appl. 32, 896–906. https://doi.org/10.1002/rra.2906 (2016).Article 

Google Scholar 
Steel, E. A., Sowder, C. & Peterson, E. E. Spatial and temporal variation of water temperature regimes on the Snoqualmie River network. J. Am. Water Resour. Assoc. 52, 769–787. https://doi.org/10.1111/1752-1688.12423 (2016).Article 

Google Scholar 
Kearney, M. R., Matzelle, A. & Helmuth, B. Biomechanics meets the ecological niche: The importance of temporal data resolution. J. Exp. Biol. 215, 922–933. https://doi.org/10.1242/jeb.059634 (2012).Article 
PubMed 

Google Scholar 
Burgmer, T., Hillebrand, H. & Pfenninger, M. Effects of climate-driven temperature changes on the diversity of freshwater macroinvertebrates. Oecologia 151, 93–103. https://doi.org/10.1007/s00442-006-0542-9 (2007).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Isaak, D. J., Young, M. K., Nagel, D. E., Horan, D. L. & Groce, M. C. The cold-water climate shield: Delineating refugia for preserving salmonid fishes through the 21st century. Glob. Change Biol. 21, 2540–2553. https://doi.org/10.1111/gcb.12879 (2015).ADS 
Article 

Google Scholar 
Steel, E. A., Beechie, T. J., Torgersen, C. E. & Fullerton, A. H. Envisioning, quantifying, and managing thermal regimes on river networks. Bioscience 67, 506–522. https://doi.org/10.1093/biosci/bix047 (2017).Article 

Google Scholar 
Budescu, D. V. Dominance analysis: A new approach to the problem of relative importance of predictors in multiple regression. Psychol. Bull. 114, 542–551. https://doi.org/10.1037/0033-2909.114.3.542 (1993).Article 

Google Scholar 
Singhal, B. B. S. & Gupta, R. P. Applied Hydrogeology of Fractured Rocks. 2 edn, 408 (Springer, 2010).Shimizu, T. Relation between scanty runoff from mountainous watershed and geology, slope and vegetation (in Japanese with English summary). Bull. Forestry Forest Prod. Res. Inst. 310, 109–128 (1980).
Google Scholar 
Iwasaki, K., Nagasaka, Y. & Nagasaka, A. Geological effects on the scaling relationships of groundwater contributions in Forested Watersheds. Water Resour. Res. 57, e2021WR029641. https://doi.org/10.1029/2021WR029641 (2021).ADS 
Article 

Google Scholar 
Ishiyama, N. et al. The role of geology in creating stream climate-change refugia along climate gradients. bioRxiv, 2022.2005.2002.490355, https://doi.org/10.1101/2022.05.02.490355 (2022).Kanno, Y., Vokoun, J. C. & Letcher, B. H. Paired stream-air temperature measurements reveal fine-scale thermal heterogeneity within headwater brook trout stream networks. River Res. Appl. 30, 745–755. https://doi.org/10.1002/rra.2677 (2014).Article 

Google Scholar 
Snyder, C. D., Hitt, N. P. & Young, J. A. Accounting for groundwater in stream fish thermal habitat responses to climate change. Ecol. Appl. 25, 1397–1419. https://doi.org/10.1890/14-1354.1 (2015).Article 
PubMed 

Google Scholar 
Carslaw, D. C. & Ropkins, K. Openair—an R package for air quality data analysis. Environ. Model. Softw. 27–28, 52–61. https://doi.org/10.1016/j.envsoft.2011.09.008 (2012).Article 

Google Scholar 
Pinheiro, J. C. & Bates, D. M. Mixed-Effects Models in S and S-PLUS. (Springer, 2000).Gelman, A. & Hill, J. Data Analysis Using Regression and Multilevel/Hierarchical Models. (Cambridge University Press, 2006).Harrison, X. A. et al. A brief introduction to mixed effects modelling and multi-model inference in ecology. PeerJ 6, e4794. https://doi.org/10.7717/peerj.4794 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Clarke, P. When can group level clustering be ignored? Multilevel models versus single-level models with sparse data. J. Epidemiol. Commun. Health 62, 752. https://doi.org/10.1136/jech.2007.060798 (2008).CAS 
Article 

Google Scholar 
Theall, K. P. et al. Impact of small group size on neighbourhood influences in multilevel models. J. Epidemiol. Commun. Health 65, 688–695. https://doi.org/10.1136/jech.2009.097956 (2011).Article 

Google Scholar 
Nakagawa, S. & Schielzeth, H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol. Evol. 4, 133–142. https://doi.org/10.1111/j.2041-210x.2012.00261.x (2013).Article 

Google Scholar 
Nakagawa, S., Johnson, P. C. D. & Schielzeth, H. The coefficient of determination R2 and intra-class correlation coefficient from generalized linear mixed-effects models revisited and expanded. J. Royal Soc. Interface 14, 20170213. https://doi.org/10.1098/rsif.2017.0213 (2017).Article 

Google Scholar 
Lüdecke, D., Ben-Shachar, M. S., Patil, I., Waggoner, P. & Makowski, D. performance: An R package for assessment, comparison and testing of statistical models. J. Open Source Softw. 6, 3139. https://doi.org/10.21105/joss.03139 (2021).ADS 
Article 

Google Scholar 
Hair, J. F., Black, W. C., Babin, B. J. & Anderson, R. E. Multivariate Data Analysis: A Global Perspective. 7 edn, (Prentice Hall, 2009).Azen, R. & Budescu, D. V. The dominance analysis approach for comparing predictors in multiple regression. Psychol. Methods 8, 129–148. https://doi.org/10.1037/1082-989x.8.2.129 (2003).Article 
PubMed 

Google Scholar 
Grömping, U. Estimators of relative importance in linear regression based on variance decomposition. Am. Stat. 61, 139–147. https://doi.org/10.1198/000313007X188252 (2007).MathSciNet 
Article 

Google Scholar 
Luo, W. & Azen, R. Determining predictor importance in hierarchical linear models using dominance analysis. J. Educ. Behav. Stat. 38, 3–31. https://doi.org/10.3102/1076998612458319 (2013).Article 

Google Scholar 
R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/. (2020).Erickson, T. R. & Stefan, H. G. Linear air/water temperature correlations for streams during open water periods. J. Hydrol. Eng. 5, 317–321. https://doi.org/10.1061/(ASCE)1084-0699(2000)5:3(317) (2000).Article 

Google Scholar 
Webb, B. W., Clack, P. D. & Walling, D. E. Water–air temperature relationships in a Devon river system and the role of flow. Hydrol. Process. 17, 3069–3084. https://doi.org/10.1002/hyp.1280 (2003).ADS 
Article 

Google Scholar 
Gu, Z., Gu, L., Eils, R., Schlesner, M. & Brors, B. Circlize implements and enhances circular visualization in R. Bioinformatics.
30, 2811–2812. https://doi.org/10.1093/bioinformatics/btu393 (2014).Sugimoto, S., Nakamura, F. & Ito, A. Heat budget and statistical analysis of the relationship between stream temperature and riparian forest in the Toikanbetsu River Basin, Northern Japan. J. For. Res. 2, 103–107. https://doi.org/10.1007/BF02348477 (1997).Article 

Google Scholar 
Dugdale, S. J., Malcolm, I. A., Kantola, K. & Hannah, D. M. Stream temperature under contrasting riparian forest cover: Understanding thermal dynamics and heat exchange processes. Sci. Total Environ. 610–611, 1375–1389. https://doi.org/10.1016/j.scitotenv.2017.08.198 (2018).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Timm, A., Ouellet, V. & Daniels, M. Riparian land cover, water temperature variability, and thermal stress for aquatic species in urban streams. Water 13, 2732. https://doi.org/10.3390/w13192732 (2021).Article 

Google Scholar 
Mitchell, S. A simple model for estimating mean monthly stream temperatures after riparian canopy removal. Environ. Manage. 24, 77–83. https://doi.org/10.1007/s002679900216 (1999).CAS 
Article 
PubMed 

Google Scholar 
Horne, J. P. & Hubbart, J. A. A spatially distributed investigation of stream water temperature in a contemporary mixed-land-use watershed. Water 12, 1756. https://doi.org/10.3390/w12061756 (2020).Article 

Google Scholar 
Graham, C. B., Barnard, H. R., Kavanagh, K. L. & McNamara, J. P. Catchment scale controls the temporal connection of transpiration and diel fluctuations in streamflow. Hydrol. Process. 27, 2541–2556. https://doi.org/10.1002/hyp.9334 (2013).ADS 
Article 

Google Scholar 
Sun, H., Kasahara, T., Otsuki, K., Saito, T. & Onda, Y. Spatio-temporal streamflow generation in a small, steep headwater catchment in Western Japan. Hydrol. Sci. J. 62, 818–829. https://doi.org/10.1080/02626667.2016.1266635 (2017).Article 

Google Scholar 
Sophocleous, M. Interactions between groundwater and surface water: The state of the science. Hydrogeol. J. 10, 52–67. https://doi.org/10.1007/s10040-001-0170-8 (2002).ADS 
CAS 
Article 

Google Scholar 
Arnott, S., Hilton, J. & Webb, B. W. The impact of geological control on flow accretion in lowland permeable catchments. Hydrol. Res. 40, 533–543. https://doi.org/10.2166/nh.2009.017 (2009).Article 

Google Scholar 
Calvache, M. L., Duque, C., Fontalva, J. M. G. & Crespo, F. Processes affecting groundwater temperature patterns in a coastal aquifer. Int. J. Environ. Sci. Technol. 8, 223–236. https://doi.org/10.1007/BF03326211 (2011).Article 

Google Scholar 
Nejadhashemi, A. P., Wardynski, B. J. & Munoz, J. D. Evaluating the impacts of land use changes on hydrologic responses in the agricultural regions of Michigan and Wisconsin. Hydrol. Earth Syst. Sci. 2011, 3421–3468, https://doi.org/10.5194/hessd-8-3421-2011 (2011).Macedo, M. N. et al. Land-use-driven stream warming in southeastern Amazonia. Philos. Trans. R Soc. Lond. B Biol. Sci. 368, 20120153–20120153. https://doi.org/10.1098/rstb.2012.0153 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Carlson, K. M. et al. Influence of watershed-climate interactions on stream temperature, sediment yield, and metabolism along a land use intensity gradient in Indonesian Borneo. J. Geophys. Res. Biogeosci. 119, 1110–1128. https://doi.org/10.1002/2013JG002516 (2014).Article 

Google Scholar 
Wang, Y. I., He, B. I. N. & Takase, K. Effects of temporal resolution on hydrological model parameters and its impact on prediction of river discharge. Hydrol. Sci. J. 54, 886–898. https://doi.org/10.1623/hysj.54.5.886 (2009).Article 

Google Scholar 
Levin, S. A. The problem of pattern and scale in ecology: The Robert H MacArthur award lecture. Ecology 73, 1943–1967. https://doi.org/10.2307/1941447 (1992).Article 

Google Scholar 
García Molinos, J. & Donohue, I. Downscaling the non-stationary effect of climate forcing on local-scale dynamics: The importance of environmental filters. Clim. Change 124, 333–346. https://doi.org/10.1007/s10584-014-1077-4 (2014).ADS 
Article 

Google Scholar 
Newman, E. A., Kennedy, M. C., Falk, D. A. & McKenzie, D. Scaling and complexity in landscape ecology. Front. Ecol. Evolution https://doi.org/10.3389/fevo.2019.00293 (2019).Article 

Google Scholar 
Atkinson, S. E., Woods, R. A. & Sivapalan, M. Climate and landscape controls on water balance model complexity over changing timescales. Water Resour. Res. 38, 50-51–50-17, https://doi.org/10.1029/2002WR001487 (2002).Engel, M. et al. Controls on spatial and temporal variability in streamflow and hydrochemistry in a glacierized catchment. Hydrol. Earth Syst. Sci. 23, 2041–2063. https://doi.org/10.5194/hess-23-2041-2019 (2019).ADS 
Article 

Google Scholar 
Karlsen, R. H. et al. Landscape controls on spatiotemporal discharge variability in a boreal catchment. Water Resour. Res. 52, 6541–6556. https://doi.org/10.1002/2016WR019186 (2016).ADS 
Article 

Google Scholar 
Parmesan, C. & Yohe, G. A globally coherent fingerprint of climate change impacts across natural systems. Nature 421, 37–42. https://doi.org/10.1038/nature01286 (2003).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Weiskopf, S. R. et al. Climate change effects on biodiversity, ecosystems, ecosystem services, and natural resource management in the United States. Sci. Total Environ. 733, 137782. https://doi.org/10.1016/j.scitotenv.2020.137782 (2020).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Radchuk, V. et al. Adaptive responses of animals to climate change are most likely insufficient. Nat. Commun. 10, 3109. https://doi.org/10.1038/s41467-019-10924-4 (2019).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Kingsford, R. T. Conservation management of rivers and wetlands under climate change—a synthesis. Mar. Freshw. Res. 62, 217–222. https://doi.org/10.1071/MF11029 (2011).CAS 
Article 

Google Scholar  More

This section is devoted to show results and to highlight eventual effects of the interplay between the nonlinearity characterizing the system dynamics and the presence of noisy fluctuations for the irradiance variable.Analysis of experimental dataThe need of taking into account noisy fluctuations of such an environmental variable is well demonstrated in Fig. 1. In the first panel (a) the experimental time behaviour of the irradiance is shown. This noisy curve is based on the experimental data (purple points) of the Boussole buoy located in the Gulf of Lion, collected over a period of nine years, precisely from 2004 to 2013. The time series of the experimental data presents quite a few gaps in time due to the malfunction of the buoy. This aspect has been remedied by merging the experimental data with those of the OASIM model validated for the Boussole site61 (yellow points). The latter is a multispectral atmospheric radiative transfer model that is in turn forced by experimental-model data based on ECMWF ERAINTERIM reanalyses which provide, for example, cloud cover data. The radiative model is partly stochastic since it considers the effects stemming from the presence of clouds, averaged along a single day (this explains why the yellow points are slightly less scattered). We see that the OASIM model accurately reproduces the profile which emerges from the experimental data. Further, we stress that the experimental data are only used in this initial analysis. In the biogeochemical simulations the irradiance signal is fully reconstructed starting from a realistic seasonal cycle combined with a range of different random fluctuations, and the information from OASIM is not used. In the second panel (b) the daily (black points) as well as the three-month (red points) running mean of the experimental series are plotted. Figure 1c shows the irradiance noisy fluctuations (INF) which have been obtained by subtracting the three-month running mean curve (3MRM, red curve in Fig. 1b) from the daily running mean one (DRM, black curve in Fig. 1b) and normalizing with respect to the mean of the 3MRM ((overline{3MRM})), namely (INF = (DRM – 3MRM) / overline{3MRM}). We see that a seasonal overall trend with higher oscillations during the winter time can be seen, implying that the characteristics of the noise may change over the year. Moreover, a slight imbalance between positive and negative values of the noisy fluctuations (that is, different values of the maximum fluctuation intensity) is present. The physical reason for the occurrence of such an aspect can be ascribed to the fact that the maximum value of solar irradiance corresponds to that measured during a sunny day. Conversely, the minimum level tends to zero corresponding to a dense darkness. While the former is close to the mean value of the solar irradiance (most of all in summer), the latter is much further away and then a natural asymmetry arises in the random fluctuations. However, it should be noted that, apart from the intense spikes, the asymmetry is not so pronounced, as proved by the mean value (red line in Fig. 1c) which is practically zero, namely (0.4%) of the (overline{3MRM}). Therefore, basing on this last observation, to model the noise affecting the irradiance dynamics, as a first approximation we consider a symmetric Gaussian autocorrelated noise as described in the next subsection.On the basis of such experimental results, we postulate the hypothesis that random fluctuations of light cannot be neglected, most of all in the study of ecological systems where light profoundly determines the system dynamics, governing fundamental processes at the basis of of the food web.Figure 1(a) Experimental data (purple points) of the stochastic solar irradiance collected by the Boussole buoy in a time-window of 9 years (2004-2013); the yellow points are the data generated by the OASIM model used to fill the gaps present in the experimental time-series due to malfunctioning of the buoy. (b) Daily (black points) and three-month (red points) running mean of the light curve in panel (a). (c) Irradiance noisy fluctuations (INF), obtained by subtracting the three-month running mean curve (3MRM) from the daily running mean one (DRM) and normalizing with respect to the mean value of 3MRM ((overline{3MRM})), namely (INF = (DRM – 3MRM) / overline{3MRM}); the red line represents the mean value of such fluctuations. Data already presented and validated in61.Full size imageSolar irradianceThe solar irradiance forcing is derived considering a deterministic seasonal oscillation combined with an Ornstein-Uhlenbeck process. The coefficient of variation (CV) of simulated light forcing, Fig. 2, (CV=sigma / mu) ((mu) and (sigma) being mean value and standard deviation calculated over both time and numerical realizations), is shown for 231 (D-tau) pairs. D and (tau) represent the intensity of a Gaussian noise source and the auto-correlation time of the fluctuations, respectively (see Eqs. (2) and (3)).Each pixel represents the mean value on time of CV calculated with respect to 1000 different stochastic realizations. Figure 2Coefficient of variation ((CV=sigma / mu)) of irradiance resulting from numerical integration of model equations for 231 (D-tau) different scenarios.Full size imageIt is easy to see the agreement between the results obtained from the numerical integration and the theoretical ones derivable from Eq. (5) by putting (text {var}{F_L(0)}=0) and (t gg 1), getting (sigma ^2_L=D / 2tau). In Fig. 2, indeed, the maximum values of (sigma) lie in the upper left part of the plot corresponding to small (high) values of (tau) (D). As it is clear the values of D have been chosen in order to obtain a relative standard deviation ranging from (5%mu) to (60%mu). We underline that, in this case, it is possible to interchangeably consider (sigma) and CV since the dependence of CV on D and (tau) does not differ from that of (sigma) (meaning that the dependence of (sigma) is not altered by dividing by (mu)) (results not shown).Effects on population dynamicsIn this section the noise-induced effects on the population dynamics are examined. The nine planktonic populations present a different qualitative behaviour of the CV, compared to that of the irradiance. In this case, the CV is characterized by a strong non-monotonic dependence on the parameter (tau). This aspect can be appreciated in Fig. 3 where different curves of CV versus the time correlation parameter are shown for different fixed values of D.Figure 3Coefficient of variation ((CV=sigma / mu)) of the nine planktonic populations resulting from numerical integration of model equations plotted versus the considered values of (tau); the different curves are related to different values of the noise intensity D.Full size imageThe existence of a maximum value for CV can be appreciated for each species. Although the qualitative behaviour is the same for all strains, particular attention has to be payed on diatoms and nanoflagellates. All the other species, indeed, present a percent variation of standard deviation between (2%) and (15%). In the case of nanoflagellates, instead, the D-dependent range is (20-90%), while diatoms reach values over the (100%) for the highest values of D. Therefore, these two species, in particular, and the whole system, in general, are extremely sensitive to the auto-correlation time which characterizes the noise.We note that the different curves related to the different selected values of D approach the horizontal axis, tending asymptotically to vanish as (tau) increases. Such a behaviour can be explained by the fact that high values of (tau) give rise to a more correlated dynamics, so that (tau rightarrow infty) implies fully correlated time-behaviours corresponding to the deterministic case. In this instance, then, all the different realizations give the same results, making the standard deviation vanish. The same happens, independently of the value of (tau), for low values of noise intensity for which the corresponding curves approach the same almost vanishing value (see orange, gray and yellow lines). Differently from the previous case, when (tau rightarrow 0) the noise tends to a delta-correlated noise, that is a white noise; for (tau ne 0), instead, the noise spectrum is not flat, being characterized by a Cauchy-Lorentz distribution. The strong nonmonotonicity of CV with respect to (tau), emerging when there are relatively high values of CV, implies a greater variability of the system biomass. Lower values of CV indicate that the system dynamics is less influenced by the presence of noise where very little or no differences with respect to the deterministic case are present. Conversely, high values of CV clearly demonstrate the remarkable signature of the presence of an impacting noise source. It is interesting to note that the noise influence on the ecosystem strongly depends on both (tau) and D, that is, just an intense noise is not enough to generate a greater response of the ecosystem. In particular, experimental data are characterized by a CV approximately equal to 0.361, which corresponds to values of D and (tau) lying on the diagonal strip in Fig. 2 ranging from ((tau ,D)=(0.5,10^4)) to ((tau ,D)=(365,10^7)). Finally we note the presence of a noise suppression effect. High values of D, indeed, can generate slight effects when the correlation time (tau) does not take on suitable values.The results shown here are an extension of the previous work by Benincà et al.56. There, the authors analyse a simpler, less realistic model of two interacting populations, whose dynamics is affected by a randomly fluctuating temperature. In that case, moreover, the deterministic oscillations of the temperature are suppressed, and the system exhibits intrinsic Lotka-Volterra oscillations whose frequency match with the characteristic one(s) of the noise. On the contrary, here, the observed maximum response (see Fig. 3) cannot be interpreted as a synchronization effect, since our model does not present intrinsic Lotka-Volterra-like oscillations and the periodic population variability is only due to the deterministic forcing(s).The nonmonotonic behaviour of the CV can be then interpreted as the signature of the intimate interplay between the ecological system and the noise. This interplay, indeed, has a pivotal role in both determining the dynamics of the populations and defining the characteristics of the ecosystem.In Fig. 3 it can be observed that the value of (tau) for which CV is maximum strongly depends on the noise intensity D. In particular, it is possible to note that the peaks in Fig. 3 move towards higher values of (tau) as the noise intensity increases. Thus, Fig. 3 demonstrates that the maximum-response effect to the random fluctuations is sensitive to the noise intensity D.However, it is important to underline that the response of the system to the noisy signal does not depend on the yearly oscillations induced by the deterministic forcings. Indeed, by considering constant the deterministic part of all external forcings (temperature, irradiance, wind and salinity), the non monotonic behaviour of CV with respect to both (tau) and D is still present, provided that the populations are not extinct (plot not shown). In this scenario indeed, besides dinoflagellates, diatoms and nanoflagellates are practically extinct as well, exhibiting thus a constant vanishing variance. All the other strains, instead, present qualitatively the same nonmonotonicity with only slight differences (shift of the peaks and different mean values of the CV curves), probably due to the extinction of diatoms and nanoflagellates which causes relevant differences in the system dynamics. More specifically, the system’s response seems to depend on both the noise intensity and the correlation time (see Fig. 3).In this scenario (absence of seasonal driving) we have studied the dependence on both parameters D and (tau) of the probability density functions (PDFs) of the non-vanishing populations. In Fig. 4, the PDFs of bacteria (B1), picophytoplankton (P3), microzooplankton (Z5) and etherotrophic nanoflagellates (Z6) are plotted for (tau =0.5) and eight different values of the parameter D.Figure 4Dependence of the probability density functions of non-vanishing populations on the parameter D for (tau =0.5). The curves are normalized within the interval taken into account. For this reason the relative peaks of the curves in the bottom panels have different values compared to those of the top panels. However, the figure aims at showing the existence of the value of the noise intensity for which the system is more sensitive as well as the generation of a stationary out-of-equilibrium state induced by the noise.Full size imageWe see that the mean value and the variance of these populations are strongly affected by the presence of random fluctuations in the irradiance. Specifically, as the noise intensity increases the mean values of picophytoplankton and bacteria concentrations exhibit a shift. In particular, the results indicate that picophytoplankton is disavantaged by the presence of a noisy component in the irradiance, which indeed tends to inhibit its ability to absorbe the solar light, slowing down its growth. As a consequence, since phytoplankton and bacteria compete for the same resources, as the former declines the latter are favoured, with a compensation mechanism which allows their predators (zooplankton populations) to be almost not affected by the noisy behaviour of the irradiance. Further, we note that for intermediate values of the noise intensity ((D = 10^4 – 10^5)) a maximum of the variance occurs (the PDFs are clearly spread on a wider range of values). Such an effect indicates that the noisy behaviour of irradiance strongly influences the whole ecosystem dynamics. Moreover, the nonmonotonic behaviour of the variance (its PDFs become larger and then tighter again as the noise intensity increases) indicates that the noise pushes the ecosystem away from equilibrium, driving it towards a non-equilibrium steady state. Finally, we note that the nonmonotonic behaviour of CV as a function of the noise intensity remains also in the presence of seasonal driving.Figure 5Coefficient of variation ((CV=sigma / mu)) of nine planktonic populations resulting from numerical integration of model equations plotted versus the considered values of D; different curves correspond to different values of the correlation time (tau).Full size imageFigure 5 shows indeed the nonmonotonic response of the ecosystem to the change of D when the deterministic seasonal cycling of the four environmental parameters (temperature, irradiance, wind and salinity) is present. It is easy to observe that also in this instance the major noise-induced effect appears in nanoflagellates and diatoms with a percent standard deviation of 50(%) and 100(%), respectively. The coalescence of different curves (related to different values of (tau)), as D decreases, is due to the fact that for (D rightarrow 0) the impact of the noise is negligible and the evolution of the system practically resembles the deterministic one. On the contrary, for higher values of D remarkable differences arise and clear peaks of CV appear in the considered range of variation.These plots show that, for a fixed value of (tau), there exists a value of the noise intensity for which the planktonic concentrations are maximally spread around their mean values (corresponding to the maximum value of CV and then of the variance). Moreover, such a nonmonotonic behaviour suggests the presence of a resonance, which can be interpreted as the effect of the interplay between the nonlinearity of the system and the environmental random fluctuations.Also in this case, the interplay between the two parameters D and (tau) in determining and characterizing the dynamics of the ecosystem transparently emerges. The value of D corresponding to the maximum value of CV, indeed, basically depends on the specific value of (tau).Finally, we point out that the different dynamic scenarios identified by the D-(tau) couples can be experienced by the system during the year, since the two parameters may seasonally vary depending on the different weather conditions. In other words, a seasonally varying noise (see Fig. 1c) may cause the nine populations explore different regions of the D-(tau) space during the year. Therefore, the results reported in this paper can highlight the detectable yearly variability of a marine ecosystem which does not stem from the deterministic seasonal variation of environmental parameters.Effects on the organic carbonIn this subsection the effects of the irradiance noise on the biogechemistry are analysed. In Fig. 6 the dependence on (tau) of both the CV [panel (a)] and the mean value concentration [panel (b)] of detritus, labile dissolved organic carbon (L-DOC), semi-labile dissolved organic carbon (SL-DOC) and gross primary production (GPP) are shown. All these biogeochemical properties are correlated with carbon cycling. Gross primary production is related to the amount of carbon entering in the ecosystem, and is related to the maximum energy available in the ecosystem progressively dissipated in the trophic web. Gross primary production is directly affected by light fluctuation and its CV shape is very similar to that of the irradiance, Fig. 2. We selected also detritus and DOC because they are important indicators for the carbon cycling dynamics and are related to the cycling of chemicals like heavy metals62. The different curves, related to different values of D, approach the same (vanishing) value for large (tau). As previously discussed for the CV [Fig. 6(a)] of biomass concentrations, this circumstance is due to the fact that, in this case, the system dynamics tends to the deterministic case, characterized by a unique possible realization implying a vanishing standard deviation. For high correlation times thus the system is insensitive to the noise intensity. On the contrary, for small values of (tau), different values of D lead to significant differences of the variance. In particular, detritus, L-DOC and SL-DOC exhibit a clear non-monotonic behaviour whose maximum value depends on the combined values of D-(tau). Only the GPP presents a decreasing monotonic behaviour.The dependence of the mean value concentration on (tau), instead, is qualitatively the same for all the four parameters. Also in this case we can note a diversification with respect to D occurring at small (tau) and a (deterministic) constant value arising for low (high) values of D ((tau)).These results manifest that not only the population dynamics, but also all the biogeochemical processes are profoundly affected by the presence of stochastic environmental variables. The values and the behaviour of the examined quantities are indeed determined by the intimate interplay between the intensity and the time correlation of the noise fluctuations.Figure 6(a) Coefficient of variation ((CV=sigma / mu)) and (b) mean value concentration ((mu)) of detritus, labile dissolved organic carbon (L-DOC), semi-labile dissolved organic carbon (SL-DOC) and gross primary production (GPP) resulting from numerical integration of model equations plotted versus the considered values of (tau); the different curves are related to different values of the correlation time D.Full size image More

Armbrust, E. V. The life of diatoms in the world’s oceans. Nature 459, 185–192 (2009).CAS 
Article 

Google Scholar 
Mann, D. G. The species concept in diatoms. Phycologia 38, 437–495 (1999).Article 

Google Scholar 
Smetacek, V. Diatoms and the ocean carbon cycle. Protist 150, 25–32 (1999).CAS 
Article 

Google Scholar 
Rynearson, T. A. et al. Major contribution of diatom resting spores to vertical flux in the sub-polar North Atlantic. Deep. Res. Part I Oceanogr. Res. Pap. 82, 60–71 (2013).CAS 
Article 

Google Scholar 
Allen, J. T. et al. Diatom carbon export enhanced by silicate upwelling in the northeast Atlantic. Nature 437, 728–732 (2005).CAS 
Article 

Google Scholar 
Boyd, P. W., Strzepek, R., Fu, F. & Hutchins, D. A. Environmental control of open-ocean phytoplankton groups: Now and in the future. Limnol. Oceanogr. 55, 1353–1376 (2010).CAS 
Article 

Google Scholar 
Hátún, H., Somavilla, R., Rey, F., Johnson, C. & Mathis, M. The subpolar gyre regulates silicate concentrations in the North Atlantic. Sci. Rep. 1–9 https://doi.org/10.1038/s41598-017-14837-4 (2017).Bopp, L. Response of diatoms distribution to global warming and potential implications: a global model study. Geophys. Res. Lett. 32, 2–5 (2005).Article 
CAS 

Google Scholar 
Warner, A. J. & Hays, G. C. Sampling by the Continuous Plankton Recorder survey. Prog. Oceanogr. 6611, 237–256 (1994).Article 

Google Scholar 
Edwards, M., Beaugrand, G., Reid, P. C., Rowden, A. A. & Jones, M. B. Ocean climate anomalies and the ecology of the North Sea. Mar. Ecol. Prog. Ser. 239, 1–10 (2002).Article 

Google Scholar 
Allen, S. et al. Interannual stability of phytoplankton community composition in the North-East atlantic. Mar. Ecol. Prog. Ser. 655, 43–57 (2020).Article 

Google Scholar 
Edwards, M. & Richardson, A. J. Impact of climate change on marine pelagic phenology and trophic mismatch. Nature 430, 881–884 (2004).CAS 
Article 

Google Scholar 
Wihsgott, J. U. et al. Observations of vertical mixing in autumn and its effect on the autumn phytoplankton bloom. Prog. Oceanogr. 177, 1157–1165 (2019).Article 

Google Scholar 
Kamykowski, D. & Zentara, S. J. Predicting plant nutrient concentrations from temperature and sigma-T in the upper kilometer of the world ocean. Deep. Res. Part A-Oceanogr. Res. Pap. 33, 89–105 (1986).CAS 
Article 

Google Scholar 
Kamykowski, D. & Zentara, S. J. Changes in world ocean nitrate availability through the 20th century. Deep. Res. Part I-Oceanogr. Res. Pap. 52, 1719–1744 (2005).Article 

Google Scholar 
Behrenfeld, M. J. et al. Climate-driven trends in contemporary ocean productivity. Nature 444, 752–755 (2006).CAS 
Article 

Google Scholar 
López-Urrutia, A., San Martin, E., Harris, R. P. & Irigoien, X. Scaling the metabolic balance of the oceans. Proc. Natl. Acad. Sci. USA 103, 8739–8744 (2006).Article 
CAS 

Google Scholar 
Richardson, A. J. & Schoeman, D. S. Climate impact on plankton ecosystems in the Northeast Atlantic. Science 305, 1609–1612 (2004).CAS 
Article 

Google Scholar 
Edwards, M., Johns, D., Leterme, S., Svendsen, E. & Richardson, A. Regional climate change and harmful algal blooms in the northeast Atlantic. Limnol. Oceanogr. 51, 820–829 (2006).Article 

Google Scholar 
Hinder, S. L. et al. Changes in marine dinoflagellate and diatom abundance under climate change. Nat. Clim. Chang. 2, 271–275 (2012).Article 

Google Scholar 
Batten, S. et al. CPR sampling: the technical background, materials and methods, consistency and comparability. Prog. Oceanogr. 58, 193–215 (2003).Article 

Google Scholar 
Reid, P. C. et al. The Continuous Plankton Recorder: concepts and history, from plankton indicator to undulating recorders. Prog. Oceanogr. 58, 117–173 (2003).Article 

Google Scholar 
Edwards, M., Beaugrand, G., Hays, G. C., Koslow, J. A. & Richardson, A. J. Multi-decadal oceanic ecological datasets and their application in marine policy and management. Trends Ecol. Evol. 25, 602–610 (2010).Article 

Google Scholar 
Edwards, M. et al. North Atlantic warming over six decades drives decreases in krill abundance with no associated range shift. Commun. Biol. 4, 1–10 (2021).Article 

Google Scholar 
Richardson, A. J. et al. Using continuous plankton recorder data. Prog. Oceanogr. 68, 27–74 (2006).Article 

Google Scholar 
Hélaouët, P., Beaugrand, G. & Reygondeau, G. Reliability of spatial and temporal patterns of C. finmarchicus inferred from the CPR survey. J. Mar. Syst. 153, 18–24 (2016).Article 

Google Scholar 
Owens, N. J. P. et al. All plankton sampling systems underestimate abundance: response to “Continuous plankton recorder underestimates zooplankton abundance” by J.W. Dippner and M. Krause. J. Mar. Syst. 128, 240–242 (2013).Article 

Google Scholar 
Jonas, T. D., Walne, A., Beaugrand, G., Gregory, L. & Hays, G. C. The volume of water filtered by a Continuous Plankton Recorder sample: the effect of ship speed. J. Plankton Res. 26, 1499–1506 (2004).Article 

Google Scholar 
O’Reilly, C. H., Zanna, L. & Woollings, T. Assessing external and internal sources of Atlantic multidecadal variability using models, proxy data, and early instrumental indices. J. Clim 32, 7727–7745 (2019).Article 

Google Scholar 
Qin, M., Dai, A. & Hua, W. Quantifying contributions of internal variability and external forcing to atlantic multidecadal variability since 1870. Geophys. Res. Lett. 47, 1–11 (2020).
Google Scholar 
Mann, M. E., Steinman, B. A. & Miller, S. K. Absence of internal multidecadal and interdecadal oscillations in climate model simulations. Nat. Commun. 11, 1–9 (2020).Article 
CAS 

Google Scholar 
Enfield, D. B., Mestas-Nuñez, A. M. & Trimble, P. J. The Atlantic Multidecadal Oscillation and its relation to rainfall and river flows in the continental U.S. Geophys. Res. Lett. 28, 2077–2080 (2001).Article 

Google Scholar 
Gray, S. T., Graumlich, L. J., Betancourt, J. L. & Pederson, G. T. A tree-ring based reconstruction of the Atlantic Multidecadal Oscillation since 1567 A.D. Geophys. Res. Lett. 31, 2–5 (2004).Article 

Google Scholar 
Edwards, M., Beaugrand, G., Helaouët, P., Alheit, J. & Coombs, S. Marine ecosystem response to the Atlantic Multidecadal Oscillation. PLoS One 8, e57212 (2013).CAS 
Article 

Google Scholar 
Jones, P. D., New, M., Parker, D. E. & Martin, S. & Rigor, I. G. Surface air temperature and its changes over the past 150 years. Rev. Geophys. 37, 173–199 (1999).Article 

Google Scholar 
Beaugrand, G. et al. Reorganization of North Atlantic marine copepod biodiversity and climate. Science296, 1692–1694 (2002).CAS 
Article 

Google Scholar 
Alvain, S., Moulin, C., Dandonneau, Y. & Bréon, F. M. Remote sensing of phytoplankton groups in case 1 waters from global SeaWiFS imagery. Deep Sea Res. Part I Oceanogr. Res. Pap. 52, 1989–2004 (2005).Article 

Google Scholar 
Kalnay, E. et al. The NCEP/NCAR 40-year reanalysis project. Bull. Am. Meteorol. Soc. 77, 437–471 (1996).Article 

Google Scholar 
Huang, B. et al. Extended reconstructed sea surface temperature, version 5 (ERSSTv5): upgrades, validations, and intercomparisons. J. Clim. 30, 8179–8205 (2017).Article 

Google Scholar 
Lam, N. S. N. Spatial interpolation methods: a review. Am. Cartogr. 10, 129–150 (1983).Article 

Google Scholar 
Beaugrand, G., McQuatters-Gollop, A., Edwards, M. & Goberville, E. Long-term responses of North Atlantic calcifying plankton to climate change. Nat. Clim. Chang. 3, 263–267 (2012).Article 
CAS 

Google Scholar 
Beaugrand, G. et al. Prediction of unprecedented biological shifts in the global ocean. Nat. Clim. Chang. 9, 237–243 (2019).Article 

Google Scholar  More

Chapron, G. et al. Recovery of large carnivores in Europe’s modern human-dominated landscapes. Science 346, 1517–1519 (2014).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Lute, M. L., Carter, N. H., López-Bao, J. V. & Linnell, J. D. C. Conservation professionals’ views on governing for coexistence with large carnivores. Biol. Cons. 248, 108668 (2020).Article 

Google Scholar 
Gantchoff, M. G. & Belant, J. L. Regional connectivity for recolonizing American black bears (Ursus americanus) in southcentral USA. Biol. Cons. 214, 66–75 (2017).Article 

Google Scholar 
Ripple, W. J. et al. Status and ecological effects of the world’s largest carnivores. Science 343, 25 (2014).Article 
CAS 

Google Scholar 
Kays, R. et al. Does hunting or hiking affect wildlife communities in protected areas?. J. Appl. Ecol. 54, 242–252 (2017).Article 

Google Scholar 
Schipper, J. et al. The status of the world’s land and marine mammals: diversity, threat, and knowledge. Science 322, 225–230 (2008).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Smith, J. A., Wang, Y. & Wilmers, C. C. Top carnivores increase their kill rates on prey as a response to human-induced fear. Proc. R. Soc. B Biol. Sci. 282, 20142711 (2015).Article 

Google Scholar 
Stillfried, M., Belant, J. L., Svoboda, N. J., Beyer, D. E. & Kramer-Schadt, S. When top predators become prey: Black bears alter movement behaviour in response to hunting pressure. Behav. Proc. 120, 30–39 (2015).Article 

Google Scholar 
Støen, O.-G. et al. Physiological evidence for a human-induced landscape of fear in brown bears (Ursus arctos). Physiol. Behav. 152, 244–248 (2015).PubMed 
Article 
CAS 

Google Scholar 
Evans, M. J., Rittenhouse, T. A. G., Hawley, J. E. & Rego, P. W. Black bear recolonization patterns in a human-dominated landscape vary based on housing: New insights from spatially explicit density models. Landsc. Urban Plan. 162, 13–24 (2017).Article 

Google Scholar 
LaRue, M. A. et al. Cougars are recolonizing the midwest: Analysis of cougar confirmations during 1990–2008. J. Wildl. Manag. 76, 1364–1369 (2012).Article 

Google Scholar 
Cove, M. V., Fergus, C., Lacher, I., Akre, T. & McShea, W. J. Projecting mammal distributions in response to future alternative landscapes in a rapidly transitioning region. Remote Sens. 11, 2482 (2019).ADS 
Article 

Google Scholar 
Frid, A. & Dill, L. Human-caused disturbance stimuli as a form of predation risk. Conserv. Ecol. 6, 25 (2002).
Google Scholar 
Clinchy, M. et al. Fear of the human “super predator” far exceeds the fear of large carnivores in a model mesocarnivore. Behav. Ecol. 27, 1826–1832 (2016).
Google Scholar 
Suraci, J. P., Clinchy, M., Zanette, L. Y. & Wilmers, C. C. Fear of humans as apex predators has landscape-scale impacts from mountain lions to mice. Ecol. Lett. 22, 1578–1586 (2019).PubMed 
Article 

Google Scholar 
Gaynor, K. M., Hojnowski, C. E., Carter, N. H. & Brashares, J. S. The influence of human disturbance on wildlife nocturnality. Science 360, 1232–1235 (2018).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Smith, J. A., Thomas, A. C., Levi, T., Wang, Y. & Wilmers, C. C. Human activity reduces niche partitioning among three widespread mesocarnivores. Oikos 127, 890–901 (2018).Article 

Google Scholar 
Tucker, M. A. et al. Moving in the Anthropocene: Global reductions in terrestrial mammalian movements. Science 359, 466–469 (2018).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Carter, N. H., Brown, D. G., Etter, D. R. & Visser, L. G. American black bear habitat selection in northern Lower Peninsula, Michigan, USA, using discrete-choice modeling. Ursus 21, 57–71 (2010).Article 

Google Scholar 
Naidoo, R. & Burton, A. C. Relative effects of recreational activities on a temperate terrestrial wildlife assemblage. Conserv. Sci. Pract. 2, e271 (2020).
Google Scholar 
Geffroy, B., Samia, D. S. M., Bessa, E. & Blumstein, D. T. How nature-based tourism might increase prey vulnerability to predators. Trends Ecol. Evol. 30, 755–765 (2015).PubMed 
Article 

Google Scholar 
Geffroy, B. et al. Evolutionary dynamics in the Anthropocene: Life history and intensity of human contact shape antipredator responses. PLoS Biol. 18, e3000818 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Beeco, J. A., Hallo, J. C. & Brownlee, M. T. J. GPS visitor tracking and recreation suitability mapping: tools for understanding and managing visitor use. Landsc. Urban Plan. 127, 136–145 (2014).Article 

Google Scholar 
Thorsen, N. H. et al. Smartphone app reveals that lynx avoid human recreationists on local scale, but not home range scale. Sci. Rep. 12, 1–13 (2022).Article 
CAS 

Google Scholar 
Evans, M. J., Hawley, J. E., Rego, P. W. & Rittenhouse, T. A. G. Hourly movement decisions indicate how a large carnivore inhabits developed landscapes. Oecologia 190, 11–23 (2019).ADS 
PubMed 
Article 

Google Scholar 
Carlos, A. W. D., Bright, A. D., Teel, T. L. & Vaske, J. J. Human-black bear conflict in urban areas: an integrated approach to management response. Hum. Dimens. Wildl. 14, 174–184 (2009).Article 

Google Scholar 
Johnson, H. E. et al. Human development and climate affect hibernation in a large carnivore with implications for human–carnivore conflicts. J. Appl. Ecol. 55, 663–672 (2018).Article 

Google Scholar 
Gould, N. P., Powell, R., Olfenbuttel, C. & DePerno, C. S. Growth and reproduction by young urban and rural black bears. J. Mammal. 102, 1165–1173 (2021).Article 

Google Scholar 
Ditmer, M. A., Noyce, K. V., Fieberg, J. R. & Garshelis, D. L. Delineating the ecological and geographic edge of an opportunist: The American black bear exploiting an agricultural landscape. Ecol. Model. 387, 205–219 (2018).Article 

Google Scholar 
McFadden-Hiller, J. E. Jr. & Belant, J. L. Spatial distribution of black bear incident reports in michigan. PLoS One 11, e0154474 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Ladle, A., Steenweg, R., Shepherd, B. & Boyce, M. S. The role of human outdoor recreation in shaping patterns of grizzly bear-black bear co-occurrence. PLoS One 13, e0191730 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Wilbur, R. C., Lischka, S. A., Young, J. R. & Johnson, H. E. Experience, attitudes, and demographic factors influence the probability of reporting human–black bear interactions. Wildl. Soc. Bull. 42, 22–31 (2018).Article 

Google Scholar 
Lustig, E. J., Lyda, S. B., Leslie, D. M., Luttbeg, B. & Fairbanks, W. S. Resource selection by recolonizing American Black Bears. J. Wildl. Manage. 85, 531–542 (2021).Article 

Google Scholar 
Sun, C. C., Fuller, A. K., Hare, M. P. & Hurst, J. E. Evaluating population expansion of black bears using spatial capture-recapture. J. Wildl. Manage. 81, 814–823 (2017).Article 

Google Scholar 
Kautz, T. M. et al. Large carnivore response to human road use suggests a landscape of coexistence. Glob. Ecol. Conserv. 30, e01772 (2021).Article 

Google Scholar 
Michigan Department of Natural Resources (MIDNR) (2021).Blount, J. D., Chynoweth, M. W., Green, A. M. & Şekercioğlu, Ç. H. Review: COVID-19 highlights the importance of camera traps for wildlife conservation research and management. Biol. Cons. 256, 108984 (2021).Article 

Google Scholar 
Weather Atlas. https://www.weather-atlas.com/enEvans, J. S. Spatial Analysis and Modelling Utilities. Package ‘spatialEco’. https://cran.r-project.org/web/packages/spatialEco/spatialEco.pdf (2021).Díaz-Ruiz, F., Caro, J., Delibes-Mateos, M., Arroyo, B. & Ferreras, P. Drivers of red fox (Vulpes vulpes) daily activity: prey availability, human disturbance or habitat structure?. J. Zool. 298, 128–138 (2016).Article 

Google Scholar 
Moore, J. F. et al. Comparison of species richness and detection between line transects, ground camera traps, and arboreal camera traps. Anim. Conserv. 23, 561–572 (2020).Article 

Google Scholar 
Parsons, A. W. et al. Urbanization focuses carnivore activity in remaining natural habitats, increasing species interactions. J. Appl. Ecol. 56, 1894–1904 (2019).Article 

Google Scholar 
Allen, M. L., Sibarani, M. C., Utoyo, L. & Krofel, M. Terrestrial mammal community richness and temporal overlap between tigers and other carnivores in Bukit Barisan Selatan National Park, Sumatra. Anim. Biodiv. Conserv. 43(1), 97–107 (2020).Article 

Google Scholar 
Tian, C. et al. Temporal niche patterns of large mammals in Wanglang National Nature Reserve, China. Glob. Ecol. Conserv. 22, e01015 (2020).Article 

Google Scholar 
Meredith, M. & Ridout, M. Estimates of coefficient of overlapping for animal activity patterns. Package ‘overlap’. https://cran.r-project.org/web/packages/overlap/overlap.pdf (2020).RStudio Team. RStudio: Integrated Development for R. RStudio, PBC, Boston, MA. http://www.rstudio.com/ (2021).Ridout, M. S. & Linkie, M. Estimating overlap of daily activity patterns from camera trap data. JABES 14, 322–337 (2009).MathSciNet 
MATH 
Article 

Google Scholar 
Lashley, M. A. et al. Estimating wildlife activity curves: comparison of methods and sample size. Sci. Rep. 8, 1–11 (2018).CAS 
Article 

Google Scholar 
Rowcliffe, M. Animal Activity Statistics. Package ‘activity’. https://cran.r-project.org/web/packages/activity/activity.pdf (2021).MacKenzie, D. I., Nichols, J. D., Hines, J. E., Knutson, M. G. & Franklin, A. B. Estimating site occupancy, colonization, and local extinction when a species is detected imperfectly. Ecology 84, 2200–2207 (2003).Article 

Google Scholar 
Wei, T., & Simko, V. Visualization of a Correlation Matrix. Package ‘corrplot’. https://cran.r-project.org/web/packages/corrplot/corrplot.pdf (2017).Norton, D. C. et al. Female American black bears do not alter space use or movements to reduce infanticide risk. PLoS One 13, e0203651 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Ditmer, M. A. et al. Behavioral and physiological responses of American black bears to landscape features within an agricultural region. Ecosphere 6, 1–21 (2015).Article 

Google Scholar 
Clark, D. et al. Using machine learning methods to predict the movement trajectories of the Louisiana black bear. SMU Data Sci. Rev. 5, 25 (2021).
Google Scholar  More

Growth conditions and plant materialsTwo experiments were conducted concurrently (sites A and B) in the same screen house in 2019 between the months of May and July at the Landmark University Greenhouse and Hydroponic Technology Center, a section of the Teaching and Research Farm of the University in Omu-Aran, Kwara State Nigeria. Experiment at site B was conducted simultaneously as A so as to validate the results of experiment A. Landmark University lies within Latitude 8° 7′ 26.21388″ and 5° 5′ 0.1788″. Both experiments (A & B) involved tomato (Solanum lycopersicum L. variety cherry) and cucumber (Cucumis sativus L. variety marketer) crops. For each crop, seeds were sown into a separate seed tray filled with coco peat (Coco peat, SRIMATHI EXPORT, INDIA). Cocopeat is the mesocarp tissue or husk after the grinding of coconut fruit. It has a lightweight and high water and nutrient holding capacities, it has an acceptable pH, electrical conductivity, and other chemical attributes27. Rice husk is the by-product of rice after milling. The rice husk used was collected from the rice processing mill of Landmark University. Rice husk is a highly porous and light weighted material with a very high specific area28.Two sets of seed trays (one for organic and another for inorganic fertilizers) were used each for tomato and cucumber crops in the nursery. Both were raised in the nursery for two weeks before transplanting. Black grow bags (30 × 17 cm) filled with a coco peat/rice husk (1:4 ratio by volume) mixture with a weight of about 10 kg were arranged in a screen house. Both the nursery and establishment of crop proper take place in a screen house. The screen house has a galvanized iron as the frame, a UV covering on top, side net for screening insect pests the floor fairly covered with granite. Temperature and relative humidity within the screen house during the period of the experiment was monitored using a Thermograph and a Barograph, and they were at an average of 31 °C and 75%, respectively.The grow bags were randomly placed in the screen house for the unbiased application of amendments. For both tomato and cucumber crops, the treatment comprised of six (6) levels of liquid organic fertilizer (5, 15, 25, 35, 45, 55 mL), in-organic fertilizer, and a control (ordinary borehole water). Levels of organic fertilizers were selected based on the recommendation of 20 mL of liquid organic fertilizer by29. The eight (8) treatments both for tomato and cucumber were arranged in a Completely Randomized Design replicated three times. One healthy plant was maintained per grow bag and four grow bags represent a treatment and there were 32 plants per block each for tomato and cucumber. For both crops, the experiment lasted for 90 days.Organic and in-organic nutrient solutionsThe liquid organic fertilizer used was obtained from the biomass of Mexican sunflower (Tithonia diversifolia). Fresh biomass (mainly leaves and stems) of the plant was collected from the Teaching and Research Farms of Landmark University, Nigeria. After rinsing, they were cut with a sterile knife into pieces of ≤ 1 cm size. A sample was taken for initial physicochemical analyses by grinding in a sterile mortal, diluted with sterile water and analyzed. The biomass was then soaked in sterile water inside a clean container, and allowed to ferment spontaneously for a period of 14 days. During the fermentation, samples were taken every 4 days for microbial analyses of the major players during the fermentation. At the end of fermentation, the mixture was separated using a sieve of mesh size ≤ 2 mm. The liquid portion was then refrigerated prior to the planting regime while another sample was taken to ascertain the physicochemical and microbial qualities of the produced liquid fertilizer. The chemical analysis is presented in Table 4. For inorganic fertilizer, Water soluble fertilizers employed in hydroponics were used (Hydroponics fertilizer, Anmol chemicals, India); calcium nitrate 650 mg L−1, potassium nitrate 450 mg L−1, magnesium 400 mg L−1, chelate 20 mg L−1, mono-ammonium phosphate 400 mg L−1. The electrical conductivity (EC) of the solution was 1.9 dS m-1.Irrigation and fertigationThe tomato and cucumber plants were fertigated morning and evening daily for one hour on each occasion according to the treatments. Preparation of the nutrient solution was with borehole water and was supplied to plants by an online pressure drip irrigation system set at 2.0 L h-1 using an arrowhead on each tomato and cucumber plant. Different tanks (250 L) were installed according to the various treatments making a total of 8 tanks. The organic fertilizer was diluted according to the various treatments equivalent to 1.25, 3.75, 6.25, 8.75, 11.25, and 13.75 L per 250 L of water respectively for 5, 15, 25, 35, 45, and 55 mL treatments. The nutrient solutions were refilled when the consumption is less than 20% of the initial volume (250 L) in the tank. One day per week, crops were irrigated with ordinary water to wash out pipes and prevent deposits of salts. The same concentration of nutrient was used from transplanting to the termination of the study for both tomato and cucumber crops, however, at the flowering of the crops, the volume of fertigation was increased to 3.0 L h-1 to be able to cope with the size of the plants.Trellising, pest and diseases controlFor both tomato and cucumber crops, plant vines were supported by twisting them around a wire that is- attached to the roof of the screen house and 2 m from the ground. Lateral outgrowths were cut off every week to ensure a sturdy single stem. Pests and diseases were scouted every day. Whiteflies, aphids, and other insects were controlled with orizon (Producer, location of producer) (active ingredient, acetamiprid, and abamectin) using 0.133% v/v. Fungi were controlled using ridomil gold (Producer, Location of producer) at 2% w/v.Determination of growth and yield of tomato and cucumberThree tomato and cucumber plants were randomly selected for each treatment for the determination of growth parameters (plant height, leaf area, number of leaves per plant, and stem diameter) at mid the flowering stage of tomato and cucumber plants.The leaf area of tomato was calculated using the model (A = KL2) developed by Lyon30, where L = Length of tomato leaf, K = constant which is 0.1551, and A = leaf area of tomato. Similarly, the leaf area of cucumber was calculated using A = 0.88LW – 4.27, where L = cucumber leaf length and W = cucumber leaf width, A = leaf area of cucumber31.Tomato fruits were ready for harvest from 65 days after transplanting, harvestings were done twice every week (Mondays and Fridays) for up to 85 days after transplanting. Similarly, harvesting of cucumber fruits started 35 days after transplanting and harvestings were also done twice a week (Mondays and Fridays), harvesting was carried out till 60 days after transplanting. Tomato and cucumber fruit yields were counted and weighed at each harvest.Analysis of tomato and cucumber leaves and fruitsAt the 50% flowering stage of tomato and cucumber plants, ten leaf samples were collected from each treatment. The leaf samples were oven-dried at 75 °C for 24 h and thereafter grounded. The grounded samples were later analyzed for nitrogen (N), phosphorous (P), potassium (K), calcium (Ca), and magnesium (Mg) content using the method of described by32. At harvest, four matured tomato and cucumber fruits of uniform size were selected per treatment, and their nutrient compositions were determined using the method of33.Statistical analysisAll data collected on the growth, yield, leaf, and fruit nutrient contents of tomato and cucumber were subjected to analysis of variance (ANOVA). The SPSS V 21.0 (New York, USA) software was used to perform ANOVA and Duncan’s multiple range test (DMRT) was used to compare means at a 5% probability level.
Ethical approvalI confirm that all the research meets ethical guidelines and adheres to the legal requirements of the study country.Compliance with international, national and/or institutional guidelinesExperimental research (either cultivated or wild), comply with relevant institutional, national, and international guidelines and legislation. Experimental studies were carried out in accordance with relevant institutional, national or international guidelines or regulation. More