fish school experiment

Schooling Behavior of Fish: General Ideas, Terms and Concepts, Prevalence, Applied Aspects

• Published: 14 December 2023
• Volume 63 , pages 1219–1233, ( 2023 )

Cite this article

• A. O. Kasumyan 1 &
• D. S. Pavlov 2  

232 Accesses

2 Citations

Explore all metrics

The article opens a thematic issue of the Journal of Ichthyology , dedicated to Dmitry Viktorovich Radakov, a researcher who made an outstanding contribution to the knowledge of schooling behavior of fish. The article deals with many terms and concepts that are widely used in scientific publications on the schooling behavior of fish, but still do not have a generally accepted definition—crowd, aggregation, shoal, school, flock. Attention is drawn to the difficulties caused by the fact that the use of these terms is most often based not on objective criteria, but on the preferences, views or beliefs of individual researchers. General ideas about the prevalence of schooling behavior in fish of different taxonomy, lifestyle, condition and age; about the importance of the visual structuring of the environment for the manifestation of schooling, about the difficulty of dividing fish into facultative and obligate schooling, about the importance of knowledge about schooling behavior for solving applied problems have been considered.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Subscribe and save.

• Get 10 units per month
• Download Article/Chapter or eBook
• 1 Unit = 1 Article or 1 Chapter
• Cancel anytime

Price includes VAT (Russian Federation)

Instant access to the full article PDF.

Rent this article via DeepDyve

Institutional subscriptions

Similar content being viewed by others

Problems and Prospects of Studying Schooling Behavior of Fish

Continuity of schooling behavior of fish.

On the Problem of the Evolutionary Origin of Schooling Behavior of Fish

Kaganovskii, A.G., Composition of stocks and behavior of the Far Eastern sardine Sardinops sagax melanosticta (Temm and Schleg) in connection with oceanographic conditions, Cand. Sci. (Biol.) Dissertation , Sverdlovsk: Moscow State Univ., 1943.

Allee, W.C., Animal Aggregations, a Study in General Sociology , Chicago: Univ. Chicago Press, 1931. https://doi.org/10.5962/bhl.title.7313

Baskin, L.M., Tolpa i stado (Crowd and Herd), St. Petersburg: Nestor-Istoriya, 2017.

Bernhardt, B., Lampert, K.P., Leese, F., et al., Are shoals of minnow Phoxinus phoxinus formed by close kin?, J. Fish. Biol. , 2012, vol. 80, no. 3, pp. 713–721. https://doi.org/10.1111/j.1095-8649.2011.03198.x

Article   CAS   PubMed   Google Scholar  

Bourke, A.F.G., Principles of Social Evolution, Oxford: Oxford Univ. Press, 2011. https://doi.org/10.1093/acprof:oso/9780199231157.001.0001

Breder, C.M., Studies on social groupings in fishes, Bull. AMNH , 1959, vol. 117, no. 1, pp. 393–482.

Google Scholar  

Breder, C.M., Fish schools as operational structures, Fish. Bull. , 1976, no. 74, pp. 471–502.

Brehmer, P., Gerlotto, F., Laurent, C., et al., Schooling behaviour of small pelagic fish: Phenotypic expression of independent stimuli, Mar. Ecol. Proc. Ser. V , 2007, vol. 334, pp. 263–272. https://doi.org/10.3354/meps334263

Article   Google Scholar  

Brown, C. and Warburton, K., Social mechanisms enhance escape responses in shoals of rainbowfish, Melanotaenia duboulayi, Environ. Biol. Fish., 1999, vol. 56, no. 4, pp. 455–459.

Budaev, S.V., “Personality” in the guppy ( Poecilia reticulata ): A correlation study of exploratory behavior and social tendency, J. Comp. Psychol. , 1997, vol. 111, no. 4, pp. 399–411. https://doi.org/10.1037/0735-7036.111.4.399

Charef, A., Ohshimo, S., Aoki, I., and Absi, N.A., Classification of fish schools based on evaluation of acoustic descriptor characteristics, Fish. Sci. , 2010, vol. 76, no. 1, pp. 1–11. https://doi.org/10.1007/s12562-009-0186-x

Article   CAS   Google Scholar  

Chen, M.J., Coss, R.G., and Goldthwaite, R.O., Timing of dispersal in juvenile jewel fish during development is unaffected by available space, Develop. Psychobiol. , 1983, vol. 16, no. 4, pp. 303–310. https://doi.org/10.1002/dev.420160405

Conçalves, A.R., Silva, M.F., Vendrell, C.L., and Almada, V.C., Agonistic behaviour and shoal composition of juvenile Diplodus sargus : First field observations, Environ. Biol. Fish., 2015, vol. 98, no. 4, pp. 1015–1021. https://doi.org/10.1007/s10641-014-0334-8

Connaughton, M.A. and Taylor, M.H., Seasonal and daily cycles in sound production associated with spawning in the weakfish, Cynoscion regalis, Environ. Biol. Fish., 1995, vol. 42, no. 3, pp. 233–240. https://doi.org/10.1007/BF00004916

Dagorn, L. and Fréon, P., Tropical tuna associated with floating objects: A simulation study of the meeting point hypothesis, Can. J. Fish. Aquat. Sci. , 1999, vol. 56, no. 6, pp. 984–993. https://doi.org/10.1139/f98-209

Dagorn, L. and Holland, K., Report of the International workshop on current status and new directions for studying schooling and aggregation behavior of pelagic fish, PFRP (Pelagic Fish. Res. Program) , 2003, vol. 8, no. 3, pp. 7–8.

Darkov, A.A. and Mochek, A.D., Behavior of the Black Sea pickarel ( Spicara smarus (L.)) during reproduction, Zool. Zh. , 1979, vol. 58, no. 7, pp. 1011–1015.

de Bie, J., Manes, C., and Kemp, P.S., Collective behaviour of fish in the presence and absence of flow, Anim. Behav. , 2020, vol. 167, pp. 151–159. https://doi.org/10.1016/j.anbehav.2020.07.003

Disler, N.N., Organy chuvstv sistemy bokovoi linii i ikh znachenie v povedenii ryb (Sensory Organs of the Lateral Line System and Their Significance in the Behavior of Fish), Moscow: Akad. Nauk SSSR, 1960.

Dmitrieva, E.N., Stages of development of the nonmigratory bream, Tr. Inst. Morfol. Zhivotnykh Akad. Nauk SSSR , 1960, no. 28, pp. 41–78.

Dugatkin, L.A. and Godin, J.-G.J., Predator inspection, shoaling and foraging under predation hazard in the Trinidadian guppy, Poecilia reticulata, Environ. Biol. Fish., 1992, vol. 34, no. 3, pp. 265–276. https://doi.org/10.1007/BF00004773

Economakis, A.E. and Lobel, P.S., Aggregation behavior of the grey reef shark, Carcharhinus amblyrhynchos , at Johnston Atoll, Central Pacific Ocean, Environ. Biol. Fish., 1998, vol. 51, no. 2, pp. 129–139. https://doi.org/10.1023/A:1007416813214

Ehrlich, P.R. and Ehrlich, A.H., Coevolution: Heterospecific schooling in Caribbean reef fishes, Am. Natur. , 1973, vol. 107, no. 953, pp. 157–160. https://doi.org/10.1086/282823

Eibl-Eibesfeldt, I., Ethology: The Biology of Behavior , New York: Holt, Rinehart and Winston, 1970.

Erisman, B.E., Konotchick, T.H., and Blum, S., Observations of spawning in the leather bass, Dermatolepis dermatolepis (Teleostei: Epinephelidae), at Cocos Island, Costa Rica, Environ. Biol. Fish. , 2009, vol. 85, no. 1, pp. 15–20. https://doi.org/10.1007/s10641-009-9463-x

Fernald, R.D. and Hirata, N.R., The ontogeny of social behavior and body coloration in the African cichlid fish Haplochromis burtoni, Z. Tierpsychol. , 1979, vol. 50, no. 2, pp. 180–187. https://doi.org/10.1111/j.1439-0310.1979.tb01025.x

Fernö, A., advances in understanding of basic behavior: consequences for fish capture studies, ICES Mar. Sci. Symp. , 1993, no. 196, pp. 5–11.

Fréon, P. and Dagorn, L., Review of fish associative behaviour: Toward a generalisation of the meeting point hypothesis, Rev. Fish Biol. Fish. , 2000, vol. 10, no. 2, pp. 183–207. https://doi.org/10.1023/A:1016666108540

Gerasimov, V.V., Feeding behavior of the Murmansk herring in and out of the school under aquatic conditions, Tr. Murmansk. Morsk. Biol. Inst. Akad. Nauk SSSR , 1962, no. 4, pp. 254–259.

Gerasimov, V.V., Specificity of imitation in fish, Tr. Murmansk. Morsk. Biol. Inst. Akad. Nauk SSSR , 1964, no. 5, pp. 177–180.

Gerasimov, V.V., Ekologo-fiziologicheskie zakonomernosti stainogo povedeniya ryb (Ecological and Physiological Patterns of Schooling Behavior of Fish) Moscow: Nauka, 1983.

Girsa, I.I., Osveshchennost’ i povedenie ryb (Illumination and Behavior of Fish), Moscow: Nauka, 1981.

Girsa, I.I. and Lapin, Yu.E., Rhythmic behavior of some fish in the Kandalaksha Bay of the White Sea, Biol. Morya , 1985, no. 5, pp. 55–57.

Griffiths, S.W. and Magurran, A.E., Familiarity in schooling fish: How long does it take to acquire?, Anim. Behav. , 1997, vol. 53, no. 5, pp. 945–949. https://doi.org/10.1006/anbe.1996.0315

Herbert-Read, J.E., Understanding how animal groups achieve coordinated movement, J. Exp. Biol. , 2016, vol. 219, no. 19, pp. 2971–2983. https://doi.org/10.1242/jeb.129411

Article   CAS   PubMed   PubMed Central   Google Scholar  

Holubová, M., Blabolil, P., Čech, M., et al., Species-specific schooling behaviour of fish in the freshwater pelagic habitat: An observational study, J. Fish. Biol. , 2020, vol. 97, no. 1, pp. 64–74. https://doi.org/10.1111/jfb.14326

Article   PubMed   Google Scholar  

In memory of D.V. Radakov, Vopr. Ikhtiol. , 1979, vol. 79, no. 3, pp. 546–548.

Kasumyan, A.O. and Ponomarev, V.Y., Study of the behavior of zebrafish, Brachydanio rerio , in response to natural chemical food stimuli, J. Ichthyol. , 1986, vol. 26, no. 5, pp. 96–105.

Kasumyan, A.O. and Pavlov, D.S., Stainoe povedenie ryb (Schooling Behavior of Fish), Moscow: KMK, 2018.

Keenleyside, M.H.A., Some aspects of the schooling behaviour of fish, Behaviour , 1955, vol. 8, no. 1, pp. 183–247. https://doi.org/10.1163/156853955X00229

Khodorevskaya, R.P., Ruban, G.I., and Pavlov, D.S., Povedenie, migratsii, raspredelenie i zapasy osetrovykh ryb Volgo-Kaspiiskogo basseina (Behavior, Migrations, Distribution and Stocks of Sturgeons in the Volga–Caspian Basin), Moscow: KMK, 2007.

Klimley, A.P., Schooling in Sphyrna lewini , a species with low risk of predation: A non-egalitarian state, Z. Tierpsychol. , 1985, vol. 70, no. 4, pp. 297–319. https://doi.org/10.1111/j.1439-0310.1985.tb00520.x

Kotrschal, A., Szorkovszky, A., Romenskyy, M., et al., Brain size does not impact shoaling dynamics in unfamiliar groups of guppies (poecilia reticulata), Behav. Processes , 2018, vol. 147, pp. 13–20. https://doi.org/10.1016/j.beproc.2017.12.006

Kowalko, J.E., Rohner, N., Rompani, S.B., et al., Loss of schooling behavior in cavefish through sight-dependent and sight-independent mechanisms, Curr. Biol. , 2013, vol. 23, no 19, pp. 1874–1883. https://doi.org/10.1016/j.cub.2013.07.056

Kratt, L.F. and Smith, R.J.F., Agonistic behaviour of age O, age I and non-breeding adult Arctic grayling Thymallus arcticus (Pallas), J. Fish. Biol. , 1979, vol. 15, no. 4, pp. 389–404. https://doi.org/10.1111/j.1095-8649.1979.tb03622.x

Krause, J., The influence of hunger on shoal size choice by three-spined sticklebacks, Gasterosteus aculeatus, J. Fish. Biol. , 1993, vol. 43, no. 5, pp. 775–780. https://doi.org/10.1111/j.1095-8649.1993.tb01154.x

Krause, J. and Godin, J.-G.J., Predator preferences for attacking particular prey group sizes: Consequences for predator hunting success and prey predation risk, Anim. Behav. , 1995, vol. 50, no. 2, pp. 465–473. https://doi.org/10.1006/anbe.1995.0260

Krause, J., Godin, J.-G.J., and Rubenstein, D., Group choice as a function of group size differences and assessment time in fish: The influence of species vulnerability to predation, Ethology , 1998, vol. 104, no. 1, pp. 68–74. https://doi.org/10.1111/j.1439-0310.1998.tb00030.x

Kronauer, D.J.C. and Levine, J.D., The ultimate and proximate underpinnings of social behavior, J. Exp. Biol. , 2017, vol. 220, no. 1, pp. 4–5. https://doi.org/10.1242/jeb.152785

Kropotkin, P.A., Vzaimopomoshch’ kak faktor evolyutsii (Mutual Aid as a Factor in Evolution), Moscow: Samoobrazovanie, 2011.

Kukhorenko, K.G., Formation of defensive reflexes on fishing gear in mackerel in the Atlantic, in Izuchenie povedeniya ryb v svyazi s sovershenstvovaniem tekhniki ikh lova (Study of Fish Behavior in Connection with the Improvement of Their Fishing Technique), Moscow: Nauka, 1977, pp. 91–97.

Lange, N.O., Development stages of roach in various ecological conditions, Tr. Inst. Morfol. Zhivotnykh Akad. Nauk SSSR , 1960a, no. 28, pp. 5–40.

Lange, N.O., Stages of development of the Kuban and Don roach Rutilus rutilus heckeli (Nordmann) and vobla Rutilus rutilus caspicus (Jakowlew), Tr. Inst. Morfol. Zhivotnykh Akad. Nauk SSSR , 1960b, no. 25, pp. 47–98.

Lawson, G.L., Kramer, D.L., and Hunte, W., Size-related habitat use and schooling behavior in two species of surgeonfish ( Acanthurus bahianus and A. coeruleus ) on a fringing reef in Barbados, West Indies, Environ. Biol. Fish. , 1999, vol. 54, no. 1, pp. 19–33. https://doi.org/10.1023/A:1007477527663

Lebedev, N.V., Elementarnye populyatsii ryb (Elementary Fish Populations), Moscow: Pishch. Prom-st, 1967.

LeBoeuf, A.C., Benton, R., and Keller, L., The molecular basis of social behavior: Models, methods and advances, Curr. Opin. Neurobiol. , 2013, vol. 23, no. 1, pp. 3–10. https://doi.org/10.1016/j.conb.2012.08.008

Leis, J.M. and Carson-Ewart, B.M., Complex behaviour by coral-reef fish larvae in open-water and near-reef pelagic environments, Environ. Biol. Fish. , 1998, vol. 53, no. 3, pp. 259–266. https://doi.org/10.1023/A:1007424719764

Luo, J., Tidal transport of the bay anchovy, Anchoa mitchilli , in darkness, J. Fish. Biol. , 1993, vol. 42, no. 4, pp. 531–539. https://doi.org/10.1111/j.1095-8649.1993.tb00357.x

Maeda, N. and Hidaka, T., Ethological function of the parr marks in a Japanese trout, Oncorhynchus masou f. Ishikawai, Zool. Mag. , 1979, vol. 88, no. 1, pp. 34–42.

Magurran, A.E. and Bendelow, J.A., Conflict and co-operation in White Cloud Mountain minnow schools, J. Fish. Biol. , 1990, vol. 37, no. 1, pp. 77–83. https://doi.org/10.1111/j.1095-8649.1990.tb05929.x

Magurran, A.E. and Seghers, B.H., Population differences in the schooling behaviour of newborn guppies, Poecilia reticulata, Ethology , 1990, vol. 84, no. 4, pp. 334–342. https://doi.org/10.1111/j.1439-0310.1990.tb00807.x

Magurran, A.E. and Seghers, B.H., Predator inspection behaviour covaries with schooling tendency amongst wild guppy, Poecilia reticulata , populations in Trinidad, Behaviour , 1994, vol. 128, nos. 1–2, pp. 121–134. https://doi.org/10.1163/156853994X00073

Magurran, A.E., Seghers, B.H., Shaw, P.W., and Carvalho, G.R., Schooling preferences for familiar fish in the guppy, Poecilia reticulata, J. Fish. Biol. , 1994, vol. 45, no. 3, pp. 401–406. https://doi.org/10.1111/j.1095-8649.1994.tb01322.x

Manteifel’, B.P., Ekologicheskie i evolyutsionnye aspekty povedeniya zhivotnykh (Ecological and Evolutionary Aspects of Animal Behavior), Moscow: Nauka, 1987.

Marti, Yu.Yu., Promyslovaya razvedka ryby (Commercial Reconnaissance of Fish), Moscow: Pishchepromizdat, 1948.

Martignac, F., Daroux, A., Bagliniere, J.-L., et al., The use of acoustic cameras in shallow waters: New hydroacoustic tools for monitoring migratory fish population. A review of DIDSON technology, Fish Fish. , 2015, vol. 16, no. 3, pp. 486–510. https://doi.org/10.1111/faf.12071

McNicol, R.E. and Noakes, D.L.G., Environmental influences on territoriality of juvenile brook char Salvelinus fontinalis in a stream environment, Environ. Biol. Fish. , 1984, vol. 10, nos. 1–2, pp. 29–42. https://doi.org/10.1007/BF00001660

Mesyatsev, I.I., The structure of gregarious fish schools, Izv. Akad. Nauk SSSR. Ser. Biol. , 1937, no. 3, pp. 745–769.

Miller, N. and Gerlai, R., From schooling to shoaling: Patterns of collective motion in zebrafish ( Danio rerio ), PLOS ONE , 2012, vol. 7, no. 11, Article e48865. https://doi.org/10.1371/journal.pone.0048865

Mochek, A.D., Etologicheskaya organizatsiya pribrezhnykh soobshchestv morskikh ryb (Ethological Organization of Coastal Marine Fish Communities), Moscow: Nauka, 1987.

Moller, P., Electric signals and schooling behavior in a weakly electric fish, Marcusenius cyprinoides L. (Mormyriformes), Science , 1976, vol. 193, no. 4254, pp. 697–699. https://doi.org/10.1126/science.948747

Moreno, G., Dagorn, L., Capello, M., et al., Fish aggregating devices (FADs) as scientific platforms, Fish. Res. , 2016, vol. 178, pp. 122–129. https://doi.org/10.1016/j.fishres.2015.09.021

Morrow, J.E., Schooling behaviour in fishes, Q. Rev. Biol. , 1948, vol. 23, no. 1, pp. 27–38. https://doi.org/10.1086/396078

Neat, F.C., Aggregating behaviour, social interactions and possible spawning in the deep-water fish Coryphaenoides rupestris , J. Fish. Biol. , 2017, vol. 91, no. 3, pp. 975–980. https://doi.org/10.1111/jfb.13386

Nikol’skii, G.V., On the biological significance of a school in fish, Tr. soveshch. po voprosam povedeniya i razvedki ryb (1953) (Proc. Meeting on the Behavior and Intelligence of Fish (1953)), Moscow: Akad. Nauk SSSR, 1955, pp. 104–107.

Nikol’skii, G.V., Ekologiya ryb (Fish Ecology), Moscow: Vyssh. Shk., 1974.

Panov, E.N., Individual’noe. Kollektivnoe. Sotsial’noe v prirode i v obshchestve. Begstvo ot odinochestva (Individual. Collective. Social in Nature and in Society. Escape from Loneliness), Moscow: LKI, 2011a.

Panov, E.N., Znaki, simvoly, yazyki. Kommunikatsiya v tsarstve zhivotnykh i v mire lyudei (Signs, Symbols, Languages. Communication in the Animal Kingdom and in the Human World), Moscow: URSS, 2011b.

Parr, A.E., A contribution to the theoretical analyses of the schooling behaviour of fishes, Occ. Pap. Bingham Oceanogr. Coll. , 1927, vol. 1, pp. 1–32.

Partridge, B.L., The effect of school size on the structure and dynamics of minnow school, Anim. Behav. , 1980, vol. 28, no. 1, pp. 68–77. https://doi.org/10.1016/S0003-3472(80)80009-1

Partridge, B.L., The structure and function of fish schools, Sci. Am. , 1982, vol. 246, no. 6, pp. 114–123. https://doi.org/10.1038/scientificamerican0682-114

Pavlov, D.S. and Kasumyan, A.O., Patterns and mechanisms of schooling behavior in fish: A review, J. Ichthyol. , 2000, vol. 40, no. Suppl. 2, pp. 163–S231.

Pavlov, D.S. and Kasumyan, A.O., Stainoe povedenie ryb (Schooling Behavior of Fish), Moscow: Mosk. Gos. Univ., 2003.

Paxton, C.G.M., Shoaling and activity levels in Corydoras, J. Fish. Biol. , 1997, vol. 51, no. 3, pp. 496–502. https://doi.org/10.1111/j.1095-8649.1997.tb01507.x

Peuhkuri, N., Size-assortative shoaling in fish: The effect of oddity on foraging behaviour, Anim. Behav. , 1997, vol. 54, no. 2, pp. 271–278. https://doi.org/10.1006/anbe.1996.0453

Peuhkuri, N., Shoal composition, body size and foraging in sticklebacks, Behav. Ecol. Sociobiol. , 1998a, vol. 43, nos. 4–5, pp. 333–337. https://doi.org/10.1007/s002650050499

Peuhkuri, N., Body size and food-patch finding in a stickleback shoal, J. Fish. Biol. , 1998b, vol. 53, no. 3, pp. 687–689. https://doi.org/10.1111/j.1095-8649.1998.tb01011.x

Pitcher, T.J., Heuristic definitions of shoaling behaviour, Animal Behav. , 1983, vol. 31, no. 2, pp. 611–613. https://doi.org/10.1016/S0003-3472(83)80087-6

Pitcher, T.J. and Parrish, B.L., Functions of shoaling behavior in teleosts, in Behaviour of Teleost Fishes , London: Chapman and Hall, 1993, pp. 262–439.

Book   Google Scholar  

Pyanov, A.I., Fish learning in response to trawl fishing, ICES Mar. Sci. Symp. , 1993, vol. 196, pp. 12–16.

Radakov, D.V., On the importance of a school for predatory fish in catching prey, in Pitanie khishchnykh ryb i ikh vzaimootnosheniya s kormovymi organizmami (Nutrition of Predatory Fish and Their Relationship with Forage Organisms), Moscow: Nauka, 1965, pp. 173–178.

Radakov, D.V., Stainost’ ryb kak ekologicheskoe yavlenie (Schooling of Fish as an Ecological Phenomenon), Moscow: Nauka, 1972.

Radakov, D.V., Schooling in the Ecology of Fish , New York: John Wiley, 1973.

Radakov, D.V. and Darkov, A.A., On the interaction of fish in schools in the zone of action of the trawl model, in Izuchenie povedeniya ryb v svyazi s sovershenstvovaniem tekhniki ikh lova (Study of Fish Behavior in Connection with the Improvement of Their Fishing Technique), Moscow: Nauka, 1977, pp. 135–140.

Rangeley, R.W. and Kramer, D.L., Tidal effects on habitat selection and aggregation by juvenile pollock Pollachius virens in the rocky intertidal zone, Mar. Ecol. Proc. Ser. , 1995, vol. 126, pp. 19–29. https://doi.org/10.3354/meps126019

Ranta, E. and Kaitala, V., School size affects individual feeding success in three-spined sticklebacks ( Gasterosteus aculeatus L.), J. Fish. Biol. , 1991, vol. 39, no. 5, pp. 733–737. https://doi.org/10.1111/j.1095-8649.1991.tb04402.x

Ranta, E. and Lindström, K., Assortative schooling in three-spined sticklebacks?, Ann. Zool. Fennici , 1990, vol. 27, no. 2, pp. 67–75.

Ranta, E., Juvonen, S.K., and Peuhkuri, N., Further evidence for the size-assortative schooling in sticklebacks, J. Fish. Biol. , 1992a, vol. 41, no. 4, pp. 627–630. https://doi.org/10.1111/j.1095-8649.1992.tb02689.x

Ranta, E., Lindström, K., and Peuhkuri, N., Size matters when three-spined sticklebacks go to school, Anim. Behav. , 1992b, vol. 43, no. 1, pp. 160–162. https://doi.org/10.1016/S0003-3472(05)80082-X

Sakakura, Y. and Tsukamoto, K., Social rank in schools of juvenile yellowtail, Seriola quinqueradiata, J. Appl. Ichthyol. , 1998, vol. 14, nos. 1–2, pp. 69–73. https://doi.org/10.1111/j.1439-0426.1998.tb00616.x

Sbikin, Yu.N., Some aspects of the social and defensive behavior of young sturgeons (Acipenseridae), Zool. Zh. , 1996, vol. 75, no. 3, pp. 383–390.

Seghers, B.H., Schooling behavior in the guppy ( Poecilia reticulata ): An evolutionary response to predation, Evolution , 1974, vol. 28, no. 3, pp. 486–489. https://doi.org/10.2307/2407174

Seghers, B.H., Facultative schooling behavior in the spottail shiner ( Notropis hudsonius ): possible costs and benefits, Environ. Biol. Fish. , 1981, vol. 6, no. 1, pp. 21–24. https://doi.org/10.1007/BF00001795

Shaw, E., some new thoughts on the schooling of fishes, FAO Fish Rep. , 1969, vol. 2, no. 62, pp. 218–231.

Shaw, E., Schooling fishes: The school, a truly egalitarian form of organization in which all members of the group are alike in influence, offers substantial benefits to its participants, Am. Sci. , 1978, vol. 66, no. 2, pp. 166–175.

Shelton, D.S., Shelton, S.G., Daniel, D.K., et al., Collective behavior in wild zebrafish, Zebrafish , 2020, vol. 17, no. 4, pp. 243–252. https://doi.org/10.1089/zeb.2019.1851

Article   PubMed   PubMed Central   Google Scholar  

Shiraishi, M., Azuma, N., and Aoki, I., Large schools of Japanese sardines, Sardinops melanostictus , mate in single pair units at night, Environ. Biol. Fish. , 1996, vol. 45, no. 4, pp. 405–409. https://doi.org/10.1007/BF00002533

Soin, S.G., Kasumyan, A.O., and Pashchenko, N.I., Ecological and morphological analysis of the development of the minnow, Phoxinus phoxinus (Cyprinidae), J. Ichthyology , 1981, vol. 21, no. 4, pp. 90–105.

Social Behaviour: Genes, Ecology and Evolution , Cambridge: Cambridge Univ. Press, 2010.

Tinbergen, N., Sotsial’noe povedenie zhivotnykh (Social Behavior of Animals), Moscow: Mir, 1993.

Van Havre, N. and Fitzgerald G.J., Shoaling and kin recognition in the threespine stickleback ( Gasterosteus aculeatus L.), Biol. Behav. , 1988, vol. 13, no. 4, pp. 190–201.

Vasnetsov, V.V., Ereemeva, E.F., Lange, N.O., et al., Stages of development of commercial semi-anadromous fishes of the Volga and Don – bream, carp, vobla, roach and pike perch, Tr. Inst. Morfol. Zhivotnykh Akad. Nauk SSSR , 1957, no. 17, pp. 7–79.

Williams, G.C., Measurements of consociation among fishes and comments on the evolution of schooling, Publ. Mich. State Univ. Mus. Biol. Ser. , 1964, vol. 2, no. 7, pp. 349–384.

Wilson, E.O., Sociobiology: The New Synthesis , Cambridge, MA: Belknap Press of Harvard Univ. Press, 1975.

Zilanov, V.K., Features of the behavior of some fish species of the North Atlantic in the zone of artificial light, in Izuchenie povedeniya ryb v svyazi s sovershenstvovaniem tekhniki ikh lova (Study of Fish Behavior in Connection with the Improvement of Their Fishing Technique), Moscow: Nauka, 1977, pp. 40–46.

Download references

ACKNOWLEDGMENTS

The authors express their sincere gratitude to A.A. Kazhlaev and L.S. Alekseeva (Moscow State University), who provided great assistance in preparing the article for publication. The authors are sincerely grateful to P.I. Kirillov (Institute of Ecology and Evolution, Russian Academy of Sciences) for careful and constructive editing of the text and illustrations, which improved the quality of the article.

The article was prepared within the framework of scientific projects of the state assignment of the Moscow State University No. 121032300100-5 and the Institute of Ecology and Evolution, Russian Academy of Sciences No. 121122300056-3 in the Unified State Information System for Accounting the Results of Civil Research, Development and Technological Works.

Author information

Authors and affiliations.

Lomonosov Moscow State University, Moscow, Russia

A. O. Kasumyan

Severtsov Institute of Problems of Ecology and Evolution, Russian Academy of Science (IPEE RAS), Moscow, Russia

D. S. Pavlov

You can also search for this author in PubMed   Google Scholar

Corresponding author

Correspondence to A. O. Kasumyan .

Ethics declarations

Ethics approval and consent to participate.

All applicable international, national, and/or institutional guidelines for the care and use of animals were followed.

CONFLICT OF INTEREST

The author of this work declares that he has no conflicts of interest.

Additional information

Translated by S. Avodkova

Publisher’s Note.

Pleiades Publishing remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Reprints and permissions

About this article

Kasumyan, A.O., Pavlov, D.S. Schooling Behavior of Fish: General Ideas, Terms and Concepts, Prevalence, Applied Aspects. J. Ichthyol. 63 , 1219–1233 (2023). https://doi.org/10.1134/S0032945223070020

Download citation

Received : 02 March 2023

Revised : 30 March 2023

Accepted : 30 March 2023

Published : 14 December 2023

Issue Date : December 2023

DOI : https://doi.org/10.1134/S0032945223070020

Share this article

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

• schooling behavior
• terminology
• aggregation
• Find a journal
• Publish with us
• Track your research

• Physics of Living Systems

Energy conservation by collective movement in schooling fish

• Yangfan Zhang author has email address

George V. Lauder

• Department of Organismic and Evolutionary Biology, Harvard University, 26 Oxford St, Cambridge, Massachusetts, 02138, USA
• https://doi.org/ 10.7554/eLife.90352.2
• Open access
• Copyright information

eLife assessment

The authors provide an important series of metabolic measurements characterizing group dynamics in fish, rationalizing that schooling behavior presents several benefits. The strength of evidence supporting this conclusion is solid , but the specific methodological and analytical approaches taken should be considered for further interpretation.

• https://doi.org/ 10.7554/eLife.90352.2.sa2

Significance of findings

important : Findings that have theoretical or practical implications beyond a single subfield

• fundamental

Strength of evidence

solid : Methods, data and analyses broadly support the claims with only minor weaknesses

• exceptional

During the peer-review process the editor and reviewers write an eLife assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife assessments

Many animals moving through fluids exhibit highly coordinated group movement that is thought to reduce the cost of locomotion. However, direct energetic measurements demonstrating the energy-saving benefits of fluid-mediated collective movements remain elusive. By characterizing both aerobic and anaerobic metabolic energy contributions in schools of giant danio ( Devario aequipinnatus ), we discovered that fish schools have a concave upward shaped metabolism– speed curve, with a minimum metabolic cost at ∼1 body length s -1 . We demonstrate that fish schools reduce total energy expenditure (TEE) per tail beat by up to 56% compared to solitary fish. When reaching their maximum sustained swimming speed, fish swimming in schools had a 44% higher maximum aerobic performance and used 65% less non-aerobic energy compared to solitary individuals, which lowered the TEE and total cost of transport by up to 53%, near the lowest recorded for any aquatic organism. Fish in schools also recovered from exercise 43% faster than solitary fish. The non-aerobic energetic savings that occur when fish in schools actively swim at high speed can considerably improve both peak and repeated performance which is likely to be beneficial for evading predators. These energetic savings may underlie the prevalence of coordinated group locomotion in fishes.

One-Sentence Summary

Fish schools showed a U -shaped metabolism-speed curve and reduced the energy use per tail beat up to 56% at high swimming speeds compared to solitary fish.

Introduction

Newton’s laws of motion underpin animal locomotion, from fine-scale maneuvers to long-distance migrations ( 1 ). As speed increases or during migration involving sustained movement over long distances, animals often move in coordinated groups: e.g. , migratory birds in V-formation ( 2 ), ducklings swimming in formation ( 3 ), cyclists in a peloton ( 4 ), elite marathon runners ( 5 ), and fish schools ( 6 ). To overcome gravitational or fluid dynamic resistance to forward motion, animals use both aerobic and glycolytic metabolism (oxidative and substrate-level phosphorylation at the cellular level) to generate energy and sustain movement in the air, on land, and in the water. In water, a fluid that is 50 times more viscous than air and contains much less O 2 per Kg than air, the need for aquatic animals to reduce fluid dynamic drag for energy conservation is even greater than for aerial or terrestrial locomotion. Hence, we use fish schooling behaviour as a model system to explore whether the fluid dynamics of collective movement (active and directional movement of animal groups along a common trajectory ( 7 )) can enable energetic savings compared to locomotion by a solitary individual. Fish schooling behaviours are some of the most prominent social and group activities exhibited by vertebrates.

Fish may school for many reasons, including foraging, reproduction, and migration, and fish schools often exhibit high-speed movement during evasion from predators. Hence, natural selection can put pressure on the formation of schools and individual interactions within a school. Besides ‘safety in numbers’ where individuals in a large school have a lower probability of becoming prey, can there be more tangible energetic benefits of high-speed swimming as a school? Since fluid drag scales as velocity squared ( 8 ), movement at higher speeds places a premium on mechanisms to conserve energy. Fast and unsustainable swimming is fueled by aerobic metabolism with a major contribution from anaerobic glycolysis ( 9 ) ( 10 ). During high-speed swimming aerobic capacity is maximized ( 10 ), and energy use is more likely to compete with other activities than in low-speed swimming which only occupies a small portion of aerobic capacity. The direct consequences of unsustainably engaging anaerobic glycolysis include both fatigue and metabolic perturbation (such as changing blood acidity) that need time to recover after bouts of intense locomotion. Animals become vulnerable during recovery because of their hindered ability to repeat peak locomotor performance in the presence of predators.

Another factor of importance in the experimental evaluation of the energetic cost of collective locomotion in fishes is that endurance and efficiency are often key to achieving lifetime fitness, especially in fish species that undergo significant migrations. Fish groups usually undergo long-distance migrations at slow speeds ranging from 0.25–1.0 BL s -1 as recorded by tags on migrating fish (Fig. S1). The common migratory speed of ∼1 BL s -1 for marine and anadromous fish could be related to the minimum aerobic cost of transport ( 11 ). However, there is currently no direct measurement of the absolute cost of teleost fish locomotion demonstrating a minimum group metabolic rate ( M O 2 ) at an optimal speed ( U opt ). Fish schools in nature swim over a wide range of speeds (Fig. S1) but no study has directly characterized the total energy expenditure (TEE), and accounted for the potentially substantial anaerobic metabolic costs involved in higher speed locomotion. As a result, the swimming TEE performance curve of fish schools from the minimum swimming speed to the critical swimming speed ( U crit ) has not yet been measured.

More broadly, despite the widespread notion that collective motion saves energy ( 6 ), very few studies have directly measured the energetic cost of collective movement compared to the cost of movement by solitary individuals. In the canonical case of V-formation flight in birds, energetic saving is exclusively inferred from indirect measurements, e.g., heart rate ( 2 ), and flapping frequency and phase ( 12 ) ( Table 1 ). There are no direct measurements of metabolic energy consumption of collective movement in birds, and indeed this would be extremely challenging to accomplish. Although some energetic aspects of schooling fishes have been studied ( 13 ) ( 14 ) ( 15 ) ( 16 ) ( 17 ), these studies have analyzed a limited range of speeds and focused on aerobic metabolism ( see Table 1 ). In fact, several field and laboratory kinematic studies suggest that both bird flocks and fish schools do not always conserve aerobic energy ( 17 ) ( 18 ) ( 19 ) and indicate that moving in a group can actually involve increased costs. Without directly measuring the energetic cost of movement, inferences solely based on kinematics do not include complex group interactions, and the interaction between collective behaviour and fluid dynamics. Consequently, we do not yet have a holistic and mechanistic view of where and how the collective movement by fish might conserve energy.

A summary of experimental studies that directly estimate the energetic saving of group movement in aerial and aquatic vertebrates that move through a fluid environment.

No energetic measurements have made for freely-moving bird flocks. Three studies measured the energetics of fish schools over a range of narrow speeds; two other studies measured the energetics of fish schooling at one speed. But no studies have quantified both the aerobic and anaerobic energetic cost of swimming as a group compared to solitary locomotion.

For swimming fish, fluid dynamic experiments have shown how collective movement can improve swimming efficiency due to interactions among deforming bodies and through interactions between moving animals and the local fluid environment (these energy-saving mechanisms are summarized in Fig. 1 ). Our current understanding of energy saving mechanisms during collective locomotion in fishes is largely based on computational fluid dynamic models with a few analyses using robotic systems ( 6 ) ( 20 ) ( 21 ) ( 22 ) ( 23 ) at low speed, but the link between such models and fish metabolic energy saving over a wide speed range is not yet established. We expect that the need to conserve energy in animals moving against air or water resistance should be greater at higher speeds. We hypothesize that fish in schools can reduce the total cost of high-speed locomotion relative to solitary movement, leading to reduced non-aerobic energetic costs and time of recovery for schooling fish swimming at speeds occupying the majority (> 50%) of their aerobic capacity. Our secondary objective is to determine the differences (if any) between the swimming performance curves (metabolism versus speed) of solitary fish and fish schools.

A summary of biomechanical principles underlying proposed hydrodynamic advantages of schooling behaviour in fish.

When fish swim into free-stream flow ( U fs ), experimental data show that fish schools are dynamic with fish changing position frequently. Regardless of fish position within a school and changes in relative position, theoretical and robotic analyses have demonstrated that at least four mechanisms (indicated by numbers) provide an advantage in the form of reduced power consumption. 1. Reduced oncoming velocity ( U 2 ) requires less thrust for a fish swimming in the wake between two leading fish ( 6 ); 2. The Knoller-Betz effect of leading edge suction reduces costs for a trailing fish due to accelerated oscillating flow at the head ( U 1 ) ( 24 ) ( 78 ); 3. Added mass “push” from follower to leader can reduce costs for the leader in front of another fish ( 79 ) ( 23 ) ( 78 ); 4. Wall effects benefit neighbouring fish where swimming next to another fish reduces swimming costs ( 80 ) ( 20 ). These principles suggest that regardless of the relative positions of the individuals within the fish school, the fish school as a collective unit should be able to save metabolic energy compared to a solitary fish swimming in U fs .

To evaluate these hypotheses, we directly measured both the aerobic and non-aerobic energy used by schooling fishes over a wide range of water velocities (0.3–8.0 body lengths s -1 ; BL s -1 ), and then compared the swimming performance curve to that of solitary fish. We equipped a high-resolution ( see Methods for our approach to enhancing the signal-to-noise ratio in our respirometer) swim-tunnel respirometer with two orthogonal high-speed cameras to quantify three-dimensional (3-D) fish kinematics. This enabled us to simultaneously measure energetics, 3-D dynamics of fish schools (n=5 replicate schools), and the kinematics of solitary fish (n=5) of a model species, giant danio ( Devario aequipinnatus ), that exhibits an active directional group swimming from near still water to maximum sustained speeds (equivalent to a Reynolds number range of 6.4•10 3 to 1.8•10 5 ; Fig. S2), using a U crit test.

We discovered that both solitary fish and fish schools have a concave upward shaped aerobic metabolic rate ( M O 2 )-speed curve over the lower portion of their speed range (0.3–3 BL s -1 , Fig. 2C, D ). The aerobic locomotor cost at 0.3–3 BL s -1 showed no statistical difference between solitary fish and fish schools ( F 1,80 ≤ 2.40 , p ≥ 0.13). Fish schools swimming at ∼1.0 BL s -1 consume less energy than at slower speeds ( F 9,40 = 24.7, p ≤ 0.0007), while swimming at 3 BL s -1 consumes a similar amount of energy to maintain position at a water velocity of 0.3 BL s -1 ( p = 0.85; Fig. 2D ). Danio schools have a minimum aerobic cost ( M O 2min = 212.9 mg O 2 kg -1 h -1 ) of 1.25 BL s -1 . M O 2min at this optimal speed ( U opt ) was lower than both the M O 2 of aggregating behaviours (328.5 mg O 2 kg -1 h -1 ) exhibited before U crit test (Table S1, Fig. S3) and for the group at rest in still water M O 2 (267.9 mg O 2 kg -1 h -1 ) (by 35% and 21% respectively; F 2,12 = 16.7, p ≤ 0.02, Fig. 2A, B ). Swimming at U opt reduces the energetic cost below that of swimming at either slower or higher speeds. Indeed, the Strouhal number of fish schools decreased from 2.1 at the lowest speed to 0.3 at the energetic minimum U opt , and then increased to 0.4 at the higher U crit .

Measurements of aerobic and anaerobic locomotor cost of fish schools and solitary fish.

( A ) Average traces of metabolic rate ( M O 2 ) of fish schools over a 40-hour experiment. Following the first 18-h quiescent state, a critical swimming speed ( U crit ) test quantifies the aerobic cost of active swimming. The ensuing 18-h measurement of excess post-exercise oxygen consumption (EPOC) quantifies the anaerobic cost. ( B ) Comparison of M O 2 for conditions of aggregating behaviour, minimum demand speed, and resting condition with minimal flow ( M O 2aggregate , M O 2min , M O 2rest ) ( C ) Comparisons of concave upward shaped M O 2 -speed curve over the entire range (0.3–8 body length s -1 , BL s -1 ) and ( D ) the concave shaped M O 2 -speed curve at the lower speeds (0.3–3 BL s -1 ). ( E, F) Percentage (%) aerobic scope used by fish schools and solitary fish during the U crit test. ( G , H ) Comparisons of EPOC and EPOC durations between fish schools and solitary fish. Statistical significance is denoted by asterisk(s). Green colour = school data (n=5 schools); blue colour = solitary fish data (n=5 individual fish); shading indicates the 95% confidence interval. Statistical details are available in the statistical analyses section.

We also characterized the kinematic features that result in elevated M O 2 at low water velocities. Danio maintaining body position in low-velocity aggregations had a higher tail beat frequency ( f TB ) than solitary fish ( F 1,80 ≥ 9.8, p ≤ 0.002). However, at U opt and the highest speed of active directional schooling, fish in schools had a lower f TB ( F 1,80 ≥ 4.6, p ≤ 0.035) than solitary fish ( Fig. 4B ). The 3-D angular heading of schooling fish transitioned from omnidirectional to pointing against water flow when water velocity is above 0.75 BL s -1 , and individual fish orient into the flow starting at 0.3 BL s -1 ( Fig. 4C ). Although the body angle and turning frequency of solitary fish and schools decreased with water velocity, schooling fish had a higher turning frequency ( F 1,80 ≥ 15.6, p < 0.002) below 0.75 BL s -1 than solitary fish ( Fig. 4D,E ).

Modeling of simultaneous aerobic and non-aerobic costs of fish schools and solitary fish for a critical swimming speed ( U crit ) test. (A) Modeling the O 2 cost of the metabolic rate ( M O 2 )-speed curve and the ensuing recovery cost (excess post-exercise oxygen consumption, EPOC) as a function of speed. After U crit , fish returned to the same resting M O 2 ( M O 2rest ) as a pre-test. ( B, C ) In addition to M O 2 (solid line & filled symbols), we modeled the total O 2 cost (dash line & half-filled symbols) for fish schools and solitary fish and when performing the U crit swimming test. The estimated partitioning of aerobic and non-aerobic contributions to swimming are denoted (red-&-bold) with respect to speed for 4–8 BL s -1 is shown below each graph. ( D, E ) Using total O 2 cost, we computed total energy expenditure (10-min period per point) and the total cost of transport (including both aerobic metabolism, high-energy phosphate, and anaerobic glycolysis) for both fish schools and solitary fish. Statistical significance is denoted by asterisk(s). Green colour = school data (n=5); blue colour = solitary fish data (n=5); shading indicates the 95% confidence interval. See methods for modeling and statistical details.

Three-dimensional characterization of swimming kinematics and fish schooling dynamics as a function of speed.

(A) Total energy expenditure (TEE) per tail beat, ( B ) Tail beat frequency, ( C ) The angle of fish to free-stream water flow, measured as the mean and the S.D. of the angles of the individuals within the school. ( D ) Three-dimensional angle of fish to the frontal plane. ( E ) Turning frequency, ( F ) Three-dimensional school length. ( G ) Three-dimensional distances among all individuals in the school and the S.D. of the distance. The upper and lower boundaries of the metrics are indicated. Statistical significance is denoted by asterisk(s). Green colour = school data (n=3-4); blue colour = solitary fish data (n=3-4); shading indicates the 95% confidence interval. See methods for details of three-dimensional reconstruction and statistics.

We discovered that, across the entire 0.3–8 BL s -1 range, M O 2 -speed curves of fish schools are concave upward shaped (reached 1053.5 mg O 2 kg -1 h -1 at 8 BL s -1 ) whereas solitary fish showed an upward concave curve that reaching a plateau of 760.2 mg O 2 kg -1 h -1 (∼10% CV in 6–8 BL s -1 , Fig. 2C ), demonstrating that schooling dynamics results in a 44% higher maximum aerobic performance ( F 1,80 = 30.0, p < 0.001). This increased maximum aerobic performance translated to a lower use of aerobic capacity compared to solitary fish over 0.3–8 BL s -1 ( see the discussion for how fish swimming in school can improve aerobic performance). Collectively, fish schools used a 36% lower proportion of their aerobic scope than solitary fish (Wilcoxon test: p = 0.0002, Fig. 2F ). Fish schools used a 38% lower proportion of aerobic scope (34 vs. 55%) at 4 BL s -1 (50% U crit & > 50% aerobic scope onward), and consistently used ∼25% lower proportion of their aerobic scope at 6 and 7 BL s -1 than solitary fish ( F 1,103 ≥ 4.8, p ≤ 0.03, Fig. 2E ).

Given that fish schools had a higher maximum aerobic performance, generally used a lower proportion of their aerobic scope, and individuals within the schools have a 14% lower f TB than the solitary fish at 8 BL s -1 (11.5 vs. 13.3 Hz, F 1,80 = 15.1, p < 0.001), we predicted that fish schools, compared with solitary fish, use less anaerobic energy to supplement aerobic energy for the high-speed movement approaching their aerobic limit. Hence, we measured post-exercise O 2 utilization, a majority of which is used to restore high-energy phosphate storage and glycolytically induced metabolic perturbations (EPOC: Excess Post-exercise O 2 Consumption). We discovered that fish schools have a 65% lower EPOC (0.69 vs. 1.95 mg O 2 ; t 8 = 4.5, p = 0.0021), and recover 43% faster than solitary fish (8 vs. 14 h; t 8 = 2.8, p = 0.025; Fig. 2G, H ).

To estimate the relative proportions of aerobic and anaerobic energy contributions for locomotion at each swimming speed, we modeled EPOC in addition to M O 2 , the aerobic cost (> 50% U crit & aerobic scope, Fig. 3A ) (Table S2). We also estimated the total O 2 cost during the entire swimming process for each school and individual and calculated total energy expenditure (TEE) ( 10 ). The TEE of fish schools was 38–53% lower than that of the solitary fish between 5– 8 BL s -1 ( F 1,103 ≥ 7.4, p ≤ 0.008; Fig. 3D ). TEE of fish schools was only 42–143 % higher than the aerobic metabolic rate at 5–8 BL s -1 ( F 1,96 ≥ 3.5, p ≤ 0.001) and anaerobic metabolic energy only accounted for 29–58 % of the TEE depending on speed ( Fig. 3B ). In contrast, TEE of solitary fish was 131–465 % higher than the aerobic metabolic rate between 5–8 BL s -1 ( F 1,112 ≥ 10.7, p ≤ 0.001), where anaerobic energy accounted for 62–81% of the TEE depending on speed ( Fig. 3C ).

Schooling dynamics reduced the total cost (aerobic plus anaerobic) of transport (TCOT) by an average of 43% compared to swimming alone ( F 1,103 = 6.9, p = 0.01), and most of this energy conservation happens at higher speeds when fish approach their aerobic limit. Schooling dynamics in danio enables an extremely shallow rate of increase in TCOT with speed compared to that of an individual ( Fig. 3E ). The TCOT of solitary fish increased by 490% (6.5 kJ km -1 kg - 1 ) at 8 BL s -1 , whereas the schooling TCOT increased by only 200% (3.6 kJ km -1 kg -1 ) at 8 BL s - 1 . Therefore, individual fish form a more energy-efficient biological entity when they collectively move as a school.

To answer the question of how schooling dynamics reduces TEE, we combine video analysis of fish tail beat kinematics with simultaneous aerobic and anaerobic measurements to compute energy expended per tail beat (TEE•beat -1 ), and compared values for fish in schools to those for solitary fish. Schooling fish reduced TEE•beat -1 by 30–56% at higher speeds compared to solitary fish ( F 1,81 ≥ 7.3, p ≤ 0.008), a substantial reduction in TEE•beat -1 consumed both by the school, and by individual fish within a school ( Fig. 4A ). Notably, the energetic benefits during active directional schooling occur when fish schools become more streamlined, as the length of fish school increased with speed and plateaued beyond 2 BL s -1 ( Fig. 4F ), while the 3-D distance among individuals stayed relatively constant at ∼1.2 BL ( Fig 4G ). These results directly demonstrate schooling dynamics benefits the swimming kinematics of individual fish within the school and results in a net outcome of up to 53% TEE reduction in fish schools compared to solitary fish.

Simultaneous characterization of energetics and kinematics enables an integrated understanding of both the physiology and physics of fish schooling behaviour. Hydrodynamic models, kinematic measurements, and robotic analyses of fish schools indicate that the cost of swimming can be reduced when fish swim beside neighbouring fish ( 20 ), behind or in front of another fish ( 24 ), or behind leading individuals ( 6 ) ( Fig. 1 ). We have demonstrated here that fish schools downshifted an overall concave upward swimming energetics performance curve at higher speeds by ∼43%. Energy use as swimming speed increases changes in a non-linear manner, with higher energy use at the lowest speeds ( Fig. 2D ), and the performance curve below 3.0 BL s -1 is U -shaped. Energy saving occurs through a substantial reduction in the non-aerobic energy contribution when fish approach aerobic limits. One of the key ecological benefits of the reduced use of glycolysis is a faster recovery time from fatigue ( Fig. 2H ). This would enable fish schools in nature to repeat high-performance movements. Considering that high-speed maneuvers are extremely common in predatory evasion ( 25 ) food-searching group motion ( 26 ), and in fish schools in the ocean (Fig. S1), energy savings by fish schooling at high speed and the faster recovery that followed could have a considerable impact on lifetime fitness ( 27 ) ( 28 ). Moreover, considerable energy savings occur even though we observed that individual danio regularly change positions within the school and do not maintain stable inter-individual locations. We expect that this is because individuals within the school can use multiple hydrodynamic mechanisms to save energy depending on their location ( Fig. 1 ). We propose that energetic savings by collective movement may not require either fixed positional arrangements among individuals or specific kinematic features (such as tail beat synchronization) compared to solitary locomotion. Based on these results, we regard fish schooling as a highly robust behaviour that provides considerable energetic benefits that are not necessarily sensitive to specific fish locations or movements. Future studies are necessary to decipher each of the specific hydrodynamic mechanisms for energy saving.

Schooling dynamics enhances aerobic performance and reduces non-aerobic energy use

We discovered that a significant amount of energy conservation for active directional swimming occurs at speeds above 3 BL s -1 . In nature, fish schools routinely exhibit active directional collective locomotion above ∼6 BL s -1 (Fig. S1, Table S3) ( 29 ), a speed that engages anaerobic glycolysis. Yet, there are no previous measurements of the anaerobic cost of schooling in fish. Hence, important and previously unrecognized benefits of active directional collective locomotion are (1) increased aerobic performance, (2) a reduced use of aerobic capacity and anaerobic energy, and (3) a resultant faster recovery from the associated metabolic perturbations and costs of swimming at higher speeds. When animals are approaching their maximum metabolic rate, the highest attainable rate of O 2 uptake limits aerobic performance during the unsteady state of high-speed locomotion and anaerobic glycolysis is engaged to support the maximum metabolic demand.

Fish swimming at higher speeds (> 50% U crit ) routinely exhibit a burst-&-glide gait mode of locomotion ( 30 , 31 ). Repeated bursting is substantially fueled by anaerobic glycolysis and the gliding phase enables peak M O 2 ( 32 ) to replenish the venous O 2 content ( 33 ). Fish within a school are known to increase gliding time and decrease burst time by 19% when they are trailing the leading fish ( 34 ), e.g, a lower f TB of fish schools at 8 BL s -1 is related to the extended gliding time ( Fig. 4B ). We reason that the higher proportion of time spent gliding likely enables more bouts of peak M O 2 and enhances the maximum aerobic performance ( 32 ) ( 33 ). This enhanced maximum aerobic performance increases the metabolic ceiling and enables fish to use a lower proportion of their aerobic scope for locomotion ( Fig. 2E ). Physiological studies of exercise metabolism and locomotion suggest that a lower proportional use of aerobic scope during movement relates to a reduced accumulation of anaerobic end-products ( 9 ). Collectively, our results showed that the increased aerobic performance ( Fig. 2C ) and the reduced use of the aerobic capacity ( Fig. 2E ) for swimming above 50% U crit in fish schools likely plays a key role in reducing the use of glycolysis and the accumulation of lactate ( e.g. EPOC, Fig. 2G ) for high-speed swimming ( 10 , 31 ).

We present direct evidence of substantial non-aerobic energy saving by demonstrating the 65% lower EPOC ( Fig. 2G ) used by fish that swam in schools. Anaerobic glycolysis is crucial in permitting continued movement when aerobic limits are reached at high swimming speeds ( 35 ). Often vertebrates at higher speeds use more fast twitch fibres that generate high-frequency contractile force in part through anaerobic glycolysis ( 36 ) ( 37 ). The higher EPOC for single fish is unlikely to have been confounded by any possible stress effect of swimming as a solitary fish. Stressed solitary fish would have elevated aerobic metabolic rates in low water velocity, which we show is not the case for our experiments ( see Fig. 2D ). Indeed, to mitigate the stress response, we acclimated fish to the solitary conditions and habituated fish to the respirometer for a quiescent day before the experiment ( see more details in experimental animal and protocols). When vertebrates move at high speed, the typical ‘stress’ neurotransmitters ( e.g. adrenaline, catecholamines & cortisol) usually increase ( 38 ) ( 39 ), regardless of whether vertebrates move as an individual or in a group. These endocrinological responses underpin the physiological mechanisms that enable higher aerobic performance ( e.g. rising heart rate) ( 40 ) ( 41 ). Moreover, after the highest swimming speeds immediately followed by an active recovery process, fish likely prioritize the recovery of the essential functions (as the levels of the neurotransmitters decline fairly rapidly) for the potential of repeated locomotion performance (speedy recovery allows fish to restore their swimming performance) ( 31 ) ( 42 ) ( 43 ).

To present a complete energetic profile showing the total energy expenditure (TEE) of locomotion for fish schooling, we model the non-aerobic cost (EPOC) on top of the aerobic swimming performance curve. The model is rooted in the use of aerobic and glycolytic energy when fish swim above 50% U crit ( 30 ) and a commonly observed inflection point of faster anaerobic end-product accumulation is at ∼50% aerobic scope ( 9 ) ( see Methods for detailed physiological bases and criteria for EPOC modeling). As a result, we observed a TEE being ∼2.5 fold higher than the aerobic metabolic energy expenditure which agrees with the theoretical estimates based on individual yearling sockeye salmon ( 44 ). This model elucidates how aerobic and anaerobic metabolic energy constitute TEE, fuel muscles, generate locomotor thrust and overcome hydrodynamic drag during swimming.

We also demonstrated that the need for energy saving in fish schools at lower speeds (< 3 BL s -1 or 38% U crit Fig. 2E , 3D, E ) is not as crucial as at higher speeds where we demonstrated substantial energy savings in the schools compared to solitary individuals. Fluid drag is exponentially less at lower speeds compared to higher speeds. Fish predominately use aerobic metabolism to support low-speed steady swimming for prolonged periods ( 10 ) which does not require lengthy recovery as swimming at higher speeds where glycolysis contributes to the energetic demands. The costs of low-speed swimming are often less than < 20% of aerobic capacity ( Fig. 2E ), which leaves the majority of aerobic capacity (>50%) for other activities. Therefore, the benefits of reducing the energetic cost of locomotion are likely not major factors underlying behaviours such as low-speed milling ( 45 ) and aggregation ( 46 ) in fishes, where other ecological drivers such as feeding, predator avoidance, and reproduction likely dominate the fitness landscape.

Schooling dynamics and energy conservation

Our 3D kinematic analyses shed light on the complex interactions between schooling and hydrodynamics that enable energy saving by fish collective movement. One of the key possible mechanisms is local interactions among individuals, where the kinematics of individual fish respond to the wakes shed by neighbours ( 47 ). Kinematic and simulation studies indicate that 2– 3 coordinated fish can save energy ( 48 ) ( 20 ) through local interactions. We directly demonstrate here that eight coordinated fish can save energy, and recent simulations show that schooling benefits can extend to at least 23 fish ( 21 ). Although these studies appear to suggest that the energetic benefits of schooling are scalable, we are still some ways from proving the exact mechanisms of how the large school size in nature is coordinated and whether or not local interactions are still one of the mechanisms for energy saving when scaling up school size. Future hydrodynamic and long-term 3-D tracking studies are needed to investigate the specific schooling formations and hydrodynamic mechanisms used by live fish for energy saving (as summarized in Fig.1 ) ( 49 ).

We show that danio in a school keep a relatively consistent 3-D distance from neighbours and have small variation in mean position even as individuals routinely change position within the school ( Fig. 4G ). Computational fluid dynamic analyses of fish schools show that a small 3D neighbouring distance ( e.g. < 0.4 BL s -1 ) increases drag and can increase swimming costs, whereas a large 3-D neighbouring distance ( e.g. > 2 BL s -1 ) reduces the hydrodynamic benefits ( 48 ). There may thus be an optimal mean distance among individuals within a school to maximize energetic savings. Whether inter-individual distance should change as swimming speed increases remains an open question and may be a function of school size. While we cannot completely preclude the possibility that elongation of danio fish school might be due to the weaker individuals falling behind, the fact that the school maintained a stable 3-D distance as speed increased and that fish continued to change position within the school suggests that the elongated school volume at higher speeds may reflect changing hydrodynamic interactions among individuals ( Fig. 1 ).

To encapsulate the complex interaction among animal physiology, kinematics and fluid dynamics, more comprehensive quantification beyond simple kinematic metrics are necessary. For example, the wing flapping frequency of birds in a V-formation can be higher than for solitary flight despite strong indirect evidence of overall energy saving by V-formation flying ( 2 ) ( 12 ) ( 19 ), and flapping phase of adjacent birds may be informative for understanding the fluid dynamic advantages of V-formation flying ( 12 ). Likewise, phase matching of body motion between neighbouring fish can help individual fish within the school to boost thrust or reduce drag ( 20 ). Fish can also adjust their body stiffness and maintain f TB while reducing the amount of muscle activity needed to generate movement ( 50 ) ( 51 ) ( 52 ). Thus, we did not observe consistent and substantial changes in the f TB of fish within the schools compared to solitary swimmers despite demonstrating substantial energetic savings by the group at high speed. Although reduced f TB has been used in past studies as an indirect indicator of energy saving in fishes ( 53 ) ( 54 ), we caution that f TB is not necessarily an indicator of energy use in group dynamics where fish constantly and dynamically interact with complex vortices and flow generated by other individuals. By simultaneously measuring kinematics and energetics, we discovered a downshifted performance curve of TEE per tail beat in fish schools at higher speeds ( Fig. 4A ), even with limited alteration in f TB at lower swimming speeds within the school. This suggests that the locomotor muscle fibres in the body musculature of fish within the school need to generate less force as indicated by the lower measured TEE. Fish can possibly fine-tune undulatory motion to harness kinetic energy from nearby vortices ( 20 ) ( 48 ), reducing biological energy contribution for thrust.

Aerobic metabolic rate–speed curve of fish schools

We discovered that the aerobic metabolic rate ( M O 2 )–speed curve at speeds less than 3 BL s -1 is U -shaped, with swimming at 0.3 BL s -1 utilizing the same aerobic energy as moving at 3 BL s -1 ( Fig. 2D ). The entire aerobic energy-speed performance curve is concave upward in an extended J -shape ( Fig. 2C ), and non-linear, in contrast to several previous analyses of energy use with speed in fish which have shown linear relationships ( 55 ) ( 56 ) ( 57 ). By quantifying energy use over a wide speed range and by measuring aerobic energy use at 14 distinct speeds, we showed that locomotion at the lowest speeds involves increased energy use relative to swimming at speeds up to 3 BL s -1 , similar to previous results from skate locomotion ( 58 ). This suggests that analyses of fish energy use as speed increases would benefit from increasing the number of tested speeds, and that efforts to extend speed measurements to both the lowest and highest speeds would enable a greater level of precision in measuring fish locomotor performance.

One key finding in this regard is that fish schools show a minimal absolute energetic cost ( M O 2 ) at a mean group swimming speed of ∼1.0 BL s -1 that is higher than the minimum swimming speeds tested. Thus, fish schools can swim at a speed with the least amount of energy use and potentially extend the distance travelled with the same amount of metabolic substrate onboard. This swimming speed of minimum cost closely matches the migratory speeds of carangiform and subcarangiform migratory fishes (often as schools) which are in the range of 0.5 to 1.5 BL s -1 as recorded by tags on migrating fish ( 59 , 60 ) (Fig. S1; Tables S1, S3). Although D. aequipinnatus is not a migratory species, our results suggest that the migratory speed of fish schools likely occurs at the speed showing the minimum aerobic metabolic rate for long-distance locomotion. Indeed, the Strouhal number when fish schools swim at U opt , showing a minimum aerobic energy cost of locomotion, is 0.3, a hallmark of efficient swimming for many fish species ( 61 ). Fish of different size swimming at 1 BL s -1 will necessarily move at different Reynolds numbers, and hence the scaling of body size to swimming speed needs to be considered in future analyses of other species that differ in size.

The exact energetic mechanisms underpinning the chosen migratory speeds of fish schools would benefit from more in-depth studies. Since both fish schools and solitary individuals have U opt of ∼1 BL s -1 , the M O 2 –speed curve of a fish school might be the average of individual curves from fish within the school. However, pushing and pulling forces fromhydrodynamic interactions of neighbouring fish could help swimmers settle into stable arrangements, a phenomenon known as the Lighthill conjecture ( 62 ). Further studies will need to explore whether individuals within the school are more frequently located at points recommended by the Lighthill conjecture which could potentially result in minimum metabolic costs of locomotion at U opt in fish schools ( 49 ). Research on fish swimming in large aquaculture circular tanks has discovered that holding a large school of salmonids at a water velocity of ∼1 BL s -1 resulted in healthy growth and good conditions ( 63 ). This phenomenon seems to suggest that the minimum metabolic costs of locomotion yield the largest available aerobic capacity for growth and other functions.

Higher locomotor costs at the lowest speeds are likely caused by higher postural costs for active stability adjustments in the near-still fluid. The direct link between higher energetic costs caused by a higher 3-D body angle and higher turning frequency at the lowest speeds in solitary danio ( Fig. 4D, E ) supports previous results in skates ( 58 ) where energetic costs are high at the very low swimming speeds due to the increased energetic cost involved in maintaining body stability and generating lift in negatively buoyant fishes. We expect more species that move in the fluid to show a higher postural cost in near-still fluid but careful measurements of both body kinematics and position along with high-resolution respirometry and accurate low-speed flow control will be required to demonstrate this phenomenon in a diversity of fish species. Elevated costs should be indicated at the lowest speeds as a direct increase in M O 2 and not just as COT. When the denominator ( i.e. speed) for deriving COT is less than 1, COT tends to skew towards the higher value at the lowest speeds. Mechanistically, fish at low fluid speeds have lower muscle and propeller efficiencies ( 44 ) with reduced fluid-assisted stability. Increased active stability adjustments such as fin movements ( 58 ) and tilting behaviours that enable the fish’s body to act as a hydrofoil ( 64 ) are used for stability control (see supplementary video). Moreover, danio do not show intermittent swimming gaits at low speed as observed in some labriform swimmers ( 65 ) ( 66 ), and hence gait-specific kinematics most likely do not play a role in the elevated locomotion cost that we observe at the lowest speeds. Also, the higher turning frequency at < 0.75 BL s -1 in fish schools can relate to the need for coordination with neighbours. Nevertheless, as fish move more rapidly over a speed range, we note that the total cost of transport of ∼1.2 kJ km -1 kg -1 at 3 BL s -1 is among the lowest recorded for aquatic organisms, and is less than the TCOT (1.5 kJ km -1 kg -1 ) for jellyfish (the most energetically efficient low-speed swimmers) ( 4 ) swimming at < 2 BL s -1 ( 67 ).

In summary, our experiments on giant danio have demonstrated substantial energy conservation resulting from schooling dynamics across a wide range of speeds in fish. Direct measurement of both aerobic and non-aerobic energy use is critical for understanding the rapid collective movement of animals. Fish schooling in the high-drag viscous aquatic medium serves as a model for understanding how group movement by animals can be a more energy-efficient biological collective than movement by isolated individuals. By increasing maximum aerobic performance, fish schools save anaerobic energy and reduce the recovery time after peak swimming performance. Furthermore, by decreasing the proportion of metabolic capacity and recovery time devoted to locomotion, animals can apportion more energy to other fitness-related activities, such as digestion, growth and reproduction. More broadly, comprehending how the collective dynamics of animal movements in the water, land, and air can modify the energy use profiles of individuals provide a better understanding of the ecological and evolutionary implications of group locomotion.

Materials and Methods

Experimental animals.

The experiments were performed on giant danio ( Devario aequipinnatus ) that were acquired from a local commercial supplier near Boston, Massachusetts USA (Table S4). Five schooling groups are randomly distributed and housed separately in five 37.9 l aquaria (n=8 per tank). The five solitary individuals are housed separately in five 9.5 l aquaria (n=1 per tank). The individual housing condition acclimated the single D. aequipinnatus to the solitary environment and helped to reduce any isolation stress that might elevate whole-organism metabolic rate. In fact, the aerobic locomotion cost of solitary individuals showed no statistical difference from (in fact, being numerically lower) that of fish schools at a very low testing speed. The flow speed is similar to some areas of the aerated home aquarium for each individual fish. This suggests that the stress of solitary fish likely does not meaningfully contribute to the higher locomotor costs ( see experimental protocol for more details on mitigating the stress). The condition factor showed no difference between solitary fish and fish schools (0.81 vs. 0.99; t 8 = 2.14, p = 0.065). All aquaria have self-contained thermal control (28 °C), an aeration system (>95 % air saturation, % sat.) and a filtration system. Water changes (up to 50% exchange ratio) were carried out weekly. Fish were fed ad libitum daily (TetraMin, Germany). Animal holding and experimental procedures were approved by the Harvard Animal Care IACUC Committee (protocol number 20-03-3).

Integrated Biomechanics & Bioenergetic Assessment System (IBAS)

The core of our Integrated Biomechanics & Bioenergetic Assessment System (IBAS) (Fig. S4) is a 9.35-l (respirometry volume plus tubing) customized Loligo® swim-tunnel respirometer (Tjele, Denmark). The respirometer has an electric motor, and a sealed shaft attached to a propeller located inside the respirometer. Regulating the revolutions per minute (RPM) of the motor controls water velocity of the tunnel. The linear regression equation between RPM and water velocity (V) is established (V = 0.06169 • RPM – 5.128, R 2 = 0.9988, p < 0.0001) by velocity field measured by particle image velocimetry (PIV). Hence, the aerobic costs during locomotion are measured through the regulation of the water velocity.

The swim-tunnel respirometer is oval-shaped. The central hollow space of the respirometer increases the turning radius of the water current. As a result, the water velocity passing the cross-section of the swimming section (80 × 80 × 225 mm) is more homogenous (validated by direct particle image velocimetry (PIV) of flow in the working section of the respirometer following the procedures in ( 68 )). Moreover, a honeycomb flow straightener (80 × 80 × 145 mm) is installed in the upstream section of the swimming section that was specifically designed to create relatively uniform flow across the working section (also confirmed by direct flow imaging). Based on our flow visualization analyses, the effective distance of slow-moving fluid due to boundary layer was <2.5mm at speeds above 2 BL s -1 . The boundary layer played an ever-diminishing role at higher speeds (> 4 BL s -1 ) when energy saving of fish schools becomes more predominant. Danio (∼10 mm wide) cannot effectively hide in the narrow boundary layer created by our flow baffle system. In addition, the convex hull volume of the fish school did not change as speed increased, suggesting that the fish school was not flattening against the wall of the swim tunnel, a typical feature when fish schools are benefiting from wall effects. In nature, fish in the centre of the school effectively swim against a ‘wall’ of surrounding fish where they can benefit from hydrodynamic interactions with neighbours.

To standardize the experimental apparatus and avoid instrument variation, solitary fish and fish schools are measured in the same swim-tunnel respirometer, and the same overall experimental design is used as in previous studies ( 14 ) ( 13 ) ( 15 ). The ratios of respirometer:individual volume ( r RI ) in our experiments was 2200 for individual fish (we used larger solitary D. aequipinnatus to increase the signal-to-noise ratio), and 693 for fish schools. The r RI is essentially a balance between giving enough space for fish to exhibit natural locomotion and reducing wall effects and generating a reliable signal-to-noise ratio to measure the decline of dissolved O 2 (DO) in water ( 69 ). The increase in the signal-to-noise ratio can come from the better technology of the O 2 probe (we used an optical O 2 probe which has a higher measurement sensitivity than the Winkler method or Electrode & Galvanic probes used in the previous studies, see ( 70 ) for review) and modifications to the swim-tunnel respirometer. We used a water loop is installed 95 cm downstream of the propeller with water returned to the respirometer 240 cm upstream of the swimming section. Flow in the water loop moves (produced by an in-line circulation pump, Universal 600, EHEIM GmbH & Co KG, Deizisau, Germany) in the same direction as the water flow in the swimming tunnel. A high-resolution fibre optic O 2 probe (Robust oxygen probe OXROB2, PyroScience GmbH, Aachen, Germany) is sealed in the loop in the downstream of the circulation pump where there is better mixing to continuously measure the DO level in the water (recording frequency ∼1 Hz, response time < 15s). The designated water sampling loop together with the water mixing by the propeller and water pumps effectively reduces the noise-to-signal ratio. As a result, the respirometry system reaches a stable signal-to-noise ratio once the sampling window is longer than 1.67 mins ( see Fig. S5), well within the duration of the velocity step to obtain a stable signal-to-noise ratio for calculating M O 2 ( 32 ).

The oxygen probe was calibrated to anoxic (0 % sat., a solution created by super-saturated sodium sulphite and bubbling nitrogen gas) and fully aerated water (100 % sat.). After fish removal, the background M O 2 in the swim-tunnel respirometer was measured for a 20-min sealed period before and after each trial to calculate the average background M O 2 (< 6% of fish M O 2 ), which was used to correct for the M O 2 of fish ( 69 ). The pre-filtered water (laboratory grade filtration system) is constantly disinfected by UV light (JUP-01, SunSun, China) located in an external water reservoir to suppress the growth of microbial elements. Water changes of 60% total volume occurred every other day and a complete disinfection by sodium hypochlorite is conducted weekly (Performance bleach, Clorox & 1000 ppm).

To simultaneously measure schooling dynamics and swimming kinematics, the customized oval-shaped swim-tunnel respirometer is located on a platform with an open window beneath the swimming section. The platform is elevated 243 mm above the base to allow a front surface mirror to be installed at a 45° angle. This mirror allows a high-speed camera (FASTCAM Mini AX50 type 170K-M-16GB, Photron Inc., United States, lens: Nikon 50mm F1.2, Japan) to record the ventral view. The second camera (FASTCAM Mini AX50 type 170K-M-16GB, Photron Inc., United States, lens: Nikon 50mm F1.2, Japan) is positioned 515 mm to the side of the swimming section to record a lateral view. Synchronized lateral and ventral video recordings were made at 125 fps, and each frame was 1024 by 1024 pixels. To avoid light refraction passing through the water and distorting the video recordings, the swim-tunnel respirometry is not submerged in a water bath. Temperature regulation of the respirometer is achieved by regulating room temperature, installing thermal insulation layers on the respirometer and replenishing the water inside the respirometer from a thermally regulated (28 °C, heater: ETH 300, Hydor, United States & chiller: AL-160, Baoshishan, China) water reservoir (insulated 37.9-l aquarium) located externally.

The aerated (100% sat., air pump: whisper AP 300, Tetra, China) reservoir water is flushed (pump: Universal 2400, EHEIM GmbH & Co KG, Deizisau, Germany) to the respirometer through an in-line computer-controlled motorized ball valve (U.S. Solid) installed at the in-flow tube. The other in-line one-way valve is installed at the out-flow tube. The out-flow tube is also equipped with a one-way valve. The valve is shut during the measurement period, a precautionary practice to eliminate the exchange of water between the respirometer and the external reservoir when the water moves at a high velocity inside the respirometer. This flushing was manually controlled to maintain DO above 80 % sat. Every time the respirometer was closed to measure M O 2 , the water temperature fluctuates no more than 0.2 °C. The water temperature inside the respirometer is measured by a needle temperature probe (Shielded dipping probe, PyroScience GmbH, Aachen, Germany) sealed through a tight rubber port of the respirometer.

To allow fish to reach the undisturbed quiescent state during the trial (another stress mitigation practice), the entire IBAS is covered by laser blackout sheets (Nylon Fabric with Polyurethane Coating; Thorlabs Inc, New Jersey, United States). The room lights are shut off and foot traffic around IBAS is restrained to the absolute minimum. Fish are orientated by dual small anterior spots of white light (lowest light intensity, Model 1177, Cambridge Instruments Inc, New York, United States) for orientation (one to the top and the other to the side) of the swimming section. The test section is illuminated by infrared light arrays to allow high-speed video recording with minimal disturbance of fish behaviour.

Experimental Protocol

The energy use of vertebrates at lower speeds is primarily aerobic, while for high-speed movement anaerobic metabolic pathways are activated to supply the additional (largely shorter-term) energy needs ( 9 ). While whole-animal aerobic metabolism is measured by oxygen (O 2 ) uptake rate ( M O 2 ), the non-aerobic O 2 cost (mostly through high-energy phosphate and anaerobic glycolysis) is measured as excess post-exercise O 2 consumption (EPOC)( 10 ). Both metabolic energy sources contribute to the total energy expenditure (TEE) required for movement. For a given workload, the higher the maximum aerobic performance of the animal, the less the need for the anaerobic energy contribution ( 9 , 35 ). The experimental protocol captures both metabolic energy contributions by measuring M O 2 during locomotion and EPOC afterwards.

We studied five replicate schools and five replicate individuals. Swimming performance test trials were conducted with Devario aequipinnatus fasted for 24 hours, a sufficient period for a small size species at 28 °C ( i.e. high resting M O 2 ) to reach an absorptive state. In fact, we observed no specific dynamic action, the amount of O 2 consumed for digestion during the first diurnal cycle (Fig. S3). Prior to the swimming performance test, testing fish were gently weighed and placed in the swim-tunnel respirometer. The fish swam at 35% U crit for 30 mins to help oxidize the inevitable but minor lactate accumulation during the prior handling and help fish become accustomed to the flow conditions in the swim-tunnel respirometer ( 43 ). After this time, the fish to be tested were habituated (>20 hours) to the respirometer environment under quiescent and undisturbed conditions. The long habituation period also helps to reduce the stress and further reduce the likelihood that the solitary fish might be more stressed than the fish schools and showed an elevated M O 2 ( 13 ). During this time, we used an automatic system to measure the resting M O 2 for at least 19 hours. Relays (Cleware GmbH, Schleswig, Germany) and software (AquaResp v.3, Denmark) were used to control the intermittent flushing of the respirometer with fresh water throughout the trial to ensure O 2 saturation of the respirometer water. M O 2 was calculated from the continuously recorded DO level (at 1 Hz) inside the respirometer chamber. The intermittent flow of water into the respirometer occurred over 930 s cycles with 30 s where water was flushed into the respirometer and 900 s where the pumps were off and the respirometer was a closed system. The first 240 s after each time the flushing pump was turned off were not used to measure M O 2 to allow O 2 levels inside the respirometer to stabilize. The remaining 660 s when the pumps were off during the cycle were used to measure M O 2 . The in-line circulation pump for water in the O 2 measurement loop stayed on throughout the trial.

We characterize the aerobic costs for the swimming performance of fish using an established incremental step-wise critical swimming speed ( U crit ) test ( 10 ). The first preliminary trial determined the U crit of this population of Devario aequipinnatus as 8 BL s -1 . Characterizing the swimming performance curve required a second preliminary trial to strategically select 9 water velocities (0.3, 0.5, 0.8, 1.0, 1.3, 1.5, 1.8, 2.3, 2.8 BL s -1 ) to bracket the hypothesized upward concave shaped aerobic metabolism-speed curve at the lower speed (< 40% U crit ). Additional five water velocities (3.8, 4.9, 5.9, 6.9, 8.0 BL s -1 ) are used to characterize the exponentially increasing curve to the maximum and sustained swimming speed, U crit ( see Fig. S6). Altogether, 14 points provide a reliable resolution to characterize the swimming performance curve. The cross-section of the danio school is ∼10% of the cross-sectional area of the swim tunnel, hence the blocking effect is negligible ( 71 ). At each water velocity, fish swam for 10 mins ( 58 ) to reach a steady state in M O 2 at low speed ( see Fig. S7). Above 40% U crit , M O 2 can become more variable ( 32 ). Hence, in this protocol, we focus on measuring the sustained aerobic energy expenditure by calculating the average M O 2 for each 10-min velocity step using Eqn 1 . At the 5 th min of each velocity step, both ventral and lateral-view cameras are triggered simultaneously to record 10-sec footage at 125 frames per second, at 1/1000 shutter speed and 1024 ×1024 pixel resolution. Thus, both data streams of M O 2 and high-speed videos are recorded simultaneously. The U crit test is terminated when 12.5% of fish in the school or a solitary individual touches the back grid of the swimming section for more than 20 secs ( 43 ). The U crit test lasted ∼140 mins and estimates the aerobic portion of TEE over the entire range of swimming performance.

To measure the contribution of non-aerobic O 2 cost, where the majority of the cost is related to substrate-level phosphorylation, and to calculate the TEE for swimming over the entire speed range, we measured EPOC after the U crit test for the ensuing 19 hours, recorded by an automatic system. Most previous measurements of EPOC have used a duration of ∼5 hours ( see review ( 42 )), but our extended measurement period ensured that longer duration recovery O 2 consumption during EPOC was measured completely as fish were exercised to U crit ( see summary table in 16 ). The intermittent flow of water into the respirometer occurred over 30 s to replenish the DO level to ∼95% sat. For the following 900 s the flushing pump remained closed, and the respirometer becomes a closed system, with the first 240 s to allow O 2 saturation inside the respirometer to stabilize. The remaining 660 s when the flushing pump was off during the cycle were used to measure M O 2 ( see Eqn. 1 ). The cycle is automated by computer software (AquaResp v.3) and provided 74 measurements of M O 2 to compute EPOC. Upon the completion of the three-day protocol, the school or individual fish are returned to the home aquarium for recovery. The fish condition was closely monitored during the first 48 hours after the experiment, during which no mortality was observed.

Bioenergetic measurement and modeling

To estimate the steady-rate whole-animal aerobic metabolic rate, M O 2 values were calculated from the sequential interval regression algorithm ( Eqn. 1 ) using the DO points continuously sampled (∼1 Hz) from the respirometer.

Where d DO / d t is the change in O 2 saturation with time, V r is the respirometer volume, V f is the fish volume (1 g body mass = 1 ml water), S o is the water solubility of O 2 (calculated by AquaResp v.3 software) at the experimental temperature, salinity and atmospheric pressure, t is a time constant of 3600 s h -1 , M f is fish mass, and a is the sampling window duration, i is the next PO 2 sample after the preceding sampling window.

To account for allometric scaling, the M O 2 values of solitary fish were transformed to match the size of the individual fish in the school ( see Table S4) using an allometric scaling exponent (b = 0.7546). The calculation of the scaling relationship [Log 10 ( M O 2 ) = b•Log10( M ) + Log10( a ), where M is the body mass & a is a constant] was performed by least squares linear regression analysis (y = 0.7546 • x + 0.2046; R 2 = 0.6727, p < 0.0001) on the 180 data points of metabolic rate and body mass from a closely related species (the best available dataset to our knowledge). (The mass scaling for tail beat frequency was not conducted because of the lack of data for D. aequipinnatus and its related species. Using the scaling exponent of distant species for mass scaling of tail beat frequency will introduce errors of unknown magnitude.) The allometrically scaled M O 2 values were used to derive other energetic metrics (listed below & such as aerobic scope) for the solitary fish. The energetic metrics of fish schools are calculated from the mass-specific M O 2 .

The resting oxygen uptake ( M O 2rest ), the minimum resting metabolic demands of a group of fish or a solitary individual, is calculated from a quantile 20% algorithm ( 74 ) using the M O 2 estimated between the 10 th –18 th hour and beyond the 32 nd –51 st hour of the trial. These are the periods of quiescent state when fish completed the EPOC from handling and swimming test.

The M O 2 for the aggregation ( Fig. 2B , M O 2aggregate ) is calculated as the average value using the M O 2 estimated between the 10 th –18 th hour without the effects of any tests.

Minimum oxygen uptake ( M O 2min ) is the lowest M O 2 value recorded in the entire trial, which always occurred at the optimal speed when a school of fish collectively reached the lowest M O 2 value.

Active oxygen uptake ( M O 2active ) is the highest average M O 2 when fish are actively swimming ( 10 ).

The aerobic scope (the metric for aerobic capacity) is the numerical difference between M O 2active and M O 2min ( i.e. M O 2active – M O 2min ) ( 10 ).

The percentage of aerobic scope (% aerobic scope) is calculated by normalizing the M O 2 value at a water velocity as a % aerobic scope: [( M O 2 – M O 2min ) / ( M O 2active – M O 2min )]∗100%. The apportioning of aerobic scope to swimming performance is computed across the entire range of swim speeds.

The excess post-exercise oxygen consumption (EPOC) is an integral area of M O 2 measured during post-exercise recovery, from the end of U crit until reached M O 2rest plus 10% or within 24 hours post-exercise, whichever endpoint occurred first ( 72 ). This approach reduces the likelihood of overestimating EPOC due to any potential spontaneous activities ( 72 ). To account for the allometric scaling effect, we used the total amount of O 2 consumed (mg O 2 ) by the standardized body mass of fish (1.66 g) for fish schools and solitary fish.

We model EPOC ( i.e. non-aerobic O 2 cost) and M O 2 during U crit to estimate a total O 2 cost over the duration of the swimming performance test. Our conceptual approach was pioneered by Brett ( 10 ) in fish ( 75 ) and is also used in sports science ( 9 ). Mathematical modeling of EPOC and M O 2 during swimming was applied to study the effects of temperature on the total cost of swimming for migratory salmon ( 43 ). We improved the mathematical modeling by applying the following physiological and physics criteria. The first criterion is that significant accumulation of glycolytic end-product occurred when fish swimming above 50% U crit ( 31 ) which corresponds to > 40% M O 2max (or ∼ 50% aerobic scope) ( 9 ). This is also when fish start unsteady-state burst-&-glide swimming gait ( 31 ). The second criterion is that the integral area for the non-aerobic O 2 cost during swimming can only differ by ≤ 0.09% when compared to EPOC. The non-aerobic O 2 cost during swimming is the area bounded by modeled M O 2 and measured M O 2 as a function of time when fish swim > 50% U crit ( see Fig. 2A & Table S2). The third criterion is that total energy expenditure is expected to increase exponentially with swimming speed (Fig. S8). Specifically, these curves were fitted by power series or polynomial models, the same models that describe the relationship between water velocity and total power and energy cost of transport (Fig. S8). Following these criteria, the non-aerobic O 2 cost at each swimming speed is computed by a percentage (%) modifier based on the aerobic O 2 cost (Table S2). The exponential curve of total O 2 cost as swimming speed of each fish school or solitary individual was derived by an iterative process until the difference between non-aerobic O 2 cost and EPOC met the 2 nd criterion. The sum of non-aerobic O 2 cost and aerobic cost gives the total O 2 cost.

Total energy expenditure (TEE) is calculated by converting total O 2 cost to kJ × kg -1 using an oxy-calorific equivalent of 3.25 cal per 1 mg O 2 ( 76 ).

Cost of transport (COT), in kJ × km -1 × kg -1 is calculated by dividing TEE by speed (in km × h -1 ) ( 58 ).

Three-dimensional kinematic data extraction from high-speed videography

We used two synchronized 10-sec high-speed videos (lateral and ventral views, at each speed) for kinematic analyses. We calibrated the field of view of the high-speed cameras using a direct linear transformation for three-dimensional (3-D) kinematic reconstruction (DLTdv8)( 77 ) by applying a stereo calibration to the swimming section of the respirometer ( see Fig. S9). We digitized the anatomical landmarks of fish ( see Fig. S10) to obtain the X, Y, and Z coordinates for each marker at the 1 st sec, 5 th sec and 10 th sec for videos recorded at each speed. These coordinates are used to calculate the following kinematic parameters. All the calculations are validated on the known length and angle of test objects inserted into the tank working section.

The 3-D distance between the tip of the nose of each fish in the school per frame is calculated in vector Eqn. 2

Where spatial coordinates of two neighbouring fish are ( X a , Y a , Z a ) and ( X b , Y b , Z b ) respectively.

The 3-D angle of each fish in the school to the frontal plane per frame is calculated by Eqn. 3

where the spatial coordinates of the caudal peduncle, nose of fish and right-angle crosshair between the peduncle and the nose are ( X 1 , Y 1 , Z 1 ), ( X 3 , Y 3 , Z 3 ) and ( X 2 , Y 2 , Z 2 ) respectively ( see Fig. S9).

The fish’s angle to water flow per frame is calculated using the arctangent function in Excel (Microsoft, United States) using the spatial coordinates of the caudal peduncle and nose of the fish.

The school length (in X -axis) is calculated as a 3-D distance ( Eqn. 2 ) between the nose of the first fish of the school to the caudal peduncle of the last fish in the school.

The values (3-D distance of the individuals, 3-D body angle, 3-D school length) calculated above were averaged among the three frames at each speed as a representative kinematic feature of the schools (or solitary fish) for each speed. The standard deviation of the angle to water flow in fish schools is calculated from the angular values of the individuals.

Tail beat frequency ( f TB ), in Hz, is calculated as the number of tail beats observed within one second. We sampled 10 tail beats from fish exhibiting steady-state swimming to obtain an average and representative f TB for each solitary fish at each speed. We sampled 10 tail beats of steady-state swimming for three individuals in the center of the fish school for an average and representative f TB for each fish school.

Turning frequency is the total number of events when fish made a 90° turn in the 10-sec video. To make the total number of turns per individual and fish schools comparable (multiple individuals inevitably have a larger absolute number of turns), we derived the average turning frequency per schooling fish by dividing the total number of turns of fish school by the number of fish in the school.

Total energy expenditure per tail beat (kJ kg -1 beat -1 ) is calculated by TEE • ( f TB • 60) -1 to estimate the total metabolic energy spent to achieve one stride length.

Additional calculations of fluid dynamic metrics were:

Strouhal number = ( f TB • tail beat amplitude) • U -1 ( 8 ). Tail beat amplitude is measured as the average distance of five peak-to-peak oscillation amplitudes of the tip of the fish’s tail. The measurement is conducted on the calibrated high-speed video in video analysis software (Phontron FASTCAM Viewer 4, Photron USA, Inc.).

Reynolds number = (water density • U • fish fork length) • water dynamic viscosity -1 ( 8 ). Water density and dynamic viscosity are given at 28 °C.

Statistical analyses

Measurement points are presented as mean ± s.e.m. For the metrics that failed normality tests, logarithm transformations were applied to meet the assumptions of normality of residuals, homoscedasticity of the residuals, and no trend in the explanatory variables. Since fish schools exhibit the features of a coherent functional unit, the statistical model treated one school as one sample size, and the experiment measured five replicates of fish schools. The majority of statistical comparisons used the General Linear Model (Fixed factors: solitary fish vs. fish schools & swimming speeds, the label of solitary fish & fish schools as a random effect) with Holm–Šídák post-hoc tests. The few comparisons that used other statistical models are listed below. The comparison of M O 2min at U opt , M O 2 of aggregating behaviours exhibited at the lowest speed ( M O 2aggregate ) and resting M O 2 ( M O 2rest ) in fish schools used one-way ANOVA with Holm– Šídák post-hoc tests. The comparison of EPOC between solitary fish and fish schools used a two-tailed Student’s t-test. The comparison of the duration of EPOC between solitary fish and fish schools used a two-tailed Student’s t-test. The overall difference in percentage aerobic scope (% aerobic scope) between fish schools and solitary individuals is compared by the Wilcoxon singed-rank test over the entire range of 14 swimming speeds. The statistical analyses were conducted in SPSS v.28 (SPSS Inc. Chicago, IL, USA). The best-fitting regression analyses were conducted using Prism v.9.4.1 (GraphPad Software, San Diego, CA, USA). 95% C.I. values were presented for all regression models as shaded areas around the regression or data points. Statistical significance is denoted by *, **, ***, **** for p -values of ≤ 0.05, ≤ 0.01, ≤ 0.001, ≤ 0.0001 respectively.

Acknowledgements

Many thanks to members of Lauder Laboratory for numerous discussions about fish schooling behaviour, for comments on the manuscript, and to Cory Hahn for fish care.

Funding provided by the National Science Foundation grant 1830881 (GVL), the Office of Naval Research grants N00014-21-1-2661 (GVL), N00014-16-1-2515 (GVL), 00014-22-1-2616 (GVL), and a Postdoctoral Fellowship of Natural Sciences and Engineering Research Council of Canada PDF - 557785 – 2021 (YZ).

Author contributions

Y.Z. and G.L. conceptualized the study. Y.Z. performed experiments and data analyses and wrote the manuscript. Y.Z. and G.L. provided manuscript edits and comments and approved the final version.

Competing interests

Authors declare that they have no competing interests.

Data and materials availability

All data generated or analysed during this study are included in this published article (and its supplementary information files).

Supplementary Materials

Introductory text

Figs. S1 to S9

Tables S1 to S4

• Biewener A.
• Weimerskirch H.
• Clerquin Y.
• Alexandre P.
• Jiraskova S.
• van Druenen T.
• Toparlar Y.
• Andrianne T.
• Shinichiro I.
• Lauder G. V.
• Laforgia J.
• Withers R. T.
• Brett J. R.
• Beamish F. W. H.
• Portugal S. J.
• Hubel T. Y.
• Wilson A. M.
• Usherwood J. R.
• Parker Jr F. R.
• Abrahams M. V.
• Colgan P. W.
• Burgerhout E.
• Tudorache C.
• Brittijn S. A.
• Palstra A. P.
• Dirks R. P.
• van den Thillart G. E. E. J. M.
• Coughlin D. J.
• Partridge B. L.
• Pitcher T. J.
• Roskilly K.
• Graving J. M.
• Bak-Coleman J.
• Couzin I. D.
• (潘宇) Pan Y.
• (董海波) Dong H.
• Schreiber C. M.
• Moored K. W.
• Bart-Smith H.
• Jones K. D.
• Dohring C. M.
• Platzer M. F.
• Stewart J. D.
• Butler R. H.
• Munns R. J.
• Magurran A. E.
• Winfield I. J.
• Farrell A. P.
• Holder P. E.
• Lawrence M. J.
• Clark T. D.
• Suski C. D.
• Weber J.-M.
• Danylchuk A. J.
• Cooke S. J.
• Misund O. A.
• Peake S. J.
• Gilbert M. J. H.
• Clutterham S. M.
• Fegely J. F.
• Xanthopoulos C. J.
• Gaesser G. A.
• Brooks G. A.
• Jayne B. C.
• Loughna P. T.
• Goldspink G.
• Milligan C. L.
• Brown H. F.
• Difrancesco D.
• Noble S. J.
• Ruggiero D. A.
• Fernandez-Pardal J.
• Saavedra J. M.
• Claireaux G.
• Jørgensen S. M.
• Hinch S. G.
• Healey M. C.
• Costanzo A.
• Hemelrijk C. K.
• Rountree R. A.
• Camazine S.
• Deneubourg J.-L.
• Franks N. R.
• Theraula G.
• Bonabeau E.
• Koumoutsakos P.
• Hover F. S.
• Triantafyllou M. S.
• Zuyev G. V.
• Belyayev V. V.
• Killen S. S.
• Lindström J.
• McKenzie D. J.
• Steffensen J. F.
• Domenici P.
• Bushnell P. G.
• Johansen K.
• Sepulveda C.
• Dickson K. A.
• Di Santo V.
• Kenaley C. P.
• Block B. A.
• Carey F. G.
• Teo S. L. H.
• Boustany A.
• Stokesbury M. J. W.
• Farwell C. J.
• Williams T. D.
• Triantafyllou G. S.
• Cochran-Carney J.
• Mivehchi A.
• Quinn D. B.
• Timmerhaus G.
• Lazado C. C.
• Cabillon N. A. R.
• Reiten B. K. M.
• Johansen L.-H.
• Wardle C. S.
• Cathcart K.
• Gellman E. D.
• Tandler T. R.
• Ellerby D. J.
• Dabiri J. O.
• Matthews D. G.
• Marquez C. A.
• Zimmerman J. F.
• Ardoña H. A. M.
• Kleber A. G.
• Parker K. K.
• Svendsen M. B. S.
• Kline R. J.
• Parkyn D. C.
• Murie D. J.
• Wootton H. F.
• Morrongiello J. R.
• Audzijonyte A.
• Svendsen J. C.
• Jordan A. D.
• Aarestrup K.
• Brafield A. E.
• Solomon D. J.
• Hedrick T. L.
• Berlinger F.
• Sheshmani A.
• Haj-Hariri H.
• Daghooghi M.
• Borazjani I.

Article and author information

Author information, yangfan zhang.

• Corresponding author. Email: [email protected]

Version history

• Sent for peer review : July 10, 2023
• Preprint posted : July 26, 2023
• Reviewed Preprint version 1 : October 27, 2023
• Reviewed Preprint version 2 : January 19, 2024
• Version of Record published : February 20, 2024

© 2023, Yangfan Zhang & George V. Lauder

This article is distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use and redistribution provided that the original author and source are credited.

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Be the first to read new articles from eLife

share this!

July 24, 2024

This article has been reviewed according to Science X's editorial process and policies . Editors have highlighted the following attributes while ensuring the content's credibility:

fact-checked

peer-reviewed publication

trusted source

Simulations reveal mechanism behind dynamic fish schooling

by Tohoku University

Schools of fish are able to pull off complex, coordinated maneuvers without ever colliding with one another. They move in unison but don't follow a leader.

To try and understand the intricacies of this collective animal behavior, researchers from Tohoku University have developed a model that simulates the group motion of fish based on visual cues . Incorporating the tendency that fish focus on nearby, quickly moving fish, the model uncovers the mechanism behind dynamic fish schooling.

The research is published in the journal PNAS Nexus .

"Fish have a wide angle of view and can detect many other fish in a school," explains Susumu Ito, "However, a recent experimental finding shows that each fish selects a single fish from a few targets and traces its motion. It is a spectacular example of selective decision making."

Attending to every single fish in the school would require an immense amount of information to be processed. Similar to how we can focus on just the words we are reading on a page of text, fish can focus on the most salient target that determines their next move. While the fish swimming directly in front may seem like the best option, it's actually fish that are slightly to the sides that tend to draw attention.

Ito and his team constructed a model that takes visual attention into account to elucidate the role of selective visual interaction in a large group of fish. It incorporates the feature of retinal ganglion cells that fire preferentially for targets that are closer and moving more quickly. Visual attention is then guided in the direction of the strongest signal. Only fish that fall within that spotlight of visual attention can influence the motion of the individual fish.

Using numerical simulations , the team of researchers found that when a fish is following three targets swimming in a row, it tends to be attracted to the leftward or rightward target, as they have a larger apparent size. A slender fish facing straight ahead viewed directly from behind will look much smaller than a fish exposing its longer profile. These results replicate the selective tracking motion seen in previous experiments.

Furthermore, the model reproduced various collective patterns of fish schools: a rotating vortex, straight, random and turning. In the turning pattern, fish repeatedly alternate between straight and rotating motion, so that the school dynamically reshapes itself.

"The selective tracking behavior is observed also in locust and fly," adds Ito. "We expect to extend the model to the group motion of various organisms in the future. A three-dimensional version of the model may also be able to explain the formation of a huge fish school known as the bait-ball."

Journal information: PNAS Nexus

Provided by Tohoku University

Explore further

Feedback to editors

Transforming marine waste and carbonated water into hydrogels via CO₂ release behavior

7 minutes ago

Cracking the code of performance degradation in solid oxide cells at the atomic level

15 minutes ago

Addressing key challenges of photoresponsive, visible light 3D printing

16 minutes ago

Volcanic caves research advances the search for life on Mars

25 minutes ago

In collective animal movements, speed matters—scientists use 'force map' to investigate fish schools

Electron imaging reveals the vibrant colors of the outermost electron layer

26 minutes ago

Viscosity of materials key to cell differentiation: New insights could improve design of biomaterials

38 minutes ago

Neat, precise and brighter than ever: New technologies improve temporal coherence of XFEL pulses

46 minutes ago

Simulations reveal black holes inherit magnetic fields from parent stars

First successful synthesis of elusive antibiotic compounds since their discovery 50 years ago

Relevant physicsforums posts, universal wing- and fin-beat frequency scaling.

9 hours ago

Exercises for loosening joints

14 hours ago

Chornobyl Dark coloured tree frogs

Nov 16, 2024

Why is it easier on the heart in a horizontal position?

Nov 15, 2024

Marburg outbreak in Rwanda

Nov 12, 2024

How would an infection prevent anesthetic from numbing my toe?

Nov 11, 2024

More from Biology and Medical

Related Stories

Model shows how fish synchronize tail fins to save energy

Nov 6, 2023

Schooling fish expend less energy in turbulent water compared to solitary swimmers, study finds

Jun 6, 2024

Hydrodynamic model of fish orientation in a channel flow

Jun 10, 2022

Study finds schools of fish can make less noise than a solitary swimmer

Apr 9, 2024

Finding the beat of collective animal motion: Scientists show reciprocity is key to driving coordinated movements

May 23, 2024

Fish in schools can take it easy

May 9, 2018

Recommended for you

Researchers characterize regulating mechanism of orderly zygotic genome activation in early embryos

A single cell's siesta: How non-moving single-celled organisms manage to avoid bright light

Mathematical model can help understand coexistence in nature

2 hours ago

Archaea shape microbiomes by using molecular spring-loaded daggers

3 hours ago

How marine worms regenerate lost body parts: Return of cells to stem cell-like state could be key

Methylmercury: How microbes create the most toxic form of mercury

Let us know if there is a problem with our content.

Use this form if you have come across a typo, inaccuracy or would like to send an edit request for the content on this page. For general inquiries, please use our contact form . For general feedback, use the public comments section below (please adhere to guidelines ).

Please select the most appropriate category to facilitate processing of your request

Thank you for taking time to provide your feedback to the editors.

Your feedback is important to us. However, we do not guarantee individual replies due to the high volume of messages.

E-mail the story

Your email address is used only to let the recipient know who sent the email. Neither your address nor the recipient's address will be used for any other purpose. The information you enter will appear in your e-mail message and is not retained by Phys.org in any form.

Newsletter sign up

Get weekly and/or daily updates delivered to your inbox. You can unsubscribe at any time and we'll never share your details to third parties.

More information Privacy policy

Donate and enjoy an ad-free experience

We keep our content available to everyone. Consider supporting Science X's mission by getting a premium account.

E-mail newsletter