Philippa Kaur

Model validationBased on the time-series validations of water temperature and DO concentration, model accuracy improved gradually, despite several discrepancies at the beginning of the simulation (Supplementary Fig. S1). The model is primarily driven by a set of boundary data, including wind speed, solar radiation, and precipitation data24,25. From this perspective, more high-quality boundary data promotes better numerical reproducibility. However, meteorological data collection was challenging due to the early observation equipment limitations and low observational accuracy compared to current data. The temporal inconsistency of accuracy in observational data has been eliminated to a large extent by fitting a regression curve24. Spatial resolution is the other issue. Possessing spatially constant values for all boundary conditions complicates the numerical reproducibility of variations on finer scales.The relationship between turnovers and the curve shape of water temperature versus DO concentration is theoretically sound27,28. In the last stage of stratification in the lake, water temperature and DO concentration near the bottom are more likely to slightly increase due to thermal diffusion and DO supplies from the upper water. If a turnover occurs, the whole column of water is mixed strongly (Supplementary Fig. S3). Bottom water temperature decreases due to surface water cooling, and DO concentration increases, due to surface water replenishment and increased oxygen solubility. If the turnover fails, only the partial column of water is mixed, causing a delay in the timing of deep-water renewal (Supplementary Fig. S3). However, the upper water in later months, like that in March, has been rapidly warmed, resulting in an increase in the bottom water temperature. For example, in 2007 and 2016, the simulated water temperature and DO concentration fluctuated within a limited range in February and then skyrocketed in March, after mixing with the warmed surface water (blue points in Supplementary Fig. S4). On the other hand, explicit definitions of turnover timing are challenging. The threshold used to judge turnover timing is reliable because the results matched the observation. The turnover timing varied by 36 days in Lake Biwa during the simulation period, which is comparable to that observed in other lakes, such as approximately 21 days in Heiligensee, Germany over a 17-year timespan29, 16 days in Lake Washington over a 40-year timespan30, and 28 days in Blelham Tarn over a 41-year timespan31.Variables affecting the mixing regimeDetermining variables that affect the mixing regime is essential to improve understanding and enable future projections16,17,18. Air temperature, wind speed, cloud cover, precipitation, water density, and lake transparency are all potential variables. We, here, compared the above variables to the turnover timing in Lake Biwa. The meteorological inputs in this study provided data for air temperature, wind speed, cloud cover, and precipitation. Water density and particulate organic carbon (POC) concentration representing lake transparency were the model’s outputs. The annual averages and cold season (November–April) values of the above variables were calculated over the simulation period (Supplementary Fig. S6). Annual averages illustrate general long-term warming trends18, while cold season values particularly determine the timing of turnover17. However, in Lake Biwa, air temperature during the cold season fluctuated greatly compared to the annual averages. A random forest analysis17 has been conducted between the turnover timing and the above two variable sets (cold season values versus annual averages) in Lake Biwa, and the cold season values better explained the turnover timing (35.39% versus 18.48%). The results agree with the conclusion drawn from the previous sensitivity tests, which indicated the relative importance of air temperature and solar radiation during winter based on 40 scenarios32.The importance of variables was estimated based on the random forest analysis using the cold season data (Fig. 4a). Wind speed dominates the timing of turnover, which is consistent with the previous studies17,25. The POC concentration, the difference in water density between the surface and bottom, and cloud cover have moderate effects on the timing of turnover. However, air temperature is less important, which is contrary to the turnover mechanism17,24,32. A re-confirmation was conducted of the relationship between turnover timing and air temperature (Fig. 4b and Supplementary Fig. S7). The cool air generally encourages an early turnover, albeit with several anomalies. The turnover timing between 1976 and 1990 remained constant independent of climate change, and the period coincidently had a substantial nutrient fluctuation (Fig. 3). As a result, it is essential to investigate the nutrient status further.Figure 4Analysis results of the relationship between potential variables and turnover timing: (a) the importance of variables importance using a random forest analysis, and (b) the relationship between the cold season air temperature and the timing of turnover. Variable importance is calculated using the percentage increase in mean square error (MSE) and the increase in node purity. Higher values illustrate the greater importance of the variable. Variables include air temperature (AT), precipitation (pptn.), cloud cover (CC), the difference in density (DD), POC, and wind speed (WS).Full size imageLake nutrient concentrationsBecause phosphorus is the limiting nutrient in Lake Biwa and DIP concentrations can be effectively limited by regulating external loadings as practiced (Fig. 3), DIP concentrations become the focus of this discussion for nutrient status. However, the DIP concentrations disproportionately responded to the external loadings of total phosphorus (TP) in Lake Biwa. Although external TP loading itself fails to determine lake phosphorus concentrations due to the hydrodynamics of lakes33, Lake Biwa exhibited insignificant changes in the inflow rate or the retention time (and see an example of the surface flow in Supplementary Fig. S8). Therefore, it can be assumed that the hydraulic loading remained constant, and the input nutrient concentrations were proportionate to the external nutrient loadings in Lake Biwa. This finding contradicts a recent meta-analysis that highlighted a deterministic relationship between input nutrient concentrations and lake nutrient concentrations, based on steady-state mass balance models6. The possible reason is the dynamics of the lake’s ecosystem22, which have been considered in this study. For example, the surface DIP concentrations were almost nonexistent regardless of the external TP loadings in Lake Biwa, supporting that phosphorus is the limiting nutrient in Lake Biwa34,35. The low DIP concentrations at the surface may be caused by the rapid recycling of phosphorus because the amount of phosphorus available for phytoplankton is easily affected by the feedback mechanism between phytoplankton photosynthesis and the phosphorus released from the water35,36.Hypoxia and strategiesThe variations in DO concentration are the public’s top concern as it relates to hypoxia, a key indicator of water quality. Lake bottom, among all water depths, is more sensitive to small changes in oxygen conditions12. In Lake Biwa, the annual minimum DO concentrations ranged from 2 to 5.5 mg/L over the last 60 years (Supplementary Fig. S9). The decrease in DO concentrations in the early period, typically till the 1980s, was mainly caused by nutrient enrichments (Fig. 3). The nutrient enrichment-induced heavy eutrophication eventually accelerates the rate of DO depletion2. After eutrophication was controlled in the 1980s, climate change became the dominant stressor23. There remains much uncertainty surrounding the relationship between climatic variables-related turnover timing and hypoxia in Lake Biwa12. We, therefore, first investigate the relationship between hypoxia and turnover timing, and then concentrate on nutrients to alleviate hypoxia.Although the relationship between turnover timing and DO concentrations is quite weak (R2 = 0.10), there is a general decrease in DO concentrations with increasing turnover timing (Fig. 5a). On the other hand, a linear relationship has been found between DIP concentrations and DO concentrations, with an R2 of 0.67 (Fig. 5b). The slope of –0.841 μgP/mgDO means an increase in DIP concentrations by approximately 0.841 μgP/L causes a decrease in DO concentrations by 1 mg/L. Note that the simulation results were compared over the whole period, and eutrophication-induced hypoxia differs theoretically from climate-induced hypoxia. Additional testing has been conducted to distinguish the effects of two stressors (eutrophication- and climate-induced hypoxia; Supplementary Fig. S10). Before 1980 when eutrophication progressed, the annual minimum DO concentrations and the DIP concentrations had a stronger linear relationship (R2 = 0.89). Although waste-water treatment has improved conditions in the lake, climate change induced alteration of turnover timing may adversely influence water quality. However, the relationship weakened dramatically with an R2 of 0.10 after 1980, when climate change dominated hypoxia. The lower R2 value indicates that climate-related hypoxia is more complex as concluded previously37,38. The two possibilities are as follows. First, there can be a legacy of hypoxia related to eutrophication. The DO recovery at the bottom of Lake Biwa was complicated by the low DO concentration in 1980 and the delayed timing of turnover; similar phenomena have been observed in the Lake of Zurich22. Second, ecosystem dynamics could help explain the difficulty in predicting hypoxia at the bottom. Phytoplankton fully exploits phosphorus at the surface, as explained above, then the death and sinking of the surface phytoplankton are accompanied by the sedimentation of phosphorus to the bottom as modeled. Bacteria break down the sinking phytoplankton, releasing phosphorus and consuming DO in the process. Additional DO consumption lowers the bottom DO concentration, which in turn encourages phosphorus release from the sediment in a low DO environment22. Such unfavorable feedback between DIP and DO concentrations are strengthened by prolonged stratification and eventually accelerates the development of hypoxia. However, future research is necessary because this numerical model simplified the relationship between water and sediment. The sinking of organic carbon into sediment is integrated in the model, and due to the decomposition of organic carbon in the sediment, nutrients are released into and oxygen is depleted in the water. Despite that, the trends between DO and DIP concentrations stay the same under climate change (Fig. 5b), and thus controlling lake phosphorus is beneficial to the Lake Biwa hypoxia.Figure 5The linear regression results of the relationship: (a) between turnover timing and annual minimum concentration of DO, (b) between the annual minimum concentration of DO and annual average concentration of DIP. The simulation results at the monitoring station were used for analysis.Full size image More

Kazarov, E. The Role of Zoos in Creating a Conservation Ethic in Visitors. SIT Digital Collections (2022). at https://digitalcollections.sit.edu/isp_collection/584.Hosey, G. How does the zoo environment affect the behaviour of captive primates?. Appl. Anim. Behav. Sci. 90, 107–129 (2005).
Google Scholar 
Morgan, K. & Tromborg, C. Sources of stress in captivity. Appl. Anim. Behav. Sci. 102, 262–302 (2007).
Google Scholar 
Sherwen, S. & Hemsworth, P. The visitor effect on zoo animals: Implications and opportunities for zoo animal welfare. Animals 9, 366 (2019).PubMed Central 

Google Scholar 
Chamove, A., Hosey, G. & Schaetzel, P. Visitors excite primates in zoos. Zoo Biol. 7, 359–369 (1988).
Google Scholar 
Tetley, C. L. & O’Hara, S. J. Ratings of animal personality as a tool for improving the breeding, management and welfare of zoo mammals. Anim. Welf. UFAW J. 21(4), 463 (2012).CAS 

Google Scholar 
Stoinski, T. S., Jaicks, H. F. & Drayton, L. A. Visitor effects on the behavior of captive western lowland gorillas: The importance of individual differences in examining welfare. Zoo Biol. 31(5), 586–599 (2012).PubMed 

Google Scholar 
Queiroz, M. B. & Young, R. J. The different physical and behavioural characteristics of zoo mammals that influence their response to visitors. Animals 8(8), 139 (2018).PubMed Central 

Google Scholar 
Fanson, K. V. & Wielebnowski, N. C. Effect of housing and husbandry practices on adrenocortical activity in captive Canada lynx (Lynx canadensis). Anim. Welf. 22, 159–165 (2013).CAS 

Google Scholar 
Pirovino, M. et al. Fecal glucocorticoid measurements and their relation to rearing, behavior, and environmental factors in the population of pileated gibbons (Hylobates pileatus) held in European zoos. Int. J. Primatol. 32(5), 1161–1178 (2011).
Google Scholar 
Williams, I., Hoppitt, W. & Grant, R. The effect of auditory enrichment, rearing method and social environment on the behavior of zoo-housed psittacines (Aves: Psittaciformes); implications for welfare. Appl. Anim. Behav. Sci. 186, 85–92 (2017).
Google Scholar 
Fernandez, E., Tamborski, M., Pickens, S. & Timberlake, W. Animal–visitor interactions in the modern zoo: Conflicts and interventions. Appl. Anim. Behav. Sci. 120, 1–8 (2009).
Google Scholar 
Hosey, G. & Skyner, L. Self-injurious behavior in zoo primates. Int. J. Primatol. 28, 1431–1437 (2007).
Google Scholar 
Mallapur, A., Sinha, A. & Waran, N. Influence of visitor presence on the behaviour of captive lion-tailed macaques (Macaca silenus) housed in Indian zoos. Appl. Anim. Behav. Sci. 94, 341–352 (2005).
Google Scholar 
Davey, G. Visitors’ Effects on the Welfare of Animals in the Zoo: A Review. J. Appl. Anim. Welf. Sci. 10, 169–183 (2007).CAS 
PubMed 

Google Scholar 
Jones, H., McGregor, P., Farmer, H. & Baker, K. The influence of visitor interaction on the behavior of captive crowned lemurs (Eulemur coronatus) and implications for welfare. Zoo Biol. 35, 222–227 (2016).CAS 
PubMed 

Google Scholar 
Cook, S. & Hosey, G. R. Interaction sequences between chimpanzees and human visitors at the zoo. Zoo Biol. 14(5), 431–440 (1995).
Google Scholar 
Baker, K. C. Benefits of positive human interaction for socially-housed chimpanzees. Anim. Welf. (South Mimms, Engl.nd) 13(2), 239 (2004).CAS 

Google Scholar 
Carder, G. & Semple, S. Visitor effects on anxiety in two captive groups of western lowland gorillas. Appl. Anim. Behav. Sci. 115, 211–220 (2008).
Google Scholar 
Wood, W. Interactions among environmental enrichment, viewing crowds, and zoo chimpanzees (Pantroglodytes). Zoo Biol. 17, 211–230 (1998).
Google Scholar 
Todd, P., Macdonald, C. & Coleman, D. Visitor-associated variation in captive Diana monkey (Cercopithecus diana diana) behaviour. Appl. Anim. Behav. Sci. 107, 162–165 (2007).
Google Scholar 
Davis, N., Schaffner, C. & Smith, T. Evidence that zoo visitors influence HPA activity in spider monkeys (Ateles geoffroyii rufiventris). Appl. Anim. Behav. Sci. 90, 131–141 (2005).
Google Scholar 
Sherwen, S. L. et al. Effects of visual contact with zoo visitors on black-capped capuchin welfare. Appl. Anim. Behav. Sci. 167, 65–73 (2015).
Google Scholar 
Choo, Y., Todd, P. & Li, D. Visitor effects on zoo orangutans in two novel, naturalistic enclosures. Appl. Anim. Behav. Sci. 133, 78–86 (2011).
Google Scholar 
Sherwen, S., Magrath, M., Butler, K., Phillips, C. & Hemsworth, P. A multi-enclosure study investigating the behavioural response of meerkats to zoo visitors. Appl. Anim. Behav. Sci. 156, 70–77 (2014).
Google Scholar 
Hosey, G. & Druck, P. The influence of zoo visitors on the behaviour of captive primates. Appl. Anim. Behav. Sci. 18, 19–29 (1987).
Google Scholar 
Mitchell, G. et al. More on the ‘influence’of zoo visitors on the behaviour of captive primates. Appl. Anim. Behav. Sci. 35(2), 189–198 (1992).
Google Scholar 
Sellinger, R. & Ha, J. The effects of visitor density and intensity on the behavior of two captive jaguars (Panthera onca). J. Appl. Anim. Welfare Sci. 8, 233–244 (2005).CAS 

Google Scholar 
Azevedo, C., Lima, M., Silva, V., Young, R. & Rodrigues, M. Visitor Influence on the Behavior of Captive Greater Rheas (Rhea americana, Rheidae Aves). J. Appl. Anim. Welfare Sci. 15, 113–125 (2012).
Google Scholar 
Das Gupta, M., Das, A., Sumy, M. C. & Islam, M. M. An explorative study on visitor’s behaviour and their effect on the behaviour of primates at Chittagong zoo. Bangladesh J. Vet. Anim. Sci. 5(2), 24–32 (2017).
Google Scholar 
Hemsworth, P. Human–animal interactions in livestock production. Appl. Anim. Behav. Sci. 81, 185–198 (2003).
Google Scholar 
Stoinski, T., Czekala, N., Lukas, K. & Maple, T. Urinary androgen and corticoid levels in captive, male Western lowland gorillas (Gorilla g. gorilla): Age- and social group-related differences. Am. J. Primatol. 56, 73–87 (2002).CAS 
PubMed 

Google Scholar 
Stoinski, T., Lukas, K., Kuhar, C. & Maple, T. Factors influencing the formation and maintenance of all-male gorilla groups in captivity. Zoo Biol. 23, 189–203 (2004).
Google Scholar 
Olsson, I. & Westlund, K. More than numbers matter: The effect of social factors on behaviour and welfare of laboratory rodents and non-human primates. Appl. Anim. Behav. Sci. 103, 229–254 (2007).
Google Scholar 
Martin, J. E. Early life experiences: Activity levels and abnormal behaviours in resocialised chimpanzees. Anim Welf. 11(4), 419–436 (2002).CAS 

Google Scholar 
Birkett, L. P. & Newton-Fisher, N. E. How abnormal is the behaviour of captive, zoo-living chimpanzees?. PLoS ONE 6(6), e20101 (2011).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ballen, C., Shine, R. & Olsson, M. Effects of early social isolation on the behaviour and performance of juvenile lizards Chamaeleo calyptratus. Anim. Behav. 88, 1–6 (2014).
Google Scholar 
Coe, C., Mendoza, S., Smotherman, W. & Levine, S. Mother-infant attachment in the squirrel monkey: Adrenal response to separation. Behav. Biol. 22, 256–263 (1978).CAS 
PubMed 

Google Scholar 
Mendoza, S., Smotherman, W., Miner, M., Kaplan, J. & Levine, S. Pituitary-adrenal response to separation in mother and infant squirrel monkeys. Dev. Psychobiol. 11, 169–175 (1978).CAS 
PubMed 

Google Scholar 
Gilbert, M. & Baker, K. Social buffering in adult male rhesus macaques (Macaca mulatta): Effects of stressful events in single vs. pair housing. J. Med. Primatol. 40, 71–78 (2010).PubMed 
PubMed Central 

Google Scholar 
Schapiro, S. Effects of social manipulations and environmental enrichment on behavior and cell-mediated immune responses in rhesus macaques. Pharmacol. Biochem. Behav. 73, 271–278 (2002).CAS 
PubMed 

Google Scholar 
Chen, W. et al. Effects of social isolation and re-socialization on cognition and ADAR1 (p110) expression in mice. PeerJ 4, e2306 (2016).PubMed 
PubMed Central 

Google Scholar 
Glatston, A., Geilvoet-Soeteman, E., Hora-Pecek, E. & Van Hooff, J. The influence of the zoo environment on social behavior of groups of cotton-topped tamarins Saguinus oedipus oedipus. Zoo Biol. 3, 241–253 (1984).
Google Scholar 
Mitchell, G. et al. Effects of visitors and cage changes on the behaviors of mangabeys. Zoo Biol. 10, 417–423 (1991).
Google Scholar 
Geissmann, T. & Orgeldinger, M. The relationship between duet songs and pair bonds in siamangs Hylobates syndactylus. Anim. Behav. 60, 805–809 (2000).CAS 
PubMed 

Google Scholar 
Palombit, R. Pair bonds in monogamous apes: A comparison of the siamang hylobates syndactylus and the white-handed gibbon hylobates lar. Behaviour 133, 321–356 (1996).
Google Scholar 
Rutberg, A. The evolution of monogamy in primates. J. Theor. Biol. 104, 93–112 (1983).ADS 
CAS 
PubMed 

Google Scholar 
Giorgi, A., Montebovi, G., Vitale, A. & Alleva, E. A behavioural case study of early social isolation of a subadult white-handed gibbon (Hylobates lar). Folia Primatol. 89, 287–294 (2018).
Google Scholar 
Skynner, L. A., Amory, J. R. & Hosey, G. The effect of visitors on the self-injurious behaviour of a male pileated gibbon (Hylobates pileatus). Zool. Garten 74(1), 38–41 (2004).
Google Scholar 
Smith, K. & Kuhar, C. Siamangs (Hylobates syndactylus) and white-cheeked gibbons (Hylobates leucogenys) show few behavioral differences related to zoo attendance. J. Appl. Anim. Welfare Sci. 13, 154–163 (2010).CAS 

Google Scholar 
Lukas, K. E. et al. Longitudinal study of delayed reproductive success in a pair of white-cheeked gibbons (Hylobates leucogenys). Zoo Biol. 21, 413–434 (2002).
Google Scholar 
Cooke, C. & Schillaci, M. Behavioral responses to the zoo environment by white handed gibbons. Appl. Anim. Behav. Sci. 106, 125–133 (2007).
Google Scholar 
Mootnick, A. & Baker, E. Masturbation in captiveHylobates (gibbons). Zoo Biol. 13, 345–353 (1994).
Google Scholar 
Geissmann, T. Reassessment of age of sexual maturity in gibbons (hylobates spp.). American Journal of Primatology 23, 11–22 (1991).Altmann, J. Observational study of behavior: Sampling methods. Behaviour 49(3–4), 227–266 (1974).CAS 
PubMed 

Google Scholar 
Pomerantz, O. & Terkel, J. Effects of positive reinforcement training techniques on the psychological welfare of zoo-housed chimpanzees (Pan troglodytes). Am. J. Primatol. 71, 687–695 (2009).PubMed 

Google Scholar 
Orgeldinger, M. Protective and territorial behavior in captive siamangs (Hylobates syndactylus). Zoo Biol. 16, 309–325 (1997).
Google Scholar 
Fox, J. et al. Package ‘car’. Vienna: R Foundation for Statistical Computing, 16 https://cran.uni-muenster.de/web/packages/car/car.pdf (2012).Magnusson, A., Skaug, H., Nielsen, A., Berg, C., Kristensen, K., Maechler, M., van Bentham, K., Bolker, B., Brooks, M. & Brooks, M. M. Package ‘glmmtmb’. R Package Version 0.2. 0 (2017).Hartig, F., & Hartig, M. F. Package ‘DHARMa’. Vienna, Austria: R Development Core Team (2017).Troisi, A. Displacement activities as a behavioral measure of stress in nonhuman primates and human subjects. Stress 5, 47–54 (2002).PubMed 

Google Scholar 
Baker, K. & Aureli, F. Behavioural indicators of anxiety: An empirical test in chimpanzees. Behaviour 134, 1031–1050 (1997).
Google Scholar 
Vick, S. J. & Paukner, A. Variation and context of yawns in captive chimpanzees (Pan troglodytes). Am. J. Primatol. Off. J. Am. Soc. Primatol. 72(3), 262–269 (2010).
Google Scholar 
Norscia, I. & Palagi, E. When play is a family business: Adult play, hierarchy, and possible stress reduction in common marmosets. Primates 52, 101–104 (2010).PubMed 

Google Scholar 
Held, S. & Špinka, M. Animal play and animal welfare. Anim. Behav. 81, 891–899 (2011).
Google Scholar 
Davey, G. Visitor behavior in zoos: A review. Anthrozoös 19, 143–157 (2006).
Google Scholar 
Nimon, A. & Dalziel, F. Cross-species interaction and communication: a study method applied to captive siamang (Hylobates syndactylus) and long-billed corella (Cacatua tenuirostris) contacts with humans. Appl. Anim. Behav. Sci. 33, 261–272 (1992).
Google Scholar 
Suomi, S. Early determinants of behaviour: Evidence from primate studies. Br. Med. Bull. 53, 170–184 (1997).CAS 
PubMed 

Google Scholar 
Anderson, J. & Chamove, A. Self-aggression and social aggression in laboratory-reared macaques. J. Abnorm. Psychol. 89, 539–550 (1980).CAS 
PubMed 

Google Scholar 
Mallapur, A. & Choudhury, B. Behavioral abnormalities in captive nonhuman primates. J. Appl. Anim. Welfare Sci. 6, 275–284 (2003).CAS 

Google Scholar 
Barlow, C., Caldwell, C. & Lee, P. Individual differences and response to visitors in zoo-housed diana monkeys (Cercopithecus diana diana). Cabdirect.org (2022). at https://www.cabdirect.org/cabdirect/abstract/20123180753.Gartner, M. & Weiss, A. Studying primate personality in zoos: Implications for the management, welfare and conservation of great apes. International Zoo Yearbook 52, 79–91 (2018).
Google Scholar 
Mitchell, G., Raymond, E., Ruppenthal, G. & Harlow, H. Long-term effects of total social isolation upon behavior of rhesus monkeys. Psychol. Rep. 18, 567–580 (1966).
Google Scholar 
Martín, O., Vinyoles, D., García-Galea, E. & Maté, C. Improving the welfare of a zoo-housed male drill (Mandrillus leucophaeus poensis) aggressive toward visitors. J. Appl. Anim. Welfare Sci. 19, 323–334 (2016).
Google Scholar 
Ross, S., Melber, L., Gillespie, K. & Lukas, K. The impact of a modern, naturalistic exhibit design on visitor behavior: A cross-facility comparison. Visitor Stud. 15, 3–15 (2012).
Google Scholar 
Quadros, S., Goulart, V., Passos, L., Vecci, M. & Young, R. Zoo visitor effect on mammal behaviour: Does noise matter?. Appl. Anim. Behav. Sci. 156, 78–84 (2014).
Google Scholar 
Bonnie, K., Ang, M. & Ross, S. Effects of crowd size on exhibit use by and behavior of chimpanzees (Pan troglodytes) and Western lowland gorillas (Gorilla gorilla) at a zoo. Appl. Anim. Behav. Sci. 178, 102–110 (2016).
Google Scholar  More

We investigated whether ricefield mosquito larval control and/or rice cultivation practices are associated with malaria vector densities through a systematic review and meta-analysis. Forty-seven experimental studies were eligible for inclusion in the qualitative analysis and thirty-three studies were eligible for the meta-analysis. It was demonstrated that the use of fish, chemical and biological larvicides in rice fields were effective in controlling larval malaria vector densities at all developmental stages. Intermittent irrigation, however, could only significantly reduce late-stage larvae. Based on a limited number of studies, meta-analyses on other forms of larval control such as monomolecular surface films (MSFs), neem, copepods and Azolla failed to demonstrate any consistent reduction in anopheline numbers. Similarly, rice cultivation practices such as plant variety and density, type of levelling and pesticide application were not generally associated with reduced malaria vectors. Nonetheless, in one study, minimal tillage was observed to reduce average numbers of larvae throughout a cropping season. In another study, herbicide application increased larval abundance over a 4-week period, as did one-time drainage in a third study.
Despite their different modes of action, the use of chemical and bacterial larvicides and MSFs were all relatively effective measures of larval control in rice fields, varying between a 57% to 76% reduction in vector abundance compared to no larviciding. Their effects were highest (often reaching 100% reduction) only shortly following application but did not persist for longer than two weeks. These larvicides mostly had short residual half-lives because they were applied to paddy water which was naturally not completely stagnant: there was a small but constant process of water loss (through drainage, evapotranspiration and percolation) and replacement through irrigation. Hence, even with a residual formulation, weekly re-application would be needed for sustained control47,40,41,50. This would be very labour- and cost-intensive to scale-up, to ensure that larvicides are evenly distributed across vast areas (even at plot/sub-plot level) throughout at least one 5-month long rice-growing season per year42,51. Aerial application (including unmanned aerial vehicles), although widely used in the US and Europe, is unlikely to be a feasible delivery system for smallholders in SSA, even in large irrigation schemes26,27,48,49. Furthermore, if synthetic organic chemicals were to be considered for riceland malaria vector control, their management in the current landscape of insecticide resistance across Africa must be considered.Biological control using fish was found to be, in general, slightly more effective than (chemical, bacterial and MSF) larviciding. The degree of effectiveness was dependent on the fish species and their feeding preferences: surface-feeding, larvivorous species provided better anopheline control than bottom-feeding selective feeders4,43. Selecting the most suitable fish for local rice fields is not straightforward; many criteria need to be considered4,52,53. Generally, fish were well-received by rice farmers, perceived to contribute to increased yield by reducing weeds and pests and providing fertiliser through excrement43,44. This was reportedly also observed in Guangxi, China, where a certain proportion of the field had to be deepened into a side-trench where the fish could take shelter when the fields were drained. Even with this reduction in rice production area, carp rearing still increased yields by 10% and farmer’s income per hectare by 70%53. Unfortunately, none of the eligible studies in this review had included yield or water use as an outcome. Future entomological studies need to measure these critical agronomic variables so that studies of vector control in rice can be understood by, and transferred to, agronomists. In SSA, irrigated rice-fish farming can be scaled up provided that an inventory of fish species suitable for specific locations is available and that water is consistently available in fields (an important limiting factor in African irrigation schemes)54. Lessons can be learnt from successful large-scale rice-fish systems in Asia, where they have served as win–win solutions for sustainable food production and malaria control16,55.Overall, there was only limited evidence that intermittent irrigation is effective at reducing late-instar anopheline larvae in rice fields. This finding contrasts with prior reviews, which found mixed results (regardless of larval stage) but emphasised that success was site-specific4,17,56. This contrast is presumably due to the inclusion criteria of our systematic review. These reviews excluded studies in various geographical settings and some older studies that reported successful anopheline control with intermittent irrigation but lacked either a contemporaneous control arm, adequate replication or adequate differentiation between culicines and anophelines16,57,50,51,52,61. It seems, from our review, that intermittent irrigation does not prevent the recruitment of early instars (and in one case, may have encouraged oviposition31) but tends to prevent their development into late-stage immatures. This important conclusion is, however, based only on four studies; more evidence is urgently needed where future trials should consider the basic principles of modern trials with adequate replication, controls and differentiation between larval instars and species.Generally, it is observed that drainage, passive or active, did not reliably reduce overall numbers of mosquito immatures. In India and Kenya, closer inspection revealed that soils were not drying sufficiently, so any stranded larvae were not killed31,46. Highlighted by van der Hoek et al.29 and Keiser et al.17, water management in rice fields is very dependent on the physical characteristics of the soil and the climate and is most suited to places that not only favour rapid drying, but also have a good control of water supply17,56. Moreover, repeated drainage, although directed against mosquitoes, can also kill their aquatic predators62. Since mosquitoes can re-establish themselves in a newly flooded rice field more quickly than their predators, intermittent irrigation with more than a week between successive drying periods can permit repeated cycles of mosquito breeding without any predation pressure. Its efficacy against malaria vectors is therefore highly reliant on the timing of the wetting and drying periods. Further site-specific research on timing, especially with regards to predator–prey interactions within the rice agroecosystem, is required to find the perfect balance.Another limitation in intermittent irrigation is that it cannot be applied during the first two to three weeks following transplanting, because rice plants must remain flooded to recover from transplanting shock. Unfortunately, this time coincides with peak vector breeding. Thus, other methods of larval control would be required to fill this gap. To agronomists, intermittent irrigation provides benefits to farmers, as it does not penalise yield but significantly reduces water consumption. Nonetheless, farmer compliance seems to be variable, especially in areas where water availability is inconsistent and intermittent irrigation would potentially require more labour31,32,39. Importantly, rice farmers doubted their ability to coordinate water distribution evenly amongst themselves, suggesting that there may be sharing issues, as in the “tragedy of the commons”63. Instead, they said that they preferred to have an agreed authority to regulate water46.No general conclusions could be made on the effect on malaria vectors of other rice cultivation practices (apart from water management) because only one study was eligible for each practice. Nevertheless, these experiments on pesticide application, tillage and weed control, as well as another study on plant spacing (not eligible since glass rods were used to simulate rice plants), do illustrate that small changes in agronomic inputs and conditions can have considerable effects on mosquito densities, not just rice yield36,38,64. Moreover, in partially- or shallowly-flooded plots, the larvae are often concentrated in depressions (usually footprints), suggesting that rice operations which leave or remove footprints (e.g. hand-weeding, drum seeders, levelling) will influence vector breeding4.Our study has some important limitations. First, in most trials, the units of intervention were replicate plots of rice, and success was measured as a reduction in larval densities within treated plots. This design focuses on the identification of effective and easy-to-implement ways of growing rice without growing mosquitoes, on the assumption that higher vector densities are harmful. However, from a public health perspective, the need for epidemiological outcomes is often, and reasonably, stressed22,65. Nonetheless, from a farmers’ perspective, it is also important to consider whether the vectors emerging from their rice fields significantly contribute to the local burden of malaria and to determine how this contribution can be minimised. There is evidence that riceland vectors do increase malaria transmission, since human biting rates are much higher in communities living next to rice schemes than their non-rice counterparts66 and that additional riceland vectors may intensify transmission and malaria prevalence in rice communities15. Hence, when investigating how rice-attributed malaria risk can be minimised, mosquito abundance as measured in the experimental rice trials is a useful indicator of potential impact on epidemiological outcomes.Second, larval density was not always separated into larval developmental stages. This can be misleading because some interventions work by reducing larval survival (but not by preventing oviposition) and development to late instars and pupae. Therefore, an intervention could completely eliminate late-stage larvae and pupae but have little effect on the total number of immatures. This was illustrated in our meta-analyses of intermittent irrigation in Table 3 and Supplementary Table 5, and could have been the case for some studies that failed to demonstrate consistent reductions in overall anopheline numbers but did not differentiate between larval instars34,45,67,60,69. We infer that when monitoring mosquito immatures in rice trials, it is important to distinguish between larval instars and pupae. Pupae should always be counted separately since its abundance is the most direct indicator of adult productivity70.Third, experimental trials rarely reported the timing of intervention application or accounted for different rice-growing phases, or “days after transplantation”, in the outcome. Both aspects are important to consider since an intervention may be suited to control larvae during certain growth phases but not others. This is illustrated by Djegbe et al.38, where, compared to deep tillage, minimal tillage could significantly reduce larvae during the early stages of rice cultivation but not during tillering and maturation38. In contrast, other interventions, such as Azolla and predatory copepods, took time to grow and accumulate, and were more effective during the later stages of a rice season45,67,71. This differentiation is important because it can identify components that could potentially form a complementary set of interventions against riceland malaria vectors, each component being effective at different parts of the season. Since rice fields, and hence the dynamics of riceland mosquito populations, vary from place to place, this set of interventions must also be robust. Special attention must be paid to the early stages of rice cultivation, particularly the first few weeks after transplanting (or sowing), since, with many vector species, a large proportion of adult mosquitoes are produced during this time.Fourth, the analysis of entomological counts is often inadequate. Many studies failed to provide the standard deviation (or any other measure of error) for larval counts and could not be included in the quantitative analysis. Often, due to the extreme (and not unexpected) variability of larval numbers, sample sizes were insufficient to calculate statistically significant differences between treatments. Fifth, a high risk of bias was found across both CTS and CITS studies, including high heterogeneity and some publication bias. Study quality was, in general, a shortcoming and limited the number of eligible studies for certain interventions, including intermittent irrigation. Moreover, there are conspicuous a priori reasons for bias in such experimental trials: trial locations are frequently chosen to maximise the probability of success.Finally, few studies were conducted in African countries, where the relationship between rice and malaria is most important because of the efficiency, and the “rice-philic” nature, of the vector An. gambiae s.l.15. In particular, there was a lack of studies on the effectiveness and scalability of biological control and rice cultivation practices. There is also very little information (particularly social science studies) on the views and perspectives of African rice farmers on mosquitoes in rice and interventions to control them72,73.In the future, as malaria declines (particularly across SSA), the contribution of rice production to increased malaria transmission is likely to become more conspicuous15. Unless this problem is addressed, rice growing will probably become an obstacle to malaria elimination. Current default methods of rice production provide near-perfect conditions for the larvae of African malaria vectors. Therefore, we need to develop modified rice-growing methods that are unfavourable to mosquitoes but still favourable for the rice. Although larviciding and biological control may be appropriate, their unsustainable costs remain the biggest barrier to uptake amongst smallholder farmers. Future investigations into riceland vector control should pay more attention to interventions that may be useful to farmers.Supported by medical entomologists, agronomists should lead the research task of identifying cultivation methods that achieve high rice productivity whilst suppressing vector productivity. Rice fields are a major global source of greenhouse gases, and agronomists have responded by successfully developing novel cultivation methods that minimise these emissions while maintaining yield. We need the same kind of response from agronomists, to achieve malaria control co-benefits within rice cultivation. At present, only a few aspects of rice cultivation have been investigated for their effects on mosquitoes, and the potential of many other practices for reducing anopheline numbers are awaiting study. Due to the spatial and temporal heterogeneity of rice agroecosystems, it is likely that no single control method can reduce mosquito numbers throughout an entire cropping season and in all soil types and irrigation methods. Thus, effective overall control is likely to come from a combination of local, site-specific set of complementary methods, each of which is active and effective during a different phase of the rice-growing season. More

The effects of extreme drought on soil biochemical propertiesAs shown in Fig. 1A, the range of SOC during the early, midterm and late extreme drought experiments, were 73.53–251.44 g kg−1, 54.75–256.16 g kg−1, and 66.37–282.16 g kg−1, respectively. Concomitantly, DOC was 171.85–323.74 mg kg−1, 158.15 – 504.62 mg kg−1, and 166.63–418.43 mg kg−1, MBC was 247.80 – 461.69 mg kg−1, 257.90–450.98 mg kg−1, and 264.10–458.15 mg kg−1, respectively (Fig. 1B, C). The variation ranges of soil TN were 3.50–16.60 g kg−1, 4.70–34.5 g kg−1, and 6.70–32.50 g kg−1, respectively (Fig. 1D). Similarly, the variation ranges of NH4+ were 5.96–12.03 g kg−1, 5.39–12.59 g kg−1, and 5.74–13.03 g kg−1, NO3− were 2.27–8.79 mg kg−1, 5.07–9.62 mg kg−1, and 5.09–9.52 mg kg−1, respectively (Fig. 1E, F). The changes of SOC and NH4+ with soil depth were significantly different in different extreme drought periods and decreased significantly with the increase of soil depth (Table 1, P  More

During recent decades, pathogens that originated in bats have become an increasing public health concern. A major challenge is to identify how those pathogens spill over into human populations to generate a pandemic threat1. Many correlational studies associate spillover with changes in land use or other anthropogenic stressors2,3, although the mechanisms underlying the observed correlations have not been identified4. One limitation is the lack of spatially and temporally explicit data on multiple spillovers, and on the connections among spillovers, reservoir host ecology and behavior, and viral dynamics. We present 25 years of data on land-use change, bat behavior, and spillover of Hendra virus from Pteropodid bats to horses in subtropical Australia. These data show that bats are responding to environmental change by persistently adopting behaviors that were previously transient responses to nutritional stress. Interactions between land-use change and climate now lead to persistent bat residency in agricultural areas, where periodic food shortages drive clusters of spillovers. Pulses of winter flowering of trees in remnant forests appeared to prevent spillover. We developed integrative Bayesian network models based on these phenomena that accurately predicted the presence or absence of clusters of spillovers in each of 25 years. Our long-term study identifies the mechanistic connections among habitat loss, climate, and increased spillover risk. It provides a framework for examining causes of bat virus spillover and for developing ecological countermeasures to prevent pandemics. More

These authors contributed equally: Yannick Garcin, Enno Schefuß, Greta C. Dargie, Simon L. LewisAix Marseille University, CNRS, IRD, INRAE, CEREGE, Aix-en-Provence, FranceYannick Garcin & Ghislain GassierInstitute of Geosciences, University of Potsdam, Potsdam, GermanyYannick GarcinMARUM—Center for Marine Environmental Sciences, University of Bremen, Bremen, GermanyEnno SchefußSchool of Geography, University of Leeds, Leeds, UKGreta C. Dargie, Bart Crezee, Dylan M. Young, Andy J. Baird, Paul J. Morris & Simon L. LewisSchool of Geography and Sustainable Development, University of St Andrews, St Andrews, UKDonna Hawthorne, Ian T. Lawson & George E. BiddulphIFP Energies Nouvelles, Earth Sciences and Environmental Technologies Division, Rueil-Malmaison, FranceDavid SebagInstitute of Earth Surface Dynamics, Geopolis, University of Lausanne, Lausanne, SwitzerlandDavid SebagFaculté des Sciences et Techniques, Université Marien Ngouabi, Brazzaville, Republic of the CongoYannick E. Bocko & Y. Emmanuel Mampouya WeninaÉcole Normale Supérieure, Université Marien Ngouabi, Brazzaville, Republic of the CongoSuspense A. IfoÉcole Normale Supérieure d’Agronomie et de Foresterie, Université Marien Ngouabi, Brazzaville, Republic of the CongoMackline MbembaFaculté de Gestion des Ressources Naturelles Renouvelables, Université de Kisangani, Kisangani, Democratic Republic of the CongoCorneille E. N. Ewango & Joseph Kanyama TabuFaculté des Sciences, Université de Kisangani, Kisangani, Democratic Republic of the CongoCorneille E. N. EwangoInstitut Supérieur Pédagogique de Mbandaka, Mbandaka, Democratic Republic of the CongoOvide Emba & Pierre BolaSchool of Geography, Geology and the Environment, University of Leicester, Leicester, UKGenevieve Tyrrell, Arnoud Boom & Susan E. PageSchool of Water, Energy and Environment, Cranfield University, Bedford, UKNicholas T. GirkinBritish Geological Survey, Centre for Environmental Geochemistry, Keyworth, UKChristopher H. VaneInstitute of Earth Sciences, University of Lausanne, Lausanne, SwitzerlandThierry AdatteNEIF Radiocarbon Laboratory, Scottish Universities Environmental Research Centre (SUERC), Glasgow, UKPauline GulliverSchool of Biosciences, University of Nottingham, Nottingham, UKSofie SjögerstenDepartment of Geography, University College London, London, UKSimon L. Lewis More

Automated tracking of four social insect speciesFifty queenright colonies were used in the tracking experiments (Table 1). Honeybee colonies (subspecies A. mellifera carnica) were housed in the campus apiary of the University of Lausanne. Colonies of L. niger were raised from single mated queens collected on campus. T. nylanderi colonies were collected from the University of Lausanne campus, and L. acervorum colonies collected from Anzeindaz, Switzerland. These four species were chosen because of their abundance and easy availability in Switzerland, and because they – or closely-related species – have previously been used as model systems for the study of spatial organisation in social insects20,21,23,27,44. The colony sizes used in our experiments (Table 1) fell within the natural range of sizes experienced by these species in nature, either as recently founded colonies (L. niger colonies are founded by a single queen and progressively grow from a few workers to mature sizes of up to 40,000 workers over the course of several years; new honeybee colonies are founded by swarms counting 2400–41,000 bees61) or as mature colonies (all colonies of T. nylanderi and L. acervorum used in our experiments were mature colonies collected whole from the field).In all species, a paper tag bearing a unique two-dimensional barcode was glued to the thorax of individuals to allow automated tracking of their movements (Fig. S1). In the ants, tagging of all individuals was performed in a single session two days before the beginning of the experiment, whilst in the bees, newly-emerged workers (one-day-old or less) were tagged every 3 days over the 21 days prior to the beginning of the experiment (Supplementary Note 1).Tagged colonies were kept in glass observation nests with a single entrance (internal nest dimensions, A. mellifera: 69 × 45 × 4 cm, L. niger: 70 × 40 × 8 mm; L. acervorum: 63 × 42 × 2 mm, T. nylanderi: 63 × 42 × 1.5 mm). The honeybee observation nests also included a 64 × 44 cm wooden frame enclosing a double-sided wax comb containing honey, pollen, and developing brood20. Bees were free to move between both sides of the comb. In all species, individuals were allowed to freely exit and enter the nest. Ants were provided with ad libitum food (Drosophila, sugar solution) and water in a foraging arena, while bees foraged on natural resources outside. Both the ant and honeybee observation nests were exposed to diurnal cycles of temperature and light (Supplementary Note 1).High resolution digital video cameras operating at two frames per second were used to identify the location and orientation of each tag across successive images22. All colonies were continuously tracked for three days, which corresponded to the inter-cohort time in the honeybee colonies. The trajectories of each worker, and the physical contacts between workers (Fig. S18 and Supplementary Note 14) were extracted using an existing software pipeline62.Building bipartite site-visit networksTo quantify the spatial preferences of individual ants and bees, the interior of the nest was discretised into a regular hexagonal lattice (Fig. 1a, b). Because the worker body lengths of our four study species span an order of magnitude (from ~ 1.5 mm for T. nylanderi to ~ 15 mm for A. mellifera), the width of the hexagonal bins were defined as 1/4 of the mean worker body-length.To characterise the spatial preferences of different individuals to different parts of the nest, we counted the number of times ({n}_{i}^{s}) that each individual i visited each hexagonal site s. A visit by individual i to site s began when i crossed the border into s, and was terminated when i crossed the border out of s, regardless of the amount of time spent inside. To prevent stationary individuals located on the border between two adjacent sites from rapidly accumulating many single-frame visits to the two sites, successive visits to a same site were only counted when at least 20s elapsed between the end of the previous visit and the start of the next.The site-visit data were used to construct a bipartite network, in which individuals (layer 1) were connected by undirected edges to the sites (layer 2) they visited (Fig. 1c, d). Because individuals typically made multiple visits to the same sites, each edge i–s was weighted according to the total number of times individual i visited site s, that is, ({n}_{i}^{s}).Partitioning site-visit networks into modulesThe extent to which the site-visit networks were partitioned into discrete ‘modules’ (i.e., set of workers with similar space-use patterns and the set of sites that they exhibit strong ties to) was assessed using the DIRTLPAwb+ algorithm for partitioning weighted bipartite networks39. This algorithm searches for the partition that maximises the number and strength of the links within modules, whilst minimizing connections between modules. The number of modules was not specified a priori by the user, but was identified by the algorithm. All site-visit networks had positive modularity (Fig. S3), indicating that they could be partitioned into a set of well-separated modules (Figs. 1e–h, S2, and S4–S5). The modules in each partition were then assigned functional labels according to the following rules. First, the module whose sites were on average closest to the nest entrance was labelled ‘forager’ module. Second, the module or modules with the greatest spatial overlap with the brood pile in the ant colonies or the broodnest(s) in the honeybee colonies were labelled ‘nurse’ module(s). After defining the forager and nurse modules, the remaining modules (if any) were labelled as follows. If there was only one module remaining after identifying the nurse and forager modules, as was typically the case in honeybee colonies, it was labelled ‘peripheral’. If there were two modules remaining, as was typically the case in ant colonies, then the module whose sites were on average closer to the nest borders (i.e., to the periphery of the nest) was labelled ‘peripheral’, and the remaining module labelled ‘intermediate’. In some cases, the DIRTLPAwb+ algorithm identified five or more modules (9.0% of all iterations across all species and colonies). In these cases, the supernumerary modules never contained more than 1 or 2 individuals, and as they could not be unequivocally assigned using our labelling scheme, they were left unclassified for these iterations.Validating network modulesAs a network constructed by a purely random process could exhibit apparent modular structure by chance, we tested whether the discovered modules represent statistically significant entities. To do so, we produced 1000 null model random networks for each observed network using an established permutation method for bipartite networks63 (Supplementary Note 2). Comparisons between the maximum modularity of the observed networks with that of the corresponding random networks showed that, in all four species, the observed modularity was significantly greater than expected by chance (Fig. S3).Constructing worker task profilesA unique labour profile for each ant and each honeybee was constructed by estimating the activity of each worker in the following four tasks:1. Entrance guarding: workers were classed as guarding when they were (i) within two body lengths of the entrance, (ii) roughly facing the entrance, i.e., with a body alignment diverging from the direct heading to the entrance by no more than π/2 radians, and (iii) ‘on station’ at the entrance, as defined by trajectory coordinates with an associated first passage time (ref. 64; time taken for the individual to pass beyond a circle centred on its current location with a radius of two body-lengths) in excess of 500s.2. Patrolling: workers were classed as patrolling65 when they were (i) active, and (ii) ‘roaming’, as defined by first passage times of More

Podos, J. Correlated evolution of morphology and vocal signal structure in Darwin’s finches. Nature 409, 185–188 (2001).ADS 
CAS 
PubMed 

Google Scholar 
Bradbury, J. W. & Vehrencamp, S. Principles of Animal Communication (Sinauer Associated, 1998).
Google Scholar 
Podos, J. & Warren, P. S. The evolution of geographic variation in birdsong. Adv. Study Behav. 37, 403–458 (2007).
Google Scholar 
Charlton, B. D., Owen, M. A. & Swaisgood, R. R. Coevolution of vocal signal characteristics and hearing sensitivity in forest mammals. Nat. Commun. 10, 1–7 (2019).CAS 

Google Scholar 
Soltis, J. Vocal communication in African elephants (Loxodonta africana). Zoo Biol. 29, 192–209 (2010).PubMed 

Google Scholar 
King, S. L. & Janik, V. M. Bottlenose dolphins can use learned vocal labels to address each other. Proc. Natl. Acad. Sci. 110, 13216–13221 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ravignani, A. et al. Ontogeny of vocal rhythms in harbor seal pups: An exploratory study. Curr. Zool. 65, 107–120 (2019).PubMed 

Google Scholar 
Rauber, R. & Manser, M. B. Effect of group size and experience on the ontogeny of sentinel calling behaviour in meerkats. Anim. Behav. 171, 129–138 (2021).
Google Scholar 
Janik, V. M. & Slater, P. J. Vocal learning in mammals. Adv. Study Behav. 26, 59–100 (1997).
Google Scholar 
Fitch, W. T. Production of vocalizations in mammals. In Encyclopedia of Language and Linguistics 2nd edn, Vol. 1 115–121 (Elsevier, 2006).
Google Scholar 
Fletcher, N. H. A frequency scaling rule in mammalian vocalization. In Handbook of Behavioral Neuroscience Vol. 19 51–56 (Elsevier, 2010).
Google Scholar 
Fant, G. Acoustic Theory of Speech Production (Mouton, 1960).
Google Scholar 
Taylor, A. & Reby, D. The contribution of source–filter theory to mammal vocal communication research. J. Zool. 280, 221–236 (2010).
Google Scholar 
Fitch, W. T., Neubauer, J. & Herzel, H. Calls out of chaos: The adaptive significance of nonlinear phenomena in mammalian vocal production. Anim. Behav. 63, 407–418 (2002).
Google Scholar 
Domning, D. P. & Hayek, L. A. C. Interspecific and intraspecific morphological variation in manatees (Sirenia: Trichechus). Mar. Mamm. Sci. 2, 87–144 (1986).
Google Scholar 
Anderson, P. K. & Barclay, R. M. R. Acoustic signals of solitary Dugongs: Physical characteristics and behavioral correlates. J. Mamm. 76, 1226–1237 (1995).
Google Scholar 
Sousa-Lima, R., Paglia, A. P. & Da Fonseca, G. Signature information and individual recognition in the isolation calls of Amazonian manatees, Trichechus inunguis (Mammalia: Sirenia). Anim. Behav. 63, 301–310 (2002).
Google Scholar 
Sousa-Lima, R. S., Paglia, A. P. & Fonseca, G. A. B. Gender, age, and identity in the isolation calls of Antillean manatees (Trichechus manatus manatus). Aquat. Mamm. 34, 109–122 (2008).
Google Scholar 
O’Shea, T. J. & Poché, L. B. Aspects of underwater sound communication in Florida manatees (Trichechus manatus latirostris). J. Mamm. 87, 1061–1071 (2006).
Google Scholar 
Rosas, F. C. W. Biology, conservation and status of the Amazonian manatee Trichechus inunguis. Mamm. Rev. 24, 49–59 (1994).
Google Scholar 
Meirelles, A. C. O. & Carvalho, V. L. Peixe-boi marinho: biologia e conservação no Brasil. Aquasis, Bambu Editora e Artes Gráficas, São Paulo (2016).Alvarez-Alemán, A., Beck, C. A. & Powell, J. A. First report of a Florida manatee (Trichechus manatus latirostris) in Cuba. Aquat. Mamm. 36, 148 (2010).
Google Scholar 
Castelblanco-Martínez, D. N. et al. First documentation of long-distance travel by a Florida manatee to the Mexican Caribbean. Ethol. Ecol. Evol. 1–12 (2021).Packard, J. M. & Wetterqvist, O. F. Evaluation of manatee habitat systems on the northwestern Florida coast. Coast. Manag. 14, 279–310 (1986).
Google Scholar 
Luna, F. D. O. et al. Genetic connectivity of the West Indian manatee in the southern range and limited evidence of hybridization with Amazonian manatees. Front. Mar. Sci. 7, 1089 (2021).
Google Scholar 
Hartman, D. Ecology and behavior of the manatee (Trichechus manatus) in Florida. Spec. Publ. Am. Soc. Mammal. 5, 153 (1979).
Google Scholar 
D’AffonsecaNeto, J. A. & Vergara-Parente, J. E. Sirenia (peixe-boi-da-Amazônia, Peixe-boi-marinho). In Tratado de Animais Selvagens: medicina veterinária (eds Cubas, Z. S. et al.) 701–714 (Roca, 2006).
Google Scholar 
Deutsch, C. J., Reid, J. P., Bonde, R. K., Easton, D. E., Kochman, H. I. & O’Shea, T. J. Seasonal movements, migratory behavior, and site fidelity of West Indian manatees along the Atlantic coast of the United States. Wildl. Monogr. 1–77 (2003).Laist, D. W. & Reynolds, J. E. III. Influence of power plants and other warm-water refuges on Florida manatees. Mar. Mamm. Sci. 21, 739–764 (2005).
Google Scholar 
Reynolds, J. E. Aspects of the social behaviour and herd structure of a semi-isolated colony of West Indian manatees, Trichechus manatus. Mammalia 45, 431–452 (1981).
Google Scholar 
Dantas, G. A. Ontogenia do padrão vocal individual do peixe-boi da Amazônia Trichechus inunguis (Sirenia, trichechidae). Dissertação (Instituto Nacional de Pesquisas da Amazônia, 2009).
Google Scholar 
Brady, B., Moore, J. & Love, K. Behavior related vocalizations of the Florida manatee (Trichechus manatus latirostris). Mar. Mamm. Sci. 1–15 (2021).Nowacek, D. P., Casper, B. M., Wells, R. S., Nowacek, S. M. & Mann, D. A. Intraspecific and geographic variation of West Indian manatee (Trichechus manatus spp.) vocalizations. J. Acoust. Soc. Am. 114, 66–69 (2003).ADS 
PubMed 

Google Scholar 
Rycyk, A. M. et al. First characterization of vocalizations and passive acoustic monitoring of the vulnerable African manatee (Trichechus senegalensis). J. Acoust. Soc. Am. 150, 3028–3037 (2021).ADS 
PubMed 

Google Scholar 
Landrau-Giovannetti, N., Mignucci-Giannoni, A. A. & Reidenberg, J. S. Acoustical and anatomical determination of sound production and transmission in West Indian (Trichechus manatus) and Amazonian (T. inunguis) manatees. Anat. Rec. 297, 1896–1907 (2014).
Google Scholar 
Morton, E. On the occurrence and significance of motivation-structural rules in some bird and mammal sounds. Am. Nat. 111, 855–869 (1977).
Google Scholar 
Borges, J. C. et al. Growth pattern differences of captive born Antillean manatee (Trichechus manatus) calves and those rescued in the Brazilian northeastern coast. J. Zoo Wildl. Med. 43, 494–500 (2012).PubMed 

Google Scholar 
Lima, D. S., Vergara-Parente, J. E., Young, R. J. & Paszkiewicz, E. Training of Antillean manatee Trichechus manatus manatus Linnaeus, 1758 as a management technique for individual welfare. Lat. Am. J. Mar. Mamm. 4, 61–68 (2005).
Google Scholar 
K. Lisa Yang Center for Conservation Bioacoustics at the Cornell Lab of Ornithology. Raven Pro: Interactive Sound Analysis Software (Version 1.5) [Computer software]. https://ravensoundsoftware.com/ (The Cornell Lab of Ornithology, 2022).Zollinger, S. A., Podos, J., Nemeth, E., Goller, F. & Brumm, H. On the relationship between, and measurement of, amplitude and frequency in birdsong. Anim. Behav. 84, e1–e9 (2012).
Google Scholar 
Sokal, R. R. & Rohlf, F. J. Biometry (W. H. Freeman and Co., 1995).MATH 

Google Scholar 
Charrier, I. & Harcourt, R. G. Individual vocal identity in mother and pup Australian sea lions (Neophoca cinerea). J. Mamm. 87, 929–938 (2006).
Google Scholar 
Green, S. & Salkind, N. J. Using SPSS for Windows and Macintosh: Analyzing and Understanding Data 4th edn. (Prentice Hall, 2003).
Google Scholar 
Charlton, B. D. et al. Cues to body size in the formant spacing of male koala (Phascolarctos cinereus) bellows: Honesty in an exaggerated trait. J. Exp. Biol. 214(20), 3414–3422 (2011).PubMed 

Google Scholar 
IBM Corp. Released 2020. IBM SPSS Statistics for Windows, Version 27.0.Best, R. C. The aquatic mammals and reptiles of the Amazon. In The Amazon 371–412 (Springer, 1984).
Google Scholar 
Gerhardt, H. C. The evolution of vocalization in frogs and toads. Annu. Rev. Ecol. Syst. 25, 293–324 (1994).
Google Scholar 
Mendoza, P. et al. Growth curve of Amazonian manatee (Trichechus inunguis) in captivity. Aquat. Mamm. 45 (2019).Schwarz, L. K. Methods and models to determine perinatal status of Florida manatee carcasses. Mar. Mamm. Sci. 24, 881–898 (2008).
Google Scholar 
Hauser, M. D., Chomsky, N. & Fitch, W. T. The faculty of language: What is it, who has it, and how did it evolve?. Science 298, 1569–1579 (2002).ADS 
CAS 
PubMed 

Google Scholar 
Marler, P. & Peters, S. Developmental overproduction and selective attrition: New processes in the epigenesis of birdsong. Dev. Psychol. J. Int. Soc. Dev. Psychol. 15.4, 369–378 (1982).
Google Scholar 
Casey, C., Reichmuth, C., Costa, D. P. & Le Boeuf, B. The rise and fall of dialects in northern elephant seals. Proc. R. Soc. B 285, 2018–2176 (2018).
Google Scholar 
Hunter, M. E. et al. Puerto Rico and Florida manatees represent genetically distinct groups. Conserv. Genet. 13, 1623–1635 (2012).
Google Scholar 
Castelblanco-Martínez, D. N. et al. Analysis of body condition indices reveals different ecotypes of the Antillean manatee. Sci. Rep. 11, 1–14 (2021).
Google Scholar 
McCracken, K. G. & Sheldon, F. H. Avian vocalizations and phylogenetic signal. Proc. Natl. Acad. Sci. 94, 3833–3836 (1997).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Moron, J. R. et al. Whistle variability of Guiana dolphins in South America: Latitudinal variation or acoustic adaptation? Mar. Mamm. Sci. 1–32 (2018)Luís, A. R. et al. Vocal universals and geographic variations in the acoustic repertoire of the common bottlenose dolphin. Sci. Rep. 11, 1–9 (2021).
Google Scholar 
Ey, E. & Fisher, J. The ‘Acoustic adaptations hypothesis’ a review of the evidence from birds, anurans and mammals. Bioacoustics 19, 21–48 (2009).
Google Scholar 
Miksis-Olds, J. L. & Tyack, P. L. Manatee (Trichechus manatus) vocalization usage in relation to environmental noise levels. J. Acoust. Soc. Am. 125, 1806–1815 (2009).ADS 
PubMed 

Google Scholar 
Sun, W., Wang, Z., Jamalabdollahi, M. & Reza Zekavat, S. A. Experimental study on the difference between acoustic communication channels in freshwater rivers/lakes and in oceans. In 2014 48th Asilomar Conference on Signals, Systems and Computers, 333–337 (2004)Rivera Chavarría, M., Castro, J. & Camacho, A. The relationship between acoustic habitat, hearing and tonal vocalizations in the Antillean manatee (Trichechus manatus manatus, Linnaeus, 1758). Biol. Open 4, 1237–1242 (2015).PubMed 
PubMed Central 

Google Scholar 
Gaspard, J. C. III. et al. Audiogram and auditory critical ratios of two Florida manatees (Trichechus manatus latirostris). J. Exp. Biol. 215, 1442–1447 (2012).PubMed 

Google Scholar 
Gerstein, E. R., Gerstein, L., Forsythe, S. E. & Blue, J. E. The underwater audiogram of the West Indian manatee (Trichechus manatus). J. Acoust. Soc. Am. 105, 3575–3583 (1999).ADS 
CAS 
PubMed 

Google Scholar 
Klishin, V., Pezo, R., Popov, V., Ya, A. & Supin. Some characteristics of hearing of the Brazilian manatee, Trichechus inunguis. Aquat. Mamm. 16 (1990).Johnson, M., de Soto, N. A. & Madsen, P. T. Studying the behaviour and sensory ecology of marine mammals using acoustic recording tags: A review. Mar. Ecol. Prog. Ser. 395, 55–73 (2009).ADS 

Google Scholar  More