Ecology

Xia, J. et al. Terrestrial carbon cycle affected by non-uniform climate warming. Nat. Geosci. 7, 173–180 (2014).CAS 
Article 

Google Scholar 
Tang, R. et al. Increasing terrestrial ecosystem carbon release in response to autumn cooling and warming. Nat. Clim. Change 12, 380–385 (2022).CAS 
Article 

Google Scholar 
Hugelius, G. et al. Large stocks of peatland carbon and nitrogen are vulnerable to permafrost thaw. Proc. Natl Acad. Sci. USA 117, 20438–20446 (2020).CAS 
Article 

Google Scholar 
Gallego-Sala, A. V. et al. Latitudinal limits to the predicted increase of the peatland carbon sink with warming. Nat. Clim. Change 8, 907–913 (2018).CAS 
Article 

Google Scholar 
Treat, C. C. et al. Widespread global peatland establishment and persistence over the last 130,000 y. Proc. Natl Acad. Sci. USA 116, 4822–4827 (2019).CAS 
Article 

Google Scholar 
Frolking, S., Roulet, N. & Fuglestvedt, J. How northern peatlands influence the Earth’s radiative budget: sustained methane emission versus sustained carbon sequestration. J. Geophys. Res. Biogeosci. 111, G01008 (2006).
Google Scholar 
Loisel, J. et al. Expert assessment of future vulnerability of the global peatland carbon sink. Nat. Clim. Change 11, 70–77 (2021).Article 

Google Scholar 
Helbig, M. et al. Direct and indirect climate change effects on carbon dioxide fluxes in a thawing boreal forest–wetland landscape. Glob. Change Biol. 23, 3231–3248 (2017).Article 

Google Scholar 
Koebsch, F. et al. Refining the role of phenology in regulating gross ecosystem productivity across European peatlands. Glob. Change Biol. 26, 876–887 (2020).Article 

Google Scholar 
Huang, Y. et al. Tradeoff of CO2 and CH4 emissions from global peatlands under water-table drawdown. Nat. Clim. Change 11, 618–622 (2021).CAS 
Article 

Google Scholar 
Evans, C. D. et al. Overriding water table control on managed peatland greenhouse gas emissions. Nature 593, 548–552 (2021).CAS 

Google Scholar 
Helfter, C. et al. Drivers of long-term variability in CO2 net ecosystem exchange in a temperate peatland. Biogeosciences 12, 1799–1811 (2015).Article 

Google Scholar 
Järveoja, J., Nilsson, M. B., Gažovič, M., Crill, P. M. & Peichl, M. Partitioning of the net CO2 exchange using an automated chamber system reveals plant phenology as key control of production and respiration fluxes in a boreal peatland. Glob. Change Biol. 24, 3436–3451 (2018).Article 

Google Scholar 
Mäkiranta, P. et al. Responses of phenology and biomass production of boreal fens to climate warming under different water-table level regimes. Glob. Change Biol. 24, 944–956 (2018).Article 

Google Scholar 
Li, Q. et al. Abiotic and biotic drivers of microbial respiration in peat and its sensitivity to temperature change. Soil Biol. Biochem. 153, 108077 (2021).CAS 
Article 

Google Scholar 
Moore, T. R. et al. Spring photosynthesis in a cool temperate bog. Glob. Change Biol. 12, 2323–2335 (2006).Article 

Google Scholar 
Korrensalo, A. et al. Species-specific temporal variation in photosynthesis as a moderator of peatland carbon sequestration. Biogeosciences 14, 257–269 (2017).CAS 
Article 

Google Scholar 
Weltzin, J. F. et al. Response of bog and fen plant communities to warming and water-table manipulations. Ecology 81, 3464–3478 (2000).Article 

Google Scholar 
Dimitrov, D. D., Grant, R. F., Lafleur, P. M. & Humphreys, E. R. Modeling the effects of hydrology on gross primary productivity and net ecosystem productivity at Mer Bleue bog. J. Geophys. Res. Biogeosci. 116, G04010 (2011).Article 
CAS 

Google Scholar 
Bubier, J., Crill, P., Mosedale, A., Frolking, S. & Linder, E. Peatland responses to varying interannual moisture conditions as measured by automatic CO2 chambers. Glob. Biogeochem. Cycles 17, 1066 (2003).Article 
CAS 

Google Scholar 
Moore, T. R. & Knowles, R. The influence of water table levels on methane and carbon dioxide emissions from peatland soils. Can. J. Soil Sci. 69, 33–38 (1989).CAS 
Article 

Google Scholar 
Nichols, D. S. Temperature of upland and peatland soils in a north central Minnesota forest. Can. J. Soil Sci. 78, 493–509 (1998).Article 

Google Scholar 
Bellisario, L. M., Moore, T. R. & Bubier, J. L. Net ecosystem CO2 exchange in a boreal peatland, northern Manitoba. Écoscience 5, 534–541 (1998).Article 

Google Scholar 
Yu, Z. et al. Peatlands and their role in the global carbon cycle. Eos 92, 97–98 (2011).Article 

Google Scholar 
Hanson, P. J. et al. Rapid net carbon loss from a whole-ecosystem warmed peatland. AGU Adv. 1, e2020AV000163 (2020).Article 

Google Scholar 
Vincent, L. A. et al. Observed trends in Canada’s climate and influence of low-frequency variability modes. J. Clim. 28, 4545–4560 (2015).Article 

Google Scholar 
Templer, P. H. et al. Climate Change Across Seasons Experiment (CCASE): a new method for simulating future climate in seasonally snow-covered ecosystems. PLoS ONE 12, e0171928 (2017).Article 
CAS 

Google Scholar 
Peichl, M. et al. A 12-year record reveals pre-growing season temperature and water table level threshold effects on the net carbon dioxide exchange in a boreal fen. Environ. Res. Lett. 9, 055006 (2014).Article 

Google Scholar 
Helbig, M., Humphreys, E. R. & Todd, A. Contrasting temperature sensitivity of CO2 exchange in peatlands of the Hudson Bay Lowlands, Canada. J. Geophys. Res. Biogeosci. 124, 2126–2143 (2019).CAS 
Article 

Google Scholar 
Griffis, T. J., Rouse, W. R. & Waddington, J. M. Interannual variability of net ecosystem CO2 exchange at a subarctic fen. Glob. Biogeochem. Cycles 14, 1109–1121 (2000).CAS 
Article 

Google Scholar 
Bubier, J. L., Crill, P. M., Moore, T. R., Savage, K. & Varner, R. K. Seasonal patterns and controls on net ecosystem CO2 exchange in a boreal peatland complex. Glob. Biogeochem. Cycles 12, 703–714 (1998).CAS 
Article 

Google Scholar 
Park, S.-B. et al. Temperature control of spring CO2 fluxes at a coniferous forest and a peat bog in Central Siberia. Atmosphere 12, 984 (2021).CAS 
Article 

Google Scholar 
Adkinson, A. C., Syed, K. H. & Flanagan, L. B. Contrasting responses of growing season ecosystem CO2 exchange to variation in temperature and water table depth in two peatlands in northern Alberta, Canada. J. Geophys. Res. Biogeosci. 116, G01004 (2011).Article 
CAS 

Google Scholar 
Heiskanen, L. et al. Carbon dioxide and methane exchange of a patterned subarctic fen during two contrasting growing seasons. Biogeosciences 18, 873–896 (2021).CAS 
Article 

Google Scholar 
Lafleur, P. M., Roulet, N. T., Bubier, J. L., Frolking, S. & Moore, T. R. Interannual variability in the peatland-atmosphere carbon dioxide exchange at an ombrotrophic bog. Glob. Biogeochem. Cycles 17, 1036 (2003).Article 
CAS 

Google Scholar 
Joiner, D. W., Lafleur, P. M., McCaughey, J. H. & Bartlett, P. A. Interannual variability in carbon dioxide exchanges at a boreal wetland in the BOREAS northern study area. J. Geophys. Res. Atmos. 104, 27663–27672 (1999).CAS 
Article 

Google Scholar 
McVeigh, P., Sottocornola, M., Foley, N., Leahy, P. & Kiely, G. Meteorological and functional response partitioning to explain interannual variability of CO2 exchange at an Irish Atlantic blanket bog. Agric. For. Meteorol. 194, 8–19 (2014).Article 

Google Scholar 
Helbig, M. et al. Increasing contribution of peatlands to boreal evapotranspiration in a warming climate. Nat. Clim. Change 10, 555–560 (2020).CAS 
Article 

Google Scholar 
Bourgault, M.-A., Larocque, M. & Garneau, M. How do hydrogeological setting and meteorological conditions influence water table depth and fluctuations in ombrotrophic peatlands? J. Hydrol. X 4, 100032 (2019).Article 

Google Scholar 
Yurova, A., Wolf, A., Sagerfors, J. & Nilsson, M. Variations in net ecosystem exchange of carbon dioxide in a boreal mire: modeling mechanisms linked to water table position. J. Geophys. Res. Biogeosci. 112, G02025 (2007).Article 
CAS 

Google Scholar 
Laine, A. M. et al. Warming impacts on boreal fen CO2 exchange under wet and dry conditions. Glob. Change Biol. 25, 1995–2008 (2019).Article 

Google Scholar 
Chivers, M. R., Turetsky, M. R., Waddington, J. M., Harden, J. W. & McGuire, A. D. Effects of experimental water table and temperature manipulations on ecosystem CO2 fluxes in an Alaskan rich fen. Ecosystems 12, 1329–1342 (2009).CAS 
Article 

Google Scholar 
Juszczak, R. et al. Ecosystem respiration in a heterogeneous temperate peatland and its sensitivity to peat temperature and water table depth. Plant Soil 366, 505–520 (2013).CAS 
Article 

Google Scholar 
Hao, D. et al. Estimating hourly land surface downward shortwave and photosynthetically active radiation from DSCOVR/EPIC observations. Remote Sens. Environ. 232, 111320 (2019).Article 

Google Scholar 
O’Donnell, J. A., Romanovsky, V. E., Harden, J. W. & McGuire, A. D. The effect of moisture content on the thermal conductivity of moss and organic soil horizons from black spruce ecosystems in interior Alaska. Soil Sci. 174, 646–651 (2009).Article 
CAS 

Google Scholar 
Nijp, J. J. et al. Rain events decrease boreal peatland net CO2 uptake through reduced light availability. Glob. Change Biol. 21, 2309–2320 (2015).Article 

Google Scholar 
Zhang, Y., Commane, R., Zhou, S., Williams, A. P. & Gentine, P. Light limitation regulates the response of autumn terrestrial carbon uptake to warming. Nat. Clim. Change 10, 739–743 (2020).CAS 
Article 

Google Scholar 
Samson, M. et al. The impact of experimental temperature and water level manipulation on carbon dioxide release in a poor fen in northern Poland. Wetlands 38, 551–563 (2018).Article 

Google Scholar 
Drever, C. R. et al. Natural climate solutions for Canada. Sci. Adv. 7, eabd6034 (2021).CAS 
Article 

Google Scholar 
Hemes, K. S., Runkle, B. R. K., Novick, K. A., Baldocchi, D. D. & Field, C. B. An ecosystem-scale flux measurement strategy to assess natural climate solutions. Environ. Sci. Technol. 55, 3494–3504 (2021).CAS 
Article 

Google Scholar 
Walker, T. W. N. et al. A systemic overreaction to years versus decades of warming in a subarctic grassland ecosystem. Nat. Ecol. Evol. 4, 101–108 (2020).Article 

Google Scholar 
Xu, B. et al. Seasonal variability of forest sensitivity to heat and drought stresses: a synthesis based on carbon fluxes from North American forest ecosystems. Glob. Change Biol. 26, 901–918 (2020).Article 

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

Google Scholar 
Joyce, P. et al. How robust Is the apparent break-down of northern high-latitude temperature control on spring carbon uptake? Geophys. Res. Lett. 48, e2020GL091601 (2021).Article 

Google Scholar 
Grant, R. F. et al. Changes in net ecosystem productivity of boreal black spruce stands in response to changes in temperature at diurnal and seasonal time scales. Tree Physiol. 29, 1–17 (2009).CAS 
Article 

Google Scholar 
Kwon, M. J. et al. Siberian 2020 heatwave increased spring CO2 uptake but not annual CO2 uptake. Environ. Res. Lett. 16, 124030 (2021).CAS 
Article 

Google Scholar 
Yu, Z., Griffis, T. J. & Baker, J. M. Warming temperatures lead to reduced summer carbon sequestration in the U.S. Corn Belt. Commun. Earth Environ. 2, 53 (2021).Article 

Google Scholar 
Wang, S. et al. Warmer spring alleviated the impacts of 2018 European summer heatwave and drought on vegetation photosynthesis. Agric. For. Meteorol. 295, 108195 (2020).Article 

Google Scholar 
Wang, T. et al. Emerging negative impact of warming on summer carbon uptake in northern ecosystems. Nat. Commun. 9, 5391 (2018).CAS 
Article 

Google Scholar 
Lin, X. et al. Siberian and temperate ecosystems shape Northern Hemisphere atmospheric CO2 seasonal amplification. Proc. Natl Acad. Sci. USA 117, 21079–21087 (2020).CAS 
Article 

Google Scholar 
Helbig, M. et al. Warming response of peatland CO2 sink is sensitive to seasonality in warming trends. Zenodo https://doi.org/10.5281/zenodo.6685222 (2022).Didan, K. MOD13Q1 MODIS/Terra Vegetation Indices 16-Day L3 Global 250 m SIN Grid V006 [Data set]. NASA EOSDIS Land Processes DAAC (2015); https://doi.org/10.5067/MODIS/MOD13Q1.006Harris, I., Osborn, T. J., Jones, P. & Lister, D. Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset. Sci. Data 7, 109 (2020).Article 

Google Scholar 
Lees, K. J. et al. Using spectral indices to estimate water content and GPP in Sphagnum moss and other peatland vegetation. IEEE Trans. Geosci. Remote Sens. 58, 4547–4557 (2020).Article 

Google Scholar 
Bennett, A. C., McDowell, N. G., Allen, C. D. & Anderson-Teixeira, K. J. Larger trees suffer most during drought in forests worldwide. Nat. Plants 1, 15139 (2015).Article 

Google Scholar 
Page, S. E. & Baird, A. J. Peatlands and global change: response and resilience. Annu. Rev. Environ. Resour. 41, 35–57 (2016).Article 

Google Scholar 
Juottonen, H. et al. Integrating decomposers, methane-cycling microbes and ecosystem carbon fluxes along a peatland successional gradient in a land uplift region. Ecosystems https://doi.org/10.1007/s10021-021-00713-w (2021). More

Database representativeness of global monoculture plantationsOur database includes 4756 measurements of carbon in aboveground live tree biomass in monoculture plantations, collected from 829 distinct sites across the globe (Fig. 1). These sites were primarily in Asia (59%), Europe (16%), North America (15%), and South America (4%), with the remainder located in Oceania and Africa (~5%). The dataset included 240 species from 96 distinct genera of tree; however, 33 genera were poorly represented (n  More

Vázquez, D. P., Gianoli, E., Morris, W. F. & Bozinovic, F. Ecological and evolutionary impacts of changing climatic variability. Biol. Rev. 92, 22–42. https://doi.org/10.1111/brv.12216 (2017).Article 
PubMed 

Google Scholar 
Marshall, K. E. & Sinclair, B. J. The impacts of repeated cold exposure on insects. J. Exp. Biol. 215, 1607–1613. https://doi.org/10.1242/jeb.059956 (2012).Article 
PubMed 

Google Scholar 
Bale, J. & Hayward, S. Insect overwintering in a changing climate. J. Exp. Biol. 213, 980–994 (2010).CAS 
PubMed 
Article 

Google Scholar 
Kingsolver, J. G. Feeding, growth, and the thermal environment of cabbage white caterpillars, Pieris rapae L. Physiol. Biochem. Zool. 73, 621–628 (2000).CAS 
PubMed 
Article 

Google Scholar 
Stange, E. E. & Ayres, M. P. Climate change impacts: Insects (JohnWiley & Sons, 2010).
Google Scholar 
Chapman, A. D. Numbers of Living Species in Australia and the World: Report for the Department of the Environment and Heritage Canberra, Australia (Department of the Environment and Heritage, 2006).
Google Scholar 
Colinet, H., Sinclair, B. J., Vernon, P. & Renault, D. Insects in fluctuating thermal environments. Annu. Rev. Entomol. 60, 123–140. https://doi.org/10.1146/annurev-ento-010814-021017 (2015).CAS 
Article 
PubMed 

Google Scholar 
Hickling, R., Roy, D. B., Hill, J. K. & Thomas, C. D. A northward shift of range margins in British Odonata. Glob. Change Biol. 11, 502–506. https://doi.org/10.1111/j.1365-2486.2005.00904.x (2005).ADS 
Article 

Google Scholar 
Rumpf, S. B., Hülber, K., Zimmermann, N. E. & Dullinger, S. Elevational rear edges shifted at least as much as leading edges over the last century. Glob. Ecol. Biogeogr. 28, 533–543. https://doi.org/10.1111/geb.12865 (2019).Article 

Google Scholar 
Halsch, C. A. et al. Insects and recent climate change. Proc. Natl. Acad. Sci. 118, e2002543117. https://doi.org/10.1073/pnas.2002543117 (2021).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
McCain, C. M. & Garfinkel, C. F. Climate change and elevational range shifts in insects. Curr. Opin. Insect Sci. 47, 111–118. https://doi.org/10.1016/j.cois.2021.06.003 (2021).Article 
PubMed 

Google Scholar 
Angilletta, M. J. Jr. & Angilletta, M. J. Thermal Adaptation: A Theoretical and Empirical Synthesis (Oxford University Press, 2009).Book 

Google Scholar 
Angilletta, M. J. & Dunham, A. E. The temperature-size rule in ectotherms: Simple evolutionary explanations may not be general. Am. Nat. 162, 332–342. https://doi.org/10.1086/377187 (2003).Article 
PubMed 

Google Scholar 
Jensen, J. L. W. V. Sur les fonctions convexes et les inégalités entre les valeurs moyennes. Acta Math. 30, 175–193 (1906).MathSciNet 
MATH 
Article 

Google Scholar 
Ruel, J. J. & Ayres, M. P. Jensen’s inequality predicts effects of environmental variation. Trends Ecol. Evol. 14, 361–366 (1999).CAS 
PubMed 
Article 

Google Scholar 
Kingsolver, J. G. & Woods, H. A. Thermal sensitivity of growth and feeding in Manduca sexta Caterpillars. Physiol. Zool. 70, 631–638. https://doi.org/10.1086/515872 (1997).CAS 
Article 
PubMed 

Google Scholar 
Gillooly, J. F., Brown, J. H., West, G. B., Savage, V. M. & Charnov, E. L. Effects of size and temperature on metabolic rate. Science 293, 2248–2251 (2001).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Bale, J. S. et al. Herbivory in global climate change research: Direct effects of rising temperature on insect herbivores. Glob. Change Biol. 8, 1–16 (2002).ADS 
Article 

Google Scholar 
Robinet, C. & Roques, A. Direct impacts of recent climate warming on insect populations. Integr. Zool. 5, 132–142 (2010).PubMed 
Article 

Google Scholar 
García-Robledo, C., Kuprewicz, E. K., Staines, C. L., Erwin, T. L. & Kress, W. J. Limited tolerance by insects to high temperatures across tropical elevational gradients and the implications of global warming for extinction. Proc. Natl. Acad. Sci. 113, 680–685. https://doi.org/10.1073/pnas.1507681113 (2016).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Vasseur, D. A. et al. Increased temperature variation poses a greater risk to species than climate warming. Proc. R. Soc. B Biol. Sci. 281, 20132612 (2014).Article 

Google Scholar 
Sandehson, D. E. The relation of temperature to the growth of insects. J. Econ. Entomol. 3, 113–140 (1910).Article 

Google Scholar 
Cook, W. C. Some Effects of Alternating Temperatures on the Growth and Metabolism of Cutworm Larvae (Oxford University Press, 1927).
Google Scholar 
Kingsolver, J. G., Ragland, G. J. & Diamond, S. E. Evolution in a constant environment: Thermal fluctuations and thermal sensitivity of laboratory and field populations of Manduca sexta. Evolution 63, 537–541. https://doi.org/10.1111/j.1558-5646.2008.00568.x (2009).Article 
PubMed 

Google Scholar 
Eldridge, W. H., Sweeney, B. W. & Law, J. M. Fish growth, physiological stress, and tissue condition in response to rate of temperature change during cool or warm diel thermal cycles. Can. J. Fish. Aquat. Sci. 72, 1527–1537 (2015).CAS 
Article 

Google Scholar 
Bernhardt, J. R., Sunday, J. M., Thompson, P. L. & O’Connor, M. I. Nonlinear averaging of thermal experience predicts population growth rates in a thermally variable environment. Proc. R. Soc. B Biol. Sci. 285, 20181076. https://doi.org/10.1098/rspb.2018.1076 (2018).Article 

Google Scholar 
Morissette, J., Swart, S., Maccormack, T. J., Currie, S. & Morash, A. J. Thermal variation near the thermal optimum does not affect the growth, metabolism or swimming performance in wild Atlantic salmon Salmo salar. J. Fish Biol. 98, 1585–1589. https://doi.org/10.1111/jfb.14348 (2021).CAS 
Article 
PubMed 

Google Scholar 
Bozinovic, F. et al. The mean and variance of environmental temperature interact to determine physiological tolerance and fitness. Physiol. Biochem. Zool. 84, 543–552 (2011).PubMed 
Article 

Google Scholar 
Boggs, C. L. The fingerprints of global climate change on insect populations. Curr. Opin. Insect Sci. 17, 69–73. https://doi.org/10.1016/j.cois.2016.07.004 (2016).Article 
PubMed 

Google Scholar 
Lemoine, N. P., Drews, W. A., Burkepile, D. E. & Parker, J. D. Increased temperature alters feeding behavior of a generalist herbivore. Oikos 122, 1669–1678. https://doi.org/10.1111/j.1600-0706.2013.00457.x (2013).Article 

Google Scholar 
Vangansbeke, D. et al. Prey consumption by phytoseiid spider mite predators as affected by diurnal temperature variations. Biocontrol 60, 595–603 (2015).Article 

Google Scholar 
Davies, C., Coetzee, M. & Lyons, C. L. Effect of stable and fluctuating temperatures on the life history traits of Anopheles arabiensis and An. quadriannulatus under conditions of inter-and intra-specific competition. Parasit. Vectors 9, 1–9 (2016).Article 

Google Scholar 
Delava, E., Fleury, F. & Gibert, P. Effects of daily fluctuating temperatures on the Drosophila-Leptopilina boulardi parasitoid association. J. Therm. Biol. 60, 95–102 (2016).PubMed 
Article 

Google Scholar 
Amarasekare, P. & Coutinho, R. M. Effects of temperature on intraspecific competition in ectotherms. Am. Nat. 184, E50–E65. https://doi.org/10.1086/677386 (2014).Article 
PubMed 

Google Scholar 
Jiang, L. & Morin, P. J. Temperature-dependent interactions explain unexpected responses to environmental warming in communities of competitors. J. Anim. Ecol. 73, 569–576 (2004).Article 

Google Scholar 
Novich, R. A., Erickson, E. K., Kalinoski, R. M. & DeLong, J. P. The temperature independence of interaction strength in a sit-and-wait predator. Ecosphere 5, 1–9 (2014).Article 

Google Scholar 
Fox, L. R. Cannibalism in natural populations. Annu. Rev. Ecol. Syst. 6, 87–106 (1975).Article 

Google Scholar 
Polis, G. A. The evolution and dynamics of intraspecific predation. Annu. Rev. Ecol. Syst. 12, 225–251 (1981).Article 

Google Scholar 
Nishimura, K. & Isoda, Y. Evolution of cannibalism: Referring to costs of cannibalism. J. Theor. Biol. 226, 293–302. https://doi.org/10.1016/j.jtbi.2003.09.007 (2004).ADS 
MathSciNet 
Article 
PubMed 
MATH 

Google Scholar 
Crumrine, P. W. Body size, temperature, and seasonal differences in size structure influence the occurrence of cannibalism in larvae of the migratory dragonfly, Anax junius. Aquat. Ecol. 44, 761–770 (2010).Article 

Google Scholar 
Reglero, P., Urtizberea, A., Torres, A. P., Alemany, F. & Fiksen, Ø. Cannibalism among size classes of larvae may be a substantial mortality component in tuna. Mar. Ecol. Prog. Ser. 433, 205–219 (2011).ADS 
Article 

Google Scholar 
Nilsson-Örtman, V., Stoks, R. & Johansson, F. Competitive interactions modify the temperature dependence of damselfly growth rates. Ecology 95, 1394–1406. https://doi.org/10.1890/13-0875.1 (2014).Article 
PubMed 

Google Scholar 
Pritchard, G. & Leggott, M. Temperature, incubation rates and the origins of dragonflies. Adv. Odonatol. 3, 121–126 (1987).
Google Scholar 
Hassall, C. & Thompson, D. J. The effects of environmental warming on Odonata: A review. Int. J. Odonatol. 11, 131–153 (2008).Article 

Google Scholar 
Johansson, F. & Crowley, P. H. Larval cannibalism and population dynamics of dragonflies. In Aquatic Insects: Challenges to Populations 36–54 (CABI, 2008).Rudolf, V. H. W. & Rasmussen, N. L. Ontogenetic functional diversity: Size structure of a keystone predator drives functioning of a complex ecosystem. Ecology 94, 1046–1056 (2013).PubMed 
Article 

Google Scholar 
Hyeun-Ji, L. & Johansson, F. Compensating for a bad start: Compensatory growth across life stages in an organism with a complex life cycle. Can. J. Zool. 94, 41–47 (2016).Article 

Google Scholar 
Sokolovska, N., Rowe, L. & Johansson, F. Fitness and body size in mature odonates. Ecol. Entomol. 25, 239–248. https://doi.org/10.1046/j.1365-2311.2000.00251.x (2000).Article 

Google Scholar 
Karl, T. R. Modern global climate change. Science 302, 1719–1723. https://doi.org/10.1126/science.1090228 (2003).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Meehl, G. A. et al. Climate Change 2007: the physical science basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge University Press, 2007).
Google Scholar 
Meehl, G. A. More intense, more frequent, and longer lasting heat waves in the 21st century. Science 305, 994–997. https://doi.org/10.1126/science.1098704 (2004).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Khelifa, R. Spatiotemporal Pattern of Phenology Across Geographic Gradients in Insects, in Chapter 1 (Geographic Gradients in Climate Change Response Explained by Non-linear Thermal-Performance Curves) (University of Zurich, 2017).
Google Scholar 
Boudot, J. P. & Kalkman, V. Atlas of the European Dragonflies and Damselflies (KNNV Publishing, 2015).
Google Scholar 
Norling, U. Growth, winter preparations and timing of emergence in temperate zone Odonata: Control by a succession of larval response patterns. Int. J. Odonatol. 24, 1–36 (2021).Article 

Google Scholar 
Sniegula, S. & Johansson, F. Photoperiod affects compensating developmental rate across latitudes in the damselfly Lestes sponsa. Ecol. Entomol. 35, 149–157. https://doi.org/10.1111/j.1365-2311.2009.01164.x (2010).Article 

Google Scholar 
Sniegula, S., Golab, M. J. & Johansson, F. Size-mediated priority and temperature effects on intra-cohort competition and cannibalism in a damselfly. J. Anim. Ecol. 88, 637–648. https://doi.org/10.1111/1365-2656.12947 (2019).Article 
PubMed 

Google Scholar 
Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675. https://doi.org/10.1038/nmeth.2089 (2012).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Benke, A. C. A method for comparing individual growth rates of aquatic insects with special reference to the Odonata. Ecology 51, 328–331 (1970).Article 

Google Scholar 
Nilsson-Örtman, V., Stoks, R., De Block, M. & Johansson, F. Generalists and specialists along a latitudinal transect: Patterns of thermal adaptation in six species of damselflies. Ecology 93, 1340–1352. https://doi.org/10.1890/11-1910.1 (2012).Article 
PubMed 

Google Scholar 
Eklund, A. et al. Sveriges Framtida Klimat: Underlag Till Dricksvattenutredningen (SMHI, 2015).
Google Scholar 
McPeek, M. A. Determination of species composition in the Enallagma damselfly assemblages of permanent lakes. Ecology 71, 83–98. https://doi.org/10.2307/1940249 (1990).Article 

Google Scholar 
Kirillin, G. et al. FLake-global: Online lake model with worldwide coverage. Environ. Model. Softw. 26, 683–684. https://doi.org/10.1016/j.envsoft.2010.12.004 (2011).Article 

Google Scholar 
SMHI. Advanced Climate Change Scenario Service. https://www.smhi.se.R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2019).Suhling, F., Suhling, I. & Richter, O. Temperature response of growth of larval dragonflies–An overview. Int. J. Odonatol. 18, 15–30 (2015).Article 

Google Scholar 
Padfield, D., O’Sullivan, H. & Pawar, S. rTPC and nls.multstart: A new pipeline to fit thermal performance curves in R. Methods Ecol. Evol. 1, 1. https://doi.org/10.1111/2041-210X.13585 (2021).Article 

Google Scholar 
Kleiber, M. Body size and metabolism. Hilgardia 6, 315–353 (1932).CAS 
Article 

Google Scholar 
Hemmingsen, A. Reports of the Steno Memorial Hospital and Nordisk Insulin Laboratorium. Energy Metab. Relat. Body Size Respir. Surf. Evol. 9, 6–110 (1960).
Google Scholar 
Dell, A. I., Pawar, S. & Savage, V. M. Systematic variation in the temperature dependence of physiological and ecological traits. Proc. Natl. Acad. Sci. 108, 10591–10596 (2011).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Logan, J. D., Wolesensky, W. & Joern, A. Temperature-dependent phenology and predation in arthropod systems. Ecol. Model. 196, 471–482. https://doi.org/10.1016/j.ecolmodel.2006.02.034 (2006).Article 

Google Scholar 
Pink, M. & Abrahams, M. V. Temperature and its impact on predation risk within aquatic ecosystems. Can. J. Fish. Aquat. Sci. 73, 869–876. https://doi.org/10.1139/cjfas-2015-0302 (2016).Article 

Google Scholar 
DeAngelis, D., Cox, D. & Coutant, C. Cannibalism and size dispersal in young-of-the-year largemouth bass: Experiment and model. Ecol. Model. 8, 133–148 (1980).Article 

Google Scholar 
Fagan, W. F. & Odell, G. M. Size-dependent cannibalism in praying mantids: Using biomass flux to model size-structured populations. Am. Nat. 147, 230–268 (1996).Article 

Google Scholar 
Dong, Q. & Deangelis, D. L. Consequences of cannibalism and competition for food in a smallmouth bass population: An individual-based modeling study. Trans. Am. Fish. Soc. 127, 174–191 (1998).Article 

Google Scholar 
Verheyen, J. & Stoks, R. Temperature variation makes an ectotherm more sensitive to global warming unless thermal evolution occurs. J. Anim. Ecol. 88, 624–636. https://doi.org/10.1111/1365-2656.12946 (2019).Article 
PubMed 

Google Scholar 
Starr, S. M. & McIntyre, N. E. Effects of water temperature under projected climate change on the development and survival of Enallagma civile (Odonata: Coenagrionidae). Environ. Entomol. 49, 230–237. https://doi.org/10.1093/ee/nvz138 (2020).CAS 
Article 
PubMed 

Google Scholar 
Culler, L. E., McPeek, M. A. & Ayres, M. P. Predation risk shapes thermal physiology of a predaceous damselfly. Oecologia 176, 653–660. https://doi.org/10.1007/s00442-014-3058-8 (2014).ADS 
Article 
PubMed 

Google Scholar 
Dokulil, M. T. et al. Increasing maximum lake surface temperature under climate change. Clim. Change https://doi.org/10.1007/s10584-021-03085-1 (2021).Article 

Google Scholar 
Merritt, R. W. & Cummins, K. W. An Introduction to the Aquatic Insects of North America 2nd edn. (Kendall/Hunt Publishing Company, 1984).
Google Scholar 
Verdonschot, R. & Peeters, E. T. Preference of larvae of Enallagma cyathigerum (Odonata: Coenagrionidae) for habitats of varying structural complexity. Eur. J. Entomol. 109, 229–234 (2012).Article 

Google Scholar 
McCarty, J. P., Wolfenbarger, L. L. & Wilson, J. A. eLS 1–13 (Wiley, 2017).Book 

Google Scholar 
Holzmann, K. L. Challenges in a Changing Climate: The Effect of Temperature Variation on Growth and Competition in Damselflies Independent thesis Advanced level (degree of Master (Two Years) thesis, Uppsala University. http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-467582 (2022). More

Data management and workflowsWe adopt best-practice approaches for open and reproducible research planning, execution, reporting, and management throughout the project (see e.g.29,30,31,32). Specifically, we use community-approved standards for experimental design and data collection29, and clean and manage the data using a fully scripted and reproducible data workflow, with data and code deposited at open repositories (Fig. 2).Fig. 2The data collection and management workflow of the FunCaB project. Reproducibility throughout the research process is assured as follows: Experimental design and data collection was based on best-practice community methods and protocols, adapted for the projects’ needs. Measurements were digitalized and the raw data stored in the project Open Science Foundation (OSF) repository before the raw data were cleaned and managed through code-based data curation, with version control secured via GitHub. The clean data are stored at the OSF repository, and a time-stamped version of the code to retrieve and clean data is provided through Zenodo. This data paper describes and documents the data collection and workflow, and describes how to access and use clean data, raw data, and code.Full size imageResearch site selection and basic information, and general study setupSite selectionOur study is conducted across the twelve calcareous grassland experimental sites in the Vestland Climate Grid (VCG), in south-western Norway (Fig. 1a). The VCG sites were chosen to fit within a climate grid reflecting a fully factorial design encompassing the major bioclimatic variation in Norway. Potential sites were identified using a combination of topographic maps, geological maps (NGU) and interpolated maps of summer temperature and annual precipitation using the 1960–1990 climate normal (100 m resolution gridded data, met.no; see33 and references therein). The three temperature levels (alpine, sub-alpine, boreal) and four levels of precipitation in the climate grid (Fig. 1b) were selected to reflect a difference in mean growing season temperature of ca. 2 °C between three temperature levels (alpine = 6.5 °C, sub-alpine = 8.5 °C, boreal = 10.5 °C mean temperature of the four warmest months of the year) and a difference in mean annual precipitation of 700 mm between four precipitation levels (precipitation levels 1 – 4 representing 700 mm, 1400 mm, 2100 mm, and 2800 mm, respectively). Climate data for the site selection was based on 100-m resolution downscaled data using the 1960–1990 climate normal from met.no. The final sites were selected from approximately 200 potential sites visited and surveyed in the summer of 2008, with selection criteria set to ensure that other factors such as grazing regime and history, bedrock, vegetation type and structure, slope and exposure were kept as constant as possible among the selected sites34. Geographical distance between sites is on average 15 km and ranges from 175 km to 650 m.Study system and experimental area selection within sitesAt each site, we selected an experimental area of ca. 75 –200 m2, targeting a homogeneous and representative part of the target grassland vegetation at large at that site. The experimental areas were placed on southerly-facing slopes, avoiding depressions and concave areas in the landscape and other features such as big rocks or formations that may affect light conditions, hydrology and/or snowdrift. The target vegetation type was forb-rich semi-natural upland grassland vegetation34,35, within the plant sociological association Potentillo-Festucetum ovinae tending towards Potentillo-Poligonium vivipara in the alpine sites and Nardo-Agrostion tenuis in some lowland sites36. The most common vascular plants across sites, based on sum of covers, are the graminoids Agrostis capillaris, Festuca rubra, Avenella flexuosa, Anthoxanthum odoratum, and Nardus stricta and the forbs Leucantemum vulgare, Hypericum maculatum, Silene acaulis, Alchemilla alpina, and Lotus corniculatus. Common bryophytes are Pleurotium schreberi, Hylocomium splendens, Polythricum spp, Racomitrium lanuginosum, R. fasciculare, and Dicranum spp. All sites were moderately grazed prior to the study by sheep, cattle, goats, reindeer, deer, moose, and/or horses; and the experimental areas were fenced for the duration of the study to prevent animal and human disturbance of the experimental infrastructure. The fenced area was lightly mowed at the end of each growing season to mimic past grazing pressure and minimize fence effects. For further description of the sites, see34 and for access to and further description of site-level data, see35.Block and experimental plot setupWithin these study areas we established four blocks, with a distance between the blocks ranging from one up to (in rare cases) 50 meters. Blocks were selectively placed in homogenous grassland vegetation, avoiding rocks, depressions, and other features as described above. Each block accommodates eight 25 × 25 cm plots, with at least 25 cm between adjacent plots. If a plot contained more than 10% bare rock, shrubs, or other non-grassland features, they were rejected or moved slightly to avoid these features. The plots were permanently marked with four aluminium 10-cm long pipes in the soil in the outer corners of all the 25 × 25 cm treatment plots, ensuring the pipes to fit the corners of a standardized vegetation analysis frame (aluminium frame demarking a 25 × 25 cm inner area, with poles fixed in the corners that fit into the aluminium tubes used for plot demarcation in the field). The upslope left corner tube was marked with a colour-coded waterproof tape. Note that in 31 out of 48 cases (12 sites × 4 blocks), the blocks were located within larger experimental blocks in the VCG sites, and control plots and various block-level data are then shared with other experiments in these larger blocks. Linking keys are described in the FunCaB data dictionaries below (see Fig. 3 and data records iii-vii below). For some datasets, additional plots within blocks were needed. These are described as needed below.Fig. 3Data structure for the FunCaB functional group removal experiment and associated Vestland Climate Grid (VCG) and FUNDER project data. Within each of the three projects, boxes represent data tables. The FunCaB project data tables include biomass of functional groups removed and forb species-level biomass (datasets i, ii), soil temperature and moisture (datasets iii, iv) plant community composition and the associated taxon table (dataset v), seedling recruitment (dataset vi), ecosystem carbon fluxes (dataset vii) and reflectance (dataset viii). Names of individual data tables are given in the coloured title area, and a selection of the main variables available within tables in the internal lists. For full sets of variables for each FunCaB dataset, see Tables 3–9. The lines linking three of the boxes exemplify links using species as keys across tables, note that all bold variables are shared between several tables and can be used as keys to join them. Keys can also be used to link to/from data from other projects in the VCG (for general VCG project keys, see top right hatched outline box, for keys between the FunCaB and FUNDER projects see the bottom right hatched outline box (both including an example value for each variable on the right). The (other) datasets* boxes refer to extensive datasets on plant community composition, cover, biomass, fitness, and reproduction available from previous projects in the VCG27 and upcoming datasets in the FUNDER project.Full size imageBackground abiotic and biotic data from the Vestland Climate GridThe Vestland Climate Grid field sites were established in 2008, and from a series of research projects within the grid over the years we have collected a broad range of datasets on the climate and environment, soils, land-use and environment, vegetation, and ecosystems, along with basic descriptive data of the 12 sites, as described in34. All these datasets are available from the previous projects through the VCG OSF (Open Science Framework) repository35, and the results are presented in associated papers, see for example34,37,38,39,40,41,42,43,44,45. The overall data structure, and the most relevant datasets from the VCG for the FunCaB project is laid out in Fig. 3, and briefly described below. Code to download and link these data to the FunCaB experimental data and sites are provided in the FunCaB github repository28 (see R/download_VCG_data).A new research project, ‘FUNDER – Direct and indirect climate impacts on the biodiversity and FUNctioning of the UNDERground ecosystem’ funded by the Norwegian Research Council KLIMAFORSK programme (project number 315249, 2021 – 2025) will augment the FunCaB experiment with data on the belowground components of the plant-soil ecosystem, including roots, mesofauna, fungi and microbes. These upcoming data will all link with the FunCaB and VCG project based on the given experimental, site and organismal keys, as indicated in Fig. 3.VCG Basic site-level attributesBasic descriptive data on the 12 sites include latitude, longitude, elevation, geology, land-use, soils, and their position within the climate grid (precipitation and temperature levels). These data are described in34,40, provided in35, and can be downloaded using28 (see R/download_VCG_data). For convenience, the climate grid information is also provided in the biomass dataset (see below).VCG Site-level climate dataTemperature was measured continuously at each of the 12 VCG sites at four heights (2 m and 30 cm above ground, at ground level, and 5 cm below ground), soil moisture was measured continuously with two replicate loggers ca. 5 cm below ground, and precipitation was measured at each site during the snow-free season. For these measurements, we used Delta T GP1 loggers (Delta T devices, Cambridge, UK) equipped with two temperature probes, two SM200 moisture sensors which were later replaced as necessary with SM300 and SM150T loggers, and an ARG 100 tipping bucket (EML LTD, North Shields, UK) from 2009 onwards. UTL-3 version 3.0 temperature loggers (GEOTEST AG, Zollikofen, Switzerland) were used for measuring the 2 m and 30 cm temperatures. Soil moisture was measured as the mean of four measurements taken along each side of the turf, several times during the growing season using a Delta T HH2 version 2.3 Moisture Meter with the same probes as for the GP1 logger (SM200, SM150T). These data are described in34,40, provided in35, and can be downloaded using28 (see folder R/download_VCG_data).VCG Soil chemical and structural dataOver the years, various soil chemical variables have been measured at the block level within each of the 12 VCG sites, including soil pH (2009) and % Loss-On-Ignition (2009, 2013), and available N, as sum of N available as NH4-N and NO3-N (available N per deployment period, 2010 & 2012). Soil pH was measured after adding 50 ml distilled water to 25 g soil and mixing for two hours. Loss-on-ignition (LOI), was measured by weighing dry soil (105 °C for 24, one hour in a desiccator), and burnt soil (six hours at 550 °C, one hour in the desiccator) and calculating LOI as the (burnt soil mass/dry soil mass) × 100. NH4-N and NO3-N were measured using in-situ ion exchange resin bags (IERBs) were used to measure the amount of plant-available nutrients in the soil. These data are partially described in34,40, and the full data are provided in35.VCG Litter decomposition dataDecomposition has been assessed at each of the 12 VCG sites using local plant litter and the Tea Bag Index method (Keuskamp et al., 2013). Local litter (dead leaves detached from live plants) was collected at each site in 2013 or 2014, with the specific timing of the collections at each site tuned to ensure that litter was present, not covered by snow, and not decomposed. In practice, this necessitated litter collection after snowmelt in spring in many sites. The litter was washed, dried, and stored in dark, dry, cool conditions. In 2016, five replicate litter bags containing 1 g of graminoid litter were buried at each site, and collected at four points in time after burial (1, 2, 3 and 12 months). Harvested litter bags were cleaned (soil and roots removed), dried for 48 h at 60 °C and weighed. The Tea Bag Index method46 was used in 2014, 2015 and 2016 to measure decomposition at all sites of the climate grid. At each site, 10 replicates of each tea type were buried pair-wise, 8 cm below ground and with at least 10 cm between the tea bags. For a couple of sites, the number of replicate tea bag pairs was higher in 2015 (12 replicates at the site Gudmedalen and 16 replicates at Låvisdalen). After collection, adhering soil particles and roots were removed and the tea bags were dried (48 h at 60 °C) and weighed. These data are partially described in47, and the full data are provided in35 and can be downloaded using28 (see folder R/download_VCG_data).VCG Species-level cover, biomass, and performance dataA variety of plant species and community composition, cover, biomass, fitness, and reproductive data exists for the sites and blocks in the VCG from 2008 to 2021. These data are described in e.g34,37,38,41,43,44,45,48,49,50, and provided in35.VCG Site-level plant functional traitsIn 2016 and 2017, we measured 11 leaf functional traits that are related to potential physiological growth rates and environmental tolerance of plants, following the standardized protocols in Pérez-Harguindeguy et al.51: leaf area (LA, cm2), leaf thickness (LT, mm), leaf dry matter content (LDMC, g/g), specific leaf area (SLA, cm2/g), carbon (C, %), nitrogen (N, %), phosphorus (P, %), carbon nitrogen ratio (C:N), nitrogen phosphorus ratio (N:P), carbon13 isotope ratio (δ13C, ‰), and nitrogen15 isotope ratio (δ15N, ‰). Trait data are available at the site level for species making up at least 80% of the vegetation cover in the control plots at each of the 12 VCG sites. The plants were collected outside of the experimental plots and within a 50 m perimeter from the blocks, and we aimed to collect up to five individuals from each species in each site. To avoid repeated sampling from a single clone, we selected individuals that were visibly separated from other ramets of that species. The sampled plant individuals were labelled, put in plastic bags with moist paper towels, and stored in darkness at 4 °C until processing, which was done as soon as possible and always within 4 days. These data are described in52, provided in35, and can be downloaded using28 (see folder R/download_VCG_data).Experimental designThe functional group removal experiment was designed to examine the impact of aboveground interactions among the major plant functional groups – graminoids, forbs and bryophytes – on the performance and functioning of other components of the vegetation and ecosystem. The experiment consists of eight 25 × 25 cm plots per site and block, with a fully factorial combination of removals of three plant functional groups, with treatments randomized within each block. The general experimental design, with the different removal treatments detailed, are provided as an insert to the timeline in Fig. 1c. The functional groups are delineated and abbreviated in the various datasets as follows: G = graminoids (including grasses, sedges and rushes), F = forbs (including herbaceous forbs, pteridophytes, dwarf-shrubs, and small individuals of trees and shrubs), B = Bryophytes (including mosses, liverworts, and hornworts). Note that all species are also coded by their respective functional group into which they were classified in the FunCaB taxon table. The experimental treatments are coded by functional group removed so that FGB = bare-ground gaps with all plants removed, FB = only graminoids remaining, GB = only forbs remaining, GF = only bryophytes remaining, B = graminoids and forbs remaining, F = bryophytes and graminoids remaining, G = bryophytes and forbs remaining, and C = intact vegetation controls with no vegetation removed. In 2016, four extra control (XC) plots were marked per site for aboveground biomass harvest and ecosystem carbon flux measurements. This sampling regime gave a total of 384 plots in the core FunCaB experiment, plus the additional 48 controls in 2016.Functional group removals were done once in 2015 (at peak growing season due to late snowmelt), twice per year in 2016 and 2017 (after the spring growth and at peak growing season) and annually from 2018 to 2021 (at peak growing season) as regrowth had declined (see below) and biannual removals were no longer necessary. At each sampling, all above-ground biomass of the relevant plant functional group was removed from each plot as follows: for each plot, all the above-ground parts of the relevant functional group(s) were removed using scissors and tweezers to cut the plants at the ground layer (i.e., the soil-vegetation interface). Roots and other below-ground parts were not removed, and non-target plant functional groups and litter were left intact.Species identification, taxonomy, and floraAll vascular plant species were identified to the species level in the field, with nomenclature following Lid and Lid53. Exceptions are sterile specimens of species that are not possible to identify without reproductive parts, and where flowers are either too rare or individuals too short-lived for comparisons of the position of individuals within the plots over years to be used to ascertain identifications (For example, Alchemilla spp. excluding A. alpina, and the annual Euphrasia spp.). Species identifications were confirmed by comparing records over time as described below. All unidentified specimens are included and flagged in the dataset, as described in Data Records below. The full taxon names are provided in the taxon table on OSF (Fig. 3).Dataset collection methodsDatasets (i–ii): Biomass and functional group removalAs described above, functional group removals were done once in 2015 at peak growing season, and twice per year in 2016 and 2017 (after the spring growth; at peak growing season) and annually at peak growing season from 2018 to 2021. For each removal plot and occasion, a picture was taken of the plot pre-removal, the biomass to be removed was collected in separate pre-marked paper bags for each functional group (graminoids, forbs and bryophytes), and a picture was taken post-removal. The collected biomass was then dried at 60 °C for 48 hours and weighed to the nearest 0.01 g (Model LPG-1002, VWR). From the four extra control (XC) plots in 2016, total above-ground biomass as well as litter (defined as dead biomass detached from live plants, see28) was collected at peak growing season. From these plots, biomass was sorted into functional groups as described above, except the forb functional group, which was sorted into species. The graminoid and bryophyte functional groups, each forb species, and litter were individually dried and weighed as described above. The data is available as (i) a biomass dataset, consisting of the removed biomass per plot, date, removal treatment, and functional group for all treatment plots, and the total biomass per functional group plus litter for the extra control plots in 2016, and (ii) a species-level forb biomass dataset from the extra control plots in 2016 (Fig. 3, Table 1).Datasets (iii-iv) – Soil microclimateWe measured soil temperature 3–5 cm below the soil surface for each plot using iButton temperature sensors (DS1922L, Manufacturer reports temperature accuracy of ±0.5 °C, Maxim Integrated INC., San Jose, CA, USA). The data are reported with a resolution of 0.0625 at 140 min intervals from June 2015 to July 2016. We measured soil moisture as volumetric soil moisture; expressed as % water volume per soil volume ((m3 water /m3 soil) × 100). These measurements were done c. 3–5 times during the growing seasons from 2015–2019, usually in connection with the flux and vegetation measurements, by taking the average of four measurements, one at each side of each plot (SM300, Manufacturer reports accuracy ±2.5% vol over 0 to 50% vol and 0–60 °C, Delta-T Devices, Cambridge, UK). The data is available as (iii) temperature and (iv) volumetric soil moisture % per plot and time point (temperature) or date (moisture) (Fig. 3, Table 1).Dataset (v): Vascular plant community composition and vegetation structureWe recorded the full vascular plant species composition of all experimental plots in 2015 (pre treatment), and the control plots plus the extra control plots in 2016. In 2017, 2018, and 2019, we recorded the community composition in controls and in the functional groups that remained in the experimental plots according to the plot’s treatment. At each analysis, each plot was sub-divided into 25 5 × 5 cm subplots, using a subplot overlay. We first recorded all species of vascular plants in the central five subplots, (i.e., the central + shaped area of each plot, Fig. 1c) noting the subplot cover of each species present in each of the five subplots (1 – 25% = 1, 26 – 50% = 2, 51 – 75% = 3, >76% = 4). Additionally, we noted if the individual was fertile (records circled if buds, flowers, or fruits were present). The five subplots were recorded and numbered (1-5) by row, and from left to right, starting from the top up-slope subplot. For the entire 25 × 25 cm plot, any additional species not present in one of the central subplots were recorded and their fertility noted. We then visually estimated the percentage cover of each vascular plant species in the whole plot to the nearest 1% and measured vegetation height in mm at four points within the plot. Note that the total coverage in each plot can exceed 100% due to layering of the vegetation. The vascular plant vegetation data is available as percentage cover and fertility status (sterile or fertile) per species per subplot and plot per sampling date, and vegetation height in mm per plot per sampling date (Fig. 3, Table 1).Other variables that were measured were percentage cover of bryophytes, litter, bare ground, and rock (measured per plot and per subplot) and moss layer depth in mm (mean of 4 measurements/plot), date of analysis, recorder/scribe (if any), and free-text comments. These data are available as % cover, depth in mm, date (year.month.day) and text strings per subplot and /or plot per sampling date (Fig. 3, Table 1).Dataset (vi): Seedling recruitmentThe total number of forb seedlings that emerged in the plots was recorded in 2018 and 2019. At peak growing season in 2018 (round 1, July-August, depending on site), all dicotyledonous seedlings were marked with wooden toothpicks and their x and y coordinates in the plot (mm, recorded from the bottom left hand-corner of the plot, Fig. 1c) and tentative species identity noted. Toward the end of the growing season (round 2, August-September, depending on site), each plot was revisited, seedling survival established, and any further seedlings marked. Survival (recorded when a seedling was present in subsequent surveys; recorded as mortality if absent) and new seedling emergence were followed up in the same manner in 2019 (rounds 3 and 4, respectively). Species identification was (re)assessed at all censuses and corrected if needed as the seedlings grew and identification uncertainty decreased. New seedlings were differentiated from emergent clonal ramets by looking for cotyledons or signs of above- or below-ground ramet connections. These data are available as talleys of seedlings, each with a status (dead or alive) and species identity (or NA when not identifiable), per subplot and /or plot per sampling round (Fig. 3, Table 1).Dataset (vii): Ecosystem carbon flux data and flux calculationsCarbon flux measurementsEcosystem CO2 fluxes were measured to estimate net ecosystem exchange (NEE), ecosystem respiration (Reco) and gross primary production (GPP). The dataset covers the years 2015, 2016 and 2017, and individual plots have multiple measurements for ecosystem carbon flux per year as detailed below. At peak growing season in 2015, a median of 2 sets of paired carbon flux measurements were measured pre-removal for all plots, where a paired set consist of a light and a dark flux measurement of an individual plot. In 2016, a median of 8 sets of paired measurements were made for all control plots, and a median of 7 for the 4 extra controls (see experimental design above). In the data files, some additional measurements exist for other experiments in the VCG sites (a median of 7 paired sets of measurements for controls (TTC) and graminoid removal plots (RTC), see42 for a presentation of this experiment and35 for technical details). In 2017, a median of 5 paired sets of measurements were made for all treated plots in nine of the sites, excluding the second wettest precipitation level (sites Gudmedalen, Rambera, and Arhelleren). These measurements were made ca. 1 week after the first round of plant functional group removals in that season.At each sampling occasion, a clear chamber (25 × 25 × 40 cm) equipped with two fans for air circulation and connected to an infrared gas analyzer (Li-840; Manufacturer reports accuracy within 1.5% of the reading value; LI-COR Biosciences, Lincoln, NE, USA) was used to measure CO2 fluxes at all plots. To prevent cutting of roots and disruption of water flow in the plots by installing collars, we instead attached a windshield to the bottom of the chamber and weighed it down on the ground by a heavy chain to prevent wind-air mixing. At each sampling occasion we made paired measurements of fluxes under light and dark conditions, covering the chamber with a fitted light-excluding cover for the dark measurements.NEE was estimated from measurements of CO2 flux under ambient light and dark conditions: NEElight = GPP – Reco, NEEdark = (-) Reco. We define NEE such that negative values reflect CO2 uptake in the ecosystem, and positive values reflect CO2 release from the ecosystem to the atmosphere. For each measurement, CO2 concentration was recorded at 5 s intervals over a period of 90–120 s. NEE was calculated from the temporal change of CO2 concentration within the closed chamber according to the following formula:$$NEE=frac{delta C{O}_{2}}{delta t}times frac{PV}{Rtimes Atimes (T+273.15)}$$where (delta frac{C{O}_{2}}{delta t}) is the slope of the CO2 concentration against time (µmol mol−1 s−1), P is the atmospheric pressure (kPa), R is the gasconstante (8.314 kPa m3 K−1 mol−1), T is the air temperature inside the chamber (°C), V is the chamber volume (m3) and A is the surface area (m2).Light intensity was measured as photosynthetically active radiation (PAR, µmol m−2 s−1) using a quantum sensor (Li-190; Manufacturer reports absolute calibration accuracy of ±5%; LI-COR Biosciences, Lincoln, NE, USA) placed inside the chamber. Temperature inside the chamber was measured using an iButton temperature logger (DS1922L, Manufacturer reports temperature accuracy of ±0.5 °C, Maxim Integrated, San Jose, CA, USA). Volumetric soil moisture content (m3 water/m3 soil) × 100 was measured by calculating the average of four measurements with a soil moisture sensor (SM300, Manufacturer reports moisture accuracy of ±2.5%, Delta-T Devices, Cambridge, UK), taken at each side of a plot.Data management and calculationsData from the LiCOR data logger and iButton was downloaded in the field and stored. The information from the field data sheets (metadata of CO2 measurements and plot soil moisture) was manually entered into digital worksheets, manually proof-read and stored. Data from the data logger (PAR and CO2) and the iButton temperature logger were linked based on information from the field data sheets. All measurements were first visually evaluated for quality and only measurements that showed a consistent linear relationship between CO2 over a time for a period of at least 60 s were used for NEE calculations. A second inclusion criterion was that this relationship had R2 ≤ 0.2 or R2 ≥ 0.8 for NEE measurement in light conditions and R2 ≥ 0.8 for NEE dark measurements (Reco). Measurements of NEE in light conditions with R2 ≤ 0.2 ensures representation of measurements with equal rates for Reco and GPP. Third, paired measurements that were more than 2 h apart were excluded. These data are available as raw fluxes and as GPP and Reco per plot per measurement (Fig. 3, Table 1).Dataset (vi): ReflectanceReflectance measures of Normalized Difference Vegetation Index (NDVI) were taken for each plot during the 2019 (post functional group removal) and 2021 (pre and post removal) field seasons (July-August), using a Trimble Greenseeker Handheld Crop Sensor (Trimble Inc., Sunnydale, CA, USA). As the sensor measures an elliptical plane, two measures perpendicular to each other were taken for each subplot (25 × 25 cm plot), with the centre of each ellipse being the centre of the subplot. Care was taken to ensure that sampling quadrat frames were not within the sensor range when conducting measurements (see methods Dataset ii). Measures of NDVI were taken at 60 cm above the surface where possible. Height was measured perpendicular to the sampled ground surface. These data are available as reflectance per plot per sampling date (Fig. 3, Table 1). More

MacLachlan, N. J. & Guthrie, A. J. Re-emergence of bluetongue, African horse sickness, and other Orbivirus diseases. Vet. Res. https://doi.org/10.1051/vetres/2010007 (2010).Article 
PubMed 
PubMed Central 

Google Scholar 
Koenraadt, C. J. M. et al. Bluetongue, Schmallenberg—What is next? Culicoides-borne viral diseases in the 21st Century. BMC Res. Notes 10, 77 (2014).
Google Scholar 
Dennis, S. J., Meyers, A. E., Hitzeroth, I. I. & Rybicki, E. P. African horse sickness: A review of current understanding and vaccine development in the. Viruses 11, 844 (2019).CAS 
PubMed Central 
Article 

Google Scholar 
Collins, Á. B., Doherty, M. L., Barrett, D. J. & Mee, J. F. Schmallenberg virus: A systematic international literature review (2011–2019) from an Irish perspective. Ir. Vet. J. 72, 1–22 (2019).Article 

Google Scholar 
Tkuwet, G. & Firesbhat, A. A review on African horse sickness. Eur. J. Appl. Sci. 7, 213–219 (2015).CAS 

Google Scholar 
Mellor, P. S. & Hamblin, C. African horse sickness. Vet. Res. 35, 445–466 (2004).PubMed 
Article 

Google Scholar 
Coetzee, P., Stokstad, M., Venter, E. H., Myrmel, M. & Van Vuuren, M. Bluetongue: A historical and epidemiological perspective with the emphasis on South Africa. Virol. J. 9, 198 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
Cagienard, A., Griot, C., Mellor, P. S., Denison, E. & Stärk, K. D. Bluetongue vector species of Culicoides in Switzerland. Med. Vet. Entomol. 20, 239–247 (2006).CAS 
PubMed 
Article 

Google Scholar 
Oluwayelu, D., Adebiyi, A. & Tomori, O. Endemic and emerging arboviral diseases of livestock in Nigeria: A review. Parasit. Vectors 11, 1–12 (2018).Article 

Google Scholar 
Sibhat, B., Ayelet, G., Gebremedhin, E. Z., Skjerve, E. & Asmare, K. Seroprevalence of Schmallenberg virus in dairy cattle in Ethiopia. Acta Trop. 178, 61–67 (2018).PubMed 
Article 

Google Scholar 
Aklilu, N. et al. African horse sickness outbreaks caused by multiple virus types in Ethiopia. Transbound. Emerg. Dis. 61, 185–192 (2014).CAS 
PubMed 
Article 

Google Scholar 
Rojas, J. M., Rodríguez-Martín, D., Martín, V. & Sevilla, N. Diagnosing bluetongue virus in domestic ruminants: Current perspectives. Vet. Med. Res. Rep. 10, 17 (2019).
Google Scholar 
Gizaw, D., Sibhat, D., Ayalew, B. & Sehal, M. Sero-prevalence study of bluetongue infection in sheep and goats in selected areas of Ethiopia. Ethiop. Vet. J. 20, 105 (2016).Article 

Google Scholar 
Abera, T. et al. Bluetongue disease in small ruminants in south western Ethiopia: Cross-sectional sero-epidemiological study. BMC Res. Notes 11, 112 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Mellor, P. S., Boorman, J. & Baylis, M. Culicoides biting midges: Their role as arbovirus vectors. Annu. Rev. Entomol. 45, 307–340 (2000).CAS 
PubMed 
Article 

Google Scholar 
Carpenter, S., Groschup, M. H., Garros, C., Felippe-Bauer, M. L. & Purse, B. V. Culicoides biting midges, arboviruses and public health in Europe. Antivir. Res. 100, 102–113 (2013).CAS 
PubMed 
Article 

Google Scholar 
Sick, F., Beer, M., Kampen, H. & Wernike, K. Culicoides biting midges—Underestimated vectors for arboviruses of public health and veterinary importance. Viruses 11, 376 (2019).CAS 
PubMed Central 
Article 

Google Scholar 
Blanda, V. et al. Geo-statistical analysis of Culicoides spp. distribution and abundance in Sicily, Italy. Parasit. Vectors 11, 78 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Vasić, A. et al. Species diversity, host preference and arbovirus detection of Culicoides (Diptera: Ceratopogonidae) in south-eastern Serbia. Parasit. Vectors 12, 1–9 (2019).Article 

Google Scholar 
Martin, E. et al. Culicoides species community composition and infection status with parasites in an urban environment of east central Texas, USA. Parasit. Vectors 12, 1–10 (2019).Article 

Google Scholar 
Gusmão, G. M. C., Brito, G. A., Moraes, L. S., Bandeira, M. D. C. A. & Rebêlo, J. M. M. Temporal variation in species abundance and richness of Culicoides (Diptera: Ceratopogonidae) in a tropical equatorial area. J. Med. Entomol. https://doi.org/10.1093/jme/tjz015 (2019).Article 
PubMed 

Google Scholar 
Sghaier, S. et al. New species of the genus Culicoides (Diptera Ceratopogonidae) for Tunisia, with detection of Bluetongue viruses in vectors. Vet. Ital. 53, 357–366 (2017).PubMed 

Google Scholar 
Gordon, S. J. G. et al. The occurrence of Culicoides species, the vectors of arboviruses, at selected trap sites in Zimbabwe. Onderstepoort J. Vet. Res. 82, e1–e8 (2015).PubMed 
Article 
CAS 

Google Scholar 
Villard, P. et al. Modeling Culicoides abundance in mainland France: Implications for surveillance. Parasit. Vectors 12, 1–10 (2019).Article 

Google Scholar 
Diarra, M. et al. Spatial distribution modelling of Culicoides (Diptera: Ceratopogonidae) biting midges, potential vectors of African horse sickness and bluetongue viruses in Senegal. Parasit. Vectors 11, 341 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Calvete, C. et al. Spatial distribution of Culicoides imicola, the main vector of bluetongue virus, Spain. Vet. Rec. 158, 130–131 (2006).CAS 
PubMed 
Article 

Google Scholar 
Purse, B. V. et al. Modelling the distributions of Culicoides bluetongue virus vectors in Sicily in relation to satellite-derived climate variables. Med. Vet. Entomol. 18, 90–101 (2004).CAS 
PubMed 
Article 

Google Scholar 
Purse, B. V. et al. Spatial and temporal distribution of bluetongue and its Culicoides vectors in Bulgaria. Med. Vet. Entomol. 20, 335–344 (2006).CAS 
PubMed 
Article 

Google Scholar 
Leta, S. et al. Modeling the global distribution of Culicoides imicola: An ensemble approach. Sci. Rep. 9, 1–9 (2019).ADS 
CAS 
Article 

Google Scholar 
Mulatu, T. & Hailu, A. The occurrence and identification of Culicoides species in the Western Ethiopia. Acad. J. Entomol. 12, 40–43 (2019).
Google Scholar 
Khamala, C. P. M. & Kettle, D. S. The Culicoides Latreille (Diptera: Ceratopogonidae) of East Africa. Trans. R. Entomol. Soc. Lond. 123, 1–95 (1971).Article 

Google Scholar 
Venter, G. J. Specie di Culicoides (Diptera: Ceratopogonidae) vettori del virus della Bluetongue in Sud Africa. Vet. Ital. 51, 325–333 (2015).PubMed 

Google Scholar 
Mathieu, B. et al. Development and validation of IIKC: An interactive identification key for Culicoides (Diptera: Ceratopogonidae) females from the Western Palaearctic region. Parasit. Vectors 5, 137 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
Thuiller, W., Lafourcade, B., Engler, R. & Araújo, M. B. BIOMOD—A platform for ensemble forecasting of species distributions. Ecography (Cop.) 32, 369–373 (2009).Article 

Google Scholar 
Baylis, M., Bouayoune, H., Touti, J. & El Hasnaoui, H. Use of climatic data and satellite imagery to model the abundance of Culicoides imicola, the vector of African horse sickness virus, in Morocco. Med. Vet. Entomol. 12, 255–266 (1998).CAS 
PubMed 
Article 

Google Scholar 
Diarra, M. et al. Modelling the abundances of two major culicoides (Diptera: Ceratopogonidae) species in the Niayes area of Senegal. PLoS One 10, e0131021 (2015).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Ramilo, D. W., Nunes, T., Madeira, S., Boinas, F. & da Fonseca, I. P. Geographical distribution of Culicoides (DIPTERA: CERATOPOGONIDAE) in mainland Portugal: Presence/absence modelling of vector and potential vector species. PLoS One 12, e0180606 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Ben Rais Lasram, F. et al. The Mediterranean Sea as a ‘cul-de-sac’ for endemic fishes facing climate change. Glob. Chang. Biol. 16, 3233–3245 (2010).ADS 
Article 

Google Scholar 
Tiffin, P. & Ross-Ibarra, J. Goal-oriented evaluation of species distribution models accuracy and precision: True Skill Statistic profile and uncertainty maps. PeerJ PrePints https://doi.org/10.7287/peerj.preprints.488v1 (2014).Article 

Google Scholar 
Graham, M. H. Confronting multicollinearity in ecological multiple regression. Ecology 84, 2809–2815 (2003).Article 

Google Scholar 
Demissie, G. H. Seroepidemiological study of African horse sickness in southern Ethiopia. Open Sci. Repos. Vet. Med. 10, e70081919 (2013).
Google Scholar 
Zeleke, A., Sori, T., Powel, K., Gebre-Ab, F. & Endebu, B. Isolation and identification of circulating serotypes of African horse sickness virus in Ethiopia. J. Appl. Res. Vet. Med. 3, 40–43 (2005).
Google Scholar 
Ayelet, G. et al. Outbreak investigation and molecular characterization of African horse sickness virus circulating in selected areas of Ethiopia. Acta Trop. 127, 91–96 (2013).PubMed 
Article 

Google Scholar 
Gulima, D. Seroepidemiological study of bluetongue in indigenous sheep in selected districts of Amhara National Regional State, north western Ethiopia. Ethiop. Vet. J. 13, 1–15 (2009).
Google Scholar 
Borkent, A. & Dominiak, P. Catalog of the biting midges of the world (Diptera: Ceratopogonidae). Zootaxa 4787, 1–377 (2020).Article 

Google Scholar 
Borkent, A. & Wirth, W. W. World species of biting midges (Diptera: Ceratopogonidae). Bull. Am. Museum Nat. Hist. 233, 5–195 (1997).
Google Scholar 
Guichard, S. et al. Worldwide niche and future potential distribution of Culicoides imicola, a major vector of bluetongue and African horse sickness viruses. PLoS One 9, e112491 (2014).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Becker, E. E. E., Venter, G. J., Labuschagne, K., Greyling, T. & van Hamburg, H. Occurrence of Culicoides species Diptera: Ceratopogonidae) in the Khomas region of Namibia during the winter months. Vet. Ital. 48, 45–54 (2012).PubMed 

Google Scholar 
Capela, R. et al. Spatial distribution of Culicoides species in Portugal in relation to the transmission of African horse sickness and bluetongue viruses. Med. Vet. Entomol. 17, 165–177 (2003).CAS 
PubMed 
Article 

Google Scholar 
Calvete, C. et al. Modelling the distributions and spatial coincidence of bluetongue vectors Culicoides imicola and the Culicoides obsoletus group throughout the Iberian peninsula. Med. Vet. Entomol. 22, 124–134 (2008).CAS 
PubMed 
Article 

Google Scholar 
Riddin, M. A., Venter, G. J., Labuschagne, K. & Villet, M. H. Culicoides species as potential vectors of African horse sickness virus in the southern regions of South Africa. Med. Vet. Entomol. 33, 498–511 (2019).CAS 
PubMed 
Article 

Google Scholar 
Foxi, C. et al. Role of different Culicoides vectors (Diptera: Ceratopogonidae) in bluetongue virus transmission and overwintering in Sardinia (Italy). Parasit. Vectors 9, 440 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Musuka, G. N., Mellor, P. S., Meiswinkel, R., Baylis, M. & Kelly, P. J. Prevalence of Culicoides imicola and other species (Diptera: Ceratopogonidae) ateight sites in Zimbabwe: To the editor. J. S. Afr. Vet. Assoc. 72, 62–63 (2001).CAS 
PubMed 
Article 

Google Scholar 
Meiswinkel, R. The 1996 outbreak of African horse sickness in South Africa—the entomological perspective. Arch. Virol. Suppl. 14, 69–83 (1998).CAS 
PubMed 

Google Scholar 
Jean Pierre, T. et al. Characteristics, classification and genesis of vertisols under seasonally contrasted climate in the Lake Chad Basin, Central Africa. J. Afr. Earth Sci. 150, 176–193 (2019).ADS 
CAS 
Article 

Google Scholar 
Elias, E. Characteristics of Nitisol profiles as affected by land use type and slope class in some Ethiopian highlands. Environ. Syst. Res. 6, 1–15 (2017).Article 

Google Scholar 
Nachtergaele, F. The classification of leptosols in the world reference base for soil resources.Veronesi, E., Venter, G. J., Labuschagne, K., Mellor, P. S. & Carpenter, S. Life-history parameters of Culicoides (Avaritia) imicola Kieffer in the laboratory at different rearing temperatures. Vet. Parasitol. 163, 370–373 (2009).CAS 
PubMed 
Article 

Google Scholar 
Verhoef, F. A. A., Venter, G. J. & Weldon, C. W. Thermal limits of two biting midges, Culicoides imicola Kieffer and C. bolitinos Meiswinkel (Diptera: Ceratopogonidae). Parasites Vectors 7, 384 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
Conte, A., Goffredo, M., Ippoliti, C. & Meiswinkel, R. Influence of biotic and abiotic factors on the distribution and abundance of Culicoides imicola and the Obsoletus Complex in Italy. Vet. Parasitol. 150, 333–344 (2007).CAS 
PubMed 
Article 

Google Scholar 
Martinez-de la Puente, J., Navarro, J., Ferraguti, M., Soriguer, R. & Figuerola, J. First molecular identification of the vertebrate hosts of Culicoides imicola in Europe and a review of its blood-feeding patterns worldwide: Implications for the transmission of bluetongue disease and African horse sickness. Med. Vet. Entomol. 31, 333–339 (2017).CAS 
PubMed 
Article 

Google Scholar 
Purse, B. V. et al. Impacts of climate, host and landscape factors on Culicoides species in Scotland. Med. Vet. Entomol. 26, 168–177 (2012).CAS 
PubMed 
Article 

Google Scholar 
Leta, S. et al. Updating the global occurrence of Culicoides imicola, a vector for emerging viral diseases. Sci. Data 6, 1–8 (2019).CAS 
Article 

Google Scholar  More

Camazine, S. et al. Self-organization in Biological Systems (Princeton University Press, Princeton, 2001).
Google Scholar 
Tschinkel, W. R. The nest architecture of the Florida harvester ant, Pogonomyrmex badius. J. Insect Sci. 4(1), 21 (2004).MathSciNet 
PubMed 
PubMed Central 
Article 

Google Scholar 
Reid, C. R. et al. Army ants dynamically adjust living bridges in response to a cost-benefit trade-off. Proc. Natl. Acad. Sci. 112(49), 15113–15118 (2015).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Grassé, P. P. Termitology. Termite anatomy-physiology-biology-systematics. Vol. II. Colony foundation-construction. Termitology. Termite anatomy-physiology-biology-systematics. Vol. II. Colony foundation-construction. Masson, Paris, (1984).Theraulaz, G., Bonabeau, E. & Deneubourg, J. L. The mechanisms and rules of coordinated building in social insects (In Information Processing in Social Insects, Birkhäuser, Basel, 1999).Hansell, M. & Hansell, M. H. Animal Architecture (Oxford University Press, Oxford, 2005).Book 

Google Scholar 
Peters, J. M., Peleg, O. & Mahadevan, L. Collective ventilation in honeybee nests. J. R. Soc. Interface 16(150), 20180561 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Grassé, P. P. La reconstruction du nid et les coordinations interindividuelles chezBellicositermes natalensis etCubitermes sp. la théorie de la stigmergie: Essai d’interprétation du comportement des termites constructeurs. Insectes Sociaux, 6(1):41–80 (1959).Theraulaz, G. & Bonabeau, E. Coordination in Distributed Building. Science 269(5224), 686–688 (1995).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Bonabeau, E., Theraulaz, G., Deneubourg, J. L. & Camazine, S. Self-organization in social insects. Trends Ecol. Evol. 12(5), 188–193 (1997).CAS 
PubMed 
Article 

Google Scholar 
Khuong, A. et al. Stigmergic construction and topochemical information shape ant nest architecture. Proc. Natl. Acad. Sci. 113(5), 1303–1308 (2016).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pénzes, Z. & Karsai, I. Round shape combs produced by Stigmergic scripts in social wasp. Proc. Eur. Conf. Artif. Life 93, 896–905 (1993).
Google Scholar 
Karsai, I. Decentralized control of construction behavior in paper wasps: an overview of the Stigmergy Approach. Artif. Life 5(2), 117–136 (1999).CAS 
PubMed 
Article 

Google Scholar 
Perna, A. & Theraulaz, G. When social behaviour is moulded in clay: On growth and form of social insect nests. J. Exp. Biol. 220(1), 83–91 (2017).PubMed 
Article 

Google Scholar 
Gallo, V. & Chittka, L. Cognitive Aspects of Comb-Building in the Honeybee?. Front. Psychol. 9, 900 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Hales, T. C. The Honeycomb Conjecture. Discrete Comput. Geom. 25(1), 1–22 (2001).MathSciNet 
MATH 
Article 

Google Scholar 
Tóth, L. F. What the bees know and what they do not know. Bull. Am. Math. Soc. 70(4), 468–481 (1964).MathSciNet 
MATH 
Article 

Google Scholar 
Jeanne, R. L. The Adaptiveness of Social Wasp Nest Architecture. Q. Rev. Biol. 50(3), 267–287 (1975).Article 

Google Scholar 
Karsai, I. & Pénzes, Z. Optimality of cell arrangement and rules of thumb of cell initiation in Polistes dominulus: A modeling approach. Behav. Ecol. 11(4), 387–395 (1999).Article 

Google Scholar 
Pirk, C., Hepburn, H., Radloff, S. & Tautz, J. Honeybee combs: construction through a liquid equilibrium process? Naturwissenschaften, 91(7) (2004).Karihaloo, B. L., Zhang, K. & Wang, J. Honeybee combs: How the circular cells transform into rounded hexagons. J. R. Soc. Interface 10(86), 20130299 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bauer, D. & Bienefeld, K. Hexagonal comb cells of honeybees are not produced via a liquid equilibrium process. Naturwissenschaften 100(1), 45–49 (2013).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Mermin, N. D. The topological theory of defects in ordered media. Rev. Mod. Phys. 51(3), 591–648 (1979).ADS 
MathSciNet 
CAS 
Article 

Google Scholar 
Bhattacharjee, S. M. Use of Topology in physical problems. In Topology and Condensed Matter Physics (eds Bhattacharjee, S. M. et al.) 171–216 (Springer, Singapore, 2017).MATH 
Chapter 

Google Scholar 
Griffin, S. M. & Spaldin, N. A. On the relationship between topological and geometric defects. J. Phys.: Condens. Matter 29(34), 343001 (2017).
Google Scholar 
Harris, W. F. Disclinations. Sci. Am. 237(6), 130–145 (1977).MathSciNet 
Article 

Google Scholar 
de Gennes, P.-G. The Physics of liquid crystals (Clarendon Press, Oxford, 1979).
Google Scholar 
Iorio, A. & Sen, S. Virus Structure: From Crick and Watson to a New Conjecture. In arXiv 0707, 3690 (2007).Lee, K. C., Yu, Q. & Erb, U. Mesostructure of Ordered Corneal Nano-nipple Arrays: The Role of 5–7 Coordination Defects. Sci. Rep. 6(1), 28342 (2016).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Stone, A. J. & Wales, D. J. Theoretical studies of icosahedral C60 and some related species. Chem. Phys. Lett. 128(5), 501–503 (1986).ADS 
CAS 
Article 

Google Scholar 
Ma, J., Alfè, D., Michaelides, A. & Wang, E. Stone-Wales defects in graphene and other planar sp2 -bonded materials. Phys. Rev. B 80(3), 033407 (2009).ADS 
Article 
CAS 

Google Scholar 
Heggie, M. I., Haffenden, G. L., Latham, C. D. & Trevethan, T. The Stone-Wales transformation: From fullerenes to graphite, from radiation damage to heat capacity. Philos.Trans. Royal Soc. A Math. Phys. Eng. Sci. 374(2076), 20150317 (2016).ADS 
Article 
CAS 

Google Scholar 
Eberhard, M. J. W. The Social Biology of Polistine Wasps. Misc. Publ. Museum Zoology Univ. Michigan 140, 110 (1969).
Google Scholar 
Jeanne, R. L. A latitudinal gradient in rates of ant predation. Ecology 60(6), 1211–1224 (1979).Article 

Google Scholar 
Seeley, T. & Heinrich, B. (1981). Regulation of temperature in the nests of social insects. John Wiley and Sons, Inc, pp. 224–234.Wenzel, J. W. Evolution of nest architecture. In The Social Biology Wasps (eds Ross, K. G. & Matthews, R. W.) 480–519 (Cornell University Press, Ithaca, New York, 1991).
Google Scholar 
Karsai, I. & Pénzes, Z. (1998). Nest shapes in paper wasps: Can the variability of forms be deduced from the same construction algorithm? Proceedings of the Royal Society of London. Series B: Biological Sciences, 265(1402):1261–1268.Carpenter, J. M. Phylogeny and biogeography of Polistes. In Natural History and Evolution of Paper-Wasps (eds Turillazzi, S. & Eberhard, M. J. W.) 18–57 (Oxford University Press, Oxford, Newyork, 1996).
Google Scholar 
Ceccolini, F. New records and distribution update of Polistes (Gyrostoma) wattii Cameron, 1900 (Hymenoptera: Vespidae: Polistinae). Caucasian Entomol. Bull. 15(2), 323–326 (2019).Article 

Google Scholar 
Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Baddeley, A., Rubak, E. & Turner, R. Spatial Point Patterns: Methodology and Applications with R (Chapman and Hall/CRC Press, London, 2015).MATH 
Book 

Google Scholar 
Steinhardt, P. J., Nelson, D. R. & Ronchetti, M. Bond-orientational order in liquids and glasses. Phys. Rev. B 28(2), 784–805 (1983).ADS 
CAS 
Article 

Google Scholar 
Schilling, T., Pronk, S., Mulder, B. & Frenkel, D. Monte Carlo study of hard pentagons. Phys. Rev. E 71(3), 036138 (2005).ADS 
Article 
CAS 

Google Scholar 
R Core Team. R: A language and environment for statistical computing. https://www.R-project.org/ (2020).Bishop, M. & Bruin, C. The pair correlation function: A probe of molecular order. Am. J. Phys. 52(12), 1106–1108 (1984).ADS 
CAS 
Article 

Google Scholar 
Fleury, P. A. Phase Transitions, Critical Phenomena, and Instabilities. Science 211, 125–131 (1981).ADS 
MathSciNet 
CAS 
PubMed 
MATH 
Article 

Google Scholar 
Wenzel, J. W. Endogenous factors, external cues, and eccentric construction in Polistes annularis (Hymenoptera: Vespidae). J. Insect Behavior 2(5), 679–699 (1989).Article 

Google Scholar 
Zsoldos. Effect of topological defects on graphene geometry and stability. Nanotechnol. Sci. Appl., p. 101 (2010).Ophus, C., Shekhawat, A., Rasool, H. & Zettl, A. Large-scale experimental and theoretical study of graphene grain boundary structures. Phys. Rev. B 92(20), 205402 (2015).ADS 
Article 
CAS 

Google Scholar 
Kosterlitz, J. M. (2016). Commentary on ‘Ordering, metastability and phase transitions in two-dimensional systems’ J M Kosterlitz and D J Thouless (1973 J. Phys. C: Solid State Phys. 6 1181-203)-the early basis of the successful Kosterlitz-Thouless theory. Journal of Physics: Condensed Matter28:481001.Hepburn, H. R. & Whiffler, L. A. Construction defects define pattern and method in comb building by honeybees. Apidologie 22(4), 381–388 (1991).Article 

Google Scholar 
Smith, M. L., Napp, N. & Petersen, K. H. Imperfect comb construction reveals the architectural abilities of honeybees. Proc. Natl. Acad. Sci. 118(31), e2103605118 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Nazzi, F. The hexagonal shape of the honeycomb cells depends on the construction behavior of bees. Sci. Rep. 6(1), 28341 (2016).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tarnai, T. Buckling patterns of shells and spherical honeycomb structures. Symmetry, pp. 639–652 (1989).Downing, H. & Jeanne, R. The regulation of complex building behaviour in the paper wasp, Polistes fuscatus (Insecta, Hymenoptera, Vespidae). Anim. Behav. 39(1), 105–124 (1990).Article 

Google Scholar  More

Naef-Daenzer, B. Patch time allocation and patch sampling by foraging great and blue tits. Anim. Behav. 59, 989–999 (2000).CAS 
PubMed 
Article 

Google Scholar 
Kotler, B. P., Brown, J. S. & Bouskila, A. Apprehension and time allocation in gerbils: The effects of predatory risk and energetic state. Ecology 85, 917–922 (2004).Article 

Google Scholar 
Wajnberg, E., Bernhard, P., Hamelin, F. & Boivin, G. Optimal patch time allocation for time-limited foragers. Behav. Ecol. Sociobiol. 60, 1–10 (2006).Article 

Google Scholar 
Embar, K., Kotler, B. P. & Mukherjee, S. Risk management in optimal foragers: The effect of sightlines and predator type on patch use, time allocation, and vigilance in gerbils. Oikos 120, 1657–1666 (2011).Article 

Google Scholar 
Lima, S. L. & Bednekoff, P. A. Temporal variation in danger drives antipredator behavior: The predation risk allocation hypothesis. Am. Nat. 153, 649–659 (1999).PubMed 
Article 

Google Scholar 
Beauchamp, G. & Ruxton, G. D. A reassessment of the predation risk allocation hypothesis: A comment on Lima and Bednekoff. Am. Nat. 177, 143–146 (2011).PubMed 
Article 

Google Scholar 
Ferrari, M. C. O., Sih, A. & Chivers, D. P. The paradox of risk allocation: A review and prospectus. Anim. Behav. 78, 579–585 (2009).Article 

Google Scholar 
Wolf, L. L. & Hainsworth, F. R. Foraging efficiencies and time budgets in nectar-feeding birds. Ecology 56, 117–128 (1975).Article 

Google Scholar 
Litzow, M. A. & Piatt, J. F. Variance in prey abundance influences time budgets of breeding seabirds: Evidence from pigeon guillemots Cepphus columba. J. Avian Biol. 34, 54–64 (2003).Article 

Google Scholar 
Rishworth, G. M., Tremblay, Y. & Green, D. B. Drivers of time-activity budget variability during breeding in a pelagic seabird. PLoS One 9, e116544 (2014).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Stephens, D. W., Brown, J. S. & Ydenberg, R. C. Foraging: Behavior and Ecology. (The University of Chicago Press, 2007).Orians, G. & Pearson, N. On the theory of central place foraging. In Analysis of Ecological Systems (eds. Horn, D., Mitchell, R. & Stairs, G.) 154–177 (The Ohio State University Press, 1979).Chaurand, T. & Weimerskirch, H. The regular alternation of short and long foraging trips in the blue petrel Halobaena caerulea: A previously undescribed strategy of food provisioning in a pelagic seabird. J. Anim. Ecol. 63, 275–282 (1994).Article 

Google Scholar 
Weimerskirch, H. et al. Alternate long and short foraging trips in pelagic seabird parents. Anim. Behav. 47, 472–476 (1994).Article 

Google Scholar 
Welcker, J., Beiersdorf, A., Varpe, Ø. & Steen, H. Mass fluctuations suggest different functions of bimodal foraging trips in a central-place forager. Behav. Ecol. 23, 1372–1378 (2012).Article 

Google Scholar 
Welcker, J. et al. Flexibility in the bimodal foraging strategy of a high Arctic alcid, the little auk Alle alle. J. Avian Biol. 40, 388–399 (2009).Article 

Google Scholar 
Jakubas, D., Wojczulanis-Jakubas, K., Iliszko, L. M. & Kidawa, D. Flexibility of little auks foraging in various oceanographic features in a changing Arctic. Sci. Rep. https://doi.org/10.1038/s41598-020-65210-x (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Shoji, A. et al. Dual foraging and pair coordination during chick provisioning by Manx shearwaters: Empirical evidence supported by a simple model. J. Exp. Biol. 218, 2116–2123 (2015).PubMed 
PubMed Central 

Google Scholar 
Phillips, R. A., Wakefield, E. D., Croxall, J. P., Fukuda, A. & Higuchi, H. Albatross foraging behaviour: No evidence for dual foraging, and limited support for anticipatory regulation of provisioning at South Georgia. Mar. Ecol. Prog. Ser. 391, 279–292 (2009).ADS 
Article 

Google Scholar 
Brown, Z. W., Welcker, J., Harding, A. M. A., Walkusz, W. & Karnovsky, N. J. Divergent diving behavior during short and long trips of a bimodal forager, the little auk Alle alle. J. Avian Biol. 43, 215–226 (2012).Article 

Google Scholar 
Baduini, C. L. & Hyrenbach, K. D. Biogeography of procellariiform foraging strategies: Does ocean productivity influence provisioning?. Mar. Ornithol. 31, 101–112 (2003).
Google Scholar 
Navarro, J. & González-Solís, J. Environmental determinants of foraging strategies in Cory’s shearwaters Calonectris diomedea. Mar. Ecol. Prog. Ser. 378, 259–267 (2009).ADS 
CAS 
Article 

Google Scholar 
Ochi, D., Oka, N. & Watanuki, Y. Foraging trip decisions by the streaked shearwater Calonectris leucomelas depend on both parental and chick state. J. Ethol. 28, 313–321 (2010).Article 

Google Scholar 
Congdon, B. C., Krockenberger, A. K. & Smithers, B. V. Dual-foraging and co-ordinated provisioning in a tropical Procellariiform, the wedge-tailed shearwater. Mar. Ecol. Prog. Ser. 301, 293–301 (2005).ADS 
Article 

Google Scholar 
Peck, D. R. & Congdon, B. C. Colony-specific foraging behaviour and co-ordinated divergence of chick development in the wedge-tailed shearwater Puffinus pacificus. Mar. Ecol. Prog. Ser. 299, 289–296 (2005).ADS 
Article 

Google Scholar 
Weimerskirch, H. How can a pelagic seabird provision its chick when relying on a distant food resource? Cyclic attendance at the colony, foraging decision and body condition in sooty shearwaters. J. Anim. Ecol. 67, 99–109 (1998).Article 

Google Scholar 
Stempniewicz, L. BWP update. Little Auk (Alle alle). J. Birds West. Palearct. 3, 175–201 (2001).
Google Scholar 
Wojczulanis-Jakubas, K. & Jakubas, D. When and why does my mother leave me? The question of brood desertion in the Dovekie (Alle Alle). Auk 129, 632–637 (2012).Article 

Google Scholar 
Harding, A. M. A., Van Pelt, T. I., Lifjeld, J. T. & Mehlum, F. Sex differences in little auk Alle alle parental care: Transition from biparental to paternal-only care. Ibis (Lond. 1859). 146, 642–651 (2004).Article 

Google Scholar 
Wojczulanis-Jakubas, K. et al. Duration of female parental care and their survival in the little auk Alle alle—Are these two traits linked ?. Behav. Ecol. Sociobiol. 74, 1–11 (2020).Article 

Google Scholar 
Wojczulanis, K., Dariusz, J. & Lech, S. The Little Auk Alle alle: An ecological indicator of a changing Arctic and a model organism. Polar Biol. https://doi.org/10.1007/s00300-021-02981-7 (2021).Article 

Google Scholar 
Steen, H., Vogedes, D., Broms, F., Falk-Petersen, S. & Berge, J. Little auks (Alle alle) breeding in a High Arctic fjord system: Bimodal foraging strategies as a response to poor food quality?. Polar Res. 26, 118–125 (2007).Article 

Google Scholar 
Wojczulanis-Jakubas, K., Jakubas, D., Karnovsky, N. J. & Walkusz, W. Foraging strategy of little auks under divergent conditions on feeding grounds. Polar Res. 29, 22–29 (2010).Article 

Google Scholar 
Jakubas, D., Wojczulanis-Jakubas, K., Iliszko, L., Darecki, M. & Stempniewicz, L. Foraging strategy of the little auk Alle alle throughout breeding season—switch from unimodal to bimodal pattern. J. Avian Biol. 45, 551–560 (2014).Article 

Google Scholar 
Jakubas, D., Iliszko, L., Wojczulanis-Jakubas, K. & Stempniewicz, L. Foraging by little auks in the distant marginal sea ice zone during the chick-rearing period. Polar Biol. 35, 73–81 (2012).Article 

Google Scholar 
Jakubas, D. et al. Intra-seasonal variation in zooplankton availability, chick diet and breeding performance of a high Arctic planktivorous seabird. Polar Biol. 391, 1547–1561 (2016).Article 

Google Scholar 
Jakubas, D. et al. Foraging closer to the colony leads to faster growth in little auks. Mar. Ecol. Prog. Ser. 489, 263–278 (2013).ADS 
Article 

Google Scholar 
Stempniewicz, L. Predator-prey interactions between Glaucous Gull Larus hyperboreus and Little Auk Alle alle in Spitsbergen. Acta Ornithol. 29, 155–170 (1995).
Google Scholar 
Wojczulanis-Jakubas, K., Jakubas, D. & Stempniewicz, L. Changes in the glaucous gull predatory pressure on little auks in Southwest Spitsbergen. Waterbirds 28, 430–435 (2005).Article 

Google Scholar 
Kharitonov, S. Methods and Theoretical Aspects of Seabird Studies. (Proc 5 All-Russian Mar Biol School, Marine Biological Institute, 2007).Wojczulanis-Jakubas, K., Jakubas, D. & Stempniewicz, L. Avifauna of Hornsund area, SW Spitsbergen: Present state and recent changes. Polish Polar Res. 29, 187–197 (2008).
Google Scholar 
Keslinka, K. L., Wojczulanis-Jakubas, K., Jakubas, D. & Neubauer, G. Determinants of the little auk (Alle alle) breeding colony location and size in W and NW coast of Spitsbergen. PLoS One 14, 1–20 (2019).
Google Scholar 
Kidawa, D., Barcikowski, M. & Palme, R. Parent-offspring interactions in a long-lived seabird, the Little Auk (Alle alle): Begging and provisioning under simulated stress. J. Ornithol. 158, 145–157 (2017).Article 

Google Scholar 
Welcker, J., Beiersdorf, A., Varpe, Ø. & Steen, H. Mass fluctuations suggest different functions of bimodal foraging trips in a central-place forager. Behav. Ecol. https://doi.org/10.1093/beheco/ars131 (2012).Article 

Google Scholar 
Jakubas, D. & Wojczulanis, K. Predicting the sex of Dovekies by discriminant analysis. Waterbirds 30, 92–96 (2007).Article 

Google Scholar 
Grissot, A. et al. Parental coordination of chick provisioning in a planktivorous arctic seabird under divergent conditions on foraging grounds. Front. Ecol. Evol. 7, 349 (2019).Article 

Google Scholar 
Stoffel, M. A., Nakagawa, S. & Schielzeth, H. rptR: Repeatability estimation and variance decomposition by generalized linear mixed-effects models. Methods Ecol. Evol. 8, 1639–1644 (2017).Article 

Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing. R. (2019).Wojczulanis-Jakubas, K., Jakubas, D. & Stempniewicz, L. Sex-specific parental care by incubating Little Auks (Alle alle). Ornis Fenn. 86, 140–148 (2009).
Google Scholar 
Welcker, J., Steen, H., Harding, A. M. A. & Gabrielsen, G. W. Sex-specific provisioning behaviour in a monomorphic seabird with a bimodal foraging strategy. Ibis (Lond. 1859). 151, 502–513 (2009).Article 

Google Scholar 
Kidawa, D. et al. Parental efforts of an Arctic seabird, the little auk Alle alle under variable foraging conditions. Mar. Biol. Res. 11, 349–360 (2015).Article 

Google Scholar 
Wickham, H. Hadley Wickham. Media 35, 211 (2009).
Google Scholar 
Karnovsky, N. J. et al. Inter-colony comparison of diving behavior of an Arctic top predator: Implications for warming in the Greenland Sea. Mar. Ecol. Prog. Ser. 440, 229–240 (2011).ADS 
Article 

Google Scholar 
Karnovsky, N. et al. Foraging distributions of little auks Alle alle across the Greenland Sea: Implications of present and future Arctic climate change. Mar. Ecol. Prog. Ser. 415, 283–293 (2010).ADS 
Article 

Google Scholar 
Gremillet, D. et al. Little auks buffer the impact of current Arctic climate change. Mar. Ecol. Prog. Ser. 454, 197–206 (2012).ADS 
Article 

Google Scholar 
Harding, A. M. A. et al. Flexibility in the parental effort of an Arctic-breeding seabird. Funct. Ecol. 23, 348–358 (2009).Article 

Google Scholar 
Jakubas, D. et al. Foraging effort does not influence body condition and stress level in little auks. Mar. Ecol. Prog. Ser. 432, 277–290 (2011).ADS 
Article 

Google Scholar 
Jakubas, D., Wojczulanis-Jakubas, K., Iliszko, L. M., Strøm, H. & Stempniewicz, L. Habitat foraging niche of a High Arctic zooplanktivorous seabird in a changing environment. Sci. Rep. 7, 1–14 (2017).CAS 
Article 

Google Scholar  More

Rotational velocityFigure 3 shows the effect of flow rate, nozzle diameter and number of nozzles on the rotational velocity of water in a circular tank. The results indicate that the rotational velocity increases with increasing flow rates and deceasing nozzle diameter. It could be seen that, the rotational velocity decreased from 10.1 to 5.0 cm s−1, when the nozzle diameter increased from 10 to 20 mm, respectively for 5 nozzles used, and it decreased from 5.1 to 4.0 cm s−1, when the nozzle diameter increased from 10 to 15 mm, respectively, for 10 nozzles used with 5 m3 h−1 flow rate. At 15 m3 h−1, the rotational velocity was decreased from 23.5 to 17.5, 12.0 to 7.5, 10.0 to 6.9, 7.6 to 4.7 and 5.9 to 4.0 cm s−1 when the nozzle diameter increased from 10 to 20 mm, respectively, for 5, 10, 15, 20 and 25 nozzles, respectively. The results also indicate that when the nozzle diameter increased from 20 to 25 mm, the rotational velocity decreased from 19.0 to 16.5, 12.0 to 10.0 and 7.1 to 5.5 cm s−1 for 3, 6 and 9 nozzles, respectively, with 15 m3 h−1 flow rate.Figure 3Effect of flow rate, nozzle diameter and number of nozzles on the rotational velocity of water in a circular tank.Full size imageAt 30 m3 h−1 flow rate, the highest value of the rotational velocity was 33.5 cm s−1 was found for 5 nozzles and 10 mm nozzle diameter. While, the lowest value of the rotational velocity was 7.3 cm s−1 was found for 25 nozzles and 25 mm nozzle diameter. At 45 m3 h−1 flow rate, the rotational velocity ranged from 11.0 to 49.9 cm s−1 for all treatments under study.At 60 m3 h−1 flow rate, the rotational velocity deceased from 61.0 to 50.1, 47.7 to 34.0, 36.3 to 23.0, 23.5 to 17.5, 21.0 to 15.0 and 17.0 to 11.5 cm s−1 when the nozzle diameter increased from 10 to 20 mm, respectively at 5, 10, 15, 20, 25 and 30 number of nozzles. The results also indicate that, when the nozzle diameter increased from 20 to 25 mm, the rotational velocity decreased from 56.0 to 47.0, 43.0 to 33.0, 27.0 to 22.0 and 19.0 to 16.5 cm s−1 at 3, 6, 9 and 12 nozzles, respectively.At 75 m3 h−1 flow rate, the rotational velocity deceased from 60.9 to 49.1, 48.4 to 38.0, 39.0 to 30, 31.8 to 23.0, 23.5 to 17.5 and 22.0 to 15.0 cm s−1 when the nozzle diameter increased from 10 to 20 mm, respectively for 5, 10, 15, 20, 25 and 30 nozzles, respectively. The results also indicate that, when the nozzle diameter increased from 20 to 25 mm, the rotational velocity decreased from 50.48 to 43.0 to 38.5, 33.0 to 27.5 and 23.5 to 22.0 cm s−1 for 3, 6, 9 and 12 nozzles, respectively.The results also indicate that the highest values of the rotational velocities were 10.1, 23.5, 33.5, 49.9, 60.9 and 61.0 cm s−1 were found for 5 nozzles and 10 mm nozzle diameter at 5, 15, 30, 45, 60 and 75 m3 h−1 flow rate, respectively. While, the lowest values of the rotational velocities were 4.0, 7.5 and 11.5 cm s−1 for 25 nozzles and 15 mm nozzle diameter at 5, 15 and 30 m3 h−1 flow rate, respectively. They were 11.5 and 15 cm s−1 were found for 30 nozzles and 15 mm nozzle diameter at 60 and 75 m3 h−1 flow rate, respectively. The velocity of water obtained seemed to be in the recommended range of safe and proper velocity for fish according to12. Due to it is effective compromise to allow heavy solids settle rapidly, yet sufficiently fast to create “good” hydraulics. Timmons and Youngs18 mentioned that the water velocity needed to maintain self-cleaning properties ranges from 3 to 40 cm s−1 varying greatly according to the physical properties of the biosolids. When fish swims at lower speed than its optimal, a large amount of energy will be used for higher spontaneous activity such as aggression. In contrast, when fish swim at higher speed than optimal, they become stressful, unstable, increase lactate production and fatigue6.Multiple regression analysis was carried out to obtain a relationship between the rotational velocity of water as dependent variable and different both of flow rate and nozzle diameter as independent variables. The best fit for this relationship with coefficient of determination of 0.95 and an error of 1.06% is in the following form:-$$ RV = 6.97 + 0.41Q – 0.19Dquad {text{R}}^{{2}} = 0.95 $$
(3)
where RV is the rotational velocity of water, cm s−1, Q is the water flow rate, m3 h−1, D is the nozzle diameter, mm.This equation could be applied in the range of 5 to 75 m3 h−1 water flow rate and from 10 to 25 mm of nozzle diameter.Impulse force of waterFigure 4 shows the effect of flow rate, diameter and number of nozzles on the impulse force of water in a circular tank. The results indicate that the impulse force of water increases with increasing flow rates and deceasing nozzle diameter and number of nozzles. It could be seen that, the impulse force of water decreased from 5.1 to 1.7 N, when the number of nozzles increased from 5 to 15, respectively at 10 nozzle diameter, and it decreased from 2.3 to 1.2 N, when the number of nozzles increased from 5 to 10, respectively, at 15 diameter nozzle with 5 m3 h−1 flow rate. At 15 m3 h−1, the impulse force of water was decreased from 84.7 to 9.4 N when the number of nozzles increased from 5 to 30, respectively 10 mm diameter nozzle. The results also indicate that when the number of nozzles increased from 5 to 25, the impulse force of water decreased from 14.8 to 1.4 N at 15 mm nozzle diameter, respectively, and it decreased from 9.5 to 1.9 and 5.3 to 1.3 N at 20 and 25 mm, respectively, when the number of nozzles increased from 3 to 9.Figure 4Effect of flow rate, nozzle diameter and number of nozzles on the impulse force of water in a circular tank.Full size imageAt 30 m3 h−1 flow rate, the impulse force of water deceased from 84.7 to 46.9, 56.9 to 14.8, 28.5 to 5.3, 14.9 to 3.0 and 11.8 to 2.2 N when the nozzle diameter increased from 10 to 15 mm, respectively at 5, 10, 15, 20 and 25 nozzles. The results also indicate that, when the nozzle diameter increased from 20 to 25 mm, the impulse force of water decreased from 21.4 to 14.9, 14.8 to 5.4, 5.3 to 2.2 and 2.3 to 1.9 N for 3, 6, 9 and 12 nozzles, respectively.At 45 m3 h−1 flow rate, the impulse force of water was ranged from 2.1 to 111.2 N for all treatments under this study. Also, at 60 m3 h−1 flow rate, the impulse force of water ranged from 5.1 to 151.3 N for all treatments under this study. At 75 m3 h−1 flow rate, the highest value of the impulse force of water 211.2 N was found for 5 numbers of nozzles and 10 mm nozzle diameter, respectively. While, the lowest value of the impulse force of water was 9.1 N was found for 12 nozzles and 25 mm nozzle diameter, respectively.The results also indicate that the highest value of the impulse force of water 211.2 N was found for 5 nozzles and 10 mm nozzle diameter at 75 m3 h−1 flow rate, respectively. While, the lowest value of the impulse force of water was 1.2 N was found for 10 nozzles and 15 mm nozzle diameter at 5 m3 h−1 flow rate, respectively.The results indicated that, the relationship between the rotational velocity and impulse force of water is linear relationship at the same treatments. When the rotational velocity increased from 10.7 to 37.6, 8.1 to 28.8, 10.2 to 36.0 and 11.0 to 31.9 cm s−1, the impulse force of water increased from 3.1 to 106.6, 1.8 to 31.1, 1.3 to 32.5 and 1.4 to 22.8 N, respectively, at the same treatments. The trend of these results agreed with those obtained by19.Multiple regression analysis was carried out to obtain a relationship between the impulse force of water as dependent variable and different both of flow rate and nozzle diameter as independent variables. The best fit for this relationship with coefficient of determination of 0.88 and an error of 2.13% is in the following form:-$$ F_{i} = 38.18 + 0.67Q – 2.35Dquad {text{R}}^{{2}} = 0.88 $$
(4)
This equation could be applied in the range of 5 to 75 m3 h−1 water flow rate and from 10 to 25 mm of nozzle diameter.Average velocity of waterFigure 5 shows the effect of flow rate, diameter and number of nozzles on the average velocity of water in a circular tank. The results indicate that the average velocity of water increases with increasing flow rates and deceasing nozzle diameter and number of nozzles. It could be seen that, the average velocity of water decreased from 3.32 to 1.59 cm s−1, when the number of nozzles increased from 5 to 15, respectively at 10 nozzle diameter, and it decreased from 1.13 to 1.07 cm s−1, when the number of nozzles increased from 5 to 10, respectively, at 15 diameter nozzle with 5 m3 h−1 flow rate. At 15 m3 h−1, the average velocity of water was decreased from 12.03 to 4.33 cm s−1 when the number of nozzles increased from 5 to 30, respectively 10 mm diameter nozzle. The results also indicate that when the number of nozzles increased from 5 to 25, the average velocity of water decreased from 6.93 to 2.89 cm s−1 at 15 mm nozzle diameter, respectively, and it decreased from 7.55 to 4.00 and 4.89 to 2.95 cm s−1 at 20 and 25 mm, respectively, when the number of nozzles increased from 3 to 9.Figure 5Effect of flow rate, nozzle diameter and number of nozzles on the average velocity of water in a circular tank.Full size imageAt 30 m3 h−1 flow rate, the highest value of the average velocity of water 18.51 cm s−1 was found for 5 nozzles and 10 mm nozzle diameter. While, the lowest value of the average velocity of water was 4.65 cm s−1 was found for 12 nozzles and 25 mm nozzle diameter. At 45 m3 h−1 flow rate, the average velocity of water ranged from 6.66 to 23.26 for all treatments under study, also, at 60 m3 h−1 flow rate, the average velocity of water ranged from 9.23 to 34.82 for all treatments under study. At 75 m3 h−1 flow rate, the average velocity of water ranged from 10.00 to 48.76 for all treatment of this study.The results also indicate that the highest value of the average velocity of water 48.76 cm s−1 was found for 5 nozzles and 10 mm nozzle diameter at 75 m3 h−1 flow rate, respectively. While, the lowest value of the average velocity of water was 1.07 cm s−1 was found for 10 nozzles and 15 mm nozzle diameter at 5 m3 h−1 flow rate, respectively. These results agreed with those obtained by18,20. Fish distribution in the circular tank is influenced by the heterogeneity of water velocity in the area between inlet flow and the center of the tank9. Fish distribution in the circular tank is mostly concentrated in the area between high and low velocity area. The high velocity area will be avoided by most fishes as it requires high swimming energy, while dead volumes (low velocity area) are unfavorable condition for fish (low DO and higher metabolites accumulation)21.Multiple regression analysis was carried out to obtain a relationship between the average velocity of water as dependent variable and different both of flow rate and nozzle diameter as independent variables. The best fit for this relationship with coefficient of determination of 0.91 and an error of 1.48% is in the following form:$$ V_{avg} = 6.53 + 0.26Q – 0.37Dquad {text{R}}^{{2}} = 0.91 $$
(5)
This equation could be applied in the range of 5 to 75 m3 h−1 water flow rate and from 10 to 25 mm of nozzle diameter. More