Philippa Kaur

Ceballos, G. et al. Accelerated modern human-induced species losses: entering the sixth mass extinction. Sci. Adv. 1, e1400253 (2015).Article 

Google Scholar 
Dereniowska, M. & Meinard, Y. The unknownness of biodiversity: its value and ethical significance for conservation action. Biol. Conserv. 260, 109199 (2021).Article 

Google Scholar 
Maron, M. et al. Towards a threat assessment framework for ecosystem services. Trends Ecol. Evol. 32, 240–248 (2017).Article 

Google Scholar 
Tilman, D. et al. Future threats to biodiversity and pathways to their prevention. Nature 546, 73–81 (2017).Article 
CAS 

Google Scholar 
Taborsky, B. et al. Towards an evolutionary theory of stress responses. Trends Ecol. Evol. 36, 39–48 (2021).Article 

Google Scholar 
van de Leemput, I. A., Dakos, V., Scheffer, M. & van Nes, E. H. Slow recovery from local disturbances as an indicator for loss of ecosystem resilience. Ecosystems 21, 141–152 (2018).Article 

Google Scholar 
Fagan, W. F. & Holmes, E. E. Quantifying the extinction vortex. Ecol. Lett. 9, 51–60 (2005).
Google Scholar 
Williams, N. F., McRae, L., Freeman, R., Capdevila, P. & Clements, C. F. Scaling the extinction vortex: body size as a predictor of population dynamics close to extinction events. Ecol. Evol. 11, 7069–7079 (2021).Article 

Google Scholar 
Clements, C. F. & Ozgul, A. Indicators of transitions in biological systems. Ecol. Lett. 21, 905–919 (2018).Article 

Google Scholar 
Shaffer, M. L. in Challenges in the Conservation of Biological Resources (eds. Decker, D. J., Krasny, M. E., Goff, G. R., Smith, C. R. & Gross, D. W.) 107–118 (Routledge, 2019).Scheffer, M. et al. Early-warning signals for critical transitions. Nature 461, 53–59 (2009).Article 
CAS 

Google Scholar 
Gardner, T. A. et al. The cost-effectiveness of biodiversity surveys in tropical forests. Ecol. Lett. 11, 139–150 (2008).Article 

Google Scholar 
Coulson, T., Mace, G. M., Hudson, E. & Possingham, H. The use and abuse of population viability analysis. Trends Ecol. Evol. 16, 219–221 (2001).Article 
CAS 

Google Scholar 
Clements, C. F., Drake, J. M., Griffiths, J. I. & Ozgul, A. Factors influencing the detectability of early warning signals of population collapse. Am. Nat. 186, 50–58 (2015).Article 

Google Scholar 
Patterson, A. C., Strang, A. G. & Abbott, K. C. When and where we can expect to see early warning signals in multispecies systems approaching tipping points: insights from theory. Am. Nat. 198, E12–E26 (2021).Article 

Google Scholar 
Vinton, A. C., Gascoigne, S. J. L., Sepil, I. & Salguero-Gómez, R. Plasticity’s role in adaptive evolution depends on environmental change components. Trends Ecol. Evol. https://doi.org/10.1016/j.tree.2022.08.008 (2022).Levin, S. A. The problem of pattern and scale in ecology: the Robert H. MacArthur Award lecture. Ecology 73, 1943–1967 (1992).Article 

Google Scholar 
Brown, J. H., Gillooly, J. F., Allen, A. P., Savage, V. M. & West, G. B. Toward a metabolic theory of ecology. Ecology 85, 1771–1789 (2004).Article 

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).Article 
CAS 

Google Scholar 
Haberle, I., Marn, N., Geček, S. & Klanjšček, T. Dynamic energy budget of endemic and critically endangered bivalve Pinna nobilis: a mechanistic model for informed conservation. Ecol. Model. 434, 109207 (2020).Article 

Google Scholar 
Gislason, H., Daan, N., Rice, J. C. & Pope, J. G. Size, growth, temperature and the natural mortality of marine fish. Fish Fish. 11, 149–158 (2010).Article 

Google Scholar 
Jennings, S. & Blanchard, J. L. Fish abundance with no fishing: predictions based on macroecological theory. J. Anim. Ecol. 73, 632–642 (2004).Article 

Google Scholar 
Valderrama, D. & Fields, K. H. Flawed evidence supporting the metabolic theory of ecology may undermine goals of ecosystem-based fishery management: the case of invasive Indo-Pacific lionfish in the western Atlantic. ICES J. Mar. Sci. 74, 1256–1267 (2017).Article 

Google Scholar 
Marshall, D. J. & McQuaid, C. D. Warming reduces metabolic rate in marine snails: adaptation to fluctuating high temperatures challenges the metabolic theory of ecology. Proc. R. Soc. B 278, 281–288 (2011).Article 

Google Scholar 
Rombouts, I., Beaugrand, G., Ibaňez, F., Chiba, S. & Legendre, L. Marine copepod diversity patterns and the metabolic theory of ecology. Oecologia 166, 349–355 (2011).Article 

Google Scholar 
Allen, A. P. & Gillooly, J. F. The mechanistic basis of the metabolic theory of ecology. Oikos 116, 1073–1077 (2022).Article 

Google Scholar 
Lawton, J. H. From physiology to population dynamics and communities. Funct. Ecol. 5, 155–161 (1991).Article 

Google Scholar 
Ames, E. M. et al. Striving for population-level conservation: integrating physiology across the biological hierarchy. Conserv. Physiol. 8, coaa019 (2020).Article 

Google Scholar 
Berger-Tal, O. et al. Integrating animal behavior and conservation biology: a conceptual framework. Behav. Ecol. 22, 236–239 (2011).Article 

Google Scholar 
Baruah, G., Clements, C. F., Guillaume, F. & Ozgul, A. When do shifts in trait dynamics precede population declines? Am. Nat. 193, 633–644 (2019).Article 

Google Scholar 
Dakos, V. et al. Methods for detecting early warnings of critical transitions in time series illustrated using simulated ecological data. PLoS ONE 7, e41010 (2012).Article 
CAS 

Google Scholar 
Ward, R. J., Griffiths, R. A., Wilkinson, J. W. & Cornish, N. Optimising monitoring efforts for secretive snakes: a comparison of occupancy and N-mixture models for assessment of population status. Sci. Rep. 7, 18074 (2017).Article 

Google Scholar 
Thompson, W. Sampling Rare or Elusive Species: Concepts, Designs, and Techniques for Estimating Population Parameters (Island Press, 2013).Clements, C. F., Blanchard, J. L., Nash, K. L., Hindell, M. A. & Ozgul, A. Body size shifts and early warning signals precede the historic collapse of whale stocks. Nat. Ecol. Evol. 1, 0188 (2017).Article 

Google Scholar 
Burant, J. B., Park, C., Betini, G. S. & Norris, D. R. Early warning indicators of population collapse in a seasonal environment. J. Anim. Ecol. 90, 1538–1549 (2021).Article 

Google Scholar 
Tuomainen, U. & Candolin, U. Behavioural responses to human-induced environmental change. Biol. Rev. 86, 640–657 (2011).Article 

Google Scholar 
Mazza, V., Dammhahn, M., Lösche, E. & Eccard, J. A. Small mammals in the big city: behavioural adjustments of non-commensal rodents to urban environments. Glob. Change Biol. 26, 6326–6337 (2020).Article 

Google Scholar 
Hendry, A. P., Farrugia, T. J. & Kinnison, M. T. Human influences on rates of phenotypic change in wild animal populations. Mol. Ecol. 17, 20–29 (2008).Article 

Google Scholar 
Speakman, J. R., Król, E. & Johnson, M. S. The functional significance of individual variation in basal metabolic rate. Physiol. Biochem. Zool. 77, 900–915 (2004).Article 

Google Scholar 
Péron, G. et al. Evidence of reduced individual heterogeneity in adult survival of long-lived species. Evolution 70, 2909–2914 (2016).Article 

Google Scholar 
Fleming, A. H., Clark, C. T., Calambokidis, J. & Barlow, J. Humpback whale diets respond to variance in ocean climate and ecosystem conditions in the California Current. Glob. Change Biol. 22, 1214–1224 (2016).Article 

Google Scholar 
Kirkwood, T. B. L., Rose, M. R., Harvey, P. H., Partridge, L. & Southwood, S. R. Evolution of senescence: late survival sacrificed for reproduction. Phil. Trans. R. Soc. Lond. B 332, 15–24 (1991).Article 
CAS 

Google Scholar 
Mallela, A. & Hastings, A. The role of stochasticity in noise-induced tipping point cascades: a master equation approach. Bull. Math. Biol. 83, 53 (2021).Article 

Google Scholar 
Burthe, S. J. et al. Do early warning indicators consistently predict nonlinear change in long-term ecological data? J. Appl. Ecol. 53, 666–676 (2016).Article 

Google Scholar 
Vucetich, J. A. & Waite, T. A. Erosion of heterozygosity in fluctuating populations. Conserv. Biol. 13, 860–868 (1999).Article 

Google Scholar 
Kramer, A. M. & Drake, J. M. Experimental demonstration of population extinction due to a predator-driven Allee effect. J. Anim. Ecol. 79, 633–639 (2010).Article 

Google Scholar 
Oram, E. & Spitze, K. Depth selection by Daphnia pulex in response to Chaoborus kairomone. Freshw. Biol. 58, 409–415 (2013).Article 

Google Scholar 
Trites, A. W. & Donnelly, C. P. The decline of Steller sea lions Eumetopias jubatus in Alaska: a review of the nutritional stress hypothesis. Mammal. Rev. 33, 3–28 (2003).Article 

Google Scholar 
Sibly, R. M., Barker, D., Hone, J. & Pagel, M. On the stability of populations of mammals, birds, fish and insects. Ecol. Lett. 10, 970–976 (2007).Article 

Google Scholar 
Dakos, V. et al. Ecosystem tipping points in an evolving world. Nat. Ecol. Evol. 3, 355–362 (2019).Article 

Google Scholar 
Dingemanse, N. J., Kazem, A. J. N., Réale, D. & Wright, J. Behavioural reaction norms: animal personality meets individual plasticity. Trends Ecol. Evol. 25, 81–89 (2010).Article 

Google Scholar 
Tanner, R. L. & Dowd, W. W. Inter-individual physiological variation in responses to environmental variation and environmental change: integrating across traits and time. Comp. Biochem. Physiol. A 238, 110577 (2019).Article 
CAS 

Google Scholar 
Patrick, S. C., Martin, J. G. A., Ummenhofer, C. C., Corbeau, A. & Weimerskirch, H. Albatrosses respond adaptively to climate variability by changing variance in a foraging trait. Glob. Change Biol. 27, 4564–4574 (2021).Article 
CAS 

Google Scholar 
Fayet, A. L., Clucas, G. V., Anker‐Nilssen, T., Syposz, M. & Hansen, E. S. Local prey shortages drive foraging costs and breeding success in a declining seabird, the Atlantic puffin. J. Anim. Ecol. https://doi.org/10.1111/1365-2656.13442 (2021).Pierce, C. L. Predator avoidance, microhabitat shift, and risk-sensitive foraging in larval dragonflies. Oecologia 77, 81–90 (1988).Article 
CAS 

Google Scholar 
Leibold, M. & Tessier, A. J. Contrasting patterns of body size for Daphnia species that segregate by habitat. Oecologia 86, 342–348 (1991).Article 

Google Scholar 
Charmantier, A. & Gienapp, P. Climate change and timing of avian breeding and migration: evolutionary versus plastic changes. Evol. Appl. 7, 15–28 (2014).Article 

Google Scholar 
Kopp, M. & Matuszewski, S. Rapid evolution of quantitative traits: theoretical perspectives. Evol. Appl. 7, 169–191 (2014).Article 

Google Scholar 
Williams, J. W., Ordonez, A. & Svenning, J.-C. A unifying framework for studying and managing climate-driven rates of ecological change. Nat. Ecol. Evol. 5, 17–26 (2021).Article 

Google Scholar 
Jaureguiberry, P. et al. The direct drivers of recent global anthropogenic biodiversity loss. Sci. Adv. 8, eabm9982 (2022).Article 

Google Scholar 
Chevin, L.-M., Collins, S. & Lefèvre, F. Phenotypic plasticity and evolutionary demographic responses to climate change: taking theory out to the field. Funct. Ecol. 27, 967–979 (2013).Article 

Google Scholar 
Ferriere, R. & Legendre, S. Eco-evolutionary feedbacks, adaptive dynamics and evolutionary rescue theory. Phil. Trans. R. Soc. B 368, 20120081 (2013).Article 

Google Scholar 
Rebecchi, L., Boschetti, C. & Nelson, D. R. Extreme-tolerance mechanisms in meiofaunal organisms: a case study with tardigrades, rotifers and nematodes. Hydrobiologia 847, 2779–2799 (2020).Article 

Google Scholar 
Hansson, B. & Westerberg, L. On the correlation between heterozygosity and fitness in natural populations. Mol. Ecol. 11, 2467–2474 (2002).Article 

Google Scholar 
Mammola, S., Carmona, C. P., Guillerme, T. & Cardoso, P. Concepts and applications in functional diversity. Funct. Ecol. 35, 1869–1885 (2021).Article 
CAS 

Google Scholar 
McClanahan, T. R. et al. Highly variable taxa-specific coral bleaching responses to thermal stresses. Mar. Ecol. Prog. Ser. 648, 135–151 (2020).Article 

Google Scholar 
Reside, A. E. et al. Beyond the model: expert knowledge improves predictions of species’ fates under climate change. Ecol. Appl. 29, e01824 (2019).Article 

Google Scholar 
Desjonquères, C., Gifford, T. & Linke, S. Passive acoustic monitoring as a potential tool to survey animal and ecosystem processes in freshwater environments. Freshw. Biol. 65, 7–19 (2020).Article 

Google Scholar 
Sequeira, A. M. M. et al. A standardisation framework for bio-logging data to advance ecological research and conservation. Methods Ecol. Evol. 12, 996–1007 (2021).Article 

Google Scholar 
Shimada, T. et al. Optimising sample sizes for animal distribution analysis using tracking data. Methods Ecol. Evol. 12, 288–297 (2021).Article 

Google Scholar 
Wauchope, H. S. et al. Evaluating impact using time-series data. Trends Ecol. Evol. 36, 196–205 (2021).Article 

Google Scholar 
Krause, D. J., Hinke, J. T., Perryman, W. L., Goebel, M. E. & LeRoi, D. J. An accurate and adaptable photogrammetric approach for estimating the mass and body condition of pinnipeds using an unmanned aerial system. PLoS ONE 12, e0187465 (2017).Article 

Google Scholar 
Besson, M. et al. Towards the fully automated monitoring of ecological communities. Ecol. Lett. https://doi.org/10.1111/ele.14123 (2022).Article 

Google Scholar 
Cavender-Bares, J. et al. Integrating remote sensing with ecology and evolution to advance biodiversity conservation. Nat. Ecol. Evol. 6, 506–519 (2022).Article 

Google Scholar 
Ingram, D. J., Ferreira, G. B., Jones, K. E. & Mace, G. M. Targeting conservation actions at species threat response thresholds. Trends Ecol. Evol. 36, 216–226 (2021).Article 

Google Scholar 
Keith, S. A. et al. Synchronous behavioural shifts in reef fishes linked to mass coral bleaching. Nat. Clim. Change 8, 986–991 (2018).Article 

Google Scholar 
Drake, J. M. & Griffen, B. D. Early warning signals of extinction in deteriorating environments. Nature 467, 456–459 (2010).Article 
CAS 

Google Scholar 
Enquist, B. J. et al. in Advances in Ecological Research Vol. 52 (eds Pawar, S. et al.) 249–318 (Academic Press, 2015).Wei, W. W. S. Multivariate Time Series Analysis and Applications (John Wiley & Sons, 2018).Holmes, E. E., Ward, E. J. & Wills, K. MARSS: multivariate autoregressive state-space models for analyzing time-series data. R J. 4, 11–19 (2012).Article 

Google Scholar 
Zhu, M., Yamakawa, T. & Sakai, T. Combined use of trawl fishery and research vessel survey data in a multivariate autoregressive state-space (MARSS) model to improve the accuracy of abundance index estimates. Fish. Sci. 84, 437–451 (2018).Article 
CAS 

Google Scholar 
Lai, G., Chang, W.-C., Yang, Y. & Liu, H. Modeling long- and short-term temporal patterns with deep neural networks. In The 41st International ACM SIGIR Conference on Research & Development in Information Retrieval 95–104, https://doi.org/10.1145/3209978.3210006 (ACM, 2018).Bury, T. M. et al. Deep learning for early warning signals of tipping points. Proc. Natl Acad. Sci. USA 118, e2106140118 (2021).Article 
CAS 

Google Scholar 
Lara-Benítez, P., Carranza-García, M. & Riquelme, J. C. An experimental review on deep learning architectures for time series forecasting. Int. J. Neural Syst. 31, 2130001 (2021).Article 

Google Scholar 
Guo, Q. et al. Application of deep learning in ecological resource research: theories, methods, and challenges. Sci. China Earth Sci. 63, 1457–1474 (2020).Article 

Google Scholar 
Rogers, T. L., Johnson, B. J. & Munch, S. B. Chaos is not rare in natural ecosystems. Nat. Ecol. Evol. 6, 1105–1111 (2022).Article 

Google Scholar 
Samplonius, J. M. et al. Phenological sensitivity to climate change is higher in resident than in migrant bird populations among European cavity breeders. Glob. Change Biol. 24, 3780–3790 (2018).Article 

Google Scholar 
Menzel, A. et al. Climate change fingerprints in recent European plant phenology. Glob. Change Biol. 26, 2599–2612 (2020).Article 

Google Scholar 
Koleček, J., Adamík, P. & Reif, J. Shifts in migration phenology under climate change: temperature vs. abundance effects in birds. Clim. Change 159, 177–194 (2020).Article 

Google Scholar 
Altermatt, F. et al. Big answers from small worlds: a user’s guide for protist microcosms as a model system in ecology and evolution. Methods Ecol. Evol. 6, 218–231 (2015).Article 

Google Scholar 
Beermann, A. J. et al. Multiple-stressor effects on stream macroinvertebrate communities: a mesocosm experiment manipulating salinity, fine sediment and flow velocity. Sci. Total Environ. 610–611, 961–971 (2018).Article 

Google Scholar 
Clements, C. F. & Ozgul, A. Including trait-based early warning signals helps predict population collapse. Nat. Commun. 7, 10984 (2016).Article 
CAS 

Google Scholar 
Jacquet, C. & Altermatt, F. The ghost of disturbance past: long-term effects of pulse disturbances on community biomass and composition. Proc. R. Soc. B 287, 20200678 (2020).Article 

Google Scholar 
Greggor, A. L. et al. Research priorities from animal behaviour for maximising conservation progress. Trends Ecol. Evol. 31, 953–964 (2016).Article 

Google Scholar 
Couvillon, M. J., Schürch, R. & Ratnieks, F. L. W. Waggle dance distances as integrative indicators of seasonal foraging challenges. PLoS ONE 9, e93495 (2014).Article 

Google Scholar 
Hamilton, C. D., Lydersen, C., Ims, R. A. & Kovacs, K. M. Predictions replaced by facts: a keystone species’ behavioural responses to declining Arctic sea-ice. Biol. Lett. 11, 20150803 (2015).Article 

Google Scholar 
Holt, R. E. & Jørgensen, C. Climate change in fish: effects of respiratory constraints on optimal life history and behaviour. Biol. Lett. 11, 20141032 (2015).Article 

Google Scholar 
Gauzens, B. et al. Adaptive foraging behaviour increases vulnerability to climate change. Preprint at https://doi.org/10.1101/2021.05.05.442768 (2021).Lenda, M., Witek, M., Skórka, P., Moroń, D. & Woyciechowski, M. Invasive alien plants affect grassland ant communities, colony size and foraging behaviour. Biol. Invasions 15, 2403–2414 (2013).Article 

Google Scholar 
Hertel, A. G. et al. Don’t poke the bear: using tracking data to quantify behavioural syndromes in elusive wildlife. Anim. Behav. 147, 91–104 (2019).Article 

Google Scholar 
Tini, M. et al. Use of space and dispersal ability of a flagship saproxylic insect: a telemetric study of the stag beetle (Lucanus cervus) in a relict lowland forest. Insect Conserv. Divers. 11, 116–129 (2018).Article 

Google Scholar 
Kunc, H. P. & Schmidt, R. Species sensitivities to a global pollutant: a meta-analysis on acoustic signals in response to anthropogenic noise. Glob. Change Biol. 27, 675–688 (2021).Article 

Google Scholar 
Anestis, A., Lazou, A., Pörtner, H. O. & Michaelidis, B. Behavioral, metabolic, and molecular stress responses of marine bivalve Mytilus galloprovincialis during long-term acclimation at increasing ambient temperature. Am. J. Physiol. 293, R911–R921 (2007).CAS 

Google Scholar 
Pacherres, C. O., Schmidt, G. M. & Richter, C. Autotrophic and heterotrophic responses of the coral Porites lutea to large amplitude internal waves. J. Exp. Biol. 216, 4365–4374 (2013).
Google Scholar 
Ban, S. S., Graham, N. A. J. & Connolly, S. R. Evidence for multiple stressor interactions and effects on coral reefs. Glob. Change Biol. 20, 681–697 (2014).Article 

Google Scholar 
Singh, R., Prathibha, P. & Jain, M. Effect of temperature on life-history traits and mating calls of a field cricket, Acanthogryllus asiaticus. J. Therm. Biol. 93, 102740 (2020).Article 

Google Scholar 
Pellegrini, A. Y., Romeu, B., Ingram, S. N. & Daura-Jorge, F. G. Boat disturbance affects the acoustic behaviour of dolphins engaged in a rare foraging cooperation with fishers. Anim. Conserv. 24, 613–625 (2021).Article 

Google Scholar 
McMahan, M. D. & Grabowski, J. H. Nonconsumptive effects of a range-expanding predator on juvenile lobster (Homarus americanus) population dynamics. Ecosphere 10, e02867 (2019).Article 

Google Scholar 
Vilhunen, S., Hirvonen, H. & Laakkonen, M. V.-M. Less is more: social learning of predator recognition requires a low demonstrator to observer ratio in Arctic charr (Salvelinus alpinus). Behav. Ecol. Sociobiol. 57, 275–282 (2005).Article 

Google Scholar 
Ortega, Z., Mencía, A. & Pérez-Mellado, V. Rapid acquisition of antipredatory responses to new predators by an insular lizard. Behav. Ecol. Sociobiol. 71, 1 (2017).Article 

Google Scholar 
Fox, R. J., Donelson, J. M., Schunter, C., Ravasi, T. & Gaitán-Espitia, J. D. Beyond buying time: the role of plasticity in phenotypic adaptation to rapid environmental change. Phil. Trans. R. Soc. B 374, 20180174 (2019).Article 

Google Scholar 
Pigeon, G., Ezard, T. H. G., Festa-Bianchet, M., Coltman, D. W. & Pelletier, F. Fluctuating effects of genetic and plastic changes in body mass on population dynamics in a large herbivore. Ecology 98, 2456–2467 (2017).Article 

Google Scholar 
Lomolino, M. V. & Perault, D. R. Body size variation of mammals in a fragmented, temperate rainforest. Conserv. Biol. 21, 1059–1069 (2007).Article 

Google Scholar 
Gardner, J. L., Peters, A., Kearney, M. R., Joseph, L. & Heinsohn, R. Declining body size: a third universal response to warming? Trends Ecol. Evol. 26, 285–291 (2011).Article 

Google Scholar 
Sheridan, J. A. & Bickford, D. Shrinking body size as an ecological response to climate change. Nat. Clim. Change 1, 401–406 (2011).Article 

Google Scholar 
Thoral, E. et al. Changes in foraging mode caused by a decline in prey size have major bioenergetic consequences for a small pelagic fish. J. Anim. Ecol. 90, 2289–2301 (2021).Article 

Google Scholar 
Stirling, I. & Derocher, A. E. Effects of climate warming on polar bears: a review of the evidence. Glob. Change Biol. 18, 2694–2706 (2012).Article 

Google Scholar 
Spanbauer, T. L. et al. Body size distributions signal a regime shift in a lake ecosystem. Proc. R. Soc. B 283, 20160249 (2016).Article 

Google Scholar 
Bjorndal, K. A. et al. Ecological regime shift drives declining growth rates of sea turtles throughout the West Atlantic. Glob. Change Biol. 23, 4556–4568 (2017).Article 

Google Scholar 
Eshun-Wilson, F., Wolf, R., Andersen, T., Hessen, D. O. & Sperfeld, E. UV radiation affects antipredatory defense traits in Daphnia pulex. Ecol. Evol. 10, 14082–14097 (2020).Article 

Google Scholar 
Zhang, H., Hollander, J. & Hansson, L.-A. Bi-directional plasticity: rotifer prey adjust spine length to different predator regimes. Sci. Rep. 7, 10254 (2017).Article 

Google Scholar 
Simbula, G., Vignoli, L., Carretero, M. A. & Kaliontzopoulou, A. Fluctuating asymmetry as biomarker of pesticides exposure in the Italian wall lizards (Podarcis siculus). Zoology 147, 125928 (2021).Article 

Google Scholar 
Leary, R. F. & Allendorf, F. W. Fluctuating asymmetry as an indicator of stress: implications for conservation biology. Trends Ecol. Evol. 4, 214–217 (1989).Article 
CAS 

Google Scholar 
Gavrilchuk, K. et al. Trophic niche partitioning among sympatric baleen whale species following the collapse of groundfish stocks in the Northwest Atlantic. Mar. Ecol. Prog. Ser. 497, 285–301 (2014).Article 

Google Scholar 
Kershaw, J. L. et al. Declining reproductive success in the Gulf of St. Lawrence’s humpback whales (Megaptera novaeangliae) reflects ecosystem shifts on their feeding grounds. Glob. Change Biol. 27, 1027–1041 (2021).Article 
CAS 

Google Scholar 
Rode, K. D., Amstrup, S. C. & Regehr, E. V. Reduced body size and cub recruitment in polar bears associated with sea ice decline. Ecol. Appl. 20, 768–782 (2010).Article 

Google Scholar 
Obbard, M. E. et al. Re-assessing abundance of Southern Hudson Bay polar bears by aerial survey: effects of climate change at the southern edge of the range. Arct. Sci. 4, 634–655 (2018).Article 

Google Scholar 
Hutchings, J. A. The cod that got away. Nature 428, 899–900 (2004).Article 
CAS 

Google Scholar 
Zhang, F. Early warning signals of population productivity regime shifts in global fisheries. Ecol. Indic. 115, 106371 (2020).Article 

Google Scholar 
Fulton, G. R. The Bramble Cay melomys: the first mammalian extinction due to human-induced climate change. Pac. Conserv. Biol. 23, 1–3 (2017).Article 

Google Scholar  More

Essl, F. et al. A conceptual framework for range-expanding species that track human-induced environmental change. BioScience 69, 908–919 (2019).Article 

Google Scholar 
Lenoir, J. et al. Species better track climate warming in the oceans than on land. Nat. Ecol. Evol. 4, 1044–1059 (2020).Article 

Google Scholar 
Pecl, G. T. et al. Biodiversity redistribution under climate change: impacts on ecosystems and human well-being. Science 355, eaai9214 (2017).Article 

Google Scholar 
Freeman, B. G., Lee-Yaw, J. A., Sunday, J. M. & Hargreaves, A. L. Expanding, shifting and shrinking: the impact of global warming on species’ elevational distributions. Glob. Ecol. Biogeogr. 27, 1268–1276 (2018).Article 

Google Scholar 
van Kleunen, M. et al. Global exchange and accumulation of non-native plants. Nature 525, 100–103 (2015).Article 

Google Scholar 
Alexander, J. M. et al. Lags in the response of mountain plant communities to climate change. Glob. Change Biol. 24, 563–579 (2018).Article 

Google Scholar 
Graae, B. J. et al. Stay or go—how topographic complexity influences alpine plant population and community responses to climate change. Perspect. Plant Ecol. Evol. Syst. 30, 41–50 (2018).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 (2019).Article 

Google Scholar 
Steinbauer, M. J. et al. Accelerated increase in plant species richness on mountain summits is linked to warming. Nature 556, 231–234 (2018).Article 
CAS 

Google Scholar 
Mamantov, M. A., Gibson-Reinemer, D. K., Linck, E. B. & Sheldon, K. S. Climate-driven range shifts of montane species vary with elevation. Glob. Ecol. Biogeogr. 30, 784–794 (2021).Article 

Google Scholar 
Alexander, J. M. et al. Plant invasions into mountains and alpine ecosystems: current status and future challenges. Alp. Bot. 126, 89–103 (2016).Article 

Google Scholar 
Pauchard, A. et al. Ain’t no mountain high enough: plant invasions reaching new elevations. Front. Ecol. Environ. 7, 479–486 (2009).Article 

Google Scholar 
Alexander, J. M., MIREN Consortium et al. Assembly of nonnative floras along elevational gradients explained by directional ecological filtering. Proc. Natl Acad. Sci. USA 108, 656–661 (2011).Article 
CAS 

Google Scholar 
Seipel, T. et al. Processes at multiple spatial scales determine non-native plant species richness and similarity in mountain regions around the world. Glob. Ecol. Biogeogr. 21, 236–246 (2012).Article 

Google Scholar 
Dainese, M. et al. Human disturbance and upward expansion of plants in a warming climate. Nat. Clim. Change 7, 577–580 (2017).Article 

Google Scholar 
McDougall, K. L. et al. Running off the road: roadside non-native plants invading mountain vegetation. Biol. Invasions 20, 3461–3473 (2018).Article 

Google Scholar 
Petitpierre, B. et al. Will climate change increase the risk of plant invasions into mountains? Ecol. Appl. 26, 530–544 (2016).Article 

Google Scholar 
Lembrechts, J. J. et al. Microclimate variability in alpine ecosystems as stepping stones for non‐native plant establishment above their current elevational limit. Ecography 41, 900–909 (2017).Article 

Google Scholar 
Haider, S. et al. Mountain roads and non-native species modify elevational patterns of plant diversity. Glob. Ecol. Biogeogr. 27, 667–678 (2018).Article 

Google Scholar 
Wolf, A., Zimmerman, N. B., Anderegg, W. R. L., Busby, P. E. & Christensen, J. Altitudinal shifts of the native and introduced flora of California in the context of 20th-century warming. Glob. Ecol. Biogeogr. 25, 418–429 (2016).Article 

Google Scholar 
Seipel, T., Alexander, J. M., Edwards, P. J. & Kueffer, C. Range limits and population dynamics of non-native plants spreading along elevation gradients. Perspect. Plant Ecol. Evol. Syst. 20, 46–55 (2016).Article 

Google Scholar 
Koide, D., Yoshida, K., Daehler, C. C. & Mueller-Dombois, D. An upward elevation shift of native and non-native vascular plants over 40 years on the island of Hawai’i. J. Vegetation Sci. 28, 939–950 (2017).Article 

Google Scholar 
Becker, T., Dietz, H., Billeter, R., Buschmann, H. & Edwards, P. J. Altitudinal distribution of alien plant species in the Swiss Alps. Perspect. Plant Ecol. Evol. Syst. 7, 173–183 (2005).Article 

Google Scholar 
Haider, S. et al. The role of bioclimatic origin, residence time and habitat context in shaping non-native plant distributions along an altitudinal gradient. Biol. Invasions 12, 4003–4018 (2010).Article 

Google Scholar 
Pyšek, P., Jarošík, V., Pergl, J. & Wild, J. Colonization of high altitudes by alien plants over the last two centuries. Proc. Natl Acad. Sci. USA 108, 439–440 (2011).Article 

Google Scholar 
Gottfried, M. et al. Continent-wide response of mountain vegetation to climate change. Nat. Clim. Change 2, 111–115 (2012).Article 

Google Scholar 
Lenoir, J. et al. Going against the flow: potential mechanisms for unexpected downslope range shifts in a warming climate. Ecography 33, 295–303 (2010).
Google Scholar 
Crimmins, S. M., Dobrowski, S. Z., Greenberg, J. A., Abatzoglou, J. T. & Mynsberge, A. R. Changes in climatic water balance drive downhill shifts in plant species’ optimum elevations. Science 331, 324–327 (2011).Article 
CAS 

Google Scholar 
Rumpf, S. B. et al. Range dynamics of mountain plants decrease with elevation. Proc. Natl Acad. Sci. USA 115, 1848–1853 (2018).Article 
CAS 

Google Scholar 
Scherrer, D. & Körner, C. Topographically controlled thermal-habitat differentiation buffers alpine plant diversity against climate warming. J. Biogeogr. 38, 406–416 (2011).Article 

Google Scholar 
Kelly, C. & Price, T. D. Correcting for regression to the mean in behavior and ecology. Am. Nat. 166, 700–707 (2005).Article 

Google Scholar 
Mazalla, L. & Diekmann, M. Regression to the mean in vegetation science. J. Vegetation Sci. 33, e13117 (2022).Article 

Google Scholar 
Colwell, R. K. & Lees, D. C. The mid-domain effect: geometric constraints on the geography of species richness. Trends Ecol. Evol. 15, 70–76 (2000).Article 
CAS 

Google Scholar 
Colwell, R. K., Brehm, G., Cardelús, C. L., Gilman, A. C. & Longino, J. T. Global warming, elevational range shifts, and lowland biotic attrition in the wet tropics. Science 322, 258–261 (2008).Article 
CAS 

Google Scholar 
Taheri, S., Naimi, B., Rahbek, C. & Araújo, M. B. Improvements in reports of species redistribution under climate change are required. Sci. Adv. 7, eabe1110 (2021).Article 

Google Scholar 
Haider, S. et al. Think globally, measure locally: the MIREN standardized protocol for monitoring plant species distributions along elevation gradients. Ecol. Evol. 12, e8590 (2022).Article 

Google Scholar 
Jacobsen, D. The dilemma of altitudinal shifts: caught between high temperature and low oxygen. Front. Ecol. Environ. 18, 211–218 (2020).Article 

Google Scholar 
Kueffer, C. et al. in Plant Invasions in Protected Areas Vol. 7 (eds Foxcroft, L. C. et al.) 89–113 (Springer, 2013).Halbritter, A. H., Alexander, J. M., Edwards, P. J. & Billeter, R. How comparable are species distributions along elevational and latitudinal climate gradients? Glob. Ecol. Biogeogr. 22, 1228–1237 (2013).Article 

Google Scholar 
Vitasse, Y. et al. Phenological and elevational shifts of plants, animals and fungi under climate change in the European Alps. Biol. Rev. 96, 1816–1835 (2021).Article 

Google Scholar 
Angert, A. L. et al. Do species’ traits predict recent shifts at expanding range edges? Ecol. Lett. 14, 677–689 (2011).Article 

Google Scholar 
Matteodo, M., Wipf, S., Stöckli, V., Rixen, C. & Vittoz, P. Elevation gradient of successful plant traits for colonizing alpine summits under climate change. Environ. Res. Lett. 8, 024043 (2013).Article 

Google Scholar 
Lembrechts, J. et al. Disturbance is the key to plant invasions in cold environments. Proc. Natl Acad. Sci. USA 113, 14061–14066 (2016).Article 
CAS 

Google Scholar 
Bates, D., Maechler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).Article 

Google Scholar 
Pebesma, E. Simple features for R: standardized support for spatial vector data. R J. 10, 439–446 (2018).Article 

Google Scholar 
Kahle, D. & Wickham, H. ggmap: spatial visualization with ggplot2. R J. 5, 144–161 (2013). http://journal.r-project.org/archive/2013-1/kahle-wickham.pdfSeipel, T., Haider, S. & MIREN consortium. MIREN survey of plant species in mountains (v2.0). Zenodo https://doi.org/10.5281/zenodo.5529072 (2022). More

Rockström, J. et al. A safe operation space for humanity. Nature 461, 472–475 (2009).Article 

Google Scholar 
Chancel, L., Piketty, T., Saez, E. & Zucman, G. World Inequality Report 2022 (Belknap Press, 2022).Ivanova, D. et al. Environmental impact assessment of household consumption. J. Ind. Ecol. 20, 526–536 (2016).Article 
CAS 

Google Scholar 
Steen-Olsen, K., Weinzettel, J., Cranston, G., Ercin, A. E. & Hertwich, E. G. Carbon, land, and water footprint accounts for the European Union: consumption, production, and displacements through international trade. Environ. Sci. Technol. 46, 10883–10891 (2012).Article 
CAS 

Google Scholar 
Tukker, A. et al. Environmental and resource footprints in a global context: Europe’s structural deficit in resource endowments. Glob. Environ. Change 40, 171–181 (2016).Article 

Google Scholar 
Bruckner, B., Hubacek, K., Shan, Y., Zhong, H. & Feng, K. Impacts of poverty alleviation on national and global carbon emissions. Nat. Sustain. 5, 311–320 (2022).Article 

Google Scholar 
Hubacek, K. et al. Global carbon inequality. Energy, Ecol. Environ. 2, 361–369 (2017).Article 

Google Scholar 
Yu, Y., Feng, K. & Hubacek, K. Tele-connecting local consumption to global land use. Glob. Environ. Change 23, 1178–1186 (2013).Article 

Google Scholar 
Wilting, H. C., Schipper, A. M., Bakkenes, M., Meijer, J. R. & Huijbregts, M. A. J. Quantifying biodiversity losses due to human consumption: a global-scale footprint analysis. Environ. Sci. Technol. 51, 3298–3306 (2017).Article 
CAS 

Google Scholar 
Lucas, P. L., Wilting, H. C., Hof, A. F. & Van Vuuren, D. P. Allocating planetary boundaries to large economies: distributional consequences of alternative perspectives on distributive fairness. Glob. Environ. Change 60, 102017 (2020).Article 

Google Scholar 
Beylot, A. et al. Assessing the environmental impacts of EU consumption at macro-scale. J. Clean. Prod. 216, 382–393 (2019).Article 

Google Scholar 
Koslowski, M., Moran, D. D., Tisserant, A., Verones, F. & Wood, R. Quantifying Europe’s biodiversity footprints and the role of urbanization and income. Glob. Sustain. 3, e1 (2020).Lutter, S., Pfister, S., Giljum, S., Wieland, H. & Mutel, C. Spatially explicit assessment of water embodied in European trade: a product-level multi-regional input-output analysis. Glob. Environ. Change 38, 171–182 (2016).Article 

Google Scholar 
Stadler, K. et al. EXIOBASE 3 (3.8.1) [Data set]. Zenodo https://doi.org/10.5281/ZENODO.4588235 (2021).Roadmap to a Resource Efficient Europe (European Commission, 2011).Steinmann, Z. J. N. et al. Headline environmental indicators revisited with the global multi-regional input–output database EXIOBASE. J. Ind. Ecol. 22, 565–573 (2018).Article 

Google Scholar 
Ivanova, D. et al. Mapping the carbon footprint of EU regions. Environ. Res. Lett. 12, 054013 (2017).Wiedmann, T. O. et al. The material footprint of nations. Proc. Natl Acad. Sci. USA 112, 6271–6276 (2015).Article 
CAS 

Google Scholar 
Lenzen, M. et al. Implementing the material footprint to measure progress towards Sustainable Development Goals 8 and 12. Nat. Sustain. 5, 157–166 (2022).Dorninger, C. et al. The effect of industrialization and globalization on domestic land-use: a global resource footprint perspective. Glob. Environ. Change 69, 102311 (2021).Article 

Google Scholar 
Mekonnen, M. M. & Gerbens-Leenes, W. The water footprint of food. Water 12, 12 (2020).Article 

Google Scholar 
Prell, C. & Feng, K. Unequal carbon exchanges: the environmental and economic impacts of iconic U.S. consumption items. J. Ind. Ecol. 20, 537–546 (2016).Article 

Google Scholar 
Prell, C., Feng, K., Sun, L., Geores, M. & Hubacek, K. The economic gains and environmental losses of US consumption: a world-systems and input-output approach. Soc. Forces 93, 405–428 (2014).Article 

Google Scholar 
Prell, C. Wealth and pollution inequalities of global trade: a network and input-output approach. Soc. Sci. J. 53, 111–121 (2016).Article 

Google Scholar 
World Economic Outlook (October 2022) (International Monetary Fund, 2022); https://www.imf.org/external/datamapper/datasets/WEOWilting, H. C., Schipper, A. M., Ivanova, O., Ivanova, D. & Huijbregts, M. A. J. Subnational greenhouse gas and land-based biodiversity footprints in the European Union. J. Ind. Ecol. 25, 79–94 (2021). https://doi.org/10.1111/jiec.13042Cabernard, L. & Pfister, S. A highly resolved MRIO database for analyzing environmental footprints and Green Economy Progress. Sci. Total Environ. 755, 142587 (2021).Jakob, M., Ward, H. & Steckel, J. C. Sharing responsibility for trade-related emissions based on economic benefits. Glob. Environ. Chang. 66, 102207 (2021).Article 

Google Scholar 
Wood, R. et al. The structure, drivers and policy implications of the European carbon footprint. Clim. Policy 20, S39–S57 (2020).Article 

Google Scholar 
Wood, R. et al. Growth in environmental footprints and environmental impacts embodied in trade: resource efficiency indicators from EXIOBASE3. J. Ind. Ecol. 22, 553–564 (2018).Article 

Google Scholar 
Hubacek, K., Chen, X., Feng, K., Wiedmann, T. & Shan, Y. Evidence of decoupling consumption-based CO2 emissions from economic growth. Adv. Appl. Energy 4, 100074 (2021).Article 

Google Scholar 
Wiedmann, T. & Lenzen, M. Environmental and social footprints of international trade. Nat. Geosci. 11, 314–321 (2018).Article 
CAS 

Google Scholar 
Dorninger, C. et al. Global patterns of ecologically unequal exchange: Implications for sustainability in the 21st century. Ecol. Econ. 179, 106824 (2021).Article 

Google Scholar 
Hickel, J., Dorninger, C., Wieland, H. & Suwandi, I. Imperialist appropriation in the world economy: drain from the global South through unequal exchange, 1990–2015. Glob. Environ. Change 73, 102467 (2022).Poore, J. & Nemecek, T. Reducing food’s environmental impacts through producers and consumers. Science 360, 987–992 (2018).Article 
CAS 

Google Scholar 
Ivanova, D. et al. Quantifying the potential for climate change mitigation of consumption options. Environ. Res. Lett. 15, 093001 (2020).Springmann, M. et al. Options for keeping the food system within environmental limits. Nature 562, 519–525 (2018).Article 
CAS 

Google Scholar 
Ivanova, D. & Wood, R. The unequal distribution of household carbon footprints in Europe and its link to sustainability. Glob. Sustain. 3, e18 (2020).Hickel, J., O’Neill, D. W., Fanning, A. L. & Zoomkawala, H. National responsibility for ecological breakdown: a fair-shares assessment of resource use, 1970–2017. Lancet Planet. Heal. 6, e342–e349 (2022).Article 

Google Scholar 
Otto, I. M., Kim, K. M., Dubrovsky, N. & Lucht, W. Shift the focus from the super-poor to the super-rich. Nat. Clim. Change 9, 82–84 (2019).Article 

Google Scholar 
Wiedmann, T., Lenzen, M., Keyßer, L. T. & Steinberger, J. K. Scientists’ warning on affluence. Nat. Commun. 11, 3107 (2020).Nielsen, K. S., Nicholas, K. A., Creutzig, F., Dietz, T. & Stern, P. C. The role of high-socioeconomic-status people in locking in or rapidly reducing energy-driven greenhouse gas emissions. Nat. Energy 6, 1011–1016 (2021).Article 

Google Scholar 
Jakob, M. Why carbon leakage matters and what can be done against it. One Earth 4, 609–614 (2021).Article 

Google Scholar 
Lave, L. B. Using input–output analysis to estimate economy-wide discharges. Environ. Sci. Technol. 29, 420A–426A (1995).Article 
CAS 

Google Scholar 
Wiedmann, T. A review of recent multi-region input–output models used for consumption-based emission and resource accounting. Ecol. Econ. 69, 211–222 (2009).Article 

Google Scholar 
Ewing, B. R. et al. Integrating ecological and water footprint accounting in a multi-regional input–output framework. Ecol. Indic. 23, 1–8 (2012).Article 

Google Scholar 
Brizga, J., Feng, K. & Hubacek, K. Household carbon footprints in the Baltic States: a global multi-regional input–output analysis from 1995 to 2011. Appl. Energy 189, 780–788 (2017).Hertwich, E. G. & Peters, G. P. Carbon footprint of nations: a global, trade-linked analysis. Environ. Sci. Technol. 43, 6414–6420 (2009).Article 
CAS 

Google Scholar 
Zhong, H., Feng, K., Sun, L., Cheng, L. & Hubacek, K. Household carbon and energy inequality in Latin American and Caribbean countries. J. Environ. Manag. 273, 110979 (2020).Article 

Google Scholar 
Stadler, K. et al. EXIOBASE 3: developing a time series of detailed environmentally extended multi-regional input–output tables. J. Ind. Ecol. 22, 502–515 (2018).Article 

Google Scholar 
Hardadi, G., Buchholz, A. & Pauliuk, S. Implications of the distribution of German household environmental footprints across income groups for integrating environmental and social policy design. J. Ind. Ecol. 25, 95–113 (2021).Zhang, Q. et al. Transboundary health impacts of transported global air pollution and international trade. Nature 543, 705–709 (2017).Article 
CAS 

Google Scholar 
Hoekstra, A. Y., Mekonnen, M. M., Chapagain, A. K., Mathews, R. E. & Richter, B. D. Global monthly water scarcity: blue water footprints versus blue water availability. PLoS ONE 7, e32688 (2012).Article 
CAS 

Google Scholar 
IPCC Climate Change 2007: The Physical Science Basis (eds Solomon, S. et al.) (Cambridge Univ. Press, 2007).Schmidt, S. et al. Understanding GHG emissions from Swedish consumption—current challenges in reaching the generational goal. J. Clean. Prod. 212, 428–437 (2019).Article 

Google Scholar 
Huijbregts, M. A. J. Priority Assessment of Toxic Substances in the Frame of LCA. Development and Application of the Multi-Media Fate, Exposure and Effect Model USES-LCA (Interfaculty Department of Envrionmental Science, 1999).Huijbregts, M. A. J. Priority Assessment of Toxic Substances in the Frame of LCA. Time Horizon Dependency in Toxicity Potentials Calculated with the Multi-Media Fate, Exposure and Effects Model USES-LCA (Institute for Biodiversity and Ecosystem Dynamics, 2000).International Reference Life Cycle Data System (ILCD) Handbook (Publications Office EU, 2011).Verones, F., Moran, D., Stadler, K., Kanemoto, K. & Wood, R. Resource footprints an d their ecosystem consequences. Sci. Rep. 7, 40743 (2017).Chaudhary, A., Pfister, S. & Hellweg, S. Spatially explicit analysis of biodiversity loss due to global agriculture, pasture and forest land use from a producer and consumer perspective. Environ. Sci. Technol. 50, 3928–3936 (2016).Article 
CAS 

Google Scholar 
Chaudhary, A., Verones, F., De Baan, L. & Hellweg, S. Quantifying land use impacts on biodiversity: combining species-area models and vulnerability indicators. Environ. Sci. Technol. 49, 9987–9995 (2015).Article 
CAS 

Google Scholar 
Marquardt, S. G. et al. Consumption-based biodiversity footprints—do different indicators yield different results? Ecol. Indic. 103, 461–470 (2019).Article 

Google Scholar 
World Development Indicators DataBank (World Bank, 2022); https://databank.worldbank.org/source/world-development-indicatorsWorld Population Prospects 2022 (United Nations, 2022); https://population.un.org/wpp/Natural Earth Vector (Natural Earth, 2022); https://www.naturalearthdata.com/Lahti, L., Huovari, J., Kainu, M. & Biecek, P. Retrieval and analysis of eurostat open data with the Eurostat package. R J. 9, 385–392 (2017).Castellani, V., Beylot, A. & Sala, S. Environmental impacts of household consumption in Europe: comparing process-based LCA and environmentally extended input-output analysis. J. Clean. Prod. 240, 117966 (2019).Article 

Google Scholar  More

Our results revealed statistically significant relationships between CTs and BA worldwide (P  More

Rahman, M. A. et al. Comparing the infiltration potentials of soils beneath the canopies of two contrasting urban tree species. Urban For. Urban Green. 38, 22–32. https://doi.org/10.1016/j.ufug.2018.11.002 (2019).Article 

Google Scholar 
Zölch, T., Henze, L., Keilholz, P. & Pauleit, S. Regulating urban surface runoff through nature-based solutions – An assessment at the micro-scale. Environ. Res. 157, 135–144. https://doi.org/10.1016/j.envres.2017.05.023 (2017).Article 
CAS 

Google Scholar 
Barron, O. V., Barr, A. D. & Donn, M. J. Effect of urbanisation on the water balance of a catchment with shallow groundwater. J. Hydrol. 485, 162–176. https://doi.org/10.1016/j.jhydrol.2012.04.027 (2013).Article 
ADS 

Google Scholar 
Rosenzweig, B. R. et al. The value of urban flood modeling. Earth’s Future 9, e2020EF001739. https://doi.org/10.1029/2020EF001739 (2021).Article 
ADS 

Google Scholar 
Pauleit, S., Fryd, O., Backhaus, A. & Jensen, M. B. In Encyclopedia of Sustainability Science and Technology (ed. Meyers, R. A.) 1–29 (Springer, 2020).
Google Scholar 
Rahman, M. A. et al. Traits of trees for cooling urban heat islands: A meta-analysis. Build. Environ. 170, 106606. https://doi.org/10.1016/j.buildenv.2019.106606 (2020).Article 

Google Scholar 
Ziter, C. D., Pedersen, E. J., Kucharik, C. J. & Turner, M. G. Scale-dependent interactions between tree canopy cover and impervious surfaces reduce daytime urban heat during summer. Proc. Natl. Acad. Sci. USA 116, 7575–7580. https://doi.org/10.1073/pnas.1817561116 (2019).Article 
ADS 
CAS 

Google Scholar 
Waldrop, M. M. News feature: The quest for the sustainable city. Proc. Natl. Acad. Sci. 116, 17134–17138. https://doi.org/10.1073/pnas.1912802116 (2019).Article 
ADS 
CAS 

Google Scholar 
Cleugh, H. A., Bui, E., Simon, D., Xu, J. & Mitchell, V. G. The Impact of Suburban Design on Water Use and Microclimate (2005).Chan, F. K. S. et al. “Sponge City” in China—A breakthrough of planning and flood risk management in the urban context. Land Use Policy 76, 772–778. https://doi.org/10.1016/j.landusepol.2018.03.005 (2018).Article 

Google Scholar 
Morgan, R. P. C. Soil Erosion and Conservation (Wiley, 2005).
Google Scholar 
Xu, C. et al. Surface runoff in urban areas: The role of residential cover and urban growth form. J. Clean. Prod. 262, 121421. https://doi.org/10.1016/j.jclepro.2020.121421 (2020).Article 

Google Scholar 
Ostoić, S. K. & van den Bosch, C. C. K. Exploring global scientific discourses on urban forestry. Urban For. Urban Green. 14, 129–138. https://doi.org/10.1016/j.ufug.2015.01.001 (2015).Article 

Google Scholar 
Rahman, M. A. et al. Tree cooling effects and human thermal comfort under contrasting species and sites. Agric. For. Meteorol. 287, 107947. https://doi.org/10.1016/j.agrformet.2020.107947 (2020).Article 
ADS 

Google Scholar 
Rötzer, T., Rahman, M. A., Moser-Reischl, A., Pauleit, S. & Pretzsch, H. Process based simulation of tree growth and ecosystem services of urban trees under present and future climate conditions. Sci. Total Environ. 676, 651–664. https://doi.org/10.1016/j.scitotenv.2019.04.235 (2019).Article 
ADS 
CAS 

Google Scholar 
Grote, R. et al. Functional traits of urban trees: Air pollution mitigation potential. Front. Ecol. Environ. 14, 543–550. https://doi.org/10.1002/fee.1426 (2016).Article 

Google Scholar 
Pace, R. et al. A single tree model to consistently simulate cooling, shading, and pollution uptake of urban trees. Int. J. Biometeorol. 65, 277–289. https://doi.org/10.1007/s00484-020-02030-8 (2021).Article 
ADS 

Google Scholar 
Kuehler, E., Hathaway, J. & Tirpak, A. Quantifying the benefits of urban forest systems as a component of the green infrastructure stormwater treatment network. Ecohydrology https://doi.org/10.1002/eco.1813 (2017).Article 

Google Scholar 
Rahman, M. A., Moser, A., Gold, A., Rötzer, T. & Pauleit, S. Vertical air temperature gradients under the shade of two contrasting urban tree species during different types of summer days. Sci. Total Environ. 633, 100–111. https://doi.org/10.1016/j.scitotenv.2018.03.168 (2018).Article 
ADS 
CAS 

Google Scholar 
Rahman, M. A., Smith, J. G., Stringer, P. & Ennos, A. R. Effect of rooting conditions on the growth and cooling ability of Pyrus calleryana. Urban For. Urban Green. 10, 185–192. https://doi.org/10.1016/j.ufug.2011.05.003 (2011).Article 

Google Scholar 
Schellekens, J., Scatena, F. N., Bruijnzeel, L. A. & Wickel, A. J. Modelling rainfall interception by a lowland tropical rain forest in northeastern Puerto Rico. J. Hydrol. 225, 168–184. https://doi.org/10.1016/S0022-1694(99)00157-2 (1999).Article 
ADS 

Google Scholar 
Guevara-Escobar, A., González-Sosa, E., Véliz-Chávez, C., Ventura-Ramos, E. & Ramos-Salinas, M. Rainfall interception and distribution patterns of gross precipitation around an isolated Ficus benjamina tree in an urban area. J. Hydrol. 333, 532–541. https://doi.org/10.1016/j.jhydrol.2006.09.017 (2007).Article 
ADS 

Google Scholar 
Xiao, Q. F. & McPherson, E. G. Surface water storage capacity of twenty tree species in Davis, California. J. Environ. Qual. 45, 188–198. https://doi.org/10.2134/jeq2015.02.0092 (2016).Article 
CAS 

Google Scholar 
Xiao, Q. F., McPherson, E. G., Ustin, S. L. & Grismer, M. E. A new approach to modeling tree rainfall interception. J. Geophys. Res. Atmos. 105, 29173–29188. https://doi.org/10.1029/2000jd900343 (2000).Article 
ADS 

Google Scholar 
Carlyle-Moses, D. E. & Gash, J. H. C. In Forest Hydrology and Biogeochemistry: Synthesis of Past Research and Future Directions (eds Levia, D. F. et al.) 407–423 (Springer, 2011).Chapter 

Google Scholar 
Hirano, T. et al. The difference in the short-term runoff characteristic between the coniferous catchment and the deciduous catchment: The effects of storm size on storm generation processes of small forested catchment. J. Jpn. Soc. Hydrol. Water Resour. 22, 24–39. https://doi.org/10.3178/jjshwr.22.24 (2009).Article 

Google Scholar 
Chandler, K. R. & Chappell, N. A. Influence of individual oak (Quercus robur) trees on saturated hydraulic conductivity. For. Ecol. Manage. 256, 1222–1229. https://doi.org/10.1016/j.foreco.2008.06.033 (2008).Article 

Google Scholar 
Stewart, I. D. A systematic review and scientific critique of methodology in modern urban heat island literature. Int. J. Climatol. 31, 200–217. https://doi.org/10.1002/joc.2141 (2011).Article 

Google Scholar 
Beck, H. E. et al. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Sci. Data 5, 180214. https://doi.org/10.1038/sdata.2018.214 (2018).Article 

Google Scholar 
Moreno-de las Heras, M., Nicolau, J. M., Merino-Martín, L. & Wilcox, B. P. Plot-scale effects on runoff and erosion along a slope degradation gradient. Water Resour. Res. 46, W04503. https://doi.org/10.1029/2009WR007875 (2010).Article 
ADS 

Google Scholar 
Wu, L., Peng, M., Qiao, S. & Ma, X.-Y. Effects of rainfall intensity and slope gradient on runoff and sediment yield characteristics of bare loess soil. Environ. Sci. Pollut. Res. 25, 3480–3487. https://doi.org/10.1007/s11356-017-0713-8 (2018).Article 

Google Scholar 
Rutter, A. J., Kershaw, K. A., Robins, P. C. & Morton, A. J. A predictive model of rainfall interception in forests, 1. Derivation of the model from observations in a plantation of Corsican pine. Agric. Meteorol. 9, 367–384. https://doi.org/10.1016/0002-1571(71)90034-3 (1971).Article 

Google Scholar 
Gash, J. H. C. An analytical model of rainfall interception by forests. Q. J. R. Meteorol. Soc. 105, 43–55. https://doi.org/10.1002/qj.49710544304 (1979).Article 
ADS 

Google Scholar 
Véliz-Chávez, C., Mastachi-Loza, C. A., Gonzalez-Sosa, E., Becerril-Pia, R. & Ramos-Salinas, N. M. Canopy storage implications on interception loss modeling. Am. J. Plant Sci. 05, 3032–3048. https://doi.org/10.4236/ajps.2014.520320 (2014).Article 

Google Scholar 
Fan, J., Oestergaard, K. T., Guyot, A. & Lockington, D. A. Measuring and modeling rainfall interception losses by a native Banksia woodland and an exotic pine plantation in subtropical coastal Australia. J. Hydrol. 515, 156–165. https://doi.org/10.1016/j.jhydrol.2014.04.066 (2014).Article 
ADS 

Google Scholar 
Ghimire, C. P., Bruijnzeel, L. A., Lubczynski, M. W. & Bonell, M. Rainfall interception by natural and planted forests in the Middle Mountains of Central Nepal. J. Hydrol. 475, 270–280. https://doi.org/10.1016/j.jhydrol.2012.09.051 (2012).Article 
ADS 

Google Scholar 
Pereira, F. L. et al. Modelling interception loss from evergreen oak Mediterranean savannas: Application of a tree-based modelling approach. Agric. For. Meteorol. 149, 680–688. https://doi.org/10.1016/j.agrformet.2008.10.014 (2009).Article 
ADS 

Google Scholar 
Pypker, T. G., Bond, B. J., Link, T. E., Marks, D. & Unsworth, M. H. The importance of canopy structure in controlling the interception loss of rainfall: Examples from a young and an old-growth Douglas-fir forest. Agric. For. Meteorol. 130, 113–129. https://doi.org/10.1016/j.agrformet.2005.03.003 (2005).Article 
ADS 

Google Scholar 
Ringgaard, R., Herbst, M. & Friborg, T. Partitioning forest evapotranspiration: Interception evaporation and the impact of canopy structure, local and regional advection. J. Hydrol. 517, 677–690. https://doi.org/10.1016/j.jhydrol.2014.06.007 (2014).Article 
ADS 

Google Scholar 
Price, A. G. & Carlyle-Moses, D. E. Measurement and modelling of growing-season canopy water fluxes in a mature mixed deciduous forest stand, southern Ontario, Canada. Agric. For. Meteorol. 119, 69–85. https://doi.org/10.1016/S0168-1923(03)00117-5 (2003).Article 
ADS 

Google Scholar 
Fathizadeh, O., Hosseini, S. M., Zimmermann, A., Keim, R. F. & Darvishi Boloorani, A. Estimating linkages between forest structural variables and rainfall interception parameters in semi-arid deciduous oak forest stands. Sci. Total Environ. 601–602, 1824–1837. https://doi.org/10.1016/j.scitotenv.2017.05.233 (2017).Article 
ADS 
CAS 

Google Scholar 
Livesley, S. J., Baudinette, B. & Glover, D. Rainfall interception and stem flow by eucalypt street trees—the impacts of canopy density and bark type. Urban For. Urban Green. 13, 192–197. https://doi.org/10.1016/j.ufug.2013.09.001 (2014).Article 

Google Scholar 
Xiao, Q. & McPherson, E. G. Rainfall interception by Santa Monica’s municipal urban forest. Urban Ecosyst. 6, 291–302. https://doi.org/10.1023/B:UECO.0000004828.05143.67 (2002).Article 

Google Scholar 
Rohatgi, A. WebPlotDigitizer (4.4), 2020).Team, R. C. (R Foundation for Statistical Computing, 2020).García-Palacios, P., Gross, N., Gaitán, J. & Maestre, F. T. Climate mediates the biodiversity–ecosystem stability relationship globally. PNAS 115, 8400–8405. https://doi.org/10.1073/pnas.1800425115 (2018).Article 
ADS 
CAS 

Google Scholar 
Le Provost, G. et al. Land-use history impacts functional diversity across multiple trophic groups. PNAS 117, 1573–1579. https://doi.org/10.1073/pnas.1910023117 (2020).Article 
ADS 
CAS 

Google Scholar 
El Kateb, H., Zhang, H., Zhang, P. & Mosandl, R. Soil erosion and surface runoff on different vegetation covers and slope gradients: A field experiment in Southern Shaanxi Province, China. CATENA 105, 1–10. https://doi.org/10.1016/j.catena.2012.12.012 (2013).Article 

Google Scholar 
Oliveira, P. T. S. et al. The water balance components of undisturbed tropical woodlands in the Brazilian cerrado. Hydrol. Earth Syst. Sci. 19, 2899–2910. https://doi.org/10.5194/hess-19-2899-2015 (2014).Article 
ADS 

Google Scholar 
Hümann, M. et al. Identification of runoff processes – The impact of different forest types and soil properties on runoff formation and floods. J. Hydrol. 409, 637–649. https://doi.org/10.1016/j.jhydrol.2011.08.067 (2011).Article 
ADS 

Google Scholar 
Sun, D. et al. Soil erosion and water retention varies with plantation type and age. For. Ecol. Manage. 422, 1–10. https://doi.org/10.1016/j.foreco.2018.03.048 (2018).Article 

Google Scholar 
Jost, G., Schume, H., Hager, H., Markart, G. & Kohl, B. A hillslope scale comparison of tree species influence on soil moisture dynamics and runoff processes during intense rainfall. J. Hydrol. 420–421, 112–124. https://doi.org/10.1016/j.jhydrol.2011.11.057 (2012).Article 

Google Scholar 
Sadeghi, S. M. M., Attarod, P., Van Stan, J. T. & Pypker, T. G. The importance of considering rainfall partitioning in afforestation initiatives in semiarid climates: A comparison of common planted tree species in Tehran, Iran. Sci. Total Environ. 568, 845–855. https://doi.org/10.1016/j.scitotenv.2016.06.048 (2016).Article 
ADS 
CAS 

Google Scholar 
Pretzsch, H. et al. Climate change accelerates growth of urban trees in metropolises worldwide. Sci. Rep. https://doi.org/10.1038/s41598-017-14831-w (2017).Article 

Google Scholar 
Rahman, M. A., Moser, A., Rötzer, T. & Pauleit, S. Microclimatic differences and their influence on transpirational cooling of Tilia cordata in two contrasting street canyons in Munich, Germany. Agric. For. Meteorol. 232, 443–456. https://doi.org/10.1016/j.agrformet.2016.10.006 (2017).Article 
ADS 

Google Scholar 
Nytch, C. J., Meléndez-Ackerman, E. J., Pérez, M. E. & Ortiz-Zayas, J. R. Rainfall interception by six urban trees in San Juan, Puerto Rico. Urban Ecosyst. 22, 103–115. https://doi.org/10.1007/s11252-018-0768-4 (2018).Article 

Google Scholar 
Rahman, M. A. et al. Comparative analysis of shade and underlying surfaces on cooling effect. Urban For. Urban Green. 63, 127223. https://doi.org/10.1016/j.ufug.2021.127223 (2021).Article 

Google Scholar 
Chen, L., Zhang, Z. & Ewers, B. E. Urban tree species show the same hydraulic response to vapor pressure deficit across varying tree size and environmental conditions. PLoS One https://doi.org/10.1371/journal.pone.0047882 (2012).Article 

Google Scholar 
Moser-Reischl, A., Rahman, M. A., Pauleit, S., Pretzsch, H. & Rötzer, T. Growth patterns and effects of urban micro-climate on two physiologically contrasting urban tree species. Landsc. Urban Plan. 183, 88–99. https://doi.org/10.1016/j.landurbplan.2018.11.004 (2019).Article 

Google Scholar 
Hao, M. et al. Impacts of changes in vegetation on saturated hydraulic conductivity of soil in subtropical forests. Sci. Rep. 9, 8372. https://doi.org/10.1038/s41598-019-44921-w (2019).Article 
ADS 
CAS 

Google Scholar 
Peters, E. B., McFadden, J. P. & Montgomery, R. A. Biological and environmental controls on tree transpiration in a suburban landscape. J. Geophys. Res. Biogeosci. https://doi.org/10.1029/2009jg001266 (2010).Article 

Google Scholar 
Komatsu, H., Kume, T. & Otsuki, K. Increasing annual runoff—broadleaf or coniferous forests?. Hydrol. Process. 25, 302–318. https://doi.org/10.1002/hyp.7898 (2011).Article 
ADS 

Google Scholar 
Li, X. et al. Process-based rainfall interception by small trees in Northern China: The effect of rainfall traits and crown structure characteristics. Agric. For. Meteorol. 218–219, 65–73. https://doi.org/10.1016/j.agrformet.2015.11.017 (2016).Article 
ADS 

Google Scholar 
Lukaszkiewicz, J. & Kosmala, M. Determining the age of streetside trees with diameter at breast height-based multifactorial model. Arboricult. Urban For. 34, 137–143. https://doi.org/10.48044/jauf.2008.018 (2008).Article 

Google Scholar 
Buttle, J. M. & Farnsworth, A. G. Measurement and modeling of canopy water partitioning in a reforested landscape: The Ganaraska Forest, southern Ontario, Canada. J. Hydrol. 466–467, 103–114. https://doi.org/10.1016/j.jhydrol.2012.08.021 (2012).Article 

Google Scholar 
Yang, B., Lee, D. K., Heo, H. K. & Biging, G. The effects of tree characteristics on rainfall interception in urban areas. Landsc. Ecol. Eng. 15, 289–296. https://doi.org/10.1007/s11355-019-00383-w (2019).Article 
CAS 

Google Scholar 
Klamerus-Iwan, A. & Witek, W. Variability in the Wettability and Water Storage Capacity of Common Oak Leaves (Quercus robur L). Water 10, 695. https://doi.org/10.3390/w10060695 (2018).Article 
CAS 

Google Scholar 
Van Stan, J. T., Siegert, C. M., Levia, D. F. & Scheick, C. E. Effects of wind-driven rainfall on stemflow generation between codominant tree species with differing crown characteristics. Agric. For. Meteorol. 151, 1277–1286. https://doi.org/10.1016/j.agrformet.2011.05.008 (2011).Article 
ADS 

Google Scholar 
Selbig, W. R. et al. Quantifying the stormwater runoff volume reduction benefits of urban street tree canopy. Sci. Total Environ. 806, 151296. https://doi.org/10.1016/j.scitotenv.2021.151296 (2022).Article 
ADS 
CAS 

Google Scholar 
Centre for Watershed Protection. Review of the Available Literature and Data on the Runoff and Pollutant Removal Capabilities of Urban Trees (Center for Watershed Protection, 2017).
Google Scholar 
Berland, A. et al. The role of trees in urban stormwater management. Landsc. Urban Plan. 162, 167–177. https://doi.org/10.1016/j.landurbplan.2017.02.017 (2017).Article 

Google Scholar 
Pauleit, S. Urban street tree plantings: Indentifying the key requirements. Proc. Inst. Civ. Eng. Municipal Eng. 156, 43–50. https://doi.org/10.1680/muen.2003.156.1.43 (2003).Article 

Google Scholar 
Weller, M. Tree Inventory Data of Central European Cities—Studies on the Composition and Structure of Urban Tree Populations and Derivation of Ecosystem Services. MSC thesis, Technical University of Munich, Germany (2021). More

Real-world dataThe dataset consisted of 2986 observations from 902 freshwater shallow lakes in Denmark and North America (data extracted from the LAGOSNE database on 22 February 2022 via R LAGOSNE package version 2.0.2)56 (Supplementary Fig. 9). The Danish lakes were sampled for one or several years from 1984 to 2020 (data extracted in October 2021 from https://odaforalle.au.dk/main.aspx) (Supplementary Fig. 10). Prerequisites for inclusion in the analysis were that lakes had been sampled for physical and chemical variables at least four times or at least three times over the growing season (May to September) for the Danish or North American lakes, respectively, had a mean depth of less than 3 m and were freshwater. Water chemistry samples were analysed using standard methods and data for total phosphorus (TP), total nitrogen (TN) and chlorophyll-a are included here57. The mean and range of TP, TN and chlorophyll-a for the combined sites is given in Table 1, along with the values for each region separately.To gain a longer-term perspective on the relationship between nutrients and chlorophyll-a, we calculated the across-year averages of the summer means of TP, TN and chlorophyll-a, sequentially increasing numbers of years included in the mean up to a total of a five-year mean, at which point there were only 99 lakes left in the dataset. In calculating the multi-year means we allowed a maximum gap of 2 years between observations (i.e. two observations could cover 3 years) to avoid including time series with too many missing years in between. Hence, only lakes with sufficient numbers of sequential data were included, resulting in a large drop in lake numbers as the length of the multi-year mean increased (Table 2).Numerical methodsDiagnostic tests or proxies of alternative equilibriaWe modelled the response of chlorophyll-a to TP and TN using generalised-linear models58 with Gamma distribution and an identity link on untransformed data for single-year and multiple-year means up to 5-year means. We used the Gamma distribution, as chlorophyll fit this significantly better than a normal or log-normal distribution. We used psuedo R2 of the model along with the patterns of residuals, and finally, we plotted the kernel density of the chlorophyll-a values as diagnostics of the presence, absence or prevalence of alternative equilibria in the simulated and real work data.To test how appropriate these diagnostics or proxies of alternative stable states in terms of how well they identify the existence of alternative stable states in randomly sampled multi-year data, we

1.

Simulated two scenarios for the main manuscript, with and without alternative stable states in the data, which were as close to the real-life data as possible. The results of these scenarios appear in the main text (please see details below in the “Data Simulation” section).

2.

We provide multiple scenarios with different degrees, or prevalence, of alternative stable states in the data, see simulations of alternative stable state scenarios. The results of these scenarios appear in Supplementary note 2.

Hierarchical bootstrap approachThere are a large number of permutations of data, both real-world and simulated, that can provide a mean of the two to five sequential years from each lake in the time series data. It was vital to have a method that selects the data for analysis that provides a valid and comparable representation of both real work and simulated data and the models’ errors. In order to provide this we used a non-parametric hierarchical bootstrap procedure38. The flowchart shows the data preparation and data analysis steps of the hierarchical bootstrap procedure (Fig. 4). In the first step (step 1 in Fig. 4), all possible longer-term means are calculated for each lake. To keep as much data as possible, we decided to allow for up to 2 years of gap in the data between years. Taking the 5-year mean data as an example, if data from a lake existed for the years 1991 and 1994−1997, a 5-year mean would be calculated for the years 1991, 1994, 1995, 1996 and 1997. However, if the time series would contain a larger gap, e.g. data would only exist for the years 1991 and 1995–1998, no 5-year mean could be calculated. After the selection procedure, all the 2-year, 3-year and 5-year means are transferred into a new table (step 2 in Fig. 4).Fig. 4Data preparation and analysis steps of the hierarchical bootstrap procedure.Full size imageThe procedure is the same for each temporal scale from 2-year means to 5-year means. For the example of 5 mean years, lakes are randomly sampled from the full 5-year mean dataset in step 2 (Fig. 4) with replacement up to the number of lakes as in the original dataset, for the 5-year means 99 (step 3a). Here, the same lake can appear multiple times or not at all. This step is common for every bootstrap procedure59. However, since we have nested data (5-year means within lakes), we need a second step, in which for every resampled lake in step 3a, one 5-year mean is chosen (step 3b in Fig. 4). Then the three GLM models are produced from the randomly selected data in step 3c (Fig. 4). These steps are then repeated 1000 times to get a good representation of the uncertainties of the model. To ensure a fair comparison between single-year data and their equivalent multi-year mean data, we repeated the bootstrap procedure with single years only using only the lakes for which we also calculated multi-year means. To take the five-year mean as an example, there were 99 lakes where we could calculate at least one 5-year mean observation. First, we ran the bootstrap procedure to calculate 5-year mean values of TP, TN and chlorophyll-a (1000 times) and then took single years’ values of TP, TN and chlorophyll-a (1000 times) from exactly the same 99 lakes. With this approach, exactly the same datasets with the same lakes and observations within lakes are used for the calculation of the multi-year means and their single-year counterparts, making for a robust analysis. The GLM models did not always converge. If either the TP, TN or TP*TN model with interaction did not converge, the iteration was not used in further analysis. The number of converging models equal for each iteration of random samples is given in the results.The described hierarchical approach is the best way to reflect the structure of the original data. A simple, non-hierarchical bootstrap would favour lakes with more five-year means over lakes with fewer five-year means, simply because these make up a larger part of the data. Furthermore, sampling without replacement at the lake level would result in five-year means from lakes with few data dominating the produced random dataset, as every lake would be sampled every time, which then would result in high model leverage of 5-year means from lakes with less data. In contrast, the hierarchical procedure ensures that every lake has the same chance to end up in the randomly sampled bootstrap, in the second step, it ensures that of each sampled lake, every 5-year mean has the same chance to end up in the random dataset. These notions are in agreement with the findings of an assessment on how to properly resample hierarchical data by non-parametric bootstrap38.Data simulationGeneral approach of simulation assumptions and proceduresWe generated random scatter for the generalised-linear model based on Gamma distributions for two different “populations” of lakes with two different intercepts and slopes. At first, we calculated the linear equations for the two populations:For each population i and j, 99 samples (equalling the number of lakes in real-life data with 5-year means, nlake) were generated with a specific number of data points depending on the scenario (nyear) each, hence nlake = 1−99 for each population of lakes, e.g. with 20 years (nyear = 20) each.We found the real nutrient data to be normally distributed, with total nitrogen (TN) having a range between 0.33 and 4.93 mg/L and a constant coefficient of variation (CV, with a mean CV of 0.35) across this range (the same is true for total phosphorus (TP) at a shorter range). Hence, for each nlake, the x for the nyear = 20 were generated based on the mean range (mean per lake of the real-life data) and CV (0.35) from the real-life TN concentration data, hence with a range of 0.33 to 4.33 mg/L. Therefore the values and random variability of x in the simulations are close to the true values of the TN concentrations. The x is then fed into the linear equations above.To the resulting yi and yj we added random noise based on the Gamma distribution (using the rgamma function in R). We used a Gamma distribution because the Chlorophyll-a concentration also follows a Gamma distribution. The variability of a Gamma distribution is expressed by the shape variable. The variability of chlorophyll-a, its shape value, equals 2.63. This shape value was used in the Gamma distribution of yi and yj. The final calculated yi and yj had therefore a random rate calculated as shape/yi or shape/yj. Hence, their variability in the y dimension was close to the true chlorophyll-a variability.The data from both lake populations were then pooled and randomly sampled using the same hierarchical bootstrap procedure with 500 iterations for the scenarios in the supplementary materials and with 1000 iterations for main text simulation scenarios, which is identical to what was done for the real-world data.Simulation scenarios based on characteristics of real-world dataThe real-world 5-year mean data consisted of 99 lakes with 5–20 years of data for each lake. For the simulation scenario in the main text, we therefore randomly sampled between 5 and 20 data points for each of the 99 simulated lakes based on the x distribution described above. Intercepts and slopes of the simulation, resembled the range of the true data (see scatter plots in Fig. 2 of the main manuscript).In the alternative stable state scenario, we chose two slopes and intercepts for different populations of lakes:

Population i: ai = 0, bi = 40

Population j: aj = 50, bj = 120

We based the slopes and intercepts of the ASS scenarios on the diagnostic combination defined by Scheffer and Carpenter7 which propose an abrupt shift in (a) the time series, (b) the multimodal distribution of states and (c) the dual relationship to a controlling factor. Here, the idea is that an ecosystem will jump from one state to the next at the same (nutrient) conditions (different intercept and/or slope, condition a within ref. 7), where any change in the nutrient will have different effects on algae or macrophytes (best represented by different slope, condition c), resulting in a multimodal distribution of the response (condition b). Hence, simulations are in line with what is predicted for ASS, but we took great lengths to also show other possibilities with the simulations in the Supplementary information to ensure we did not overlook any occasional occurrence of alternative equilibria.Here, the appearance of alternative stable states in the data could happen at any point in the time series of a single lake, or the entire time series could include only one of the two alternative stable states. To accommodate these alternative stable state constellations (for each of which we made a separate simulation scenario, (see Supplementary Note 2, “Simulations of alternative stable state constellations”), we forced the alternative stable state scenario to be constructed of 1/3 of data with one state, 1/3 of data with the second state and 1/3 of data where both alternative states could occur. In the latter case, the alternative stable state appeared after the first 20% but before the last 20% of the time series. Since the variability and range of x (nutrient) and y (chlorophyll -a response) is simulated as close as possible to the real-world data in all scenarios, the measures taken here (variable time series and combination of different alternative stable state scenario constellations) produce a simulation as close to the real-world data as possible. Specifically, we found the real-world nutrient data to be normally distributed, with total nitrogen (TN) having a range between 0.33 and 4.93 mg/L and a constant coefficient of variation (CV, with a mean CV of 0.35) across this range (the same is true for total phosphorus (TP) at a shorter range). Hence, for each simulated lake, the x were generated based on this mean range and CV. Furthermore, the resulting yi and yj were randomised by using a Gamma distribution (using the rgamma function in R). We used a Gamma distribution because the chlorophyll-a concentration also follows a Gamma distribution. The variability of a Gamma distribution is expressed by the shape variable. The variability of chlorophyll-a, its shape value, equals 2.63. This shape value was used in the Gamma distribution of yi and y. The final calculated yi and yj had, therefore a random rate calculated as shape/yi or shape/y. Hence, their variability in the y dimension was close to the true chlorophyll-a variability.For the scenario without alternative stable states, both populations of data had the same intercept and slope:

Population i: ai = 0, bi = 40

Population j: aj = 0, bj = 40.

Please see Supplementary Note 2 for further simulations of different potential constellations of alternative states. There we show that our approach finds alternative stable states in response to nutrient concentration, even if they appear in time series from different lakes.Assessment of diagnostic tests or proxies of alternative equilibriaWe modelled the response of chlorophyll-a to TP and TN using generalised-linear models3 with Gamma distribution and an identity link on untransformed data for single-year and multiple-year means up to 5-year means. We used the Gamma distribution, as chlorophyll fit this significantly better than a normal or log-normal distribution. We used R2 of the model along with the patterns of residuals, and finally, we plotted the kernel density of the chlorophyll-a values as diagnostics of the presence, absence or prevalence of alternative equilibria in the simulated and real work data.The comparison of how the diagnostics/proxies of alternative stable states respond to the variation in the prevalence of alternative equilibria in the simulated datasets provides a robust assessment of their ability to identify both the presence and absence of alternative equilibria. It is the response of these diagnostic tests over time, with the increase in the temporal perspective as more years are added to the mean values of TP, TN and chlorophyll-a, that are key to the identification of the presence and or absence of alternative equilibria in a given dataset. The simulations show that a dataset which contains alternative equilibria will show (1) no improvement in R2 as the temporal perspective of the data increases (more years in the multi-year mean); (2) an increased bimodality in the residuals of the models of nutrients predicting chlorophyll-a will increase as more years are added to the multi-year mean and (3) the kernel density function of chlorophyll-a will display increasingly bimodality as more years are added to the mean. In the absence of alternative equilibria, the patterns differ with an R2, and increase in unimodality of residuals and a consistent unimodal pattern in the kernel density function. Thus, the diagnostic tests provide a robust test of both the presence and absence of alternative equilibria in a given dataset.Alternative stable state assessment for real data with limited data rangeIt could be the case that alternative stable states do not appear in the full dataset but only in a limited TN and TP concentration range. We filtered and re-analyzed the data, only keeping data points within the following two ranges: – TN concentration = 0.5−2 mg/L–TP concentration = 0.05−0.4 mg/L. In the filtered data, 1329 out of the original 2876 single-year data points, 289 out of 1028 3-year mean data points and 212 out of the 864 five-mean year data points remained. The filtered data consisted of data points from 550, 48 and 27 lakes for the single-year data, 3-year means and 5-year means, respectively. The smaller range resulted in lower R² of the models, yet the pattern that multi-year means result in higher R² compared to single-year data was largely consistent, apart from the 5-year mean TN models for which both, the single-year and mean data resulted in very low R² (Supplementary Fig. 6). Furthermore, due to the lower number of samples, the errors of all proxies are higher, making conclusions more difficult than for the full data. Still, we do not see any clear indication of alternative stable states in the scatter plots (two groups of dots are not appearing (Supplementary Fig. 5), the Kernel density plots (or model residuals (Supplementary Fig. 6)). i.e. no signs of bimodality in residuals or Kernel density plots. Please see details on this analysis in the supplementary material.Details and the R code for the steps for the random multi-year sampling can be found in the supplementary materials.Reporting summaryFurther information on research design is available in the Nature Portfolio Reporting Summary linked to this article. More

Speakman, J. R. The cost of living: Field metabolic rates of small mammals. Adv. Ecol. Res. 30, 177–297 (1999).Article 

Google Scholar 
Brown, J. H., Gillooly, J. F., Allen, A. P., Savage, V. M. & West, G. B. Toward a metaboolic theory of ecology. Ecology 85(7), 1771–1789. https://doi.org/10.1890/03-9000 (2004).Article 

Google Scholar 
Larivée, M. L., Boutin, S., Speakman, J. R., McAdam, A. G. & Humphries, M. M. Associations between over-winter survival and resting metabolic rate in juvenile North American red squirrels. Funct. Ecol. 24(3), 597–607. https://doi.org/10.1111/j.1365-2435.2009.01680.x (2010).Article 

Google Scholar 
Corp, N., Gorman, M. L. & Speakman, J. R. Seasonal variation in the resting metabolic rate of male wood mice Apodemus sylvaticus from two contrasting habitats 15 km apart. J. Comp. Physiol. B 167(3), 229–239. https://doi.org/10.1007/s003600050069 (1997).Article 
CAS 

Google Scholar 
Lehto Hürlimann, M., Martin, J. G. A. & Bize, P. Evidence of phenotypic correlation between exploration activity and resting metabolic rate among populations across an elevation gradient in a small rodent species. Behav. Ecol. Sociobiol. 73(9), 131. https://doi.org/10.1007/s00265-019-2740-6 (2019).Article 

Google Scholar 
Reher, S., Rabarison, H., Montero, B. K., Turner, J. M. & Dausmann, K. H. Disparate roost sites drive intraspecific physiological variation in a Malagasy bat. Oecologia 198(1), 35–52. https://doi.org/10.1007/s00442-021-05088-2 (2022).Article 
ADS 

Google Scholar 
McDonald, R. I. et al. Research gaps in knowledge of the impact of urban growth on biodiversity. Nat. Sustain. https://doi.org/10.1038/s41893-019-0436-6 (2019).Article 

Google Scholar 
Shochat, E., Warren, P. S., Faeth, S. H., McIntyre, N. E. & Hope, D. From patterns to emerging processes in mechanistic urban ecology. Trends Ecol. Evol. 21(4), 186–191. https://doi.org/10.1016/j.tree.2005.11.019 (2006).Article 

Google Scholar 
United Nations, Department of Economic and Social Affairs, Population Division. World Urbanization Prospects 2018: Highlights. https://population.un.org/wup/Publications/ (2018).Alberti, M. et al. The complexity of urban eco-evolutionary dynamics. Bioscience 70(9), 772–793. https://doi.org/10.1093/biosci/biaa079 (2020).Article 

Google Scholar 
Birnie-Gauvin, K., Peiman, K. S., Gallagher, A. J., de Bruijn, R. & Cooke, S. J. Sublethal consequences of urban life for wild vertebrates. Environ. Rev. 24(4), 416–425. https://doi.org/10.1139/er-2016-0029 (2016).Article 

Google Scholar 
Diamond, S. E. & Martin, R. A. Physiological adaptation to cities as a proxy to forecast global-scale responses to climate change. J. Exp. Biol. 224((Suppl_1)), jeb22336. https://doi.org/10.1242/jeb.229336 (2021).Article 

Google Scholar 
Grimm, N. B. et al. Global change and the ecology of cities. Science 319(5864), 756–760. https://doi.org/10.1126/science.1150195 (2008).Article 
ADS 
CAS 

Google Scholar 
McDonnell, M. J. & Pickett, S. T. Ecosystem structure and function along urban-rural gradients: An unexploited opportunity for ecology. Ecology 71(4), 1232–1237. https://doi.org/10.2307/1938259 (1990).Article 

Google Scholar 
Francis, R. A. & Chadwick, M. A. What makes a species synurbic?. Appl. Geogr. 32(2), 514–521. https://doi.org/10.1016/j.apgeog.2011.06.013 (2012).Article 

Google Scholar 
Luniak, M. Synurbization–adaptation of animal wildlife to urban development. In Proc. 4th Int. Symposium Urban Wildl. Conserv (Tucson, University of Arizona, 2004).Coogan, S. C. P., Raubenheimer, D., Zantis, S. P. & Machovsky-Capuska, G. E. Multidimensional nutritional ecology and urban birds. Ecosphere 9(4), e02177. https://doi.org/10.1002/ecs2.2177 (2018).Article 

Google Scholar 
Lowry, H., Lill, A. & Wong, B. B. Behavioural responses of wildlife to urban environments. Biol. Rev. Camb. Philos. Soc. 88(3), 537–549. https://doi.org/10.1111/brv.12012 (2013).Article 

Google Scholar 
Łopucki, R., Klich, D., Ścibior, A. & Gołębiowska, D. Hormonal adjustments to urban conditions: Stress hormone levels in urban and rural populations of Apodemus agrarius. Urban Ecosyst. 22(3), 435–442. https://doi.org/10.1007/s11252-019-0832-8 (2019).Article 

Google Scholar 
McCleery, R. in Urban mammals in Urban Ecosystem Ecology (eds. Aitkenhead-Peterson, J., Volder, A.) 87–102 (American Society of Agronomy, 2010). https://doi.org/10.2134/agronmonogr55.c52010Uchida, K., Suzuki, K., Shimamoto, T., Yanagawa, H. & Koizumi, I. Seasonal variation of flight initiation distance in Eurasian red squirrels in urban versus rural habitat. J. Zool. 298(3), 225–231. https://doi.org/10.1111/jzo.12306 (2016).Article 

Google Scholar 
Kleerekoper, L., van Esch, M. & Salcedo, T. B. How to make a city climate-proof, addressing the urban heat island effect. Resour. Conserv. Recyl. 64, 30–38. https://doi.org/10.1016/j.resconrec.2011.06.004 (2012).Article 

Google Scholar 
Pickett, S. T. et al. Urban ecological systems: Scientific foundations and a decade of progress. J. Environ. Manag. 92(3), 331–362. https://doi.org/10.1016/j.jenvman.2010.08.022 (2011).Article 
CAS 

Google Scholar 
Rizwan, A. M., Dennis, L. Y. & Chunho, L. A review on the generation, determination and mitigation of Urban Heat Island. J. Environ. Sci. 20(1), 120–128 (2008).Article 
CAS 

Google Scholar 
Isaksson, C. Urban ecophysiology: Beyond costs, stress and biomarkers. J. Exp. Biol. 223(22), jeb203794. https://doi.org/10.1242/jeb.203794 (2020).Article 

Google Scholar 
Miles, L. S., Carlen, E. J., Winchell, K. M. & Johnson, M. T. J. Urban evolution comes into its own: Emerging themes and future directions of a burgeoning field. Evol. Appl. 14(1), 3–11. https://doi.org/10.1111/eva.13165 (2020).Article 

Google Scholar 
Gavett, A. P. & Wakeley, J. S. Blood constituents and their relation to diet in urban and rural house sparrows. Condor 88(3), 279–284. https://doi.org/10.2307/1368873 (1986).Article 

Google Scholar 
Murray, M. et al. Greater consumption of protein-poor anthropogenic food by urban relative to rural coyotes increases diet breadth and potential for human-wildlife conflict. Ecography 38(12), 1235–1242. https://doi.org/10.1111/ecog.01128 (2015).Article 

Google Scholar 
Pollock, C. J., Capilla-Lasheras, P., McGill, R. A. R., Helm, B. & Dominoni, D. M. Integrated behavioural and stable isotope data reveal altered diet linked to low breeding success in urban-dwelling blue tits (Cyanistes caeruleus). Sci. Rep. 7(1), 5014. https://doi.org/10.1038/s41598-017-04575-y (2017).Article 
ADS 
CAS 

Google Scholar 
Schulte-Hostedde, A. I., Mazal, Z., Jardine, C. M. & Gagnon, J. Enhanced access to anthropogenic food waste is related to hyperglycemia in raccoons (Procyon lotor). Conserv. Physiol. 6(1), coy026. https://doi.org/10.1093/conphys/coy026 (2018).Article 
CAS 

Google Scholar 
Fingland, K., Ward, S. J., Bates, A. J. & Bremner-Harrison, S. A systematic review into the suitability of urban refugia for the Eurasian red squirrel Sciurus vulgaris. Mamm. Rev. 52(1), 26–38. https://doi.org/10.1111/mam.12264 (2021).Article 

Google Scholar 
Jokimäki, J., Selonen, V., Lehikoinen, A. & Kaisanlahti-Jokimäki, M.-L. The role of urban habitats in the abundance of red squirrels (Sciurus vulgaris, L.) in Finland. Urban For. Urban Green. 27, 100–108. https://doi.org/10.1016/j.ufug.2017.06.021 (2017).Article 

Google Scholar 
Dausmann, K. H., Wein, J., Turner, J. M. & Glos, J. Absence of heterothermy in the European red squirrel (Sciurus vulgaris). Mammal. Biol. 78(5), 332–335. https://doi.org/10.1016/j.mambio.2013.01.004 (2013).Article 

Google Scholar 
Turner, J. M., Reher, S., Warnecke, L. & Dausmann, K. H. Eurasian red squirrels show little seasonal variation in metabolism in food-enriched habitat. Physiol. Biochem. Zool. 90(6), 655–662. https://doi.org/10.1086/694847 (2017).Article 

Google Scholar 
McNab, B. K. On the comparative ecological and evolutionary significance of total and mass-specific rates of metabolism. Physiol. Biochem. Zool. 72(5), 642–644 (1999).Article 
CAS 

Google Scholar 
Menzies, A. K. et al. Body temperature, heart rate, and activity patterns of two boreal homeotherms in winter: Homeostasis, allostasis, and ecological coexistence. Funct. Ecol. 34(11), 2292–2301. https://doi.org/10.1111/1365-2435.13640 (2020).Article 

Google Scholar 
Wauters, L. & Dhondt, A. Activity budget and foraging behaviour of the red squirrel (Sciurus vulgaris Linnaeus, 1758) in a coniferous habitat. Z. Säugetierkd. 52(6), 341–353 (1987).
Google Scholar 
Wauters, L., Swinnen, C. & Dhondt, A. A. Activity budget and foraging behaviour of red squirrels (Sciurus vulgaris) in coniferous and deciduous habitats. J. Zool. 227(1), 71–86. https://doi.org/10.1111/j.1469-7998.1992.tb04345.x (1992).Article 

Google Scholar 
Reher, S., Dausmann, K. H., Warnecke, L. & Turner, J. M. Food availability affects habitat use of Eurasian red squirrels (Sciurus vulgaris) in a semi-urban environment. J. Mammal. 97(6), 1543–1554. https://doi.org/10.1093/jmammal/gyw105 (2016).Article 

Google Scholar 
Moller, H. Foods and foraging behavior of red (Sciurus vulgaris) and grey (Sciurus carolinensis) squirrels. Mammal. Rev. 13(2–4), 81–98. https://doi.org/10.1111/j.1365-2907.1983.tb00270.x (1983).Article 

Google Scholar 
Krauze-Gryz, D. & Gryz, J. in A review of the diet of the red squirrel (Sciurus vulgaris) in different types of habitats in Red squirrels: Ecology, conservation & management in Europe (eds. Shuttleworth, C. M., Lurz, P. W. W., Hayward, M. W.) 39–50 (European Squirrel Initiative, London, 2015)Shuttleworth, C. M. in The effect of supplemental feeding on the red squirrel (Sciurus vulgaris), Doctoral dissertation (University of London, London, 1996).Birnie-Gauvin, K., Peiman, K. S., Raubenheimer, D. & Cooke, S. J. Nutritional physiology and ecology of wildlife in a changing world. Conserv. Physiol. https://doi.org/10.1093/conphys/cox030 (2017).Article 

Google Scholar 
Wist, B., Stolter, C. & Dausmann, K. H. Sugar addicted in the city: Impact of urbanisation on food choice and diet composition of the Eurasian red squirrel (Sciurus vulgaris). J. Urban Ecol. 8(1), juac012. https://doi.org/10.1093/jue/juac012 (2022).Article 

Google Scholar 
Burton, T., Killen, S. S., Armstrong, J. D. & Metcalfe, N. B. What causes intraspecific variation in resting metabolic rate and what are its ecological consequences?. Proc. Biol. Sci. 278(1724), 3465–3473. https://doi.org/10.1098/rspb.2011.1778 (2011).Article 
CAS 

Google Scholar 
Clarke, A. Costs and consequences of evolutionary temperature adaptation. Trends Ecol. Evol. 18(11), 573–581. https://doi.org/10.1016/j.tree.2003.08.007 (2003).Article 

Google Scholar 
Lovegrove, B. G. The influence of climate on the basal metabolic rate of small mammals: A slow-fast metabolic continuum. J. Comp. Physiol. B 173(2), 87–112. https://doi.org/10.1007/s00360-002-0309-5 (2003).Article 
CAS 

Google Scholar 
McNab, B. K. The energetics of endotherms. Ohio J. Sci. 74(6), 370–380 (1974).
Google Scholar 
Tattersall, G. J. et al. Coping with thermal challenges: Physiological adaptations to environmental temperatures. Compr. Physiol. 2(3), 2151–2202 (2012).Article 

Google Scholar 
Broggi, J. et al. Sources of variation in winter basal metabolic rate in the great tit. Funct. Ecol. 21(3), 528–533. https://doi.org/10.1111/j.1365-2435.2007.01255.x (2007).Article 

Google Scholar 
Schlünzen, K. H., Hoffmann, P., Rosenhagen, G. & Riecke, W. Long-term changes and regional differences in temperature and precipitation in the metropolitan area of Hamburg. Int. J. Climatol. 30(8), 1121–1136. https://doi.org/10.1002/joc.1968 (2010).Article 

Google Scholar 
Reher, S. & Dausmann, K. H. Tropical bats counter heat by combining torpor with adaptive hyperthermia. Proc. R. Soc. B Biol. Sci. 288(1942), 20202059. https://doi.org/10.1098/rspb.2020.2059 (2021).Article 

Google Scholar 
Rezende, E. L. & Bacigalupe, L. D. Thermoregulation in endotherms: Physiological principles and ecological consequences. J. Comp. Physiol. B 185(7), 709–727. https://doi.org/10.1007/s00360-015-0909-5 (2015).Article 
CAS 

Google Scholar 
Scholander, P. F., Hock, R., Walters, V., Johnson, F. & Irving, L. Heat regulation in some arctic and tropical mammals and birds. Biol. Bull. 99(2), 237–258. https://doi.org/10.2307/1538741 (1950).Article 
CAS 

Google Scholar 
Terblanche, J. S., Clusella-Trullas, S., Deere, J. A., Van Vuuren, B. J. & Chown, S. L. Directional evolution of the slope of the metabolic rate-temperature relationship is correlated with climate. Physiol. Biochem. Zool. 82(5), 495–503. https://doi.org/10.1086/605361 (2009).Article 

Google Scholar 
Gallo, K. P., Easterling, D. R. & Peterson, T. C. The influence of land use/land cover on climatological values of the diurnal temperature range. J. Clim. 9(11), 2941–2944. https://doi.org/10.1175/1520-0442(1996)009%3c2941:TIOLUC%3e2.0.CO;2 (1996).Article 
ADS 

Google Scholar 
Wang, K. et al. Urbanization effect on the diurnal temperature range: Different roles under solar dimming and brightening. J. Clim. 25(3), 1022–1027. https://doi.org/10.1175/jcli-d-10-05030.1 (2012).Article 
ADS 

Google Scholar 
Fristoe, T. S. et al. Metabolic heat production and thermal conductance are mass-independent adaptations to thermal environment in birds and mammals. Proc. Natl. Acad. Sci. USA 112(52), 15934–15939. https://doi.org/10.1073/pnas.1521662112 (2015).Article 
ADS 
CAS 

Google Scholar 
Sándor, K. et al. Urban nestlings have reduced number of feathers in Great Tits (Parus major). Ibis 163(4), 1369–1378. https://doi.org/10.1111/ibi.12948 (2021).Article 

Google Scholar 
Beliniak, A., Krauze-Gryz, D., Jasińska, K., Jankowska, K. & Gryz, J. Contrast in daily activity patterns of red squirrels inhabiting urban park and urban forest. Hystrix https://doi.org/10.4404/hystrix-00476-2021 (2021).Article 

Google Scholar 
Thomas, L. S., Teich, E., Dausmann, K., Reher, S. & Turner, J. M. Degree of urbanisation affects Eurasian red squirrel activity patterns. Hystrix 29(2), 175–180. https://doi.org/10.4404/hystrix-00065-2018 (2018).Article 

Google Scholar 
Krauze-Gryz, D., Gryz, J. & Brach, M. Spatial organization, behaviour and feeding habits of red squirrels: Differences between an urban park and an urban forest. J. Zool. 315(1), 69–78. https://doi.org/10.1111/jzo.12905 (2021).Article 

Google Scholar 
Jarman, T. E., Gartrell, B. D. & Battley, P. F. Differences in body composition between urban and rural mallards Anas platyrhynchos. J. Urban Ecol. 6(1), juaa011. https://doi.org/10.1093/jue/juaa011 (2020).Article 

Google Scholar 
Cruz-Neto, A. P. & Bozinovic, F. The relationship between diet quality and basal metabolic rate in endotherms: Insights from intraspecific analysis. Physiol. Biochem. Zool. 77(6), 877–889 (2004).Article 

Google Scholar 
Geluso, K. & Hayes, J. P. Effects of dietary quality on basal metabolic rate and internal morphology of European starlings (Sturnus vulgaris). Physiol. Biochem. Zool. 72(2), 189–197 (1999).Article 
CAS 

Google Scholar 
Seebacher, F. Is endothermy an evolutionary by-product?. Trends Ecol. Evol. 35(6), 503–511. https://doi.org/10.1016/j.tree.2020.02.006 (2020).Article 

Google Scholar 
Perissinotti, P. P., Antenucci, C. D., Zenuto, R. & Luna, F. Effect of diet quality and soil hardness on metabolic rate in the subterranean rodent Ctenomys talarum. Comp. Biochem. Physiol. Mol. Integr. Physiol. 154(3), 298–307. https://doi.org/10.1016/j.cbpa.2009.05.013 (2009).Article 
CAS 

Google Scholar 
Thorp, C. R., Ram, P. K. & Florant, G. L. Diet alters metabolic rate in the yellow-bellied marmot (Marmota flaviventris) during hibernation. Physiol. Zool. 67(5), 1213–1229. https://doi.org/10.1086/physzool.67.5.30163890 (1994).Article 

Google Scholar 
Silva, S. I., Jaksic, F. M. & Bozinovic, F. Interplay between metabolic rate and diet quality in the South American fox Pseudalopex culpaeus. Comp. Biochem. Physiol. Mol Integr. Physiol. 137(1), 33–38. https://doi.org/10.1016/j.cbpb.2003.09.022 (2004).Article 
CAS 

Google Scholar 
Rewkiewicz-Dziarska, A., Wielopolska, A. & Gill, J. Hematological indices of Apodemus agrarius (Pallas, 1771) from different urban environments. Bull. Acad. Polon. Sci. Ser. Sci. Biol. 25(4), 261–268 (1977).CAS 

Google Scholar 
Ohrnberger, S. A., Hambly, C., Speakman, J. R. & Valencak, T. G. Limits to sustained energy intake XXXII: Hot again: Dorsal shaving increases energy intake and milk output in golden hamsters (Mesocricetus auratus). J Exp. Biol. https://doi.org/10.1242/jeb.230383 (2020).Article 

Google Scholar 
Speakman, J. R. & Król, E. The heat dissipation limit theory and evolution of life histories in endotherms—Time to dispose of the disposable soma theory?. Integr. Comp. Biol. 50(5), 793–807. https://doi.org/10.1093/icb/icq049 (2010).Article 

Google Scholar 
Diamond, S. E., Chick, L. D., Perez, A., Strickler, S. A. & Martin, R. A. Evolution of thermal tolerance and its fitness consequences: Parallel and non-parallel responses to urban heat islands across three cities. Proc. R. Soc. B Biol. Sci. 285(1882), 20180036. https://doi.org/10.1098/rspb.2018.0036 (2018).Article 

Google Scholar 
Isaksson, C. & Hahs, A. Urbanization, oxidative stress and inflammation: A question of evolving, acclimatizing or coping with urban environmental stress. Funct. Ecol. 29(7), 913–923. https://doi.org/10.1111/1365-2435.12477 (2015).Article 

Google Scholar 
Sokolova, I. M. & Lannig, G. Interactive effects of metal pollution and temperature on metabolism in aquatic ectotherms: Implications of global climate change. Clim. Res. 37(2–3), 181–201 (2008).Article 

Google Scholar 
Carey, H. V., Andrews, M. T. & Martin, S. L. Mammalian hibernation: Cellular and molecular responses to depressed metabolism and low temperature. Physiol. Rev. 83(4), 1153–1181 (2003).Article 
CAS 

Google Scholar 
Pereira, M. E., Aines, J. & Scheckter, J. L. Tactics of heterothermy in eastern gray squirrels (Sciurus carolinensis). J. Mammal. 83(2), 467–477 (2002).Article 

Google Scholar 
Breuner, C. W., Wingfield, J. C. & Romero, L. M. Diel rhythms of basal and stress-induced corticosterone in a wild, seasonal vertebrate. Gambel’s white-crowned sparrow. J Exp. Zool. 284(3), 334–342. https://doi.org/10.1002/(SICI)1097-010X(19990801)284:3%3c334::AID-JEZ11%3e3.0.CO;2-# (1999).Article 
CAS 

Google Scholar 
Careau, V., Thomas, D., Humphries, M. M. & Réale, D. Energy metabolism and animal personality. Oikos 117(5), 641–653. https://doi.org/10.1111/j.0030-1299.2008.16513.x (2008).Article 

Google Scholar 
Fletcher, Q. E. et al. Seasonal stage differences overwhelm environmental and individual factors as determinants of energy expenditure in free-ranging red squirrels. Funct. Ecol. 26(3), 677–687. https://doi.org/10.1111/j.1365-2435.2012.01975.x (2012).Article 

Google Scholar 
Barthel, L. & Berger, A. Unexpected gene-flow in urban environments: The example of the European Hedgehog. Animals 10(12), 2315. https://doi.org/10.3390/ani10122315 (2020).Article 

Google Scholar 
Fusco, N. A., Carlen, E. J. & Munshi-South, J. Urban landscape genetics: are biologists keeping up with the pace of urbanization?. Current Landsc. Ecol. Rep. 6(2), 35–45. https://doi.org/10.1007/s40823-021-00062-3 (2021).Article 

Google Scholar 
Ziege, M. et al. Population genetics of the European rabbit along a rural-to-urban gradient. Sci. Rep. 10(1), 2448. https://doi.org/10.1038/s41598-020-57962-3 (2020).Article 
ADS 
CAS 

Google Scholar 
Morash, A. J., Neufeld, C., MacCormack, T. J. & Currie, S. The importance of incorporating natural thermal variation when evaluating physiological performance in wild species. J. Exp. Biol. 221(14), jeb164673. https://doi.org/10.1242/jeb.164673 (2018).Article 

Google Scholar 
Pörtner, H.-O., et al. Climate change 2022: Impacts, adaptation and vulnerability. IPCC Sixth Assessment Report (2022).Anderies, J. M., Katti, M. & Shochat, E. Living in the city: Resource availability, predation, and bird population dynamics in urban areas. J. Theor. Biol. 247(1), 36–49. https://doi.org/10.1016/j.jtbi.2007.01.030 (2007).Article 
ADS 
MATH 

Google Scholar 
Shochat, E. Credit or debit? Resource input changes population dynamics of city-slicker birds. Oikos 106(3), 622–626. https://doi.org/10.1111/j.0030-1299.2004.13159.x (2004).Article 

Google Scholar 
Koprowski, J. L. Handling tree squirrels with a safe and efficient restraint. Wildl. Soc. B 30(1), 101–103. https://doi.org/10.2307/3784642 (2002).Article 

Google Scholar 
Magris, L. & Gurnell, J. Population ecology of the red squirrel (Sciurus vulgaris) in a fragmented woodland ecosystem on the Island of Jersey Channel Islands. J. Zool. 256(1), 99–112. https://doi.org/10.1017/s0952836902000134 (2002).Article 

Google Scholar 
Bethge, J., Wist, B., Stalenberg, E. & Dausmann, K. Seasonal adaptations in energy budgeting in the primate Lepilemur leucopus. J Comp. Physiol. B 187(5–6), 827–834. https://doi.org/10.1007/s00360-017-1082-9 (2017).Article 

Google Scholar 
Dausmann, K. H., Glos, J. & Heldmaier, G. Energetics of tropical hibernation. J Comp. Physiol. B 179(3), 345–357. https://doi.org/10.1007/s00360-008-0318-0 (2009).Article 
CAS 

Google Scholar 
Kobbe, S., Nowack, J. & Dausmann, K. H. Torpor is not the only option: Seasonal variations of the thermoneutral zone in a small primate. J. Comp. Physiol. B 184(6), 789–797. https://doi.org/10.1007/s00360-014-0834-z (2014).Article 

Google Scholar 
Lighton, J. R. Measuring Metabolic Rates: A Manual for Scientists (Oxford University Press, 2018).Book 

Google Scholar 
Bethge, J., Razafimampiandra, J. C., Wulff, A. & Dausmann, K. H. Sportive lemurs elevate their metabolic rate during challenging seasons and do not enter regular heterothermy. Conserv. Physiol. 9(1), coab075. https://doi.org/10.1093/conphys/coab075 (2021).Article 

Google Scholar 
Reher, S., Ehlers, J., Rabarison, H. & Dausmann, K. H. Short and hyperthermic torpor responses in the Malagasy bat Macronycteris commersoni reveal a broader hypometabolic scope in heterotherms. J. Comp. Physiol. B 188(6), 1015–1027. https://doi.org/10.1007/s00360-018-1171-4 (2018).Article 
CAS 

Google Scholar 
Grolemund, G. & Wickham, H. Dates and times made easy with lubridate. J Stat. Softw. 40(3), 1–25 (2011).Article 

Google Scholar 
Wickham, H., François, R., Henry, L. & Müller, K. RStudio. dplyr: A Grammar of Data Manipulation (1.0. 7) (2021).Zeileis, A. & Grothendieck, G. zoo: S3 infrastructure for regular and irregular time series. J. Stat. Softw. 14(6), 1–27. https://doi.org/10.18637/jss.v014.i06 (2005).Article 

Google Scholar 
Sarkar, D. Lattice: Multivariate Data Visualization with R (Springer Science & Business Media, New York, 2008).Book 
MATH 

Google Scholar 
Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. https://doi.org/10.18637/jss.v067.i01 (2015).Article 

Google Scholar 
Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. lmerTest package: Tests in linear mixed effects models. J. Stat. Softw. 82(13), 1–26 (2017).Article 

Google Scholar 
Wickham, H. ggplot2: Elegant graphics for data analysis (Springer, 2016).Book 
MATH 

Google Scholar 
Fox, J. Effect displays in R for generalised linear models. J. Stat. Softw. 8(15), 1–27 (2003).Article 

Google Scholar 
Garamszegi, L. Z. et al. Changing philosophies and tools for statistical inferences in behavioral ecology. Behav. Ecol. 20(6), 1363–1375. https://doi.org/10.1093/beheco/arp137 (2009).Article 

Google Scholar 
Symonds, M. R. E. & Moussalli, A. A brief guide to model selection, multimodel inference and model averaging in behavioural ecology using Akaike’s information criterion. Behav. Ecol. Sociobiol. 65(1), 13–21. https://doi.org/10.1007/s00265-010-1037-6 (2010).Article 

Google Scholar 
Whittingham, M. J., Stephens, P. A., Bradbury, R. B. & Freckleton, R. P. Why do we still use stepwise modelling in ecology and behaviour?. J. Anim. Ecol. 75(5), 1182–1189. https://doi.org/10.1111/j.1365-2656.2006.01141.x (2006).Article 

Google Scholar 
Barton, K. & Barton, M. K. MuMIn: Multi-Model Inference. R package version 1.43.17; https://CRAN.R-project.org/package=MuMIn (2020).Zuur, A., Ieno, E. N., Walker, N., Saveliev, A. A. & Smith, G. M. Mixed effects models and extensions in ecology with R ( Springer Science & Business Media 2009).Burnham, K. P. & Anderson, D. R. Multimodel inference: Understanding AIC and BIC in model selection. Soc. Method. Res. 33(2), 261–304 (2004).Article 

Google Scholar 
Johnson, J. B. & Omland, K. S. Model selection in ecology and evolution. Trends Ecol. Evol. 19(2), 101–108. https://doi.org/10.1016/j.tree.2003.10.013 (2004).Article 

Google Scholar 
Lorah, J. Effect size measures for multilevel models: Definition, interpretation, and TIMSS example. Large-scale Assess. Educ. 6(1), 8. https://doi.org/10.1186/s40536-018-0061-2 (2018).Article 

Google Scholar 
Selya, A. S., Rose, J. S., Dierker, L. C., Hedeker, D. & Mermelstein, R. J. A practical guide to calculating cohen’s f2, a measure of local effect size, from PROC MIXED. Front. Psychol. 3, 111–111. https://doi.org/10.3389/fpsyg.2012.00111 (2012).Article 

Google Scholar 
Lüdecke, D. sjPlot: Data visualization for statistics in social science. R package version 2.8.5 2020; https://CRAN.R-project.org/package=sjPlot (2020).Nakagawa, S., Johnson, P. C. & Schielzeth, H. The coefficient of determination R2 and intra-class correlation coefficient from generalized linear mixed-effects models revisited and expanded. J. R. Soc. Interface 14(134), 20170213 (2017).Article 

Google Scholar 
Nakagawa, S. & Schielzeth, H. Repeatability for Gaussian and non-Gaussian data: A practical guide for biologists. Biol. Rev. Camb. Philos. Soc. 85(4), 935–956. https://doi.org/10.1111/j.1469-185X.2010.00141.x (2010).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(11), 1639–1644. https://doi.org/10.1111/2041-210X.12797 (2017).Article 

Google Scholar  More

To calculate total root biomass C colonized by AM and EcM fungi, we developed a workflow that combines multiple publicly available datasets to ultimately link fine root stocks to mycorrhizal colonization estimates (Fig. 1). These estimates were individually derived for 881 different spatial units that were constructed by combining 28 different ecoregions, 15 land cover types and six continents. In a given spatial unit, the relationship between the proportion of AM- and EcM-plants aboveground biomass and the proportion of AM- and EcM-associated root biomass depends on the prevalence of distinct growth forms. Therefore, to increase the accuracy of our estimates, calculations were made separately for woody and herbaceous vegetation and combined in the final step and subsequently mapped. Below we detail the specific methodologies we followed within the workflow and the main assumptions and uncertainties associated.Fig. 1Workflow used to create maps of mycorrhizal fine root biomass carbon. The workflow consists of two main steps: (1) Estimation of total fine root stock capable to form mycorrhizal associations with AM and EcM fungi and (2) estimation of the proportion of fine roots colonized by AM and EcM fungi.Full size imageDefinition of spatial unitsAs a basis for mapping mycorrhizal root abundances at a global scale, we defined spatial units based on a coarse division of Bailey’s ecoregions23 After removing regions of permanent ice and water bodies, we included 28 ecoregions defined according to differences in climatic regimes and elevation (deposited at Dryad-Table S1). A map of Bailey’s ecoregions was provided by the Oak Ridge National Laboratory Distributed Active Archive Center24 at 10 arcmin spatial resolution. Due to potential considerable differences in plant species identities, ecoregions that extended across multiple continents were split for each continent. The continent division was based upon the FAO Global Administrative Unit Layers (http://www.fao.org/geonetwork/srv/en/). Finally, each ecoregion-continent combination was further divided according to differences in land cover types using the 2015 Land Cover Initiative map developed by the European Space Agency at 300 m spatial resolution (https://www.esa-landcover-cci.org/). To ensure reliability, non-natural areas (croplands and urban areas), bare areas and water bodies were discarded (Table 1). In summary, a combination of 28 ecoregions, 15 land cover types and six continents were combined to define a total of 881 different spatial units (deposited at Dryad-Table S2). The use of ecoregion/land cover/continent combination provided a much greater resolution than using a traditional biome classification and allowed to account for human-driven transformations of vegetation, the latter based on the land cover data.Table 1 List of land cover categories within the ESA CCI Land Cover dataset, used to assemble maps of mycorrhizal root biomass.Full size tableMycorrhizal fine root stocksTotal root C stocksEstimation of the total root C stock in each of the spatial units was obtained from the harmonized belowground biomass C density maps of Spawn et al.20. These maps are based on continental-to-global scale remote sensing data of aboveground biomass C density and land cover-specific root-to-shoot relationships to generate matching belowground biomass C maps. This product is the best up-to-date estimation of live root stock available. For subsequent steps in our workflow, we distinguished woody and herbaceous belowground biomass C as provided by Spawn et al.20. As the tundra belowground biomass C map was provided without growth form distinction, it was assessed following a slightly different workflow (see Section 2.2.3 for more details). To match the resolution of other input maps in the workflow, all three belowground biomass C maps were scaled up from the original spatial resolution of 10-arc seconds (approximately 300 m at the equator) to 10 arc‐minutes resolution (approximately 18.5 km at the equator) using the mean location of the raster cells as aggregation criterion.As the root biomass C maps do not distinguish between fine and coarse roots and mycorrhizal fungi colonize only the fine fractions of the roots, we considered the fine root fraction to be 88,5% and 14,1% of the total root biomass for herbaceous and woody plants, respectively. These constants represent the mean value of coarse/fine root mass ratios of herbaceous and woody plants provided by the Fine-Root Ecology Database (FRED) (https://roots.ornl.gov/)25 (deposited at Dryad-Table S3). Due to the non-normality of coarse/fine root mass ratios, mean values were obtained from log-transformed data and then back-transformed for inclusion into the workflow.Finally, the belowground biomass C maps consider the whole root system, but mycorrhizal colonization occurs mainly in the upper 30 cm of the soil18. Therefore, we estimated the total fine root stocks in the upper 30 cm by applying the asymptotic equation of vertical root distribution developed by Gale & Grigal26:$$y=1-{beta }^{d}$$where y is the cumulative root fraction from the soil surface to depth d (cm), and β is the fitted coefficient of extension. β values of trees (β = 0.970), shrubs (β = 0.978) and herbs (β = 0.952) were obtained from Jackson et al.27. A mean value was then calculated for trees and shrubs to obtain a woody vegetation β value of 0.974. As a result, we estimated that 54.6% of the total live root of woody vegetation and 77.1% of herbaceous vegetation is stored in the upper 30 cm of the soil. In combination, this allowed deriving fine root C stocks in the upper 30 cm of woody and herbaceous vegetation.The proportion of root stocks colonized by AM and EcMThe proportion of root stock that forms associations with AM or EcM fungi was obtained from the global maps of aboveground biomass distribution of dominant mycorrhizal types published by Soudzilovskaia et al.14. These maps provide the relative abundance of EcM and AM plants based on information about the biomass of grass, shrub and tree vegetation at 10arcmin resolution. To match with belowground root woody plants biomass data, proportions of AM trees and shrubs underlying the maps of Soudzilovskaia et al.14 were summed up to obtain the proportion of AM woody vegetation. The same was done for EcM trees and shrubs.Our calculations are subjected to the main assumption that, within each growth form, the proportion of aboveground biomass associated with AM and EcM fungi reflects the proportional association of AM and EM fungi to belowground biomass. We tested whether root:shoot ratios were significantly different between AM and EcM woody plants (the number of EcM herbaceous plants is extremely small17). Genera were linked to growth form based on the TRY database (https://www.try-db.org/)19 and the mycorrhizal type association based on the FungalRoots database17. Subsequently, it was tested whether root:shoot ratios of genera from the TRY database (https://www.try-db.org/)19 were significantly different for AM vs EcM woody plants. No statistically significant differences (ANOVA-tests p-value = 0.595) were found (Fig. 2).Fig. 2Mean and standar error of root to shoot ratios of AM and EcM woody plant species.Full size imageEstimation of mycorrhizal fine root stocksWe calculated the total biomass C of fine roots that can potentially be colonized by AM or EcM fungi by multiplying the total woody and herbaceous fine root C biomass in the upper 30 cm of the soil by the proportion of AM and EcM of woody and herbaceous vegetation. In the case of tundra vegetation, fine root C stocks were multiplied by the relative abundance of AM and EcM vegetation without distinction of growth forms (for simplicity, this path was not included in Fig. 1, but can be seen in Fig. 3. As tundra vegetation consists mainly of herbs and small shrubs, the distinction between woody and herbaceous vegetation is not essential in this case.Fig. 3Workflow used to create mycorrhizal fine root biomass C maps specific for tundra areas.Full size imageFinally, we obtained the mean value of mycorrhiza growth form fine root C stocks in each of the defined spatial units. These resulted in six independent estimations: AM woody, AM herbaceous, EcM woody, EcM herbaceous, AM tundra and EcM tundra total fine root biomass C (Fig. 4).Fig. 4Fine root biomass stocks capable to form association with AM (a) and EcM (b) fungi for woody, herbaceous and tundra vegetation. Final AM and EcM stock result from the sum of the growth form individual maps. There were no records of fine root biomass of EcM herbaceous vegetation.Full size imageThe intensity of root colonization by mycorrhizal fungiColonization databaseThe FungalRoot database is the largest up-to-date compilation of intensity of root colonization data, providing 36303 species observations for 14870 plant species. Colonization data was filtered to remove occurrences from non-natural conditions (i.e., from plantations, nurseries, greenhouses, pots, etc.) and data collected outside growing seasons. Records without explicit information about habitat naturalness and growing season were maintained as colonization intensity is generally recorded in the growing season of natural habitats. When the intensity of colonization occurrences was expressed in categorical levels, they were converted to percentages following the transformation methods stated in the original publications. Finally, plant species were distinguished between woody and herbaceous species using the publicly available data from TRY (https://www.trydb.org/)19. As a result, 9905 AM colonization observations of 4494 species and 521 EcM colonization observations of 201 species were used for the final calculations (Fig. 5).Fig. 5Number of AM (a) and EcM (b) herbaceous and woody plant species and total observations obtained from FungalRoot database.Full size imageThe use of the mean of mycorrhizal colonization intensity per plant species is based on two main assumptions:

1)

The intensity of root colonization is a plant trait: It is known that the intensity of mycorrhizal infections of a given plant species varies under different climatic and soil conditions28,29, plant age30 and the identity of colonizing fungal species31. However, Soudzilovskaia et al.9 showed that under natural growth conditions the intraspecific variation of root mycorrhizal colonization is lower than interspecific variation, and is within the range of variations in other plant eco-physiological traits. Moreover, recent literature reported a positive correlation between root morphological traits and mycorrhizal colonization, with a strong phylogenetic signature of these correlations32,33. These findings provide support for the use of mycorrhizal root colonization of plants grown in natural conditions as a species-specific trait.

2)

The percentage of root length or root tips colonized can be translated to the percentage of biomass colonized: intensity of root colonization is generally expressed as the proportion of root length colonized by AM fungi or proportion of root tips colonized by EcM fungi (as EcM infection is restricted to fine root tips). Coupling this data with total root biomass C stocks requires assuming that the proportion of root length or proportion of root tips colonized is equivalent to the proportion of root biomass colonized. While for AM colonization this equivalence can be straightforward, EcM colonization can be more problematic as the number of root tips varies between tree species. However, given that root tips represent the terminal ends of a root network34, the proportion of root tips colonized by EcM fungi can be seen as a measurement of mycorrhizal infection of the root system and translated to biomass independently of the number of root tips of each individual. Yet, it is important to stress that estimations of fine root biomass colonized by AM and EcM as provided in this paper might not be directly comparable.

sPlot databaseThe sPlotOpen database21 holds information about the relative abundance of vascular plant species in 95104 different vegetation plots spanning 114 countries. In addition, sPlotOpen provides three partially overlapping resampled subset of 50000 plots each that has been geographically and environmentally balanced to cover the highest plant species variability while avoiding rare communities. From these three available subsets, we selected the one that maximizes the number of spatial units that have at least one vegetation plot. We further checked if any empty spatial unit could be filled by including sPlot data from other resampling subsets.Plant species in the selected subset were classified as AM and EcM according to genus-based mycorrhizal types assignments, provided in the FungalRoot database17. Plant species that could not be assigned to any mycorrhizal type were excluded. Facultative AM species were not distinguished from obligated AM species, and all were considered AM species. The relative abundance of species with dual colonization was treated as 50% AM and 50% ECM. Plant species were further classified into woody and herbaceous species using the TRY database.Estimation of the intensity of mycorrhizal colonizationThe percentage of AM and EcM root biomass colonized per plant species was spatially upscaled by inferring the relative abundance of AM and EcM plant species in each plot. For each mycorrhizal-growth form and each vegetation plot, the relative abundance of plant species was determined to include only the plant species for which information on the intensity of root colonization was available. Then, a weighted mean intensity of colonization per mycorrhizal-growth form was calculated according to the relative abundance of the species featuring that mycorrhizal-growth form in the vegetation plot. Lastly, the final intensity of colonization per spatial unit was calculated by taking the mean value of colonization across all plots within that spatial unit. These calculations are based on 38127 vegetation plots that hold colonization information, spanning 384 spatial units.The use of vegetation plots as the main entity to estimate the relative abundance of AM and EcM plant species in each spatial unit assumes that the plant species occurrences and their relative abundances in the selected plots are representative of the total spatial unit. This is likely to be true for spatial units that are represented by a high number of plots. However, in those spatial units where the number of plots is low, certain vegetation types or plant species may be misrepresented. We addressed this issue in our uncertainty analysis. Details are provided in the Quality index maps section.Final calculation and maps assemblyThe fraction of total fine root C stocks that is colonized by AM and EcM fungi was estimated by multiplying fine root C stocks by the mean root colonization intensity in each spatial unit. This calculation was made separately for tundra, woody and herbaceous vegetation.To generate raster maps based on the resulting AM and EcM fine root biomass C data, we first created a 10 arcmin raster map of the spatial units. To do this, we overlaid the raster map of Bailey ecoregions (10 arcmin resolution)24, the raster of ESA CCI land cover data at 300 m resolution aggregated to 10 arcmin using a nearest neighbour approach (https://www.esa-landcover-cci.org/) and the FAO polygon map of continents (http://www.fao.org/geonetwork/srv/en/), rasterized at 10 arcmin. Finally, we assigned to each pixel the corresponding biomass of fine root colonized by mycorrhiza, considering the prevailing spatial unit. Those spatial units that remained empty due to lack of vegetation plots or colonization data were filled with the mean value of the ecoregion x continent combination.Quality index mapsAs our workflow comprises many different data sources and the extracted data acts in distinct hierarchical levels (i.e plant species, plots or spatial unit level), providing a unified uncertainty estimation for our maps is particularly challenging. Estimates of mycorrhizal fine root C stocks are related mainly to belowground biomass C density maps and mycorrhizal aboveground biomass maps, which have associated uncertainties maps provided by the original publications. In contrast, estimates of the intensity of root colonization in each spatial unit have been associated with three main sources of uncertainties:

1)

The number of observations in the FungalRoot database. The mean species-level intensity of mycorrhizal colonization in the vegetation plots has been associated with a number of independent observations of root colonization for each plant species. We calculated the mean number of observations of each plant species for each of the vegetation plots and, subsequently the mean number of observations (per plant species) from all vegetation plots in each spatial unit. These spatial unit averaged number of observations ranged from 1 to 14 in AM and from 1 to 26 in EcM. A higher number of observations would indicate that the intraspecific variation in the intensity of colonization is better captured and, therefore, the species-specific colonization estimates are more robust.

2)

The relative plant coverage that was associated with colonization data. From the selected vegetation plots, only a certain proportion of plant species could be associated with the intensity of colonization data in FungalRoot database. The relative abundance of the plant species with colonization data was summed up in each vegetation plot. Then, we calculated the average values for each spatial unit. Mean abundance values ranged from 0.3 to 100% in both AM and EcM spatial units. A high number indicates that the dominant plant species of the vegetation plots have colonization data associated and, consequently, the community-averaged intensity of colonization estimates are more robust.

3)

The number of vegetation plots in each spatial unit. Each of the spatial units differs in the number of plots used to calculate the mean intensity of colonization, ranging from 1 to 1583 and from 1 to 768 plots in AM and EcM estimations, respectively. A higher number of plots is associated with a better representation of the vegetation variability in the spatial units, although this will ultimately depend on plot size and intrinsic heterogeneity (i.e., a big but homogeneous spatial unit may need fewer vegetation plots for a good representation than a small but very heterogeneous spatial unit).

We provide independent quality index maps of the spatial unit average of these three sources of uncertainty. These quality index maps can be used to locate areas where our estimates have higher or lower robustness. More