Philippa Kaur

Brondizio, E. S., Settele, J., Díaz, S. & Ngo, H. T. (eds). Global assessment report on biodiversity and ecosystem services of the intergovernmental science-policy platform on biodiversity and ecosystem services. https://doi.org/10.5281/zenodo.3831673 (IPBES Secretariat, 2019).UNEP/FAO. The UN Decade on Ecosystem Restoration 2021-2030 “Prevent, halt and reverse the degradation of ecosystems worldwide.” https://www.decadeonrestoration.org/ (2020).Fischer, J., Riechers, M., Loos, J., Martin-Lopez, B. & Temperton, V. M. Making the UN decade on ecosystem restoration a social-ecological endeavour. Trends Ecol. Evol. 36, 1 (2021).
Google Scholar 
Tolvanen, A. & Aronson, J. Ecological reastoration, ecosystem services, and land use: a European perspective. Ecol. Soc. 21, 47 (2016).
Google Scholar 
Strassburg, B. B. N. et al. Global priority areas for ecosystem restoration. Nature 586, 724–729 (2020).CAS 
PubMed 

Google Scholar 
Temperton, V. M. et al. Step back from the forest and step up to the Bonn Challenge: how a broad ecological perspective can promote successful landscape restoration. Restor. Ecol. 27, 705–719 (2019).
Google Scholar 
Prach, K. & Hobbs, R. J. Spontaneous succession versus technical reclamation in the restoration of disturbed sites. Restor. Ecol. 16, 363–366 (2008).
Google Scholar 
Prach, K., Šebelíková, L., Řehounková, K. & del Moral, R. Possibilities and limitations of passive restoration of heavily disturbed sites. Landsc. Res. 45, 247–253 (2019).
Google Scholar 
Gilby, B. L. et al. Applying systematic conservation planning to improve the allocation of restoration actions at multiple spatial scales. Restor. Ecol. 29, e13403 (2021).
Google Scholar 
Erdős, L. et al. The edge of two worlds: a new review and synthesis on Eurasian forest-steppes. Appl. Veg. Sci. 21, 345–362 (2018).
Google Scholar 
Poschlod, P. & WallisDeVries, M. F. The historical and socioeconomic perspective of calcareous grasslands. Lessons learnt from the distant and recent past. Biol. Conserv. 104, 361–376 (2022).
Google Scholar 
Wesche, K. et al. The Palaearctic steppe biome: a new synthesis. Biodivers. Conserv. 25, 2197–2231 (2016).
Google Scholar 
Butaye, J., Dries, A. & Honnay, O. Conservation and restoration of calcareous grasslands: a concise review of the effects of fragmentation and management on plant species. Biotechnol. Agron. Soc. Environ. 9, 111–118 (2005).
Google Scholar 
Strassburg, B. B. N. et al. Strategic approaches to restoring ecosystems can triple conservation gains and halve costs. Nat. Ecol. Evol. 3, 62–70 (2019).PubMed 

Google Scholar 
Knight, M. L. & Overbeck, G. E. How much does is cost to restore a grassland? Restor. Ecol. 29, e13463 (2021).
Google Scholar 
Albert, Á.-J. et al. Trait-based analysis of spontaneous grassland recovery in sandy old-fields. Appl. Veg. Sci. 17, 214–224 (2014).
Google Scholar 
Crouzeilles, R. et al. Achieving cost-effective landscape-scale forest restoration through targeted natural regeneration. Conserv. Lett. 13, e12709 (2020).
Google Scholar 
Seregélyes, T., Molnár, Z. S., Csomós, Á. & Bölöni, J. Regeneration potential of the Hungarian (semi)-natural habitats I. Concepts and basic data of the MÉTA database. Acta Bot. Hung. 50, 229–248 (2008).
Google Scholar 
Käyhkö, N. & Skånes, H. Change trajectories and key biotopes – Assessing landscape dynamics and sustainability. Landsc. Urban Plan 75, 300–321 (2006).
Google Scholar 
Käyhkö, N. & Skånes, H. Retrospective land cover/land use change trajectories as drivers behind the local distribution and abundance patterns of oaks in south-western Finland. Landsc. Urban Plan 88, 12–22 (2008).
Google Scholar 
Swetnam, R. D. Rural land use in England and Wales between 1930 and 1998: Mapping trajectories of change with a high resolution spatio-temporal dataset. Landsc. Urban Plan 81, 91–103 (2007).
Google Scholar 
Ruiz, J. & Domon, G. 2009. Analysis of landscape pattern change trajectories within areas of intensive agricultural use: case study in a watershed of southern Québec, Canada. Landsc. Ecol. 24, 419–432 (2009).
Google Scholar 
Eremiášová, R. & Skokanová, H. Land use changes (recorded in old maps) and delimitation of the most stable areas from the perspective of land use in the Kašperské Hory region. Landsc. Ecol. 88, 20–34 (2009).
Google Scholar 
Frondoni, R. B. M. & Capotorti, G. A landscape analysis of land cover change in the Municipality of Rome (Italy): spatio-temporal characteristics and ecological implications of land cover transitions from 1954 to 2001. Landsc. Urban Plan 100, 117–128 (2011).
Google Scholar 
Biró, M., Szitár, K., Horváth, F., Bagi, I. & Molnár, Z. S. Detection of long-term landscape changes and trajectories in a Pannonian sand region: comparing land-cover and habitat-based approaches at two spatial scales. Community Ecol. 14, 219–230 (2013).
Google Scholar 
Molnár, Z. S, Biró, M., Bartha, S. & Fekete, G. in Eurasian Steppes. Ecological Problems and Livelihoods in a Changing World (eds Werger, M. J. A. & van Staalduinen, M. A.) Ch. 7 (Springer, 2012).Mezősi, G. in The Physical Geography of Hungary. Geography of the Physical Environment (ed. Mezősi, G) Ch. 11 (Springer, 2017).Biró, M., Bölöni, J. & Molnár, Z. Use of long-term data to evaluate loss and endangerment status of Natura 2000 habitats and effects of protected areas. Conserv. Biol. 32, 660–671 (2018).PubMed 

Google Scholar 
Pe’er, G. et al. Action needed for the EU Common Agricultural Policy to address sustainability challenges. People Nat. 2, 305–316 (2020).
Google Scholar 
Benton, T. G., Bieg, C., Harwatt, H., Pudasaini, R. & Wellesley, L. Food system impacts on biodiversity loss. Three levers for food system transformation in support of nature. Chatham House, the Royal Institute of International Affairs. ISBN: 978 1 78413 433 4 (2021).Kuemmerle, T. et al. Cross-border comparison of post-socialist farmland abandonment in the Carpathians. Ecosystems 11, 614 (2008).
Google Scholar 
Feranec, J. et al. Inventory of major landscape changes in the Czech Republic, Hungary, Romania and Slovak Republic 1970s – 1990s. Int. J. Appl. Earth Observ. Geoinf. 2, 129–139 (2000).
Google Scholar 
Pyšek, P. et al. Scientists’ warning on invasive alien species. Biol. Rev. 95, 1511–1534 (2020).PubMed 

Google Scholar 
Csákvári, E. et al. Conservation biology research priorities for 2050: a Central-Eastern European perspective. Biol. Conserv. 264, 109396 (2021).
Google Scholar 
Molnár, Z. S., Bölöni, J. & Horváth, F. Threatening factors encountered: actual endangerment of the Hungarian (semi-)natural habitats. Acta Bot. Hung. 50, 199–217 (2008).
Google Scholar 
Király, G., Molnár, ZS., Bölöni, J., Csiky, J. & Vojtkó, A. Magyarország földrajzi kistájainak növényzete (in Hungarian). MTA ÖBKI, Vácrátót, 248 (2008).Botta-Dukát, Z. Invasion of alien species to Hungarian (semi-)natural habitats. Acta Bot. Hung. 50, 219–227 (2008).
Google Scholar 
Csákvári, E., Bede-Fazekas, Á., Horváth, F., Molnár, Z. & Halassy, M. Do environmental predictors affect the regeneration capacity of sandy habitats? A country-wide survey from Hungary. Glob. Ecol. Conserv. 27, e01547 (2021).
Google Scholar 
Somodi, I. et al. Implementation and application of multiple potential natural vegetation models–a case study of Hungary. J. Veg. Sci. 28, 1260–1269 (2017).
Google Scholar 
Bölöni, J., Molnár, Zs. & Kun, A. (Eds.), Magyarország élőhelyei. A hazai vegetációtípusok leírása és határozója (in Hungarian) (Habitats – Description and Identification of Vegetation Types of Hungary, ÁNÉR 2011). MTA Ökológiai és Botanikai Kutatóintézet, Vácrátót, pp. 439. ISBN 978-963-8391-51 (2011).Choi, Y. D. et al. Ecological restoration for future sustainability in a changing environment. Ecoscience 15, 53–64 (2008).CAS 

Google Scholar 
Valkó, O. et al. Abandonment of croplands: problem or chance for grassland restoration? Case studies from Hungary. Ecosyst. Health Sustain. 2, e01208 (2016).
Google Scholar 
Csecserits, A. et al. Tree plantations are hot-spots of plant invasion in a landscape with heterogeneous land-use. Agric. Ecosyst. Environ. 226, 88–98 (2016).
Google Scholar 
Pyšek P. & Richardson D. M. in Biological Invasions. Ecological Studies (Analysis and Synthesis) (ed. Nentwig, W) Ch. 7 (Springer, 2008).Reis, B. P. et al. The long-term effect of initial restoration intervention, landscape composition, and time on the progress of Pannonic sand grassland restoration. Landsc. Ecol. Eng. https://doi.org/10.1007/s11355-022-00512-y (2022).Article 

Google Scholar 
Ruprecht, E. Successfully recovered grassland: a promising example from Romanian old‐fields. Restor. Ecol. 14, 473–480 (2006).
Google Scholar 
Török, P. et al. Restoring grassland biodiversity: sowing low-diversity seed mixtures can lead to rapid favourable changes. Biol. Conserv. 143, 3 (2010).
Google Scholar 
Török, P., Vida, E., Deák, B., Lengyel, S. & Tóthmérész, B. Grassland restoration on former croplands in Europe: an assessment of applicability of techniques and costs. Biodivers. Conserv. 20, 2311–2332 (2011).
Google Scholar 
Prach, K., Jongepierová, I., Řehounková, K. & Fajmon, K. Restoration of grasslands on ex-arable land using regional and commercial seed mixtures and spontaneous succession: successional trajectories and changes in species richness. Agric. Ecosyst. Environ. 182, 131–136 (2014).
Google Scholar 
Prach, K., Chenoweth, J. & del Moral, R. Spontaneous and assisted restoration of vegetation on the bottom of a former water reservoir, the Elwha River, Olympic National Park, WA, USA. Restor. Ecol. 27, 592–599 (2019).
Google Scholar 
Török, P., Helm, A., Kiehl, K., Buisson, E. & Valkó, O. Beyond the species pool: modification of species dispersal, establishment, and assembly by habitat restoration. Restor. Ecol. 26, S65–S72 (2018).
Google Scholar 
Török, P., Bullock James M, J. M., Jiménez‐Alfaro, B. & Sonkoly, J. The importance of dispersal and species establishment in vegetation dynamics and resilience. J. Veg. Sci. 31, 935–942 (2020).
Google Scholar 
Saura, S., Bodin, Ö. & Fortin, M. J. Stepping stones are crucial for species’ long-distance dispersal and range expansion through habitat networks. J. Appl. Ecol. 51, 171–182 (2014).
Google Scholar 
Kirmer, A., Baasch, A. & Tischew, S. Sowing of low and high diversity seed mixtures in ecological restoration of surface mined-land. Appl. Veg. Sci. 15, 198–207 (2012).
Google Scholar 
Llumiquinga, Y. B. et al. Long-term results of initial seeding, mowing and carbon amendment on the restoration of Pannonian sand grassland on old fields. Tuxenia 41, 361–379 (2021).
Google Scholar 
Edwards, A. R. et al. Hay strewing, brush harvesting of seed and soil disturbance as tools for the enhancement of botanical diversity in grasslands. Biol. Conserv. 134, 372–382 (2007).
Google Scholar 
Veldman, J. W. et al. Where tree planting and forest expansion are bad for biodiversity and ecosystem services. BioScience 65, 1011–1018 (2015).
Google Scholar 
Bussion, E., Archibald, S., Fidelis, A. & Sudling, K. N. Ancient grasslands guide ambitious goals in grassland restoration. Science 377, 594–598 (2022).
Google Scholar 
Csecserits, A. et al. Regeneration of sandy old-field in the forest steppe region of Hungary. Plant Biosyst. 145, 715–726 (2011).
Google Scholar 
Szitár, K. et al. Az országos zöldinfrastruktúrahálózat kijelölésének módszertana többszempontú állapotértékelés alapján. (in Hungarian) (Methodology for designating the national green infrastructure network based on multi-criteria assessment). Term.észetvédelmi K.özlemények 27, 145–157 (2021).
Google Scholar 
Szalai, S., Szinell, C. S. & Zoboki, J. Early warning systems for drought preparedness and drought management. In Proc. Expert Group Meeting (eds Wilhite, D. A., Sivakumar, M. V. K. & Wood, D. A.) (World Meteorological Organization, 2000).Szilassi, P. et al. The link between landscape pattern and vegetation naturalness on a regional scale. Ecol. Indic. 81, 252–259 (2017).
Google Scholar 
Demeter, I., Makádi, M., Végső, B., Aranyos, T. J. & Posta, K. The effect of recycled plant residues on the microbial activity of typical sandy soil of the Nyírség region. In Abstract Book, 18th Alps-Adria Scientific Workshop https://doi.org/10.34116/NTI.2019.AA.13 (2019).Borhidi, A. Social behaviour types, the naturalness and relative ecological indicator values of the higher plants in the Hungarian Flora. Acta Bot. Hung. 39, 97–181 (1995).
Google Scholar 
Horváth, F. et al. Flóra adatbázis 1.2. Taxonlista és attribútum-állomány (Flora database 1.2. Taxon list and attribute file). MTA Ökológiai és Botanikai Kutatóintézet, Vácrátót, ISBN 9638391197 (1995).Király, G. Új Magyar Füvészkönyv. Magyarország hajtásos növényei (New Herbal Guide to the Hungarian Flora). Aggteleki Nemzeti Park Igazgatóság, Jósvafő, Hungary, 628p. (2009).Máté, A. 6260 pannon homoki gyepek. In: Haraszthy, L. (Eds.), Natura 2000 fajok és élőhelyek Magyarországon. (in Hungarian) Pro Vértes Közalapítvány, Csákvár, Hungary, pp. 817-823. ISBN: 9789630888530 (2014).Molnár, Z. S. et al. Magyarországi Élőhelytérképezési Adatbázisának (MÉTA) térképezési módszertani és Adatlapkitöltési Útmutatója (AL-KÚ) 3.3 Kézirat, (Guide on the methods of MÉTA and on the completion of the MÉTA datasheets). MTA ÖBKI, Vácrátót, Hungary, 54 pp. (2003).Molnár, Z. S. et al. A grid-based, satellite-image supported multi-attributed vegetation mapping method (MÉTA). Folia Geobotanica 42, 225–247 (2007).
Google Scholar 
Horváth, F. et al. Fact sheet of the MÉTA database 1.2. Acta Bot. Hung. 50, 11–34 (2008).
Google Scholar 
Bölöni, J., Kun, A. & Molnár, Z. S. Élőhely-ismereti Útmutató (Habitat guide). MTA ÖBKI, Vácrátót, Hungary (2003).European Environment Agency. Corine Land Cover 2006 seamless vector data (Version 17). https://www.eea.europa.eu/data-and-maps/data/clc-2006-vector-data-version-3 (2013).European Environment Agency. CLC2006 Technical Guidelines. Report No. 17/2007, ISNN 1725-2237 (2017).ESRI ArcGIS Vers. 10.2. (Environmental System Research Institute Inc., 2013).Pásztor, L. et al. Compilation of novel and renewed, goal oriented digital soil maps using geostatistical and data mining tools. Hungarian Geogr. Bull. 64, 49–64 (2015).
Google Scholar 
Hijmans, R. J. raster: geographic data analysis and modeling. R package version 2.4-20, https://cran.r-project.org/web/packages/raster/index.html (2015).R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing. https://www.r-project.org/ (2019).USGS. Shuttle Radar Topography Mission, 3 Arc Second scene SRTM_u03_n045e016-SRTM_ff03_n048e022, Unfilled Unfinished 2.0, Global Land Cover Facility, February 2000. College Park, MD, USA, University of Maryland (2004).SRTM. SRTM Mission Summary. URL: lta.cr.usgs.gov/srtm/mission_summary (2015). [Last accesed: 2016.04.22.].Szalai, S. et al. Climate of the Greater Carpathian Region. Final Technical Report. http://www.carpatclim-eu.org/ (2013).Liaw, A. & Wiener, M. Classification and regression by randomForest. R. N. 2, 18–22, https://CRAN.R-project.org/doc/Rnews/ (2002).
Google Scholar 
Breiman, L., Friedman, J., Stone, C. J. & Olshen, R. A. Classification and Regression Trees (CRC Press, 1984).Sarica, A., Cerasa, A. & Quattrone, A. Random Forest algorithm for the classification of neuroimaging data in Alzheimer’s disease: a systematic review. Front. Aging Neurosci. 6, 329 (2017).
Google Scholar 
Hothorn, T., Hornik, K. & Zeileis, A. Unbiased recursive partitioning: a conditional inference framework. J. Comput. Graph Stat. 15, 651–674 (2006).
Google Scholar 
Pebesma, E. Simple features for R: standardized support for spatial vector. Data. R. J. 10, 439–446 (2018).
Google Scholar 
Bivand, R. S., Pebesma, E. & Gomez-Rubio, V. Applied Spatial Data Analysis with R 2nd ed. (Springer, 2013).Bivand, R. S. & Wong, D. W. S. Comparing implementations of global and local indicators of spatial association. TEST 27, 716–748 (2018).
Google Scholar 
Bölöni, J., Molnár, Z. S., Horváth, F. & Illyés, E. Naturalness-based habitat quality of the Hungarian (semi-)natural habitats. Acta Bot. Hung. 50, 149–159 (2008).
Google Scholar 
Czúcz, B., Molnár, Z. S., Horváth, F. & Botta-Dukát, Z. The natural capital index of Hungary. Acta Bot. Hung. 50, 161–177 (2008).
Google Scholar  More

Callaghan, T. V. et al. Multiple effects of changes in arctic snow cover. Ambio 40, 32–45 (2011).
Google Scholar 
Cooper, E. J. Warmer shorter winters disrupt arctic terrestrial ecosystems. Annu. Rev. Ecol. Evol. Syst. 45, 271 (2014).
Google Scholar 
IPCC. IPCC, 2013: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [Stocker, T.F., D. Qin, G.-K. Plattner, M. Tignor, S.K. Allen, J. Boschung, A. Nauels, Y. Xia, V. Bex and P.M. Midgley (eds.)]. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 1535 pp. https://doi.org/10.1017/CBO9781107415324. (Cambridge University Press, 2013).Seastedt, T. R. The role of microarthropods in decomposition and mineralization processes. Annu. Rev. Entomol. 29, 25–46 (1984).
Google Scholar 
Osler, G. H. & Sommerkorn, M. Toward a complete soil C and N cycle: Incorporating the soil fauna. Ecology 88, 1611–1621 (2007).PubMed 

Google Scholar 
Coulson, S. J. et al. The terrestrial and freshwater invertebrate biodiversity of the archipelagoes of the Barents Sea, Svalbard, Franz Josef Land and Novaya Zemlya. Soil Biol. Biochem. 68, 440–470 (2014).CAS 

Google Scholar 
Hodkinson, I. D. Terrestrial and freshwater invertebrates. In Arctic Biodiversity Assessment (ed. Barry, T.) 246–274 (Arctic Council, 2013).
Google Scholar 
Strathdee, A. T. & Bale, J. S. Life on the edge: Insect ecology in arctic environments. Annu. Rev. Entomol. 43, 85–106 (1998).CAS 
PubMed 

Google Scholar 
Templer, P. H. et al. Impact of a reduced winter snowpack on litter arthropod abundance and diversity in a northern hardwood forest ecosystem. Biol. Fertil. Soils 48, 413–424 (2012).
Google Scholar 
Bokhorst, S., Metcalfe, D. B. & Wardle, D. A. Reduction in snow depth negatively affects decomposers but impact on decomposition rates is substrate dependent. Soil Biol. Biochem. 62, 157–164 (2013).CAS 

Google Scholar 
Slatyer, R. A., Nash, M. A. & Hoffmann, A. A. Measuring the effects of reduced snow cover on Australia’s alpine arthropods. Austral Ecol. 42, 844–857 (2017).
Google Scholar 
Bokhorst, S. et al. Extreme winter warming events more negatively impact small rather than large soil fauna: Shift in community composition explained by traits not taxa. Glob. Change Biol. 18, 1152–1162 (2012).ADS 

Google Scholar 
Sulkava, P. & Huhta, V. Effects of hard frost and freeze-thaw cycles on decomposer communities and N mineralisation in boreal forest soil. Appl. Soil Ecol. 22, 225–239 (2003).
Google Scholar 
Konestabo, H. S., Michelsen, A. & Holmstrup, M. Responses of springtail and mite populations to prolonged periods of soil freeze-thaw cycles in a sub-arctic ecosystem. Appl. Soil Ecol. 36, 136–146 (2007).
Google Scholar 
Coulson, S. J., Leinaas, H. P., Ims, R. A. & Søvik, G. Experimental manipulation of the winter surface ice layer: The effects on a high arctic soil microarthropod community. Ecography 23, 299–306 (2000).
Google Scholar 
Dollery, R., Hodkinson, I. D. & Jonsdottir, I. S. Impact of warming and timing of snow melt on soil microarthropod assemblages associated with Dryas-dominated plant communities on Svalbard. Ecography 29, 111–119 (2006).
Google Scholar 
Ávila-Jimenez, M. L., Coulson, S. J., Solhoy, T. & Sjoblom, A. Overwintering of terrestrial Arctic arthropods: The fauna of Svalbard now and in the future. Polar Res. 29, 127–137 (2010).
Google Scholar 
Makkonen, M. et al. Traits explain the responses of a sub-arctic Collembola community to climate manipulation. Soil Biol. Biochem. 43, 377–384 (2011).CAS 

Google Scholar 
Lindo, Z. Warming favours small-bodied organisms through enhanced reproduction and compositional shifts in belowground systems. Soil Biol. Biochem. 91, 271–278 (2015).CAS 

Google Scholar 
Hågvar, S. A review of Fennoscandian arthropods living on and in snow. Eur. J. Entomol. 107, 281–298 (2010).
Google Scholar 
Hao, C., Chen, T.-W., Wu, Y., Chang, L. & Wu, D. Snow microhabitats provide food resources for winter-active Collembola. Soil Biol. Biochem. 143, 107731 (2020).CAS 

Google Scholar 
Christenson, L. et al. Winter climate change influences on soil faunal distribution and abundance: Implications for decomposition in the Northern Forest. Northeast. Nat. 24, B209–B234 (2017).
Google Scholar 
Convey, P. et al. Survival of rapidly fluctuating natural low winter temperatures by High Arctic soil invertebrates. J. Therm. Biol. 54, 111–117 (2015).PubMed 

Google Scholar 
Krab, E. J., Monteux, S., Weedon, J. T. & Dorrepaal, E. Plant expansion drives bacteria and collembola communities under winter climate change in frost-affected tundra. Soil Biol. Biochem. 138, 107569 (2019).CAS 

Google Scholar 
Sörensen, J. G. & Holmstrup, M. Cryoprotective dehydration is widespread in Arctic springtails. J. Insect Physiol. 57, 1147–1153 (2011).PubMed 

Google Scholar 
Convey, P., Coulson, S. J., Worland, M. R. & Sjöblom, A. The importance of understanding annual and shorter-term temperature patterns and variation in the surface levels of polar soils for terrestrial biota. Polar Biol. 41, 1587–1605 (2018).
Google Scholar 
Birkemoe, T. & Leinaas, H. P. Reproductive biology of the arctic collembolan Hypogastrura tullbergi. Ecography 22, 31–39 (1999).
Google Scholar 
Birkemoe, T. & Leinaas, H. P. Effects of temperature on the development of an arctic Collembola (Hypogastrura tullbergi). Funct. Ecol. 14, 693–700 (2001).
Google Scholar 
Kankaanpää, T. et al. Spatiotemporal snowmelt patterns within a high Arctic landscape, with implications for flora and fauna. Arct. Antarct. Alp. Res. 50, e1415624 (2018).
Google Scholar 
Cooper, E. J., Dullinger, S. & Semenchuk, P. Late snowmelt delays plant development and results in lower reproductive success in the high arctic. Plant Sci. 180, 157–167 (2011).CAS 
PubMed 

Google Scholar 
Krab, E. J. et al. Winter warming effects on tundra shrub performance are species-specific and dependent on spring conditions. J. Ecol. 106, 599–612 (2018).CAS 

Google Scholar 
Wheeler, H. C., Hoye, T. T., Schmidt, N. M., Svenning, J.-C. & Forchhammer, M. C. Phenological mismatch with abiotic conditions-implications for flowering in Arctic plants. Ecology 96, 775–787 (2015).PubMed 

Google Scholar 
Wheeler, J. A. et al. The snow and the willows: Earlier spring snowmelt reduces performance in the low-lying alpine shrub Salix herbacea. J. Ecol. 104, 1041–1050 (2016).CAS 

Google Scholar 
Pollierer, M. M., Langel, R., Körner, C., Maraun, M. & Scheu, S. The underestimated importance of belowground carbon input for forest soil animal food webs. Ecol. Lett. 10, 729–736 (2007).PubMed 

Google Scholar 
Coulson, S. J., Hodkinson, I. D. & Webb, N. R. Microscale distribution patterns in high Arctic soil microarthropod communities: The influence of plant species within the vegetation mosaic. Ecography 26, 801–809 (2003).
Google Scholar 
Hodkinson, I. D. et al. Global change and Arctic ecosystems: Conclusions and predictions from experiments with terrestrial invertebrates on Spitsbergen. Arct. Alp. Res. 30, 306–313 (1998).
Google Scholar 
Førland, E. J., Benestad, R., Hanssen-Bauer, I., Haugen, J. E. & Skaugen, T. E. Temperature and precipitation development at Svalbard 1900–2100. Adv. Meteorol. 2011, 893790 (2011).
Google Scholar 
Alatalo, J. M., Jagerbrand, A. K. & Cuchta, P. Collembola at three alpine subarctic sites resistant to twenty years of experimental warming. Sci. Rep. 5, 18161 (2015).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Coulson, S. J. et al. Effects of experimental temperature elevation on high-arctic soil microarthropod populations. Polar Biol. 16, 147–153 (1996).
Google Scholar 
Decker, K. L. M., Wang, D., Waite, C. & Scherbatskoy, T. Snow removal and ambient air temperature effects on forest soil temperatures in Northern Vermont. Soil Sci. Soc. Am. J. 67, 1234–1242 (2003).ADS 
CAS 

Google Scholar 
van Pelt, W. J. J. et al. Multidecadal climate and seasonal snow conditions in Svalbard. J. Geophys. Res. Earth Surf. 121, 2100–2117 (2016).ADS 

Google Scholar 
Semenchuk, P. R. et al. Deeper snow alters soil nutrient availability and leaf nutrient status in high Arctic tundra. Biogeochemistry 124, 81–94 (2015).
Google Scholar 
Sjursen, H., Michelsen, A. & Jonasson, S. Effects of long-term soil warming and fertilisation on microarthropod abundances in three sub-arctic ecosystems. Appl. Soil Ecol. 30, 148–161 (2005).
Google Scholar 
Meehan, M. L. et al. Response of soil fauna to simulated global change factors depends on ambient climate conditions. Pedobiologia 83, 150672 (2020).
Google Scholar 
Harte, J., Rawa, A. & Price, V. Effects of manipulated soil microclimate on mesofaunal biomass and diversity. Soil Biol. Biochem. 28, 313–322 (1996).CAS 

Google Scholar 
Siepel, H. Life history tactics of soil microarthropods. Biol. Fertil. Soils 18, 263–278 (1994).
Google Scholar 
Chernova, N. M., Potapov, M. B., Savenkova, Y. Y. & Bokova, A. I. Ecological significance of parthenogenesis in Collembola. Zool. Zhurnal 88, 1455–1470 (2009).
Google Scholar 
Birkemoe, T. & Somme, L. Population dynamics of two collembolan species in an Arctic tundra. Pedobiologia 42, 131–145 (1998).
Google Scholar 
Bokhorst, S. et al. Contrasting responses of springtails and mites to elevation and vegetation type in the sub-Arctic. Pedobiologia 67, 57–64 (2018).
Google Scholar 
Widenfalk, L. A., Malmstrom, A., Berg, M. P. & Bengtsson, J. Small-scale Collembola community composition in a pine forest soil—Overdispersion in functional traits indicates the importance of species interactions. Soil Biol. Biochem. 103, 52–62 (2016).CAS 

Google Scholar 
Morgner, E. The importance of winter in annual ecosystem respiration in the High Arctic: Effects of snow depth in two vegetation types. Polar Res. 29, 474–474 (2010).
Google Scholar 
Green, K. & Slatyer, R. Arthropod community composition along snowmelt gradients in snowbeds in the Snowy Mountains of south-eastern Australia. Austral Ecol. 45, 144–157 (2020).
Google Scholar 
Ayres, E. et al. Experimentally increased snow accumulation alters soil moisture and animal community structure in a polar desert. Polar Biol. 33, 897–907 (2010).
Google Scholar 
Semenchuk, P. R., Elberling, B. & Cooper, E. J. Snow cover and extreme winter warming events control flower abundance of some, but not all species in high arctic Svalbard. Ecol. Evol. 3, 2586–2599 (2013).PubMed 
PubMed Central 

Google Scholar 
Morsdorf, M. A. et al. Deepened winter snow significantly influences the availability and forms of nitrogen taken up by plants in High Arctic tundra. Soil Biol. Biochem. 135, 222–234 (2019).CAS 

Google Scholar 
Cooper, E. J., Little, C. J., Pilsbacher, A. K. & Morsdorf, M. A. Disappearing green: Shrubs decline and bryophytes increase with nine years of increased snow accumulation in the High Arctic. J. Veg. Sci. 30, 857–867 (2019).
Google Scholar 
Mundra, S. et al. Ectomycorrhizal and saprotrophic fungi respond differently to long-term experimentally increased snow depth in the High Arctic. Microbiologyopen 5, 856–869 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Schneider, K. & Maraun, M. Feeding preferences among dark pigmented fungal taxa (“Dematiacea”) indicate limited trophic niche differentiation of oribatid mites (Oribatida, Acari). Pedobiologia 49, 61–67 (2005).
Google Scholar 
Krab, E. J., Berg, M. P., Aerts, R., van Logtestijn, R. S. P. & Cornelissen, J. H. C. Vascular plant litter input in subarctic peat bogs changes Collembola diets and decomposition patterns. Soil Biol. Biochem. 63, 106–115 (2013).CAS 

Google Scholar 
Jucevica, E. & Melecis, V. Global warming affect Collembola community: A long-term study. Pedobiologia 50, 177–184 (2006).
Google Scholar 
Krab, E. J., Oorsprong, H., Berg, M. P. & Cornelissen, J. H. C. Turning northern peatlands upside down: Disentangling microclimate and substrate quality effects on vertical distribution of Collembola. Funct. Ecol. 24, 1362–1369 (2010).
Google Scholar 
Zettel, J. Alpine Collembola—Adaptations and strategies for survival in harsh environments. Zool. Anal. Complex Syst. 102, 73–89 (2000).
Google Scholar 
Block, W. Terrestrial arthropods and low-temperature. Cryobiology 18, 436–444 (1981).CAS 
PubMed 

Google Scholar 
Semenchuk, P. R., Christiansen, C. T., Grogan, P., Elberling, B. & Cooper, E. J. Long-term experimentally deepened snow decreases growing-season respiration in a low- and high-arctic tundra ecosystem. J. Geophys. Res. Biogeosci. 121, 1236–1248 (2016).
Google Scholar 
Semenchuk, P. R. et al. Soil organic carbon depletion and degradation in surface soil after long-term non-growing season warming in High Arctic Svalbard. Sci. Total Environ. 646, 158–167 (2019).ADS 
CAS 
PubMed 

Google Scholar 
Gillespie, M. A. K. et al. Status and trends of terrestrial arthropod abundance and diversity in the North Atlantic region of the Arctic. Ambio 49, 718–731 (2020).PubMed 

Google Scholar 
Andriuzzi, W. S., Adams, B. J., Barrett, J. E., Virginia, R. A. & Wall, D. H. Observed trends of soil fauna in the Antarctic Dry Valleys: Early signs of shifts predicted under climate change. Ecology 99, 312–321 (2018).CAS 
PubMed 

Google Scholar 
Staub, B. & Delaloye, R. Using near-surface ground temperature data to derive snow insulation and melt indices for mountain permafrost applications. Permafr. Periglac. Process. 28, 237–248 (2017).
Google Scholar 
Rendos, M. et al. Organic carbon content and temperature as substantial factors affecting diversity and vertical distribution of Collembola on forested scree slopes. Eur. J. Soil Biol. 75, 180–187 (2016).
Google Scholar 
Fjellberg, A. The Collembola of the Norwegian Arctic Islands (Norsk Polarinstitutt, 1994).
Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing. https://www.R-project.org/. (R Foundation for Statistical Computing, 2020). Accessed 06 June 2020. More

Raschke, K., Monteith, J. L. & Weatherley, P. E. How stomata resolve the dilemma of opposing priorities. Phil. Trans. R. Soc. Lond. B 273, 551–560 (1976).Article 
CAS 

Google Scholar 
Brodribb, T. J. & Cochard, H. Hydraulic failure defines the recovery and point of death in water-stressed conifers. Plant Physiol. 149, 575–584 (2009).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Brodribb, T. J., Bowman, D. J. M. S., Nichols, S., Delzon, S. & Burlett, R. Xylem function and growth rate interact to determine recovery rates after exposure to extreme water deficit. New Phytol. 188, 533–542 (2010).Article 
PubMed 

Google Scholar 
Choat, B. et al. Triggers of tree mortality under drought. Nature 558, 531–539 (2018).Article 
CAS 
PubMed 

Google Scholar 
Keeling, R. F. et al. Atmospheric evidence for a global secular increase in carbon isotopic discrimination of land photosynthesis. Proc. Natl Acad. Sci. USA https://doi.org/10.1073/pnas.1619240114 (2017).Guerrieri, R. et al. Disentangling the role of photosynthesis and stomatal conductance on rising forest water-use efficiency. Proc. Natl Acad. Sci. USA 116, 16909–16914 (2019).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Grossiord, C. et al. Plant responses to rising vapor pressure deficit. New Phytol. 226, 1550–1566 (2020).Article 
PubMed 

Google Scholar 
McDowell, N. G. & Allen, C. D. Darcy’s law predicts widespread forest mortality under climate warming. Nat. Clim. Change 5, 669–672 (2015).Article 

Google Scholar 
Brienen, R. J. W. et al. Long-term decline of the Amazon carbon sink. Nature 519, 344–348 (2015).Article 
CAS 
PubMed 

Google Scholar 
McDowell, N. G. et al. Pervasive shifts in forest dynamics in a changing world. Science 368, eaaz9463 (2020).Article 
CAS 
PubMed 

Google Scholar 
Novick, K. A. et al. The increasing importance of atmospheric demand for ecosystem water and carbon fluxes. Nat. Clim. Change 6, 1023–1027 (2016).Article 
CAS 

Google Scholar 
Damour, G., Simonneau, T., Cochard, H. & Urban, L. An overview of models of stomatal conductance at the leaf level. Plant Cell Environ. 33, 1419–1438 (2010).PubMed 

Google Scholar 
Wang, Y., Sperry, J. S., Anderegg, W. R. L., Venturas, M. D. & Trugman, A. T. A theoretical and empirical assessment of stomatal optimization modeling. New Phytol. 227, 311–325 (2020).Article 
CAS 
PubMed 

Google Scholar 
Anderegg, W. R. L. et al. Woody plants optimise stomatal behaviour relative to hydraulic risk. Ecol. Lett. 21, 968–977 (2018).Article 
PubMed 

Google Scholar 
Venturas, M. D. et al. A stomatal control model based on optimization of carbon gain versus hydraulic risk predicts aspen sapling responses to drought. New Phytol. 220, 836–850 (2018).Article 
CAS 
PubMed 

Google Scholar 
Sabot, M. E. B. et al. Plant profit maximization improves predictions of European forest responses to drought. New Phytol. 226, 1638–1655 (2020).Article 
PubMed 

Google Scholar 
Eller, C. B. et al. Stomatal optimization based on xylem hydraulics (SOX) improves land surface model simulation of vegetation responses to climate. New Phytol. 226, 1622–1637 (2020).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Hickler, T., Prentice, I. C., Smith, B., Sykes, M. T. & Zaehle, S. Implementing plant hydraulic architecture within the LPJ Dynamic Global Vegetation Model. Glob. Ecol. Biogeogr. 15, 567–577 (2006).Article 

Google Scholar 
Bonan, G. B., Williams, M., Fisher, R. A. & Oleson, K. W. Modeling stomatal conductance in the earth system: linking leaf water-use efficiency and water transport along the soil–plant–atmosphere continuum. Geosci. Model Dev. 7, 2193–2222 (2014).Article 

Google Scholar 
Christoffersen, B. O. et al. Linking hydraulic traits to tropical forest function in a size-structured and trait-driven model (TFS v.1-Hydro). Geosci. Model Dev. 9, 4227–4255 (2016).Article 

Google Scholar 
Kennedy, D. et al. Implementing plant hydraulics in the Community Land Model, Version 5. J. Adv. Model. Earth Syst. 11, 485–513 (2019).Article 

Google Scholar 
Cowan, I. R. & Farquhar, G. D. Stomatal function in relation to leaf metabolism and environment. Symp. Soc. Exp. Biol. 31, 471–505 (1977).CAS 
PubMed 

Google Scholar 
Wolf, A., Anderegg, W. R. L. & Pacala, S. W. Optimal stomatal behavior with competition for water and risk of hydraulic impairment. Proc. Natl Acad. Sci. USA 113, E7222–E7230 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Sperry, J. S. et al. Predicting stomatal responses to the environment from the optimization of photosynthetic gain and hydraulic cost. Plant Cell Environ. 40, 816–830 (2017).Article 
CAS 
PubMed 

Google Scholar 
Bartlett, M. K., Detto, M. & Pacala, S. W. Predicting shifts in the functional composition of tropical forests under increased drought and CO2 from trade-offs among plant hydraulic traits. Ecol. Lett. 22, 67–77 (2019).Article 
PubMed 

Google Scholar 
Prentice, I. C., Dong, N., Gleason, S. M., Maire, V. & Wright, I. J. Balancing the costs of carbon gain and water transport: testing a new theoretical framework for plant functional ecology. Ecol. Lett. 17, 82–91 (2014).Article 
PubMed 

Google Scholar 
Wright, I. J., Reich, P. B. & Westoby, M. Least‐cost input mixtures of water and nitrogen for photosynthesis. Am. Nat.161, 98–111 (2003).Article 
PubMed 

Google Scholar 
Wang, H. et al. Towards a universal model for carbon dioxide uptake by plants. Nat. Plants 3, 734–741 (2017).Article 
CAS 
PubMed 

Google Scholar 
Maire, V. et al. The coordination of leaf photosynthesis links C and N fluxes in C3 plant species. PLoS ONE 7, e38345 (2012).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Stocker, B. D. et al. Quantifying soil moisture impacts on light use efficiency across biomes. New Phytol. 218, 1430–1449 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Stocker, B. D. et al. P-model v1.0: an optimality-based light use efficiency model for simulating ecosystem gross primary production. Geosci. Model Dev. 13, 1545–1581 (2020).Article 

Google Scholar 
Lavergne, A. et al. Historical changes in the stomatal limitation of photosynthesis: empirical support for an optimality principle. New Phytol. 225, 2484–2497 (2020).Article 
CAS 
PubMed 

Google Scholar 
Sperry, J. S. & Love, D. M. What plant hydraulics can tell us about responses to climate-change droughts. New Phytol. 207, 14–27 (2015).Article 
CAS 
PubMed 

Google Scholar 
Farquhar, G. D., von Caemmerer, S. & Berry, J. A. A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta 149, 78–90 (1980).Article 
CAS 
PubMed 

Google Scholar 
Chen, J.-L., Reynolds, J. F., Harley, P. C. & Tenhunen, J. D. Coordination theory of leaf nitrogen distribution in a canopy. Oecologia 93, 63–69 (1993).Article 
PubMed 

Google Scholar 
Buckley, T. N., John, G. P., Scoffoni, C. & Sack, L. How does leaf anatomy influence water transport outside the xylem? Plant Physiol. 168, 1616–1635 (2015).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Scoffoni, C. et al. Outside-xylem vulnerability, not xylem embolism, controls leaf hydraulic decline during dehydration. Plant Physiol. 173, 1197–1210 (2017).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Carminati, A. & Javaux, M. Soil rather than xylem vulnerability controls stomatal response to drought. Trends Plant Sci. 25, 868–880 (2020).Article 
CAS 
PubMed 

Google Scholar 
Klein, T. et al. Xylem embolism refilling and resilience against drought-induced mortality in woody plants: processes and trade-offs. Ecol. Res. 33, 839–855 (2018).CAS 

Google Scholar 
Rodriguez-Dominguez, C. M. & Brodribb, T. J. Declining root water transport drives stomatal closure in olive under moderate water stress. New Phytol. 225, 126–134 (2020).Article 
CAS 
PubMed 

Google Scholar 
Sack, L. & Holbrook, N. M. Leaf hydraulics. Annu. Rev. Plant Biol. 57, 361–381 (2006).Article 
CAS 
PubMed 

Google Scholar 
Bourbia, I., Pritzkow, C. & Brodribb, T. J. Herb and conifer roots show similar high sensitivity to water deficit. Plant Physiol. 186, 1908–1918 (2021).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Zhou, S., Duursma, R. A., Medlyn, B. E., Kelly, J. W. G. & Prentice, I. C. How should we model plant responses to drought? An analysis of stomatal and non-stomatal responses to water stress. Agric. Meteorol. 182–183, 204–214 (2013).Article 

Google Scholar 
Kanechi, M., Uchida, N., Yasuda, T. & Yamaguchi, T. Non-stomatal inhibition associated with inactivation of rubisco in dehydrated coffee leaves under unshaded and shaded conditions. Plant Cell Physiol. 37, 455–460 (1996).Article 
CAS 

Google Scholar 
Salmon, Y. et al. Leaf carbon and water status control stomatal and nonstomatal limitations of photosynthesis in trees. New Phytol. 226, 690–703 (2020).Article 
CAS 
PubMed 

Google Scholar 
Dong, N. et al. Components of leaf-trait variation along environmental gradients. New Phytol. 228, 82–94 (2020).Article 
CAS 
PubMed 

Google Scholar 
Martínez‐Vilalta, J., Poyatos, R., Aguadé, D., Retana, J. & Mencuccini, M. A new look at water transport regulation in plants. New Phytol. 204, 105–115 (2014).Article 
PubMed 

Google Scholar 
Bartlett, M. K., Klein, T., Jansen, S., Choat, B. & Sack, L. The correlations and sequence of plant stomatal, hydraulic, and wilting responses to drought. Proc. Natl Acad. Sci. USA 113, 13098–13103 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Brodribb, T. J., Holbrook, N. M., Edwards, E. J. & Gutiérrez, M. V. Relations between stomatal closure, leaf turgor and xylem vulnerability in eight tropical dry forest trees. Plant Cell Environ. 26, 443–450 (2003).Article 

Google Scholar 
Martin‐StPaul, N., Delzon, S. & Cochard, H. Plant resistance to drought depends on timely stomatal closure. Ecol. Lett. 20, 1437–1447 (2017).Article 
PubMed 

Google Scholar 
Skelton, R. P. et al. Low vulnerability to xylem embolism in leaves and stems of North American oaks. Plant Physiol. 177, 1066–1077 (2018).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Choat, B. et al. Global convergence in the vulnerability of forests to drought. Nature 491, 752–755 (2012).Article 
CAS 
PubMed 

Google Scholar 
Dewar, R. et al. New insights into the covariation of stomatal, mesophyll and hydraulic conductances from optimization models incorporating nonstomatal limitations to photosynthesis. New Phytol. 217, 571–585 (2018).Article 
CAS 
PubMed 

Google Scholar 
Hölttä, T., Lintunen, A., Chan, T., Mäkelä, A. & Nikinmaa, E. A steady-state stomatal model of balanced leaf gas exchange, hydraulics and maximal source–sink flux. Tree Physiol. 37, 851–868 (2017).Article 
PubMed 

Google Scholar 
Pivovaroff, A. L., Sack, L. & Santiago, L. S. Coordination of stem and leaf hydraulic conductance in southern California shrubs: a test of the hydraulic segmentation hypothesis. New Phytol. 203, 842–850 (2014).Article 
PubMed 

Google Scholar 
Boyer, J. S., Wong, S. C. & Farquhar, G. D. CO2 and water vapor exchange across leaf cuticle (epidermis) at various water potentials. Plant Physiol. 114, 185–191 (1997).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Deans, R. M., Brodribb, T. J., Busch, F. A. & Farquhar, G. D. Optimization can provide the fundamental link between leaf photosynthesis, gas exchange and water relations. Nat. Plants 6, 1116–1125 (2020).Article 
CAS 
PubMed 

Google Scholar 
Zhou, S.-X., Medlyn, B. E. & Prentice, I. C. Long-term water stress leads to acclimation of drought sensitivity of photosynthetic capacity in xeric but not riparian Eucalyptus species. Ann. Bot. 117, 133–144 (2016).Article 
CAS 
PubMed 

Google Scholar 
Rungwattana, K. et al. Trait evolution in tropical rubber (Hevea brasiliensis) trees is related to dry season intensity. Funct. Ecol. 32, 2638–2651 (2018).Article 

Google Scholar 
Dybzinski, R., Farrior, C., Wolf, A., Reich, P. B. & Pacala, S. W. Evolutionarily stable strategy carbon allocation to foliage, wood, and fine roots in trees competing for light and nitrogen: an analytically tractable, individual-based model and quantitative comparisons to data. Am. Nat. 177, 153–166 (2011).Article 
PubMed 

Google Scholar 
Hikosaka, K. & Anten, N. P. R. An evolutionary game of leaf dynamics and its consequences for canopy structure. Funct. Ecol. 26, 1024–1032 (2012).Article 

Google Scholar 
Franklin, O. et al. Organizing principles for vegetation dynamics. Nat. Plants 6, 444–453 (2020).Article 
PubMed 

Google Scholar 
Le Quéré, C. et al. Global carbon budget 2017. Earth Syst. Sci. Data 10, 405–448 (2018).Article 

Google Scholar 
Jasechko, S. et al. Terrestrial water fluxes dominated by transpiration. Nature 496, 347–350 (2013).Article 
CAS 
PubMed 

Google Scholar 
Xu, X., Medvigy, D., Powers, J. S., Becknell, J. M. & Guan, K. Diversity in plant hydraulic traits explains seasonal and inter-annual variations of vegetation dynamics in seasonally dry tropical forests. New Phytol. 212, 80–95 (2016).Article 
PubMed 

Google Scholar 
Wang, H. et al. Acclimation of leaf respiration consistent with optimal photosynthetic capacity. Glob. Change Biol. 26, 2573–2583 (2020).Article 

Google Scholar 
Papastefanou, P. et al. A dynamic model for strategies and dynamics of plant water-potential regulation under drought conditions. Front. Plant Sci. 11, 373 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Grieu, P., Guehl, J. M. & Aussenac, G. The effects of soil and atmospheric drought on photosynthesis and stomatal control of gas exchange in three coniferous species. Physiol. Plant. 73, 97–104 (1988).Article 

Google Scholar 
Liu, F., Andersen, M. N., Jacobsen, S.-E. & Jensen, C. R. Stomatal control and water use efficiency of soybean (Glycine max L. Merr.) during progressive soil drying. Environ. Exp. Bot. 54, 33–40 (2005).Article 
CAS 

Google Scholar 
Tezara, W., Driscoll, S. & Lawlor, D. W. Partitioning of photosynthetic electron flow between CO2 assimilation and O2 reduction in sunflower plants under water deficit. Photosynthetica 46, 127–134 (2008).Article 
CAS 

Google Scholar 
Liu, C.-C. et al. Influence of drought intensity on the response of six woody karst species subjected to successive cycles of drought and rewatering. Physiol. Plant. 139, 39–54 (2010).Article 
CAS 
PubMed 

Google Scholar 
Posch, S. & Bennett, L. T. Photosynthesis, photochemistry and antioxidative defence in response to two drought severities and with re-watering in Allocasuarina luehmannii. Plant Biol. 11, 83–93 (2009).Article 
CAS 
PubMed 

Google Scholar 
Jiang, M., Kelly, J. W. G., Atwell, B. J., Tissue, D. T. & Medlyn, B. E. Drought by CO2 interactions in trees: a test of the water savings mechanism. New Phytol. 230, 1421–1434 (2021).Article 
CAS 
PubMed 

Google Scholar 
Ennajeh, M., Tounekti, T., Vadel, A. M., Khemira, H. & Cochard, H. Water relations and drought-induced embolism in olive (Olea europaea) varieties ‘Meski’ and ‘Chemlali’ during severe drought. Tree Physiol. 28, 971–976 (2008).Article 
PubMed 

Google Scholar 
Peguero-Pina, J. J., Sancho-Knapik, D., Morales, F., Flexas, J. & Gil-Pelegrín, E. Differential photosynthetic performance and photoprotection mechanisms of three Mediterranean evergreen oaks under severe drought stress. Funct. Plant Biol. 36, 453–462 (2009).Article 
PubMed 

Google Scholar 
Liu, C.-C. et al. Exploitation of patchy soil water resources by the clonal vine Ficus tikoua in karst habitats of southwestern China. Acta Physiol. Plant. 33, 93–102 (2011).Article 

Google Scholar 
Leuning, R. A critical appraisal of a combined stomatal-photosynthesis model for C3 plants. Plant Cell Environ. 18, 339–355 (1995).Article 
CAS 

Google Scholar 
Medlyn, B. E. et al. Reconciling the optimal and empirical approaches to modelling stomatal conductance. Glob. Change Biol. 17, 2134–2144 (2011).Article 

Google Scholar 
Brodribb, T. et al. Linking xylem network failure with leaf tissue death. New Phytol. 232, 68–79 (2021).Article 
PubMed 

Google Scholar 
Klein, T. The variability of stomatal sensitivity to leaf water potential across tree species indicates a continuum between isohydric and anisohydric behaviours. Funct. Ecol. 28, 1313–1320 (2014).Article 

Google Scholar  More

Molecular Phylogeny of Paramecium bursaria and Identification of its EndosymbiontsThe SSU and ITS rDNA of the nuclear ribosomal operon were sequenced to infer the genetic variability of the investigated strains. The SSU and ITS rDNA sequences were aligned according to their secondary structure (examples are presented for the strain SAG 27.96; Fig. 1 and Supplementary Fig. 1). Additional sequences acquired from GenBank were incorporated into a dataset, which included all syngens also from references known for P. bursaria. The phylogenetic analyses revealed five highly supported lineages among the P. bursaria strains, which corresponded to their syngen assignment. As demonstrated in Fig. 2, all investigated strains belonging to the syngens R1, R2 and R5 originated from Europe, whereas the others of the syngens R3-R4 showed a worldwide distribution. The three known green algal endosymbionts, i.e., Chlorella variabilis (Cvar), Chlorella vulgaris (Cvul) and Micractinium conductrix (Mcon) showed no or only little affiliation to specific syngens.Figure 1ITS‐1 (A) and ITS-2 (B) secondary structures of Paramecium protobursaria, SAG 27.96 (syngen R1).Full size imageFigure 2Molecular phylogeny of the Paramecium bursaria species complex based on SSU and ITS rDNA sequence comparisons. The phylogenetic tree shown was inferred using the maximum likelihood method based on the datasets (2197 aligned positions of 19 taxa) using the computer program PAUP 4.0a169. For the analyses, the best model was calculated by PAUP 4.0a169. The setting of the best model was given as follows: TVM + I (base frequencies: A 0.2983, C 0.1840, G 0.2271, T 0.2906; rate matrix A–C 2.6501, A–G 8.6851, A–U 5.3270, C–G 0.91732, C–U 8.6851, G–U 1.0000) with the proportion of invariable sites (I = 0.9544). The branches in bold are highly supported in all bootstrap analyses (bootstrap values  > 50% calculated with PAUP using the maximum likelihood, neighbour—joining, and maximum parsimony). The clades are named after the syngens (color‐coded) proposed by Greczek‐Stachura et al.10 and Bomford9 in brackets. The accession numbers are given after the strain numbers. The endosymbiotic green algae identified are highlighted (Mcon—Micractinium conductrix, Cvar—Chlorella variabilis and Cvul—Chlorella vulgaris) after the origin of the P. bursaria strains. The reference strain of each syngen is marked with an asterisk. The strains used for morphological comparisons are marked with a green dot next to the strain number.Full size imageSynapomorphies of the Paramecium bursaria SyngensAs demonstrated in Fig. 2, the subdivision of the P. bursaria strains into syngens is supported by the phylogenetic analyses of the SSU and ITS rDNA sequences. To figure out if these splits were also supported by characteristic molecular signatures, we studied the secondary structures of both SSU and ITS of all available sequences. We discovered 30, respectively 23 variable positions among the SSU and ITS sequences (numbers of these positions in the respective alignments are given in Fig. 3). All syngens showed characteristic patterns among the SSU and ITS. Only the syngens R1 and R2 could not be distinguished using the SSU only, however, in combination with the ITS, each syngen is characterized by unique synapomorphies as highlighted in yellow (Fig. 3). In addition, few variable base positions within syngens (marked in blue in Fig. 3) have been recognized in the ITS regions. For comparison with literature data, we also analyzed all available sequences of the mitochondrial COI gene to find synapomorphies for the five syngens. Within this gene, only 18 variable positions at the amino acid level could be discovered of which 13 are diagnostic for the five syngens (Fig. 3).Figure 3Variable base positions among the SSU, ITS rRNA, and COI sequences of the five syngens among the Paramecium bursaria species complex. The unique synapomorphies are highlighted in yellow, variable positions marked in blue.Full size imageThe synapomorphies discovered above were used to get insights into the geographical distribution of each P. bursaria syngen. Despite the complete SSU and ITS rDNA sequences included in the phylogeny presented in Fig. 2, records of the partial SSU or ITS rDNA sequences are available in GenBank (BLASTn search; 100% identity;13). Considering the metadata of our investigated strains and of the entries in GenBank (Supplementary Table 1), we constructed three haplotype networks using the Templeton-Crandall-Sing (TCS) approach. The SSU haplotype network (Fig. 4) containing 84 records showed that the syngens R1, R2 and R5 were only found in Europe, whereas the other three syngens have been discovered around the world. A similar distribution pattern occurred when using the ITS (101 entries in GenBank). Records of syngens R1 and R5 have only been found in Europe, whereas all other syngens were distributed around the world. The 132 COI records found in GenBank by the BLASTn search were used for the haplotype network, which also showed the similar pattern (Fig. 4).Figure 4TCS haplotype networks of the five syngens inferred from SSU, ITS rRNA, and COI sequences of the Paramecium bursaria species complex. This network was inferred using the algorithm described by Clement et al.40,41. Sequence nodes corresponding to samples collected from different geographical regions.Full size imageCiliate TaxonomyConsidering all our findings, P. bursaria is morphologically highly variable, and obviously represents a cryptic species complex (Figs. 5, 6; Supplementary Table 2). The known five syngens most likely represent biological species according to Mayr14 and can be attributed to the cryptic species described by Greczek-Stachura et al.11. As mentioned above, the assignments of these cryptic species by Greczek-Stachura et al.11 have not been validly described according to the ICZN. In addition, the naming using a mixture of Latin prefix and Greek suffix is also not appropriate (the epithet bursa derived from the Greek word byrsa). Therefore, we describe the five syngens as new species as follows. The general morphological features of these species are summarized in Table 1.Figure 5Ventral views of Paramecium bursaria morphotypes in vivo: P. protobursaria (syngen R1), i.e., strains SAG 2645 (A) and PB-25 (B); P. deuterobursaria (syngen R2), i.e., strains CCAP 1660/36 (C) and CCAP 1660/34 (D); P. tritobursaria (syngen R3), i.e., strains CCAP 1660/28 (E), CCAP 1660/26 (F) and CCAP 1660/31 (G); P. tetratobursaria (syngen R4), i.e., strains CCAP 1660/25 (H) and CCAP 1660/33 (I); P. pentobursaria (syngen R5), i.e., strain CCAP 1660/30 (J). Scale bar 20 µm.Full size imageFigure 6Morphological details of the Paramecium bursaria species complex from specimens of strains PB-25 (A), CCAP 1660/30 (B), SAG 2645 (C, F, G, I, L–N), CCAP 1660/36 (D), CCAP 1660/26 (E, H), CCAP 1660/30 (J, O), CCAP 1660/16 (K) in vivo (A–F, H–O) and after silver nitrate staining (G). Adoral membranelles (A, B), endosymbiotic algae Micractinium conductrix (C), caudal and somatic cilia (D), arrows denote excretory pores of the contractile vacuoles: extruded extrusomes are shown and caudal cilia (E), ventral views showing the preoral suture and the oral opening (F), the ciliary pattern (G), arrows denote excretory pores of the contractile vacuoles (H), trichocysts and symbiotic algae underneath the pellicula (I, J), cell size variations (K), radial collecting channels (white arrows) and excretory pores (black arrows) of contractile vacuoles (L), macro- and micronucleus (M), cytopyge and characteristic rectangular pellicular pattern (N), pattern of the pellicula (O). AS anterior suture, CC caudal cilia, CP cytopyge (cell after), CV contractile vacuole, EP excretory pore of a contractile vacuole, EX extrusomes, M1–M3 membranelles 1–3, MA macronucleus, MI micronucleus, OO oral opening, S symbiotic algae, SC somatic cilia, SK somatic kineties, UM undulating membrane. Scale bars 10 µm (A, I), 20 µm (B, D–H, J, L–O), 50 µm (K).Full size imageTable 1 Main morphometric and morphological characteristics of the Paramecium bursaria syngens (min and max values).Full size table
Paramecium protobursaria sp. nov.Synonym: Paramecium primabursaria nom. inval.Description: The strains SAG 27.96 and PB-25 belong to syngen R1 according to Greczek-Stachura et al.10,11 and differ from other syngens by their SSU and ITS rDNA sequences (MT231333). From morphology, the cells are ellipsoidal to broadly ellipsoidal and dorso-ventrally flattened in vivo. The cells measure 70–164 × 44–65 µm; the single macronucleus is located around mid-cell and measures 25–38 × 11–22 µm; the adjacent single compact micronucleus measures 11–20 × 5–8 µm; the usually two (rarely one) contractile vacuoles, one in the anterior and one in the posterior cell portion have radial collecting channels and 1–3 excretory pores each; the number of ciliary rows/20 µm is 14–22; the length of the caudal cilia is 9–19 µm; the numerous trichocysts located in the cell cortex are 4–6 µm in length. The symbiotic algae belong to M. conductrix; the larger algae measure 4–7 × 4–7 µm; the smaller algal cells measure 2–5 × 2–5 µm.Geographic distribution: The investigated strains of syngen R1 were found in Europe: Göttingen, Germany; Lake Mondsee, Austria. In addition, this species has been reported from different places in Europe, Asia and North America (see details in Supplementary Table 1).Reference material: Strain SAG 27.96 and the clonal strain SAG 2645 derived from SAG 27.96 are available at the Culture Collection of Algae (SAG), University of Göttingen, Germany.Holotype: Two slides (one holotype, one paratype) with protargol-impregnated specimens from the clonal culture SAG 2645, which derived from the reference material SAG 27.96, isolated from the pond of the Old Botanical Garden of the University of Göttingen (Germany), have been deposited in the Oberösterreichisches Landesmuseum at Linz (LI, Austria).Zoobank Registration LSID: AFD967ED-BC2A-43FD-847E-5DF588BB025C.
Paramecium deuterobursaria sp. nov.Synonym: Paramecium bibursaria nom. inval.Description: The strains CCAP 1660/34 and CCAP 1660/36 belong to syngen R2 according to Greczek-Stachura et al.10,11 and differ from other syngens by their SSU and ITS rDNA sequences (OK318487). From morphology, the cells are ellipsoidal to broadly ellipsoidal and dorso-ventrally flattened in vivo. The cells measure 81–167 × 35–83 µm; the single macronucleus is located around mid-cell and measures 24–46 × 10–32 µm; the adjacent single compact micronucleus measures 10–18 × 5–9 µm, no micronucleus seen in live cells of strain CCAP 1660/34; the usually two (rarely one or three) contractile vacuoles, one in the anterior and one in the posterior cell portion have radial collecting channels and 1–3 excretory pores each; the number of ciliary rows/20 µm is 13–22; the length of the caudal cilia is 11–20 µm; the numerous trichocysts located in the cell cortex are 4–6 µm in length. The symbiotic algae belong to M. conductrix; the larger algae measure 5–7 × 4–7 µm; the smaller algal cells measure 3–5 × 2–5 µm.Geographic distribution: The investigated strains of syngen R2 were found in Europe: Zurich, Switzerland; Lake Piburg, Austria. In addition, this species has been reported from different places in Europe, Asia and Australia (see details in Supplementary Table 1).Reference material: Strain CCAP 1660/36 is available at the Culture Collection of Algae and Protozoa (CCAP) at the Scottish Association for Marine Science, Oban, Scotland.Holotype: Two slides (one holotype, one paratype) with protargol-impregnated specimens from the reference material CCAP 1660/36, isolated from Lake Piburg (Tyrol, Austria), have been deposited in the Oberösterreichisches Landesmuseum at Linz (LI, Austria).Zoobank Registration LSID: D1C20BE6-9A15-4A3D-A7E5-DFC31FF04679.
Paramecium tritobursaria sp. nov.Synonym: Paramecium tribursaria nom. inval.Description: The strains CCAP 1660/26, CCAP 1660/28 and CCAP 1660/31 belong to syngen R3 according to Greczek-Stachura et al.10,11 and differ from other syngens by their SSU and ITS rDNA sequences (MT231339). From morphology, the cells are ellipsoidal to broadly ellipsoidal and dorso-ventrally flattened in vivo. The cells measure 80–153 × 49–73 µm; the single macronucleus is located around mid-cell and measures 21–53 × 12–31 µm; the adjacent single compact micronucleus measures 9–17 × 3–6 µm; no micronucleus seen in live cells of strain CCAP 1660/28; the usually two (rarely one or three) contractile vacuoles, one in the anterior and one in the posterior cell portion have radial collecting channels and 1–3 excretory pores each; the number of ciliary rows/20 µm is 12–20; the length of the caudal cilia is 8–19 µm; the numerous trichocysts located in the cell cortex are 4–6 µm in length. The symbiotic algae belong to C. variabilis; the larger algae measure 4–7 × 3–6 µm; the smaller algal cells measure 3–5 × 2–4 µm.Geographic distribution: The investigated strains of syngen R3 were found in Europe and Asia: Lake Piburg, Austria; Tokyo, Japan; Khabarovsk region, Amur River, Russia. In addition, this species has been reported from different places in Europe, Asia, North and South America as well as in Australia (see details in Supplementary Table 1).Reference material: Strain CCAP 1660/26 is available at the Culture Collection of Algae and Protozoa (CCAP) at the Scottish Association for Marine Science, Oban, Scotland.Holotype: Two slides (one holotype, one paratype) with protargol-impregnated specimens from the reference material CCAP 1660/26, isolated from Japan, have been deposited in the Oberösterreichisches Landesmuseum at Linz (LI, Austria).Zoobank Registration LSID: CC0FBA7E-9E3A-4C37-B424-C9BFF2018EC0.
Paramecium tetratobursaria sp. nov.Synonym: Paramecium tetrabursaria nom. inval.Description: The strains CCAP 1660/25 and CCAP 1660/33 belong to syngen R4 according to Greczek-Stachura et al.10,11 and differ from other syngens by their SSU and ITS rDNA sequences (MT231347). From morphology, the cells are ellipsoidal to broadly ellipsoidal and dorso-ventrally flattened in vivo. The cells measure 65–179 × 37–79 µm; the single macronucleus is located around mid-cell and measures 18–53 × 10–29 µm; the adjacent single compact micronucleus measures 8–18 × 4–10 µm; the usually two (rarely one or three) contractile vacuoles, one in the anterior and one in the posterior cell portion have radial collecting channels and 1–3 excretory pores each; the number of ciliary rows/20 µm is 14–19; the length of the caudal cilia is 12–20 µm; the numerous trichocysts located in the cell cortex are 4–7 µm in length. The symbiotic algae belong to C. variabilis (CCAP 1660/25) and M. conductrix (CCAP 1660/33); the larger algae measure 3–6 × 3–6 µm; the smaller algal cells measure 2–5 × 1–4 µm.Geographic distribution: The investigated strains of syngen R4 are found in North- and South America: Burlington, North Carolina, USA; San Pedro de la Paz, Laguna Grande, Chile. In addition, this species has been reported from Europe (see details in Supplementary Table 1).Reference material: Strain CCAP 1660/25 is available at the Culture Collection of Algae and Protozoa (CCAP) at the Scottish Association for Marine Science, Oban, Scotland.Holotype: Two slides (one holotype, one paratype) with protargol-impregnated specimens from the reference material CCAP 1660/25, isolated from a pond in Burlington (North Carolina, USA), have been deposited in the Oberösterreichisches Landesmuseum at Linz (LI, Austria).Zoobank Registration LSID: 78BA9923-07A9-4918-AD7C-9E5E15CC9CDB.
Paramecium pentobursaria sp. nov.Synonym: Paramecium pentabursaria nom. inval.Description: The strain CCAP 1660/30 belongs to syngen R5 according to Greczek-Stachura et al.10,11 and differs from other syngens by their SSU and ITS rDNA sequences (MT231348). From morphology, the cells are ellipsoidal to broadly ellipsoidal and dorso-ventrally flattened in vivo. The cells measure 161–194 × 76–99 µm; the single macronucleus is located around mid-cell and measures 24–47 × 19–31 µm; the adjacent single compact micronucleus measures 13–20 × 4–9 µm; the usually two (rarely one or three) contractile vacuoles, one in the anterior and one in the posterior cell portion have radial collecting channels and 1–4 excretory pores each; the number of ciliary rows/20 µm is 13–19; the length of the caudal cilia is 14–25 µm; the numerous trichocysts located in the cell cortex are 5–7 µm in length. The symbiotic algae belong to C. variabilis; the larger algae measure 5–6 × 5–6 µm; the smaller algal cells measure 4–5 × 3–4 µm.Geographic distribution: The investigated strain of Syngen R5 was found in Europe: Astrakhan Nature Reserve, Russia.Reference material: Strain CCAP 1660/30 is available at the Culture Collection of Algae and Protozoa (CCAP) at the Scottish Association for Marine Science, Oban, Scotland.Holotype: Two slides (one holotype, one paratype) with protargol-impregnated specimens from the reference material CCAP 1660/30, isolated from Astrakhan Nature Reserve (Russia), have been deposited in the Oberösterreichisches Landesmuseum at Linz (LI, Austria).Zoobank Registration LSID: 6629FA71-E00F-48C6-83AB-61C0CA4823B6.Syngen Affiliation related to Ciliate Morphology, Endosymbionts and Geographic DistributionPearson-correlations of morphometric, syngen-specific and endosymbiont datasets of the P. bursaria strains revealed four significant positive correlations (p  r  > 0.75) between ciliate cell length (BLEN) and width (BWID), BWID and macronucleus width (MACWID), as well as length and width of large symbiotic algae (LSALEN and LSAWID; Fig. 7).Figure 7Pearson-correlations of morphometric, symbiont and syngen data of Paramecium strains under study. Colored dots indicate the strength of correlation, and the size of dots represent p-values. Bold squares highlight significant correlations, with − 0.75  > r  > 0.75 and p  1, accounting for 73.1% variation in total (Supplementary Table 3). Principal component axis 1 (PC1) appears to be most negatively weighted by syngen (SYN) and width of the macronucleus (MACWID), separating CCAP 1660/30 and CCAP 1660/33 from the other strains. Principal component axis 2 (PC2) is primarily positively influenced by symbiotic algae characteristics (LSALEN, LSAWID, small symbiotic algal length (SSALEN) and width (SSAWID)) and, ciliate cell length (BLEN) and width (BWID; Supplementary Table 4), partitioning strain PB-25, CCAP 1660/26 and CCAP 1660/36 from CCAP 1660/31 and SAG 27.96 (Fig. 8).Figure 8PCA of morphometric data of Paramecium bursaria strains. Only the top eight contributing variables are shown.Full size imageThe redundancy analysis (RDA; Fig. 9) revealed a large difference between morphometric features and the tested set of explanatory variables (i.e., algal species (ALSPEC), LSAWID, SSALEN, SYN and GEO) as only 26.9% of the total variation could be explained.Figure 9Ordination diagram for redundancy analysis (RDA) of morphometric data and shown syngen (SYN), geographic region (GEO), and algal features (ALSPEC, LSAWID and SSALEN) as explanatory features.Full size image More

LaBarbera, M. The evolution and ecology of body size. In Patterns and Processes in the History of Life (eds Raup, D. M. & Jablonski, D.) 69–98 (Springer, 1986).
Google Scholar 
Peters, R. H. & Peters, R. H. The Ecological Implications of Body Size Vol. 2 (Cambridge University Press, 1986).
Google Scholar 
Klingenberg, C. P. & Spence, J. On the role of body size for life-history evolution. Ecol. Entomol. 22(1), 55–68 (1997).
Google Scholar 
Blanckenhorn, W. U. The evolution of body size: What keeps organisms small?. Q. Rev. Biol. 75(4), 385–407 (2000).CAS 
PubMed 

Google Scholar 
Sibly, R. M. & Brown, J. H. Effects of body size and lifestyle on evolution of mammal life histories. Proc. Natl. Acad. Sci. USA 104(45), 17707–17712 (2007).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Schmidt-Nielsen, K. Scaling in biology: The consequences of size. J. Exp. Zool. 194(1), 287–307 (1975).CAS 
PubMed 

Google Scholar 
Calder, W. A. Size, Function, and Life History (Courier Corporation, 1996).
Google Scholar 
Gould, S. J. Allometry and size in ontogeny and phylogeny. Biol. Rev. 41(4), 587–638 (1966).CAS 
PubMed 

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

Google Scholar 
Gearty, W. & Payne, J. L. Physiological constraints on body size distributions in Crocodyliformes. Evolution 74(2), 245–255 (2020).PubMed 

Google Scholar 
Maurer, B. A., Brown, J. H. & Rusler, R. D. The micro and macro in body size evolution. Evolution 46(4), 939–953 (1992).PubMed 

Google Scholar 
Hone, D. W. & Benton, M. J. The evolution of large size: How does Cope’s Rule work?. Trends Ecol. Evol. 20(1), 4–6 (2005).PubMed 

Google Scholar 
Reeve, J. P. & Fairbairn, D. J. Predicting the evolution of sexual size dimorphism. J. Evol. Biol. 14(2), 244–254 (2001).
Google Scholar 
Blanckenhorn, W. U. Behavioral causes and consequences of sexual size dimorphism. Ethology 111(11), 977–1016 (2005).
Google Scholar 
Wu, H., Jiang, T., Huang, X. & Feng, J. Patterns of sexual size dimorphism in horseshoe bats: Testing Rensch’s rule and potential causes. Sci. Rep. 8(1), 1–13 (2018).ADS 

Google Scholar 
Cox, R. M., Butler, M. A. & John-Alder, H. B. The evolution of sexual size dimorphism in reptiles. In Sex, Size and Gender Roles: Evolutionary Studies of Sexual Size Dimorphism (eds Fairbairn, D. J. et al.) 38–49 (Oxford University Press, 2007).
Google Scholar 
Stillwell, R. C., Blanckenhorn, W. U., Teder, T., Davidowitz, G. & Fox, C. W. Sex differences in phenotypic plasticity affect variation in sexual size dimorphism in insects: From physiology to evolution. Annu. Rev. Entomol. 55, 227–245 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
Rensch, B. Die Abhängigkeit der relativen Sexualdifferenz von der Körpergrösse. Bonn. Zool. Beitr. 1, 58–69 (1950).
Google Scholar 
Rensch, B. Evolution Above the Species Level (Columbia University Press, 1960).
Google Scholar 
Shine, R. Ecological causes for the evolution of sexual dimorphism: A review of the evidence. Q. Rev. Biol. 64(4), 419–461 (1989).CAS 
PubMed 

Google Scholar 
Portik, D. M., Blackburn, D. C. & McGuire, J. A. Macroevolutionary patterns of sexual size dimorphism among African tree frogs (Family: Hyperoliidae). J. Hered. 111(4), 379–391 (2020).PubMed 

Google Scholar 
Ceballos, C. P., Adams, D. C., Iverson, J. B. & Valenzuela, N. Phylogenetic patterns of sexual size dimorphism in turtles and their implications for Rensch’s rule. Evol. Biol. 40(2), 194–208 (2013).
Google Scholar 
Amado, T. F., Martinez, P. A., Pincheira-Donoso, D. & Olalla-Tárraga, M. Á. Body size distributions of anurans are explained by diversification rates and the environment. Glob. Ecol. Biogeogr. 30(1), 154–164 (2021).
Google Scholar 
Starostová, Z., Kubička, L. & Kratochvíl, L. Macroevolutionary pattern of sexual size dimorphism in geckos corresponds to intraspecific temperature-induced variation. J. Evol. Biol. 23(4), 670–677 (2010).PubMed 

Google Scholar 
Herczeg, G., Gonda, A. & Merilä, J. Rensch’s rule inverted–female-driven gigantism in nine-spined stickleback Pungitius pungitius. J. Anim. Ecol. 79(3), 581–588 (2010).PubMed 

Google Scholar 
Liao, W. B. & Chen, W. Inverse Rensch’s rule in a frog with female-biased sexual size dimorphism. Naturwissenschaften 99(5), 427–431 (2012).ADS 
CAS 
PubMed 

Google Scholar 
Cooper, M. I. Sexual size dimorphism and the rejection of Rensch’s rule in Diplopoda (Arthropoda). J. Entomol. Zool. Stud. 6(1), 1582–1587 (2018).
Google Scholar 
Cheng, R. C. & Kuntner, M. Phylogeny suggests nondirectional and isometric evolution of sexual size dimorphism in argiopine spiders. Evolution 68(10), 2861–2872 (2014).PubMed 

Google Scholar 
Webb, T. J. & Freckleton, R. P. Only half right: Species with female-biased sexual size dimorphism consistently break Rensch’s rule. PLoS ONE 2(9), e897 (2007).ADS 
PubMed 
PubMed Central 

Google Scholar 
Gaston, K. J., Chown, S. L. & Evans, K. L. Ecogeographical rules: Elements of a synthesis. J. Biogeogr. 35(3), 483–500 (2008).
Google Scholar 
Olalla-Tárraga, M. Á. & Rodríguez, M. Á. Energy and interspecific body size patterns of amphibian faunas in Europe and North America: Anurans follow Bergmann’s rule, urodeles its converse. Glob. Ecol. Biogeogr. 16(5), 606–617 (2007).
Google Scholar 
Olalla-Tárraga, M. Á., Diniz-Filho, J. A. F., Bastos, R. P. & Rodríguez, M. Á. Geographic body size gradients in tropical regions: Water deficit and anuran body size in the Brazilian Cerrado. Ecography 32(4), 581–590 (2009).
Google Scholar 
Gouveia, S. F. & Correia, I. Geographical clines of body size in terrestrial amphibians: Water conservation hypothesis revisited. J. Biogeogr. 43(10), 2075–2084 (2016).
Google Scholar 
Pincheira-Donoso, D., Meiri, S., Jara, M., Olalla-Tárraga, M. Á. & Hodgson, D. J. Global patterns of body size evolution are driven by precipitation in legless amphibians. Ecography 42(10), 1682–1690 (2019).
Google Scholar 
Nevo, E. Adaptive color polymorphism in cricket frogs. Evolution 27(3), 353–367 (1973).PubMed 

Google Scholar 
Ashton, K. G. Do amphibians follow Bergmann’s rule?. Can. J. Zool. 80(4), 708–716 (2002).MathSciNet 

Google Scholar 
Bergmann, C. Ueber die Verhältnisse der Wärmeökonomie der Thiere zu ihrer Grösse. Göttinger Studien. 1, 595–708 (1847).
Google Scholar 
Olalla-Tárraga, M. Á., Rodríguez, M. Á. & Hawkins, B. A. Broad-scale patterns of body size in squamate reptiles of Europe and North America. J. Biogeogr. 33(5), 781–793 (2006).
Google Scholar 
Trullas, S. C., van Wyk, J. H. & Spotila, J. R. Thermal melanism in ectotherms. J. Therm. Biol. 32(5), 235–245 (2007).
Google Scholar 
Rodríguez, M. Á., López-Sañudo, I. L. & Hawkins, B. A. The geographic distribution of mammal body size in Europe. Glob. Ecol. Biogeogr. 15(2), 173–181 (2006).
Google Scholar 
Olalla-Tárraga, M. Á., Diniz-Filho, J. A. F., Bastos, R. P. & Rodriguez, M. A. Geographic body size gradients in tropical regions: Water deficit and anuran body size in the Brazilian Cerrado. Ecography 32(4), 581–590 (2009).
Google Scholar 
Womack, M. C. & Bell, R. C. Two-hundred million years of anuran body-size evolution in relation to geography, ecology and life history. J. Evol. Biol. 33(10), 1417–1432 (2020).PubMed 

Google Scholar 
Frost, D. R. Amphibian Species of the World: An online reference, version 6. http://research.amnh.org/herpetology/amphibia/index.php. Accessed 12 July 2021 (2021).
Acevedo, A. A., Armesto, O. & Palma, R. E. Two new species of Pristimantis (Anura: Craugastoridae) with notes on the distribution of the genus in northeastern Colombia. Zootaxa 4750(4), 499–523 (2020).
Google Scholar 
Heinicke, M. P., Duellman, W. E. & Hedges, S. B. Major Caribbean and Central American frog faunas originated by ancient oceanic dispersal. Proc. Natl. Acad. Sci. USA 104(24), 10092–10097 (2007).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Pinto-Sánchez, N. R. et al. The great American biotic interchange in frogs: Multiple and early colonization of Central America by the South American genus Pristimantis (Anura: Craugastoridae). Mol. Phylogenet. Evol. 62(3), 954–972 (2012).PubMed 

Google Scholar 
Zumel, D., Buckley, D. & Ron, S. R. The Pristimantis trachyblepharis species group, a clade of miniaturized frogs: Description of four new species and insights into the evolution of body size in the genus. Zool. J. Linn. Soc. zlab044 (2021).Pincheira-Donoso, D. et al. The multiple origins of sexual size dimorphism in global amphibians. Glob. Ecol. Biogeogr. 30(2), 443–458 (2021).
Google Scholar 
Woolbright, L. L. Sexual selection and size dimorphism in anuran amphibia. Am. Nat. 121(1), 110–119 (1983).
Google Scholar 
Nali, R. C., Zamudio, K. R., Haddad, C. F. & Prado, C. P. Size-dependent selective mechanisms on males and females and the evolution of sexual size dimorphism in frogs. Am. Nat. 184(6), 727–740 (2014).PubMed 

Google Scholar 
Hill, R. et al. Herpetological husbandry observations on the captive reproduction of gaige’s rain frog Pristimantis gaigeae (Dunn 1931). Herpetol. Rev. 41(4), 465 (2010).
Google Scholar 
Rojas-Rivera, A., Cortés-Bedoya, S., Gutiérrez-Cárdenas, P. D. A. & Castellanos, J. M. Pristimantis achatinus (Cachabi robber frog). Parental care and clutch size. Herpetol. Rev. 42, 588–589 (2011).
Google Scholar 
Granados-Pérez, Y. & Ramirez-Pinilla, M. P. Reproductive phenology of three species of Pristimantis in an Andean cloud forest. Revista Acad. Colomb. Ci. Exact. 44(173), 1083–1098 (2020).
Google Scholar 
Levy, D. L. & Heald, R. Biological scaling problems and solutions in amphibians. Cold Spring Harb. Perspect. Biol. 8(1), a019166 (2016).PubMed Central 

Google Scholar 
O’Donnell, M. S. & Ignizio, D. A. Bioclimatic predictors for supporting ecological applications in the conterminous United States. US Geol. Survey Data Series. 691(10), 4–9 (2012).
Google Scholar 
Valenzuela-Sánchez, A., Cunningham, A. A. & Soto-Azat, C. Geographic body size variation in ectotherms: Effects of seasonality on an anuran from the southern temperate forest. Front. Zool. 12(1), 1–10 (2015).
Google Scholar 
Parsons, J. J. The northern Andean environment. Mt. Res. Dev. 2(3), 253–264 (1982).
Google Scholar 
Navas, C. A., Carvajalino-Fernández, J. M., Saboyá-Acosta, L. P., Rueda-Solano, L. A. & Carvajalino-Fernández, M. A. The body temperature of active amphibians along a tropical elevation gradient: Patterns of mean and variance and inference from environmental data. Funct. Ecol. 27(5), 1145–1154 (2013).
Google Scholar 
Swemmer, A. M., Knapp, A. K. & Snyman, H. A. Intra-seasonal precipitation patterns and above-ground productivity in three perennial grasslands. J. Ecol. 95(4), 780–788 (2007).
Google Scholar 
Losos, J. B. Lizards in an Evolutionary Tree: Ecology and Adaptive Radiation of Anoles (Univ. of California Press, 2011).
Google Scholar 
Pincheira-Donoso, D. & Hunt, J. Fecundity selection theory: Concepts and evidence. Biol. Rev. 92(1), 341–356 (2017).PubMed 

Google Scholar 
Morrison, C. & Hero, J. M. Geographic variation in life-history characteristics of amphibians: A review. J. Anim. Ecol. 72(2), 270–279 (2003).
Google Scholar 
Morrow, C. B., Ernest, S. M. & Kerkhoff, A. J. Macroevolution of dimensionless life-history metrics in tetrapods. Proc. Royal Soc. B. 288, 20210200 (2021).
Google Scholar 
Revell, L. J., Harmon, L. J. & Collar, D. C. Phylogenetic signal, evolutionary process, and rate. Syst. Biol. 57(4), 591–601 (2008).PubMed 

Google Scholar 
Kamilar, J. M. & Cooper, N. Phylogenetic signal in primate behaviour, ecology and life history. Philos. Trans. R. Soc. Lond. B Biol. Sci. 368(1618), 20120341 (2013).PubMed 
PubMed Central 

Google Scholar 
Meyer, A. L. & Wiens, J. J. Estimating diversification rates for higher taxa: BAMM can give problematic estimates of rates and rate shifts. Evolution 72(1), 39–53 (2018).PubMed 

Google Scholar 
Rabosky, D. L. et al. Rates of speciation and morphological evolution are correlated across the largest vertebrate radiation. Nat. Commun. 4(1), 1–8 (2013).
Google Scholar 
Mendoza, A. M., Ospina, O. E., Cárdenas-Henao, H. & García-R, J. C. A likelihood inference of historical biogeography in the world’s most diverse terrestrial vertebrate genus: Diversification of direct-developing frogs (Craugastoridae: Pristimantis) across the Neotropics. Mol. Phylogenet. Evol. 85, 50–58 (2015).PubMed 

Google Scholar 
Baker, J., Meade, A., Pagel, M. & Venditti, C. Adaptive evolution toward larger size in mammals. Proc. Natl. Acad. Sci. USA 112(16), 5093–5098 (2015).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Hariharan, I. K., Wake, D. B. & Wake, M. H. Indeterminate growth: Could it represent the ancestral condition?. Cold Spring Harb. Perspect. Biol. 8(2), a019174 (2016).PubMed Central 

Google Scholar 
Amado, T. F., Bidau, C. J. & Olalla-Tárraga, M. Á. Geographic variation of body size in New World anurans: Energy and water in a balance. Ecography 42(3), 456–466 (2019).
Google Scholar 
Watters, J. L., Cummings, S. T., Flanagan, R. L. & Siler, C. D. Review of morphometric measurements used in anuran species descriptions and recommendations for a standardized approach. Zootaxa 4072, 477–495 (2016).PubMed 

Google Scholar 
Lovich, J. E. & Gibbons, J. W. A review of techniques for quantifying sexual size dimorphism. Growth Dev. Aging. 56, 269–269 (1992).CAS 
PubMed 

Google Scholar 
Lanfear, R., Calcott, B., Ho, S. Y. & Guindon, S. PartitionFinder: Combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol. Biol. Evol. 29(6), 1695–1701 (2012).CAS 
PubMed 

Google Scholar 
Drummond, A. J. & Rambaut, A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol. 7(1), 1–8 (2007).
Google Scholar 
Drummond, A. J., Ho, S. Y., Phillips, M. J. & Rambaut, A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 4(5), e88 (2006).PubMed 
PubMed Central 

Google Scholar 
Rambaut, A. FigTree, A Graphical Viewer of Phylogenetic Trees. (2007)Olalla-Tárraga, M. A., Bini, L. M., Diniz-Filho, J. A. & Rodríguez, M. Á. Cross-species and assemblage-based approaches to Bergmann’s rule and the biogeography of body size in Plethodon salamanders of eastern North America. Ecography 33(2), 362–368 (2010).
Google Scholar 
QGIS.org. QGIS Geographic Information System. QGIS Association. http://www.qgis.org. Accessed 10 July 2021 (2022).Wei, T. et al. Package ‘corrplot’. Statistician. 56(316), e24 (2017).
Google Scholar 
James, F. C. Geographic size variation in birds and its relationship to climate. Ecology 51(3), 365–390 (1970).
Google Scholar 
Hawkins, B. A. & Felizola Diniz-Filho, J. A. Beyond Rapoport’s rule: Evaluating range size patterns of New World birds in a two-dimensional framework. Glob. Ecol. Biogeogr. 15(5), 461–469 (2006).
Google Scholar 
Eager, C. standardize: Tools for standardizing variables for regression in R. R package version 0.21 (2017).Meireles, J. E., O’Meara, B. & Cavender-Bares, J. Linking leaf spectra to the plant tree of life. In Remote Sensing of Plant Biodiversity (eds Cavender-Bares, J. et al.) 155–172 (Springer, 2010).
Google Scholar 
Pagel, M. Inferring the historical patterns of biological evolution. Nature 401(6756), 877–884 (1999).ADS 
CAS 
PubMed 

Google Scholar 
Pagel, M. The maximum likelihood approach to reconstructing ancestral character states of discrete characters on phylogenies. Syst. Biol. 48(3), 612–622 (1999).
Google Scholar 
Revell, L. J. Phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol. 3(2), 217–223 (2012).
Google Scholar 
Revell, L. J. Two new graphical methods for mapping trait evolution on phylogenies. Methods Ecol. Evol. 4(8), 754–759 (2013).
Google Scholar 
Rabosky, D. L. et al. BAMM tools: An R package for the analysis of evolutionary dynamics on phylogenetic trees. Methods Ecol. Evol. 5(7), 701–707 (2014).
Google Scholar 
Rabosky, D. L. Automatic detection of key innovations, rate shifts, and diversity-dependence on phylogenetic trees. PLoS ONE 9(2), e89543 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
Rabosky, D. L., Mitchell, J. S. & Chang, J. Is BAMM flawed? Theoretical and practical concerns in the analysis of multi-rate diversification models. Syst. Biol. 66(4), 477–498 (2017).PubMed 
PubMed Central 

Google Scholar 
Plummer, M., Best, N., Cowles, K. & Vines, K. CODA: Convergence diagnosis and output analysis for MCMC. R News. 6(1), 7–11 (2006).
Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna, 2021). https://www.R-project.org. Accessed 1 June 2021 (2021).Fairbairn, D. J. Allometry for sexual size dimorphism: Pattern and process in the coevolution of body size in males and females. Annu. Rev. Ecol. Evol. Syst. 28(1), 659–687 (1997).
Google Scholar 
Fairbairn, D. J. Allometry for sexual size dimorphism: Testing two hypotheses for Rensch’s rule in the water strider Aquarius remigis. Am. Nat. 166(S4), S69–S84 (2005).PubMed 

Google Scholar 
Visser, A. G., Beevers, L. & Patidar, S. Complexity in hydroecological modelling: A comparison of stepwise selection and information theory. River Res. Appl. 34(8), 1045–1056 (2018).
Google Scholar 
Calcagno, V. & de Mazancourt, C. glmulti: An R package for easy automated model selection with (generalized) linear models. J. Stat. Softw. 34(12), 1–29 (2010).
Google Scholar 
Callaghan, S., Guilyardi, E., Steenman-Clark, L. & Morgan, M. The METAFOR project. in Earth System Modelling-Volume 1 (Springer, 2013).Garamszegi, L. Z. & Mundry, R. Multimodel-inference in comparative analyses. In Modern Phylogenetic Comparative Methods and THEIR Application in Evolutionary Biology (ed. Garamszegi, L. Z.) 305–331 (Springer Berlin, 2014).
Google Scholar  More

Ault, T. R. On the essentials of drought in a changing climate. Science 368, 256–260 (2020).CAS 

Google Scholar 
Rodell, M. et al. Emerging trends in global freshwater availability. Nature 557, 651–659 (2018).CAS 

Google Scholar 
Green, J. K. et al. Large influence of soil moisture on long-term terrestrial carbon uptake. Nature 565, 476–479 (2019).CAS 

Google Scholar 
Reichstein, M. et al. Climate extremes and the carbon cycle. Nature 500, 287–295 (2013).CAS 

Google Scholar 
Goulden, M. L. & Bales, R. C. California forest die-off linked to multi-year deep soil drying in 2012–2015 drought. Nat. Geosci. 12, 632–637 (2019).CAS 

Google Scholar 
Brodribb, T. J., Powers, J., Cochard, H. & Choat, B. Hanging by a thread? Forests and drought. Science 368, 261–266 (2020).CAS 

Google Scholar 
Seneviratne, S. I. et al. Investigating soil moisture–climate interactions in a changing climate: a review. Earth Sci. Rev. 99, 125–161 (2010).CAS 

Google Scholar 
Short Gianotti, D. J., Rigden, A. J., Salvucci, G. D. & Entekhabi, D. Satellite and station observations demonstrate water availability’s effect on continental-scale evaporative and photosynthetic land surface dynamics. Water Resour. Res. 55, 540–554 (2019).
Google Scholar 
Anderegg, W. R. L., Trugman, A. T., Bowling, D. R., Salvucci, G. & Tuttle, S. E. Plant functional traits and climate influence drought intensification and land–atmosphere feedbacks. Proc. Natl Acad. Sci. USA 116, 14071–14076 (2019).CAS 

Google Scholar 
Jung, M. et al. Recent decline in the global land evapotranspiration trend due to limited moisture supply. Nature 467, 951–954 (2010).CAS 

Google Scholar 
Novick, K. A. et al. The increasing importance of atmospheric demand for ecosystem water and carbon fluxes. Nat. Clim. Chang. 6, 1023–1027 (2016).CAS 

Google Scholar 
Liu, Y., Kumar, M., Katul, G. G., Feng, X. & Konings, A. G. Plant hydraulics accentuates the effect of atmospheric moisture stress on transpiration. Nat. Clim. Chang. 10, 691–695 (2020).CAS 

Google Scholar 
Gupta, A., Rico-Medina, A. & Caño-Delgado, A. I. The physiology of plant responses to drought. Science 368, 266–269 (2020).CAS 

Google Scholar 
Wolf, S. et al. Warm spring reduced carbon cycle impact of the 2012 US summer drought. Proc. Natl Acad. Sci. USA 113, 5880–5885 (2016).CAS 

Google Scholar 
Teuling, A. J. et al. Evapotranspiration amplifies European summer drought. Geophys. Res. Lett. 40, 2071–2075 (2013).
Google Scholar 
Mastrotheodoros, T. et al. More green and less blue water in the Alps during warmer summers. Nat. Clim. Chang. 10, 155–161 (2020).
Google Scholar 
Peterson, T. J., Saft, M., Peel, M. C. & John, A. Watersheds may not recover from drought. Science 372, 745–749 (2021).CAS 

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

Google Scholar 
Massmann, A., Gentine, P. & Lin, C. When does vapor pressure deficit drive or reduce evapotranspiration? J. Adv. Model. Earth Syst. 11, 3305–3320 (2019).
Google Scholar 
Grossiord, C. et al. Plant responses to rising vapor pressure deficit. New Phytol. 226, 1550–1566 (2020).
Google Scholar 
Orth, R. & Destouni, G. Drought reduces blue-water fluxes more strongly than green-water fluxes in Europe. Nat. Commun. 9, 3602 (2018).
Google Scholar 
Pendergrass, A. G. et al. Flash droughts present a new challenge for subseasonal-to-seasonal prediction. Nat. Clim. Chang. 10, 191–199 (2020).
Google Scholar 
Chu, H., Baldocchi, D. D., John, R., Wolf, S. & Reichstein, M. Fluxes all of the time? A primer on the temporal representativeness of FLUXNET. J. Geophys. Res. Biogeosci. 122, 289–307 (2017).
Google Scholar 
Ukkola, A. M. et al. Land surface models systematically overestimate the intensity, duration and magnitude of seasonal-scale evaporative droughts. Environ. Res. Lett. 11, 104012 (2016).
Google Scholar 
Trugman, A. T., Medvigy, D., Mankin, J. S. & Anderegg, W. R. L. Soil moisture stress as a major driver of carbon cycle uncertainty. Geophys. Res. Lett. 45, 6495–6503 (2018).
Google Scholar 
De Kauwe, M. G. et al. Do land surface models need to include differential plant species responses to drought? Examining model predictions across a mesic–xeric gradient in Europe. Biogeosciences 12, 7503–7518 (2015).
Google Scholar 
Dong, J., Lei, F. & Crow, W. T. Land transpiration–evaporation partitioning errors responsible for modeled summertime warm bias in the central United States. Nat. Commun. 13, 336 (2022).CAS 

Google Scholar 
Kennedy, D. et al. Implementing plant hydraulics in the Community Land Model, version 5. J. Adv. Model. Earth Syst. 11, 485–513 (2019).
Google Scholar 
Novick, K. A. et al. Confronting the water potential information gap. Nat. Geosci. 15, 158–164 (2022).CAS 

Google Scholar 
Liu, Y., Holtzman, N. M. & Konings, A. G. Global ecosystem-scale plant hydraulic traits retrieved using model–data fusion. Hydrol. Earth Syst. Sci. 25, 2399–2417 (2021).CAS 

Google Scholar 
Lin, Y.-S. et al. Optimal stomatal behaviour around the world. Nat. Clim. Chang. 5, 459–464 (2015).CAS 

Google Scholar 
Fan, Y., Miguez-Macho, G., Jobbágy, E. G., Jackson, R. B. & Otero-Casal, C. Hydrologic regulation of plant rooting depth. Proc. Natl Acad. Sci. USA 114, 10572–10577 (2017).CAS 

Google Scholar 
Anderegg, W. R. L. et al. Hydraulic diversity of forests regulates ecosystem resilience during drought. Nature 561, 538–541 (2018).CAS 

Google Scholar 
Lehmann, P., Merlin, O., Gentine, P. & Or, D. Soil texture effects on surface resistance to bare-soil evaporation. Geophys. Res. Lett. 45, 10398–10405 (2018).
Google Scholar 
Fatichi, S. et al. Soil structure is an important omission in Earth System Models. Nat. Commun. 11, 522 (2020).CAS 

Google Scholar 
McCormick, E. L. et al. Widespread woody plant use of water stored in bedrock. Nature 597, 225–229 (2021).CAS 

Google Scholar 
Baldocchi, D., Ma, S. & Verfaillie, J. On the inter- and intra-annual variability of ecosystem evapotranspiration and water use efficiency of an oak savanna and annual grassland subjected to booms and busts in rainfall. Glob. Chang. Biol. 27, 359–375 (2021).CAS 

Google Scholar 
Condon, L. E., Atchley, A. L. & Maxwell, R. M. Evapotranspiration depletes groundwater under warming over the contiguous United States. Nat. Commun. 11, 873 (2020).CAS 

Google Scholar 
Maxwell, R. M. & Condon, L. E. Connections between groundwater flow and transpiration partitioning. Science 353, 377–380 (2016).CAS 

Google Scholar 
Humphrey, V. et al. Sensitivity of atmospheric CO2 growth rate to observed changes in terrestrial water storage. Nature 560, 628–631 (2018).CAS 

Google Scholar 
Mankin, J. S., Seager, R., Smerdon, J. E., Cook, B. I. & Williams, A. P. Mid-latitude freshwater availability reduced by projected vegetation responses to climate change. Nat. Geosci. 12, 983–988 (2019).CAS 

Google Scholar 
Strassburg, B. B. N. et al. Global priority areas for ecosystem restoration. Nature 586, 724–729 (2020).CAS 

Google Scholar 
Zhao, M. et al. Ecological restoration impact on total terrestrial water storage. Nat. Sustain. 4, 56–62 (2021).
Google Scholar 
Feng, X. et al. Revegetation in China’s Loess Plateau is approaching sustainable water resource limits. Nat. Clim. Chang. 6, 1019–1022 (2016).
Google Scholar 
Guan, K. et al. Photosynthetic seasonality of global tropical forests constrained by hydroclimate. Nat. Geosci. 8, 284–289 (2015).CAS 

Google Scholar 
Fan, Y., Li, H. & Miguez-Macho, G. Global patterns of groundwater table depth. Science 339, 940–943 (2013).CAS 

Google Scholar 
Zhao, M., Geruo, A., Velicogna, I. & Kimball, J. S. A global gridded dataset of GRACE drought severity index for 2002–14: comparison with PDSI and SPEI and a case study of the Australia Millennium Drought. J. Hydrometeorol. 18, 2117–2129 (2017).
Google Scholar 
Watkins, M. M., Wiese, D. N., Yuan, D.-N., Boening, C. & Landerer, F. W. Improved methods for observing Earth’s time variable mass distribution with GRACE using spherical cap mascons. J. Geophys. Res. Solid Earth 120, 2648–2671 (2015).
Google Scholar 
Wiese, D. N., Landerer, F. W. & Watkins, M. M. Quantifying and reducing leakage errors in the JPL RL05M GRACE mascon solution. Water Resour. Res. 52, 7490–7502 (2016).
Google Scholar 
Adler, R. F. et al. The Version-2 Global Precipitation Climatology Project (GPCP) monthly precipitation analysis (1979–present). J. Hydrometeorol. 4, 1147–1167 (2003).
Google Scholar 
Sun, Q. et al. A review of global precipitation data sets: data sources, estimation, and intercomparisons. Rev. Geophys. 56, 79–107 (2018).
Google Scholar 
Gebremichael, M. et al. Error uncertainty analysis of GPCP monthly rainfall products: a data-based simulation study. J. Appl. Meteorol. 42, 1837–1848 (2003).
Google Scholar 
Rodell, M. et al. Basin scale estimates of evapotranspiration using GRACE and other observations. Geophys. Res. Lett. 31, L20504 (2004).
Google Scholar 
Major River Basins of the World (Global Runoff Data Centre, 2020).Pascolini-Campbell, M. A., Reager, J. T. & Fisher, J. B. GRACE-based mass conservation as a validation target for basin-scale evapotranspiration in the contiguous United States. Water Resour. Res. 56, e2019WR026594 (2020).
Google Scholar 
Fekete, B. M., Vörösmarty, C. J. & Grabs, W. High-resolution fields of global runoff combining observed river discharge and simulated water balances. Global Biogeochem. Cycles 16, 15-1–15-10 (2002).
Google Scholar 
Breiman, L. Random forests. Mach. Learn. 45, 5–32 (2001).
Google Scholar 
Myneni, R., Knyazikhin, Y. & Park, T (ed. NASA EOSDIS Land Processes DAAC) (2021).Hersbach, H. et al. The ERA5 global reanalysis. Q. J. R. Meteorol. Soc. 146, 1999–2049 (2020).
Google Scholar 
Sulla-Menashe, D. & Friedl, M. A. User Guide to Collection 6 MODIS Land Cover (MCD12Q1 and MCD12C1) Product (US Geological Survey, 2018).Zhao, M., Aa, G., Liu, Y. & Konings, A. Evapotranspiration frequently increases during droughts. Zenodo https://doi.org/10.5281/zenodo.6842054 (2022). More

Seneviratne, S. I. et al. in Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (eds Masson-Delmotte, V. et al.) 1513–1766 (Cambridge Univ. Press, 2021).
Google Scholar 
Hesketh, A. V. & Harley, C. D. G. Glob. Change Biol. https://doi.org/10.1111/gcb.16390 (2022).Article 

Google Scholar 
Jørgensen, L. B., Ørsted, M., Malte, H., Wang, T. & Overgaard, J. Nature https://doi.org/10.1038/s41586-022-05334-4 (2022).Article 

Google Scholar 
Sunday, J. M., Bates, A. E. & Dulvy, N. K. Nature Clim. Change 2, 686–690 (2012).Article 

Google Scholar 
Cossins, A. R. & Bowler, K. Temperature Biology of Animals (Chapman & Hall, 1987).
Google Scholar 
Dell, A. I., Pawar, S. & Savage, V. M. Proc. Natl Acad. Sci. USA 108, 10591–10596 (2011).Article 
PubMed 

Google Scholar 
Dillon, M. E. et al. Integr. Comp. Biol. 56, 14–30 (2016).Article 
PubMed 

Google Scholar 
Stillman, J. H. Physiology 34, 86–100 (2019).Article 
PubMed 

Google Scholar 
MacMillan, H. A. J. Exp. Biol. 222, jeb191593 (2019).Article 
PubMed 

Google Scholar 
Kingsolver, J. G. & Umbanhowar, J. J. Exp. Biol. 221, jeb167858 (2018).Article 
PubMed 

Google Scholar  More

OverviewWe compiled systematically sampled empirical taxa occurrence across the landscape, and inferentially assembled respective blue and green local food webs by combining these data with a metaweb approach. We quantified key properties of the inferred food webs, then analysed with GIS-derived environmental information how focal food-web metrics change along elevation and among different land-use types in blue versus green systems. Details are given below.Assemble food webs using a metaweb approachWe applied a metaweb method to obtain the composition and structure of multiple local food webs across a landscape spatial scale10. A metaweb is an accumulation of all interactions (here, trophic relationships) among the focal taxa. In this study, we built our metaweb based on known trophic interactions derived from literature and published datasets, which themselves were all based on primary empirical natural history observations. We further complemented or refined the trophic interactions in the metaweb based on expert knowledge of primary observations that are not yet published or only accessible in grey literature. The expert knowledge covers authors and collaborators who have specific natural history knowledge on Central European plants, herbivorous insects, birds, fish, and aquatic invertebrates. Importantly, these observations were all based on empirical observations and/or unpublished data accumulated over considerable field research experience. The respective literature we referred, as well as the metaweb itself with information source of each trophic link (online repository), are provided in Supplementary Methods. By assuming that any interaction in the metaweb will realise if the interacting taxa co-occur, the metaweb approach allows an inference of local food webs if taxa occurrence is known. Such an assumption of fixed diets may lead to an over-estimation of the locally realised trophic links32, as it essentially ignores the possible intraspecific diet variation caused by resource availability61,62, predation risk63, temperature64, ontogenetic shift65, or other genetic and environmental sources66. Therefore, the food webs we inferred systematically using this method capture trophic relationships driven by community composition (species presence versus absence) but not the above-mentioned processes. Nonetheless, since the trophic interactions were based on empirical observations, the fixed diets can be seen as collapsing all intraspecific variations of diet-determining traits (or trait-matching) at species level, within which we know realisable interactions surely exist. This, together with co-occurrence as a pre-requisite, gives realistic boundaries for the potential interaction realisation, which is plausible and non-biased when applying to localised sites. With this approach, we were addressing a systematic comparison among potential local food webs between the blue and green systems and across the selected gradients. For sensitivity analyses considering the potential inaccuracy of the metaweb approach mentioned here, please see further below Food-web metrics and analyses and Supplementary Discussion.We compiled taxa occurrence of four terrestrial and two aquatic broad taxonomic groups (“focal groups”) to assemble local green and blue communities, respectively and independently, based on the well-resolved data available. Each focal group referred to a distinct taxonomic group, and the within- and among-group trophic relationships captured most of the realised interactions. These focal groups were vascular plants, butterflies, grasshoppers, and birds in the green biome, and stream invertebrates and fishes in the blue biome. Notably, with “butterflies” we refer to their larval stage and accordingly their mostly-herbivorous trophic interactions throughout this study. Larval interactions were also the predominant interaction assessed for stream invertebrates (i.e., all interactions of stream invertebrates focussed on their aquatic stage, which is predominant larval). The occurrence data of these focal groups were compiled from highly standardised multiple-year empirical surveys of various authorities, all conducted by trained biologists with fixed protocols (Supplementary Methods). The information across sites should thus be representative and can be up-scaled to the landscape. The occurrences of plants, butterflies, birds, and stream invertebrates were from the Biodiversity Monitoring Switzerland programme (BDM Coordination Office67) managed by the Swiss Federal Office for the Environment (BAFU/FOEN). The occurrences of grasshoppers and fishes were from the Swiss database on faunistic records, info fauna (CSCF), where we further complemented fish occurrence from the data of Progetto Fiumi Project (Eawag). In terms of biological resolution, taxa were resolved to species level in most cases, while the plant and butterfly groups included some multi-species complexes. Insects of the order Ephemeroptera, Plecoptera, and Trichoptera were resolved to species, while all other stream invertebrates were resolved to family level. These were each treated as a node later in our food-web assembly, and referred to as “species”, as the species within such complexes and families mostly share the same trophic role. Spatially, the occurrence datasets adopted coordinates resolved to 1 × 1 km2. The species that were recorded in the same 1 × 1 km2 grid were considered to co-occurred. We took the co-occurring four/two focal groups to form local green/blue local communities, respectively. To obtain better co-occurrence across group-specific data from different sources (e.g., BDM and info fauna), we intentionally coarsened the grasshopper and fish occurrence to 5 × 5 km2 coordinates. This is arguably a biologically acceptable approximation considering the high mobility of these two groups. Also, we only included known stream-borne fishes and dropped pure lake-borne ones to match our stream-only invertebrate occurrence data. Across all 462 green and 465 blue communities we assembled, we covered 2016 plant, 191 butterfly, 109 grasshopper, 155 bird, 248 stream invertebrate, and 78 stream fish species. Unlike the knowledge of plant occurrence in green communities, we did not have detailed occurrence information of the basal components (e.g., primary producers) in blue ones. Therefore, we assumed three mega nodes—namely plant (including all alive or dead plant materials), plankton (including zooplankton, phytoplankton, and other algae), and detritus—as the basal nodes occurring in all blue communities, without further discrimination of identities or biology within. These adding to our focal groups thus cover major taxonomic groups as well as trophic roles from producers to top consumers in both blue and green systems.Taking the above-assembled local communities then drawing trophic links among species (nodes) according to the metaweb yielded the local food webs (illustrated in Fig. 1), representatively covering the whole Swiss area. Notably, although our understanding of trophic interactions indeed encompassed some links across the blue and green taxa (e.g., between piscivorous birds and fishes), our occurrence datasets did not present sufficient spatial grids where these taxa co-occur. We, therefore, did not include such links, nor assembled blue-green interconnected food webs, but the blue and green food webs separately instead (but see Supplementary Discussion). Also, we dropped isolated nodes, i.e., basal nodes without any co-occurring consumer and consumer nodes without any co-occurring resource, from the inferred food webs. These could possibly be passing-by species that were recorded but had no trophic interaction locally, or those that interact with non-focal taxa whose occurrence information was unknown to us. We thus had to exclude them to focus on evidence-supported occurrences and trophic interactions. Nonetheless, across all cases, isolated nodes were rather rare (averaged less than 3% of species occurred in either blue or green communities).Environmental dataWe acquired environmental data across all of Switzerland (42,000 km2) on a 1 × 1 km2 grid basis (i.e., values are averaged over the grid) from GIS databases, with which we mapped environmental conditions to the grids where we assembled food webs. These included: topographical information from DHM25 (Swisstopo, FOT), land-cover information from CLC (EEA), and climate information (averaged over the decade of 2005–2015) from CHELSA. Among environmental variables, elevation and temperature are essentially highly correlated. In this study, we took elevation as the focal environmental gradient throughout, as after accounting for the main effects of elevation on temperature, the residual temperature was not a good predictor of the food-web metrics we looked at (see next section, and Supplementary Table 4). In other words, by analysing along the elevation gradient, we already captured most of the temperature influences on food webs. Based on the labels provided by the GIS databases, we categorised the originally detailed land cover into the five major land-use types that we used in this study, namely forest, scrubland, open space, farmland, and urban area. Forest includes broad-leaved, coniferous, and mixed forests. Scrub includes bushy and herbaceous vegetation, heathlands, and natural grasslands. Open space encompasses sparsely vegetated areas, such as dunes, bare rocks, glaciers and perpetual snow. Farmland include any form of arable, pastures, and agro-forestry areas. Finally, urban area is where artificial constructions and infrastructure prevail. As each grid could contain multiple land-use types, we then defined the dominant land-use type of the grid as any of the five above that occupied more than 50% of the grid’s area. Analyses separated by land-use types with subsetted food webs (land-use-specific analyses) were based on the grids’ dominant land-use type. There were a few grids where the dominant land-use type did not belong to the focal major five, e.g., wetlands or water bodies, and a few where no single land-use type covered more than 50% of the area. Food webs of these grids were still included in the overall analyses but excluded from any land-use-specific analyses (as revealed in the difference in sample sizes between all versus land-use type subsetted food webs in Fig. 2; analyses details below).Food-web metrics and analysesWe quantified five metrics as the measures of the food webs’ structural and ecological properties. For the fundamental structure of the food webs, the number of nodes (“No. Nodes”) reflects the size of the web, meanwhile represents local species richness (though the few isolated nodes were excluded as above-mentioned). Connectance is the proportion of realised links among all potential ones (thus bounded 0–1), reflecting how connected the web is. We also derived holistic topological measures, namely nestedness and modularity. Nestedness of a food web, on the one hand, describes the tendency that some nodes’ narrower diets being subsets of other’s broader diets. We adopted a recently developed UNODF index68 (bounded 0–1) that is especially suitable for quantifying such a feature in our unipartite food webs. On the other hand, modularity (bounded 0–1 with our index) reflects the tendency of a food web to form modules, where nodes are highly connected within but only loosely connected between. Nestedness and modularity are two commonly investigated structures in ecological networks and have been considered relevant to species feeding ecology24 and the stability of the system69. Finally, we measured the level of consumers’ diet niche overlap of the food webs (Horn’s index70, bounded 0–1), which essentially depends on the arrangement of trophic relationships (thus the structure of the webs), and could have strong ecological implications as niche partitioning has been recognised to be a key mechanism that drives species coexistence71,72. We selected these fundamental and holistic properties as they are potentially more relevant to the processes that may have shaped food webs across a landscape scale (e.g., community assembly), in comparison to some node- or link-centric properties. Also, addressing similar metrics as in the literature13,69 would facilitate potential cross-study comparison or validation.To first gain a glimpse of the structure of the blue and green food webs, we performed a principal component analysis (PCA; Fig. 3a) on the inferred food webs (n = 462 and 465 in green and blue, respectively) taking the four structural metrics (number of nodes, connectance, nestedness, and modularity) as the explaining variables of blue versus green system types. We then confirmed that system type, elevation, and land-use type were all important predictors of food-web metrics (whereas the residual temperature after accounting elevation effects was not) by conducting general linear model analyses, taking the former as interactive predictors while the latter response variables (Supplementary Tables 3, 4). To check how elevation influences food-web properties in blue and green systems separately, and how food-web properties depend on each other, we ran a series of piecewise structural equation modelling (SEM)73 analyses on inferred food webs (Fig. 3b, c) whose dominant land use can be defined (n = 421 and 430 in green and blue, respectively). This was also conducted on subsetted webs of each of the five major land-use types (Supplementary Figs. 1 and 2). The SEM relationships were derived from linear mixed model analyses with dominant land-use type as a random effect (assumption tests see Supplementary Figs. 12–17). The SEM structure of direct effects was set according to the literature13,69 and is illustrated in Fig. 3b. In short, this structure tests the dependencies from elevation (an environmental predictor) to food-web metrics (ecological responses). The further dependencies among food-web metrics themselves were assigned with the principle of pointing from relative lower-level properties to higher-level ones. That is, from number of nodes (purely determined by nodes) to connectance (determined by numbers of nodes and links), further to nestedness and modularity (holistic topologies, determined further by the arrangement of links), then to diet niche overlap (ecological functional outcome).Finally, to check and visualise the exact changing patterns of food webs, we applied generalised additive models (GAMs) to reveal the relationships between food-web metrics and the whole-ranged elevation (Figs. 4 and 5), as well as a particular comparison between food webs in forests and farmlands below 1500 m a.s.l. (Supplementary Fig. 5), as this elevation segment covered most of the sites belonged to these two land-use types. We also performed a series of linear models (LMs) and least-squared slope comparisons based on land-use-specific subsets of food webs (Figs. 4 and 5; Supplementary Figs. 3 and 4), to investigate whether food-web elevational patterns are different among land-use types (assumption tests see Supplementary Tables 5 and 6). In the GAMs analyses, specifically, we simulated two sets of randomised webs, i.e., “keep-group” and “fully”, as the null models to compare with the inferred ones74. Both randomisations generated ten independently simulated webs from each input inferred local food web, keeping the same number of nodes and connectance as of the latter. On the one hand, the keep-group randomisation shuffled trophic links from an input local web but only allowed them to realised fulfilling some pre-set within- and among-group relationships. That is, in green communities, birds can feed on all groups, grasshoppers on any groups but birds, while butterflies only on plants; in blue communities, fishes can feed on all groups, while invertebrates on themselves and the basal resources. These pre-set group-wide relationships captured the majority of realistic trophic interactions compiled in our metaweb. On the other hand, the fully randomised webs shuffled trophic links disregarding the biological identity of nodes. The GAMs of nestedness, modularity, and niche overlap illustrated the patterns of these randomised webs (Fig. 5). Comparing among the three types of webs, the patterns exhibited already by fully randomised webs should be those contributed by variations in web size and connectance, while the difference between keep-group and fully randomised webs by the focal-group composition of local communities, and the difference between inferred and keep-group randomised webs further by the realistic species-specific diets. In addition, we also applied the same GAMs and LMs approach to analyse node richness, as well as both realised and potential diet generality (vulnerability for plants) of each focal group (Supplementary Figs. 6–11). These analyses provided hints about the changes in community composition and species diet breadths along elevation and among land-use types, which helped explain the detected food-web responses in mechanistic ways.In addition, to check if our findings were shaped or strongly influenced by the potential inaccuracy of using the metaweb, we repeated the above PCA, SEM, and GAM analyses as a series of sensitivity analyses. We generated food webs based on our locally inferred ones (i.e., the observations) but with random 10% link removal. This procedure mimics the effect of potential intraspecific diet variation (mentioned earlier) so that some trophic interactions in the metaweb do not realise locally. Overall, these analyses with link removal showed that our conclusions are qualitatively and quantitatively highly robust, and only very minorly affected by the such potential inaccuracy of metawebs, which is also in accordance to other food-web studies (see e.g., Pearse & Altermatt 201575). All details and outcomes of these additional analyses are given in Supplementary discussion.All metric quantification and analyses were performed under R version 4.0.3 (R Core Team76). All applied packages and functions were described in Supplementary Methods, while the R scripts performing these tasks can be accessed at the online repository provided.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More