Ecology

Brondizio, E. S., Settele, J., Díaz, S. & Ngo, H. T. IPBES (2019): Global Assessment Report on Biodiversity and Ecosystem Services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES Secretariat, 2019).Urban, M. C. Accelerating extinction risk from climate change. Science 348, 571–573 (2015).CAS 
PubMed 

Google Scholar 
Warren, R., Price, J., Graham, E., Forstenhaeusler, N. & VanDerWal, J. The projected effect on insects, vertebrates, and plants of limiting global warming to 1.5 °C rather than 2 °C. Science 360, 791–795 (2018).CAS 
PubMed 

Google Scholar 
Schloss, C. A., Nuñez, T. A. & Lawler, J. J. Dispersal will limit ability of mammals to track climate change in the Western Hemisphere. Proc. Natl Acad. Sci. USA 109, 8606–8611 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Travis, J. M. J. Climate change and habitat destruction: a deadly anthropogenic cocktail. Proc. Biol. Sci. 270, 467–473 (2003).CAS 
PubMed 
PubMed Central 

Google Scholar 
Hill, J. K. et al. Impacts of landscape structure on butterfly range expansion. Ecol. Lett. 4, 313–321 (2001).
Google Scholar 
Guo, F., Lenoir, J. & Bonebrake, T. C. Land-use change interacts with climate to determine elevational species redistribution. Nat. Commun. 9, 1315 (2018).PubMed 
PubMed Central 

Google Scholar 
McLaughlin, J. F., Hellmann, J. J., Boggs, C. L. & Ehrlich, P. R. Climate change hastens population extinctions. Proc. Natl Acad. Sci. USA 99, 6070–6074 (2002).CAS 
PubMed 
PubMed Central 

Google Scholar 
Jetz, W., Wilcove, D. S. & Dobson, A. P. Projected impacts of climate and land-use change on the global diversity of birds. PLoS Biol. 5, e157 (2007).PubMed 
PubMed Central 

Google Scholar 
Mantyka-Pringle, C. S. et al. Climate change modifies risk of global biodiversity loss due to land-cover change. Biol. Conserv. 187, 103–111 (2015).
Google Scholar 
Conradie, S. R., Woodborne, S. M., Cunningham, S. J. & McKechnie, A. E. Chronic, sublethal effects of high temperatures will cause severe declines in southern African arid-zone birds during the 21st century. Proc. Natl Acad. Sci. USA 116, 14065–14070 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Forister, M. L. et al. Compounded effects of climate change and habitat alteration shift patterns of butterfly diversity. Proc. Natl Acad. Sci. USA 107, 2088–2092 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
Oliver, T. H. & Morecroft, M. D. Interactions between climate change and land use change on biodiversity: attribution problems, risks, and opportunities. Wiley Interdiscip. Rev. Clim. Change 5, 317–335 (2014).
Google Scholar 
MacLean, S. A. & Beissinger, S. R. Species’ traits as predictors of range shifts under contemporary climate change: a review and meta-analysis. Glob. Change Biol. 23, 4094–4105 (2017).
Google Scholar 
Pacifici, M. et al. Species’ traits influenced their response to recent climate change. Nat. Clim. Change 7, 205–208 (2017).
Google Scholar 
Root, T. Energy constraints on avian distributions and abundances. Ecology 69, 330–339 (1988).
Google Scholar 
Whitfield, M. C., Smit, B., McKechnie, A. E. & Wolf, B. O. Avian thermoregulation in the heat: scaling of heat tolerance and evaporative cooling capacity in three southern African arid-zone passerines. J. Exp. Biol. 218, 1705–1714 (2015).PubMed 

Google Scholar 
McKechnie, A. E. et al. Avian thermoregulation in the heat: evaporative cooling in five Australian passerines reveals within-order biogeographic variation in heat tolerance. J. Exp. Biol. 220, 2436–2444 (2017).PubMed 

Google Scholar 
Platts, P. J. et al. Habitat availability explains variation in climate-driven range shifts across multiple taxonomic groups. Sci. Rep. 9, 15039 (2019).PubMed 
PubMed Central 

Google Scholar 
Pearson, R. G. Climate change and the migration capacity of species. Trends Ecol. Evol. 21, 111–113 (2006).PubMed 

Google Scholar 
Partners in Flight. Avian Conservation Assessment Database Version 2021 (accessed 5 February 2021); http://pif.birdconservancy.org/ACADHill, M. J. & Guerschman, J. P. The MODIS global vegetation fractional cover product 2001–2018: characteristics of vegetation fractional cover in grasslands and savanna woodlands. Remote Sens. (Basel) 12, 406 (2020).
Google Scholar 
Wiebe, K. L. & Gerstmar, H. Influence of spring temperatures and individual traits on reproductive timing and success in a migratory woodpecker. Auk 127, 917–925 (2010).
Google Scholar 
Viana, D. S. & Chase, J. M. Ecological traits underlying interspecific variation in climate matching of birds. Glob. Ecol. Biogeogr. 31, 1021–1034 (2022).
Google Scholar 
Kellermann, V., Van Heerwaarden, B., Sgrò, C. M. & Hoffmann, A. A. Fundamental evolutionary limits in ecological traits drive Drosophila species distributions. Science 325, 1244–1246 (2009).CAS 
PubMed 

Google Scholar 
Devictor, V. et al. Differences in the climatic debts of birds and butterflies at a continental scale. Nat. Clim. Change 2, 121–124 (2012).
Google Scholar 
Mason, L. R. et al. Population responses of bird populations to climate change on two continents vary with species’ ecological traits but not with direction of change in climate suitability. Clim. Change 157, 337–354 (2019).
Google Scholar 
Coyle, J. R., Hurlbert, A. H. & White, E. P. Opposing mechanisms drive richness patterns of core and transient bird species. Am. Nat. 181, E83–E90 (2013).PubMed 

Google Scholar 
Valiela, I. & Martinetto, P. Changes in bird abundance in eastern North America: urban sprawl and global footprint? BioScience 57, 360–370 (2007).
Google Scholar 
Smith, S. J., Edmonds, J., Hartin, C. A., Mundra, A. & Calvin, K. Near-term acceleration in the rate of temperature change. Nat. Clim. Change 5, 333–336 (2015).
Google Scholar 
Winkler, K., Fuchs, R., Rounsevell, M. & Herold, M. Global land use changes are four times greater than previously estimated. Nat. Commun. 12, 2501 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Rosenberg, K. V. et al. Decline of the North American avifauna. Science 366, 120–124 (2019).CAS 
PubMed 

Google Scholar 
Currie, D. J. & Venne, S. Climate change is not a major driver of shifts in the geographical distributions of North American birds. Glob. Ecol. Biogeogr. 26, 333–346 (2017).
Google Scholar 
Socolar, J. B., Epanchin, P. N., Beissinger, S. R. & Tingley, M. W. Phenological shifts conserve thermal niches in North American birds and reshape expectations for climate-driven range shifts. Proc. Natl Acad. Sci. USA 114, 12976–12981 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
Barnagaud, J.-Y. et al. Relating habitat and climatic niches in birds. PLoS ONE 7, e32819 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Ponti, R., Arcones, A., Ferrer, X. & Vieites, D. R. Seasonal climatic niches diverge in migratory birds. Ibis 162, 318–330 (2020).
Google Scholar 
Zurell, D., Gallien, L., Graham, C. H. & Zimmermann, N. E. Do long-distance migratory birds track their niche through seasons? J. Biogeogr. 45, 1459–1468 (2018).
Google Scholar 
Stephens, P. A. et al. Consistent response of bird populations to climate change on two continents. Science 352, 84–87 (2016).CAS 
PubMed 

Google Scholar 
Ralston, J., DeLuca, W. V., Feldman, R. E. & King, D. I. Population trends influence species ability to track climate change. Glob. Change Biol. 23, 1390–1399 (2017).
Google Scholar 
Magurran, A. E. et al. Long-term datasets in biodiversity research and monitoring: assessing change in ecological communities through time. Trends Ecol. Evol. 25, 574–582 (2010).PubMed 

Google Scholar 
Jarzyna, M. A. & Jetz, W. A near half-century of temporal change in different facets of avian diversity. Glob. Change Biol. 23, 2999–3011 (2017).
Google Scholar 
van der Bolt, B., van Nes, E. H., Bathiany, S., Vollebregt, M. E. & Scheffer, M. Climate reddening increases the chance of critical transitions. Nat. Clim. Change 8, 478–484 (2018).
Google Scholar 
Bowler, D. E., Heldbjerg, H., Fox, A. D., O’Hara, R. B. & Böhning-Gaese, K. Disentangling the effects of multiple environmental drivers on population changes within communities. J. Anim. Ecol. 87, 1034–1045 (2018).PubMed 

Google Scholar 
Zurell, D., Graham, C. H., Gallien, L., Thuiller, W. & Zimmermann, N. E. Long-distance migratory birds threatened by multiple independent risks from global change. Nat. Clim. Change 8, 992–996 (2018).
Google Scholar 
Northrup, J. M., Rivers, J. W., Yang, Z. & Betts, M. G. Synergistic effects of climate and land-use change influence broad-scale avian population declines. Glob. Change Biol. 25, 1561–1575 (2019).
Google Scholar 
Guisan, A. et al. Predicting species distributions for conservation decisions. Ecol. Lett. 16, 1424–1435 (2013).PubMed 
PubMed Central 

Google Scholar 
Pardieck, K. L., Ziolkowski, D. J. Jr, Lutmerding, M., Aponte, V. & Hudson, M.-A. R. North American Breeding Bird Survey Dataset 1966–2018 Version 2018.0. (US Geological Survey, 2019); https://www.sciencebase.gov/catalog/item/5d65256ae4b09b198a26c1d7Harris, D. J., Taylor, S. D. & White, E. P. Forecasting biodiversity in breeding birds using best practices. PeerJ 6, e4278 (2018).PubMed 
PubMed Central 

Google Scholar 
Wickham, H., Francois, R., Henry, L. & Müller, K. dplyr: a Grammar of Data Manipulation. R package version 1.0.0 https://cran.r-project.org/web/packages/dplyr/index.html (2020).Wickham, H. & Henry, L. tidyr: Tidy Messy Data. R package version 1.1.0 https://cran.r-project.org/web/packages/tidyr/index.html (2020).Hijmans, R. J. raster: Geographic Data Analysis and Modeling. R package version 3.0-12 https://cran.r-project.org/web/packages/raster/index.html (2015).Bivand, R., Pebesma, E. J. & Gómez-Rubio, V. Applied Spatial Data Analysis with R (Springer, 2013).Hijmans, R. J. geosphere: Spherical Trigonometry. R package version 1.5–10 https://cran.r-project.org/web/packages/geosphere/index.html (2019).Hart, E. M. & Bell, K. prism. R package version 0.0.6 https://github.com/ropensci/prism (2015).Senyondo, H. et al. rdataretriever: R interface to the data retriever. J. Open Source Softw. 6, 2800 (2021).
Google Scholar 
Morris, B. D. & White, E. P. The EcoData retriever: improving access to existing ecological data. PLoS ONE 8, e65848 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Senyondo, H. et al. Retriever: data retrieval tool. J. Open Source Softw. 2, 451 (2017).
Google Scholar 
Hurlbert, A. H. & White, E. P. Disparity between range map- and survey-based analyses of species richness: patterns, processes and implications. Ecol. Lett. 8, 319–327 (2005).
Google Scholar 
Harris, D. J. Generating realistic assemblages with a joint species distribution model. Methods Ecol. Evol. 6, 465–473 (2015).
Google Scholar 
Sheard, C. et al. Ecological drivers of global gradients in avian dispersal inferred from wing morphology. Nat. Commun. 11, 2463 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Eyres, A., Böhning-Gaese, K. & Fritz, S. A. Quantification of climatic niches in birds: adding the temporal dimension. J. Avian Biol. 48, 1517–1531 (2017).
Google Scholar 
Martin, A. E. & Fahrig, L. Habitat specialist birds disperse farther and are more migratory than habitat generalist birds. Ecology 99, 2058–2066 (2018).PubMed 

Google Scholar 
Sauer, J. R. & Link, W. A. Analysis of the North American Breeding Bird Survey using hierarchical models. Auk 128, 87–98 (2011).
Google Scholar 
García Molinos, J., Schoeman, D. S., Brown, C. J. & Burrows, M. T. VoCC: an R package for calculating the velocity of climate change and related climatic metrics. Methods Ecol. Evol. 10, 2195–2202 (2019).
Google Scholar 
Krenek, S., Berendonk, T. U. & Petzoldt, T. Thermal performance curves of Paramecium caudatum: a model selection approach. Eur. J. Protistol. 47, 124–137 (2011).PubMed 

Google Scholar 
Bahn, V. & McGill, B. J. Can niche-based distribution models outperform spatial interpolation? Glob. Ecol. Biogeogr. 16, 733–742 (2007).
Google Scholar 
Dobson, L. L., La Sorte, F. A., Manne, L. L. & Hawkins, B. A. The diversity and abundance of North American bird assemblages fail to track changing productivity. Ecology 96, 1105–1114 (2015).PubMed 

Google Scholar 
Roberts, D. R. et al. Cross-validation strategies for data with temporal, spatial, hierarchical, or phylogenetic structure. Ecography 40, 913–929 (2017).
Google Scholar 
Tikhonov, G. et al. Joint species distribution modelling with the R-package HMSC. Methods Ecol. Evol. 11, 442–447 (2020).PubMed 
PubMed Central 

Google Scholar 
Greenwell, B., Boehmke, B., Cunningham, J. & GBM Developers. gbm: Generalized boosted regression models. R package version 2.1.5 https://cran.r-project.org/web/packages/gbm/index.html (2019).Wood, S. N. Generalized Additive Models: an Introduction with R (CRC Press/Taylor & Francis Group, 2017).Jetz, W., Thomas, G. H., Joy, J. B., Hartmann, K. & Mooers, A. O. The global diversity of birds in space and time. Nature 491, 444–448 (2012).CAS 
PubMed 

Google Scholar 
Revell, L. J. phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol. 3, 217–223 (2012).
Google Scholar 
Paradis, E. & Schliep, K. ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 35, 526–528 (2019).CAS 
PubMed 

Google Scholar 
Bürkner, P.-C. brms: an R package for Bayesian multilevel models using Stan. J. Stat. Softw. 80, 1–28 (2017).
Google Scholar 
Stan Development Team. Stan Modeling Language Users Guide and Reference Manual (2020); https://mc-stan.org/users/documentation/ More

Wu, A. Social buffering of stress – Physiological and ethological perspectives. Appl. Anim. Behav. Sci. 239, 105325 (2021).
Google Scholar 
Hennessy, M. B., Kaiser, S. & Sachser, N. Social buffering of the stress response: diversity, mechanisms, and functions. Front. Neuroendocrinol. 30, 470–482 (2009).CAS 
PubMed 

Google Scholar 
Young, C., Majolo, B., Heistermann, M., Schülke, O. & Ostner, J. Responses to social and environmental stress are attenuated by strong male bonds in wild macaques. Proc. Natl Acad. Sci. USA 111, 18195–18200 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
Stanton, M. E., Patterson, J. M. & Levine, S. Social influences on conditioned cortisol secretion in the squirrel monkey. Psychoneuroendocrinology 10, 125–134 (1985).CAS 
PubMed 

Google Scholar 
Caldji, C., Diorio, J. & Meaney, M. J. Variations in maternal care in infancy regulate the development of stress reactivity. Biol. Psychiatry 48, 1164–1174 (2000).CAS 
PubMed 

Google Scholar 
Novak, M. A., Hamel, A. F., Kelly, B. J., Dettmer, A. M. & Meyer, J. S. Stress, the HPA axis, and nonhuman primate well-being: a review. Appl. Anim. Behav. Sci. 143, 135–149 (2013).PubMed 
PubMed Central 

Google Scholar 
Sapolsky, R. M., Romero, L. M. & Munck, A. U. How do glucocorticoids influence stress responses? Integrating permissive, suppressive, stimulatory, and preparative actions. Endocr. Rev. 21, 55–89 (2000).CAS 
PubMed 

Google Scholar 
Liu, D. et al. Maternal Care, hippocampal glucocorticoid receptors, and hypothalamic-pituitary-adrenal responses to stress. Sci. Ment. Heal. Stress Brain 9, 75–78 (1997).
Google Scholar 
Gjerstad, J. K., Lightman, S. L. & Spiga, F. Role of glucocorticoid negative feedback in the regulation of HPA axis pulsatility. Stress 21, 403–416 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
Spiga, F., Walker, J. J., Terry, J. R. & Lightman, S. L. HPA axis-rhythms. Compr. Physiol. 4, 1273–1298 (2014).PubMed 

Google Scholar 
Sapolsky, R. M. Why Zebras Don’t Get Ulcers (Henry Holt and Company, LLC, 2004).Campos, F. A. et al. Glucocorticoid exposure predicts survival in female baboons. Sci. Adv. 7, 1–10 (2021).
Google Scholar 
Banerjee, S. B., Arterbery, A. S., Fergus, D. J. & Adkins-Regan, E. Deprivation of maternal care has long-lasting consequences for the hypothalamic-pituitary-adrenal axis of zebra finches. Proc. R. Soc. B Biol. Sci. 279, 759–766 (2012).
Google Scholar 
Hennessy, M. B., Nigh, C. K., Sims, M. L. & Long, S. J. Plasma cortisol and vocalization responses of postweaning age guinea pigs to maternal and sibling separation: evidence for filial attachment after weaning. Dev. Psychobiol. 28, 103–115 (1995).CAS 
PubMed 

Google Scholar 
Hennessy, M. B., O’Leary, S. K., Hawke, J. L. & Wilson, S. E. Social influences on cortisol and behavioral responses of preweaning, periadolescent, and adult guinea pigs. Physiol. Behav. 76, 305–314 (2002).CAS 
PubMed 

Google Scholar 
Wiener, S. G., Johnson, D. F. & Levine, S. Influence of postnatal rearing conditions on the response of squirrel monkey infants to brief perturbations in mother-infant relationships. Physiol. Behav. 39, 21–26 (1987).CAS 
PubMed 

Google Scholar 
Girard-Buttoz, C. et al. Early maternal loss leads to short-but not long-term effects on diurnal cortisol slopes in wild chimpanzees. Elife 10, e64134 (2021).Rosenbaum, S. et al. Social bonds do not mediate the relationship between early adversity and adult glucocorticoids in wild baboons. Proc. Natl Acad. Sci. USA 117, 20052–20062 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Moss, C. Elephant Memories: Thirteen Years in the Life of an Elephant Family (Univ. Chicago Press, 1988).Douglas-Hamilton, I., Bhalla, S., Wittemyer, G. & Vollrath, F. Behavioural reactions of elephants towards a dying and deceased matriarch. Appl. Anim. Behav. Sci. 100, 87–102 (2006).
Google Scholar 
Shoshani, J., Kupsky, W. J. & Marchant, G. H. Elephant brain. Part I: gross morphology, functions, comparative anatomy, and evolution. Brain Res. Bull. 70, 124–157 (2006).PubMed 

Google Scholar 
Goldenberg, S. Z. & Wittemyer, G. Orphaned female elephant social bonds reflect lack of access to mature adults. Sci. Rep. 7, 14408 (2017).PubMed 
PubMed Central 

Google Scholar 
Goldenberg, S. Z. & Wittemyer, G. Orphaning and natal group dispersal are associated with social costs in female elephants. Anim. Behav. 143, 1–8 (2018).
Google Scholar 
Lee, P. C. Allomothering among African elephants. Anim. Behav. 35, 278–291 (1987).
Google Scholar 
Parker, J. M. et al. Poaching of African elephants indirectly decreases population growth through lowered orphan survival. Curr. Biol. 31, 4156–4162.e5 (2021).Wittemyer, G. et al. Where sociality and relatedness diverge: the genetic basis for hierarchical social organization in African elephants. Proc. R. Soc. B Biol. Sci. 276, 3513–3521 (2009).
Google Scholar 
Goldenberg, S. Z., Douglas-Hamilton, I. & Wittemyer, G. Vertical transmission of social roles drives resilience to poaching in elephant metworks. Curr. Biol. 26, 75–79 (2016).CAS 
PubMed 

Google Scholar 
Gobush, K. S., Mutayoba, B. M. & Wasser, S. K. Long-term impacts of poaching on relatedness, stress physiology, and reproductive output of adult female African elephants. Conserv. Biol. 22, 1590–1599 (2008).CAS 
PubMed 

Google Scholar 
Gobush, K. S. et al. Loxodonta africana (African Savanna Elephant). Loxodonta africana: the IUCN red list of threatened species 2021 e.T181008073A181022663 https://doi.org/10.2305/IUCN.UK.2021-1.RLTS.T181008073A181022663.en (2021).Wittemyer, G. et al. Illegal killing for ivory drives global decline in African elephants. Proc. Natl Acad. Sci. USA 111, 13117–13121 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
Wittemyer, G., Daballen, D. & Douglas-Hamilton, I. Comparative Demography of an At-Risk African Elephant Population. PLoS ONE 8, e53726 (2013).McCormick, S. D. & Romero, L. M. Conservation endocrinology. Bioscience 67, 429–442 (2017).
Google Scholar 
Wittemyer, G. The elephant population of Samburu and Buffalo Springs National Reserves, Kenya. Afr. J. Ecol. 39, 357–369 (2001).
Google Scholar 
Cockrem, J. F. Individual variation in glucocorticoid stress responses in animals. Gen. Comp. Endocrinol. 181, 45–58 (2013).CAS 
PubMed 

Google Scholar 
Taff, C. C., Schoenle, L. A. & Vitousek, M. N. The repeatability of glucocorticoids: a review and meta-analysis. Gen. Comp. Endocrinol. 260, 136–145 (2018).CAS 
PubMed 

Google Scholar 
Hooten, M. B. & Hobbs, N. T. A guide to Bayesian model selection for ecologists. Ecol. Monogr. 85, 3–28 (2015).
Google Scholar 
Wittemyer, G. & Getz, W. M. Hierarchical dominance structure and social organization in African elephants, Loxodonta africana. Anim. Behav. 73, 671–681 (2007).
Google Scholar 
Heim, C., Ehlert, U. & Hellhammer, D. H. The potential role of hypocortisolism in the pathophysiology of stress-related bodily disorders. Psychoneuroendocrinology 25, 1–35 (2000).CAS 
PubMed 

Google Scholar 
Dickens, M. J. & Romero, L. M. A consensus endocrine profile for chronically stressed wild animals does not exist. Gen. Comp. Endocrinol. 191, 177–189 (2013).CAS 
PubMed 

Google Scholar 
Ma, D., Serbin, L. A. & Stack, D. M. How children’s anxiety symptoms impact the functioning of the hypothalamus–pituitary–adrenal axis over time: a cross-lagged panel approach using hierarchical linear modeling. Dev. Psychopathol. 31, 1–15 (2018).
Google Scholar 
Blas, J., Bortolotti, G. R., Tella, J. L., Baos, R. & Marchant, T. A. Stress response during development predicts fitness in a wild, long lived vertebrate. Proc. Natl Acad. Sci. USA 104, 8880–8884 (2007).CAS 
PubMed 
PubMed Central 

Google Scholar 
Boonstra, R. Reality as the leading cause of stress: rethinking the impact of chronic stress in nature. Funct. Ecol. 27, 11–23 (2013).
Google Scholar 
Gunnar, M. R. & Vazquez, D. M. Low cortisol and a flattening of expected daytime rhythm: Potential indices of risk in human development. Dev. Psychopathol. 13, 515–538 (2001).CAS 
PubMed 

Google Scholar 
Perry, R. E. et al. Corticosterone administration targeting a hypo-reactive HPA axis rescues a socially-avoidant phenotype in scarcity-adversity reared rats. Dev. Cogn. Neurosci. 40, 100716 (2019).Fries, E., Hesse, J., Hellhammer, J. & Hellhammer, D. H. A new view on hypocortisolism. Psychoneuroendocrinology 30, 1010–1016 (2005).CAS 
PubMed 

Google Scholar 
Dorsey, C., Dennis, P., Guagnano, G., Wood, T. & Brown, J. L. Decreased baseline fecal glucocorticoid concentrations associated with skin and oral lesions in black rhinoceros (Diceros bicornis). J. Zoo. Wildl. Med. 41, 616–625 (2010).PubMed 

Google Scholar 
Pawluski, J. et al. Low plasma cortisol and fecal cortisol metabolite measures as indicators of compromised welfare in domestic horses (Equus caballus). PLoS ONE 12, 1–18 (2017).
Google Scholar 
Feng, X. et al. Maternal separation produces lasting changes in cortisol and behavior in rhesus monkeys. Proc. Natl Acad. Sci. USA 108, 14312–14317 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
González Ramírez, C. et al. The NR3C1 gene expression is a potential surrogate biomarker for risk and diagnosis of posttraumatic stress disorder. Psychiatry Res. 284, 112797 (2020).PubMed 

Google Scholar 
Cluver, L., Fincham, D. S. & Seedat, S. Posttraumatic stress in AIDS-orphaned children exposed to high levels of trauma: the protective role of perceived social support. J. Trauma. Stress 22, 106–112 (2009).PubMed 

Google Scholar 
Bastille-Rousseau, G. et al. Landscape-scale habitat response of African elephants shows strong selection for foraging opportunities in a human dominated ecosystem. Ecography 43, 149–160 (2020).
Google Scholar 
Foley, C. A. H., Papageorge, S. & Wasser, S. K. Noninvasive stress and reproductive measures of social and ecological pressures in free-ranging African elephants. Conserv. Biol. 15, 1134–1142 (2001).
Google Scholar 
Wittemyer, G., Getz, W. M., Vollrath, F. & Douglas-Hamilton, I. Social dominance, seasonal movements, and spatial segregation in African elephants: a contribution to conservation behavior. Behav. Ecol. Sociobiol. 61, 1919–1931 (2007).
Google Scholar 
Wittemyer, G., Daballen, D. & Douglas‐Hamilton, I. Differential influence of human impacts on age‐specific demography underpins trends in an African elephant population. Ecosphere 12, e03720 (2021).Brown, J. L. et al. Individual and environmental risk factors associated with fecal glucocorticoid metabolite concentrations in zoo-housed Asian and African elephants. PLoS ONE 14, 1–18 (2019).
Google Scholar 
Goldenberg, S. Z. et al. Increasing conservation translocation success by building social functionality in released populations. Glob. Ecol. Conserv. 18, e00604 (2019).Dantzer, B., Fletcher, Q. E., Boonstra, R. & Sheriff, M. J. Measures of physiological stress: a transparent or opaque window into the status, management and conservation of species? Conserv. Physiol. 2, 1–18 (2014).
Google Scholar 
Kaisin, O., Fuzessy, L., Poncin, P., Brotcorne, F. & Culot, L. A meta-analysis of anthropogenic impacts on physiological stress in wild primates. Conserv. Biol. 0, 1–14 (2020).CAS 

Google Scholar 
Ganswindt, A., Rasmussen, H. B., Heistermann, M. & Hodges, J. K. The sexually active states of free-ranging male African elephants (Loxodonta africana): defining musth and non-musth using endocrinology, physical signals, and behavior. Horm. Behav. 47, 83–91 (2005).CAS 
PubMed 

Google Scholar 
Santymire, R. M. et al. Using ACTH challenges to validate techniques for adrenocortical activity analysis in various African wildlife species. Int. J. Anim. Vet. Adv. 4, 99–108 (2012).CAS 

Google Scholar 
Watson, R. et al. Development of a versatile enzyme immunoassay for non-invasive assessment of glucocorticoid metabolites in a diversity of taxonomic species. Gen. Comp. Endocrinol. 186, 16–24 (2013).CAS 
PubMed 

Google Scholar 
Oduor, S. et al. Differing physiological and behavioral responses to anthropogenic factors between resident and non-resident African elephants at Mpala Ranch, Laikipia County, Kenya. PeerJ 8, e10010 (2020).Brown, J. L., Kersey, D. C., Freeman, E. W. & Wagener, T. Assessment of diurnal urinary cortisol excretion in Asian and African elephants using different endocrine methods. Zoo. Biol. 29, 274–283 (2010).PubMed 

Google Scholar 
Justice, C. O. et al. The moderate resolution imaging spectroradiometer (MODIS): land remote sensing for global change research. IEEE Trans. Geosci. Remote Sens. 36, 1228–1249 (1998).
Google Scholar 
Lafferty, D. J. R., Zimova, M., Clontz, L., Hackländer, K. & Mills, L. S. Noninvasive measures of physiological stress are confounded by exposure. Sci. Rep. 9, 1–6 (2019).
Google Scholar 
O’Dwyer, K., Dargent, F., Forbes, M. R. & Koprivnikar, J. Parasite infection leads to widespread glucocorticoid hormone increases in vertebrate hosts: a meta-analysis. J. Anim. Ecol. 89, 519–529 (2020).PubMed 

Google Scholar 
Parker, J. M., Goldenberg, S. Z., Letitiya, D. & Wittemyer, G. Strongylid infection varies with age, sex, movement and social factors in wild African elephants. Parasitology 147, 348–359 (2020).PubMed 

Google Scholar 
Gibbons, L., Jacobs, D. E., Fox, M. T. & Hansen, J. The RVC/FAO guide to veterinary diagnostic parasitology. McMaster egg-counting technique. http://www.rvc.ac.uk/review/Parasitology/EggCount/Purpose.htm (2004)R Core Team. A language and environment for statistical computing. https://www.r-project.org/. (2020).Rstudio Team. RStudio: integrated development for R. http://www.rstudio.com/ (2020).Plummer, M. rjags: Bayesian graphical models using MCMC. https://cran.r-project.org/package=rjags (2019).Brooks, S. P. & Gelman, A. General methods for monitoring convergence of iterative simulations. J. Comput. Graph. Stat. 7, 434–455 (1998).
Google Scholar 
Gelman, A. & Rubin, D. B. Inference from iterative simulation using multiple sequences. Stat. Sci. 7, 457–511 (1992).
Google Scholar 
Wickham, H. ggplot2: elegant graphics for data analysis. https://ggplot2.tidyverse.org (2016).Youngflesh, C. MCMCvis: tools to visualize, manipulate, and summarize MCMC output. J. Open Source Softw. 3, 640 (2018).
Google Scholar 
Parker, J. M. The Physiological Condition of Orphaned African Elephants (Loxodonta africana). Doctoral dissertation, Colorado State University. (2021). More

Cremers, G. Presence of 10 models of plant architecture in Euphorbes-Malgaches. Comptes Rendus Hebd. des. Seances de. L Academie des. Sci. Ser. D. 281, 1575–1578 (1975).
Google Scholar 
Balduzzi, M. et al. Reshaping plant biology: qualitative and quantitative descriptors for plant morphology. Front. Plant Sci. 8, 117 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Albert, C. H. et al. A multi-trait approach reveals the structure and the relative importance of intra- vs. interspecific variability in plant traits. Funct. Ecol. 24, 1192–1201 (2010).Article 

Google Scholar 
Farnsworth, K. D. & Niklas, K. J. Theories of optimization, form and function in branching architecture in plants. Funct. Ecol. 9, 355–363 (1995).Article 

Google Scholar 
Enquist, B. J. et al. in Advances in Ecological Research (eds Pawar, S.et al.), 249–318 (Academic Press, 2015).Niklas, K. J. & Spatz, H. C. Allometric theory and the mechanical stability of large trees: proof and conjecture. Am. J. Bot. 93, 824–828 (2006).PubMed 
Article 

Google Scholar 
Price, C. A. et al. The metabolic theory of ecology: prospects and challenges for plant biology. N. Phytol. 188, 696–710 (2010).Article 

Google Scholar 
Martone, P. T. et al. Mechanics without muscle: biomechanical inspiration from the plant world. Integr. Comp. Biol. 50, 888–907 (2010).PubMed 
Article 

Google Scholar 
West, G. B. & Brown, J. H. The origin of allometric scaling laws in biology from genomes to ecosystems: towards a quantitative unifying theory of biological structure and organization. J. Exp. Biol. 208, 1575–1592 (2005).PubMed 
Article 

Google Scholar 
Enquist, B. J. Universal scaling in tree and vascular plant allometry: toward a general quantitative theory linking plant form and function from cells to ecosystems. Tree Physiol. 22, 1045–1064 (2002).PubMed 
Article 

Google Scholar 
Anfodillo, T. et al. An allometry-based approach for understanding forest structure, predicting tree-size distribution and assessing the degree of disturbance. Proc. R. Soc. Lond. B Biol. Sci. 280, 20122375 (2013).
Google Scholar 
Duncanson, L. I., Dubayah, R. O. & Enquist, B. J. Assessing the general patterns of forest structure: quantifying tree and forest allometric scaling relationships in the United States. Glob. Ecol. Biogeogr. 24, 1465–1475 (2015).Article 

Google Scholar 
West, G. B., Brown, J. H. & Enquist, B. J. The fourth dimension of life: Fractal geometry and allometric scaling of organisms. Science 284, 1677–1679 (1999).CAS 
PubMed 
Article 

Google Scholar 
Winter, C. L. & Tartakovsky, D. M. Theoretical foundation for conductivity scaling. Geophys. Res. Lett. 28, 4367–4369 (2001).Article 

Google Scholar 
Reich, P. B. et al. Universal scaling of respiratory metabolism, size and nitrogen in plants. Nature 439, 457–461 (2006).CAS 
PubMed 
Article 

Google Scholar 
Choi, S. et al. Application of the metabolic scaling theory and water–energy balance equation to model large‐scale patterns of maximum forest canopy height. Glob. Ecol. Biogeogr. 25, 1428–1442 (2016).Article 

Google Scholar 
Osler, G. H. R., West, P. W. & Downes, G. M. Effects of bending stress on taper and growth of stems of young Eucalyptus regnans trees. Trees 10, 239–246 (1996).
Google Scholar 
Berthier, S. et al. Irregular heartwood formation in maritime pine (Pinus pinaster Ait): consequences for biomechanical and hydraulic tree functioning. Ann. Bot. 87, 19–25 (2001).Article 

Google Scholar 
Fournier, M. et al. Integrative biomechanics for tree ecology: beyond wood density and strength. J. Exp. Bot. 64, 4793–4815 (2013).CAS 
PubMed 
Article 

Google Scholar 
Sone, K., Noguchi, K. & Terashima, I. Dependency of branch diameter growth in young Acer trees on light availability and shoot elongation. Tree Physiol. 25, 39–48 (2005).PubMed 
Article 

Google Scholar 
Anten, N. P. & Schieving, F. The role of wood mass density and mechanical constraints in the economy of tree architecture. Am. Nat. 175, 250–260 (2010).PubMed 
Article 

Google Scholar 
Jelonek, T. et al. The biomechanical formation of trees (Prace Naukowe, Doniesienia, Komunikaty, 2019).Muller‐Landau, H. C. et al. Testing metabolic ecology theory for allometric scaling of tree size, growth and mortality in tropical forests. Ecol. Lett. 9, 575–588 (2006).PubMed 
Article 

Google Scholar 
McMahon, T. A. & Kronauer, R. E. Tree structures: deducing the principle of mechanical design. J. Theor. Biol. 59, 443–466 (1976).CAS 
PubMed 
Article 

Google Scholar 
Alméras, T. & Fournier, M. Biomechanical design and long-term stability of trees: morphological and wood traits involved in the balance between weight increase and the gravitropic reaction. J. Theor. Biol. 256, 370–381 (2009).PubMed 
Article 

Google Scholar 
West, G. B., Brown, J. H. & Enquist, B. J. A general model for the origin of allometric scaling laws in biology. Science 276, 122–126 (1997).CAS 
PubMed 
Article 

Google Scholar 
Mäkelä, A. & Valentine, H. T. Crown ratio influences allometric scaling in trees. Ecol 87, 2967–2972 (2006).Article 

Google Scholar 
Duursma, R. A. et al. Self‐shading affects allometric scaling in trees. Funct. Ecol. 24, 723–730 (2010).Article 

Google Scholar 
Pretzsch, H. & Dieler, J. Evidence of variant intra-and interspecific scaling of tree crown structure and relevance for allometric theory. Oecologia 169, 637–649 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
Lin, Y. et al. Plant interactions alter the predictions of metabolic scaling theory. PloS One 8, e57612 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Cheng, D. et al. Scaling relationship between tree respiration rates and biomass. Biol. Lett. 6, 715–717 (2010).PubMed 
PubMed Central 
Article 

Google Scholar 
Ogawa, K. Scaling relations based on the geometric and metabolic theories in woody plant species: A review. Perspect. Plant Ecol. Evol. Syst. 40, 125480 (2019).Article 

Google Scholar 
Risto, S. et al. Functional–structural plant models: a growing paradigm for plant studies. Ann. Bot. 114, 599–603 (2014).Article 

Google Scholar 
Jackson, T. et al. Finite element analysis of trees in the wind based on terrestrial laser scanning data. Agric. Meteorol. 265, 137–144 (2019).Article 

Google Scholar 
Disney, M. Terrestrial LiDAR: a three‐dimensional revolution in how we look at trees. N. Phytol. 222, 1736–1741 (2019).Article 

Google Scholar 
Malhi, Y. et al. New perspectives on the ecology of tree structure and tree communities through terrestrial laser scanning. Interface Focus 8, 20170052 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Bayer, D., Seifert, S. & Pretzsch, H. Structural crown properties of Norway spruce (Picea abies [L.] Karst.) and European beech (Fagus sylvatica [L.]) in mixed versus pure stands revealed by terrestrial laser scanning. Trees 27, 1035–1047 (2013).Article 

Google Scholar 
Lin, Y. & Herold, M. Tree species classification based on explicit tree structure feature parameters derived from static terrestrial laser scanning data. Agric. Meteorol. 216, 105–114 (2016).Article 

Google Scholar 
Tanago, J. G. et al. Estimation of above‐ground biomass of large tropical trees with terrestrial LiDAR. Methods Ecol. Evol. 9, 223–234 (2018).Article 

Google Scholar 
Takoudjou, S. M. et al. Using terrestrial laser scanning data to estimate large tropical trees biomass and calibrate allometric models: A comparison with traditional destructive approach. Methods Ecol. Evol. 9, 905–916 (2018).Article 

Google Scholar 
Dassot, M., Fournier, M. & Deleuze, C. Assessing the scaling of the tree branch diameters frequency distribution with terrestrial laser scanning: methodological framework and issues. Ann. Sci. 76, 66 (2019).Article 

Google Scholar 
Klockow, P. A. et al. Allometry and structural volume change of standing dead southern pine trees using non-destructive terrestrial LiDAR. Remote Sens. Environ. 241, 111729 (2020).Article 

Google Scholar 
Stovall, A. E., Anderson-Teixeira, K. J. & Shugart, H. H. Assessing terrestrial laser scanning for developing non-destructive biomass allometry. Ecol. Manag. 427, 217–229 (2018).Article 

Google Scholar 
Dai, J. et al. Drought-modulated allometric patterns of trees in semi-arid forests. Commun. Biol. 3, 1–8 (2020).Article 

Google Scholar 
Ogawa, K., Hagihara, A. & Hozumi, K. Growth analysis of a seedling community of Chamaecyparis obtusa. VI. Estimation of individual gross primary production by the summation method. In Transactions of the 30th Meeting of Chubu Branch of Japanese Forestry Society, 179–181 (Honda Kiyoshi, 1985).Yokota, T. & Hagihara, A. Dependence of the aboveground CO2 exchange rate on tree size in field-grown hinoki cypress (Chamaecyparis obtusa). J. Plant Res. 109, 177–184 (1996).Article 

Google Scholar 
Enquist, B. J. et al. Biological scaling: does the exception prove the rule? Nature 445, E9–E10 (2007).CAS 
PubMed 
Article 

Google Scholar 
Lau, A. et al. Estimating architecture-based metabolic scaling exponents of tropical trees using terrestrial LiDAR and 3D modelling. Ecol. Manag. 439, 132–145 (2019).Article 

Google Scholar 
Li, Y. et al. Retrieval of tree branch architecture attributes from terrestrial laser scan data using a Laplacian algorithm. Agric. Meteorol. 284, 107874 (2020).Article 

Google Scholar 
Noyer, E. et al. Biomechanical control of beech pole verticality (Fagus sylvatica) before and after thinning: theoretical modelling and ground‐truth data using terrestrial LiDAR. Am. J. Bot. 106, 187–198 (2019).PubMed 
Article 

Google Scholar 
Jackson, T. et al. A new architectural perspective on wind damage in a natural forest. Front. Glob. Chang. 1, 13 (2019).Article 

Google Scholar 
Jackson, T. Strain Measurements on 21 Trees in Wytham Woods, UK. NERC Environmental Information Data Centre. https://doi.org/10.5285/533d87d3-48c1-4c6e-9f2f-fda273ab45bc (2018).Kozłowski, J. & Konarzewski, M. Is West, Brown and Enquist’s model of allometric scaling mathematically correct and biologically relevant? Funct. Ecol. 18, 283–289 (2004).Article 

Google Scholar 
Kleiber, M. Body size and metabolic rate. Physiol. Rev. 27, 511–541 (1947).CAS 
PubMed 
Article 

Google Scholar 
Hay, M. J. M. et al. Branching responses of a plagiotropic clonal herb to localised incidence of light simulating that reflected from vegetation. Oecologia 127, 185–190 (2001).CAS 
PubMed 
Article 

Google Scholar 
Cordero, R. A., Fetcher, N. & Voltzow, J. Effects of wind on the allometry of two species of plants in an elfin cloud forest. Biotropica 39, 177–185 (2010).Article 

Google Scholar 
Anfodillo, T. et al. Allometric trajectories and “stress”: a quantitative approach. Front. Plant Sci. 7, 1681 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Louarn, G. & Song, Y. Two decades of functional-structural plant modelling: now addressing fundamental questions in systems biology and predictive ecology. Ann. Bot. 126, 501–509 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Poorter, H. & Sack, L. Pitfalls and possibilities in the analysis of biomass allocation patterns in plants. Front. Plant Sci. 3, 259 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
Thomas, S. C. Reproductive allometry in Malaysian rain forest trees: biomechanics versus optimal allocation. Evol. Ecol. 10, 517–530 (1996).Article 

Google Scholar 
Kempes, C. P. et al. Predicting maximum tree heights and other traits from allometric scaling and resource limitations. PLoS One 6, e20551 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Blanchard, E. et al. Contrasted allometries between stem diameter, crown area, and tree height in five tropical biogeographic areas. Trees 30, 1953–1968 (2016).Article 

Google Scholar 
Swetnam, T. L., O’Connor, C. D. & Lynch, A. M. Tree morphologic plasticity explains deviation from metabolic scaling theory in semi-arid conifer forests, southwestern USA. PLoS One 11, e0157582 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Loehle, C. Biomechanical constraints on tree architecture. Trees 30, 2061–2070 (2016).Article 

Google Scholar 
Guillon, T., Dumont, Y. & Fourcaud, T. Numerical methods for the biomechanics of growing trees. Comput. Math. Appl. 64, 289–309 (2012).Article 

Google Scholar 
Olson, M. E., Rosell, J. A., Muñoz, S. Z. & Castorena, M. Carbon limitation, stem growth rate and the biomechanical cause of Corner’s rules. Ann. Bot. 122, 583–592 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
West, G. B., Enquist, B. J. & Brown, J. H. A general quantitative theory of forest structure and dynamics. Proc. Natl Acad. Sci. USA 106, 7040–7045 (2009).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar  More

Mouillot, D. et al. Rare species support vulnerable functions in high-diversity ecosystems. PLoS Biol. 11, e1001569 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Leitão, R. P. et al. Rare species contribute disproportionately to the functional structure of species assemblages. Proc. R Soc. B 283, 20160084 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Enquist, B. J. et al. The commonness of rarity: Global and future distribution of rarity across land plants. Sci. Adv. 5, eaaz0414 (2019).Bevill, R. L. & Louda, S. M. Comparisons of related rare and common species in the study of plant rarity. Conserv. Biol. 13, 493–498 (1999).Article 

Google Scholar 
Murray, B. R., Thrall, P. H., Gill, A. M. & Nicotra, A. B. How plant life-history and ecological traits relate to species rarity and commonness at varying spatial scales. Austral Ecol. 27, 291–310 (2002).Article 

Google Scholar 
Gaston, K. J. Common ecology. Bioscience 61, 354–362 (2011).Article 

Google Scholar 
Kraft, N. J. et al. Community assembly, coexistence and the environmental filtering metaphor. Funct. Ecol. 29, 592–599 (2015).Article 

Google Scholar 
Gaston, K. J. What is rarity? in Rarity 1–21 (Springer, 1994).Rabinowitz, D. Seven forms of rarity. in The biological aspects of rare plant conservation (ed. Synge, H.) 205–217 (John Wiley and Sons: Chichester, UK, 1981).Sykes, L., Santini, L., Etard, A. & Newbold, T. Effects of rarity form on species’ responses to land use. Conserv. Biol. 34, 688–696 (2019).PubMed 
Article 

Google Scholar 
Patykowski, J. et al. The effect of prescribed burning on plant rarity in a temperate forest. Ecol. Evol. 8, 1714–1725 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Ames, G. M., Wall, W. A., Hohmann, M. G. & Wright, J. P. Trait space of rare plants in a fire-dependent ecosystem. Conserv. Biol. 31, 903–911 (2017).PubMed 
Article 

Google Scholar 
Foster, C. N. et al. Effects of fire regime on plant species richness and composition differ among forest, woodland and heath vegetation. Appl. Veg. Sci. 21, 132–143 (2018).Article 

Google Scholar 
Fernández-García, V. et al. Fire regimes shape diversity and traits of vegetation under different climatic conditions. Sci. Total Environ. 716, 137137 (2020).ADS 
PubMed 
Article 
CAS 

Google Scholar 
Bassett, M., Leonard, S. W. J., Chia, E. K., Clarke, M. F. & Bennett, A. F. Interacting effects of fire severity, time since fire and topography on vegetation structure after wildfire. For. Ecol. Manag. 396, 26–34 (2017).Article 

Google Scholar 
Miller, B. P., Symons, D. R. & Barrett, M. D. Persistence of rare species depends on rare events: Demography, fire response and phenology of two plant species endemic to a semiarid Banded Iron Formation range. Aust. J. Bot. 67, 268–280 (2019).Article 

Google Scholar 
Etchells, H., O’Donnell, A. J., Lachlan McCaw, W. & Grierson, P. F. Fire severity impacts on tree mortality and post-fire recruitment in tall eucalypt forests of southwest Australia. For. Ecol. Manag. 459, 117850 (2020).Article 

Google Scholar 
Bradstock, R. A., Tozer, M. G. & Keith, D. A. Effects of high frequency fire on floristic composition and abundance in a fire-prone heathland near Sydney. Aust. J. Bot. 45, 641–655 (1997).Article 

Google Scholar 
Penman, T. D., Binns, D. L., Brassil, T. E., Shiels, R. J. & Allen, R. M. Long-term changes in understorey vegetation in the absence of wildfire in south-east dry sclerophyll forests. Aust. J. Bot. 57, 533–540 (2010).Article 

Google Scholar 
Ooi, M. K. The importance of fire season when managing threatened plant species: A long-term case-study of a rare Leucopogon species (Ericaceae). J. Environ. Manage. 236, 17–24 (2019).PubMed 
Article 

Google Scholar 
Pausas, J. G., Bradstock, R. A., Keith, D. A. & Keeley, J. E. Plant functional traits in relation to fire in crown-fire ecosystems. Ecology 85, 1085–1100 (2004).Article 

Google Scholar 
Australian Bureau of Meteorology. Climate Data Online. www.bom.gov.au (2019).Abell, R. S. Geoscience map of Jervis Bay Territory and Beecroft peninsula (1:25000 scale). Australian Geological Survey Organisation (1992).Taws, N. Vegetation survey and mapping of Jervis Bay Territory. (Taws Botanical Research, 1997).Taylor, G., Abell, R. & Paterson, I. Geology, geomorphology, soils and earth resources. in Jervis Bay (eds. Cho Arthur, G., Georges, Stoutjesdikj Richard, R., & Longmore) .-. (Australian Nature Conservation Agency, 1995).Keith, D. A. Ocean shores to desert dunes: the native vegetation of NSW and the ACT (Selected Extracts). (Department of Environment and Conservation (NSW), 2004).Keith, D. A. & Tozer, M. G. Vegetation dynamics in coastal heathlands of the Sydney basin. in Proceedings of the Linnean Society of New South Wales vol. 134 (2012).Lindenmayer, D. B. et al. Contrasting mammal responses to vegetation type and fire. Wildl. Res. 35, 395–408 (2008).Article 

Google Scholar 
Bradstock, R. A. & Kenny, B. J. An application of plant functional types to fire management in a conservation reserve in southeastern Australia. J. Veg. Sci. 14, 345–354 (2003).Article 

Google Scholar 
Bowd, E. J., Banks, S. C., Bissett, A., May, T. W. & Lindenmayer, D. B. Direct and indirect disturbance impacts in forests. Ecol. Lett. https://doi.org/10.1111/ele.13741 (2021).Article 
PubMed 

Google Scholar 
Thompson, C. G., Kim, R. S., Aloe, A. M. & Becker, B. J. Extracting the variance inflation factor and other multicollinearity diagnostics from typical regression results. Basic Appl. Soc. Psychol. 39, 81–90 (2017).Article 

Google Scholar 
Fox, J. & Weisberg (Sage, 2019).
Google Scholar 
Venables, W. N. R., B. D. Modern Applied Statistics with S. Fourth Edition. (Springer, 2002).Morrison, D. A. et al. Effects of fire frequency on plant species composition of sandstone communities in the Sydney region: Inter-fire interval and time-since-fire. Aust. J. Ecol. 20, 239–247 (1995).ADS 
Article 

Google Scholar 
Burnham, K. P., Anderson, D. R. & Huyvaert, K. P. AIC model selection and multimodel inference in behavioral ecology: Some background, observations, and comparisons. Behav. Ecol. Sociobiol. 65, 23–35 (2011).Article 

Google Scholar 
Hartig, F. DHARMa: Residual diagnostics for hierarchical (multi-Level / mixed) regression models. (2020).Falster, D. et al. AusTraits, a curated plant trait database for the Australian flora. Sci. Data 8, 254 (2021).PubMed 
PubMed Central 
Article 

Google Scholar 
Tozer, M. G. & Bradstock, R. A. Fire-mediated effects of overstorey on plant species diversity and abundance in an eastern Australian heath. Plant Ecol. 164, 213–223 (2003).Article 

Google Scholar 
Gosper, C. R., Yates, C. J., Prober, S. M. & Parsons, B. C. Contrasting changes in vegetation structure and diversity with time since fire in two Australian Mediterranean-climate plant communities. Austral Ecol. 37, 164–174 (2012).Article 

Google Scholar 
Foster, C., Barton, P., Robinson, N., MacGregor, C. & Lindenmayer, D. B. Effects of a large wildfire on vegetation structure in a variable fire mosaic. Ecol. Appl. 27, 2369–2381 (2017).CAS 
PubMed 
Article 

Google Scholar 
Preston, F. W. The commonness, and rarity, of species. Ecology 29, 254–283 (1948).Article 

Google Scholar 
McGill, B. J. A renaissance in the study of abundance. Science 314, 770–772 (2006).CAS 
PubMed 
Article 

Google Scholar 
Silvertown, J. Plant coexistence and the niche. Trends Ecol. Evol. 19, 605–611 (2004).Article 

Google Scholar 
Lyons, K. G. & Schwartz, M. W. Rare species loss alters ecosystem function—invasion resistance. Ecol. Lett. 4, 358–365 (2001).Article 

Google Scholar 
Dee, L. E. et al. When do ecosystem services depend on rare species?. Trends Ecol. Evol. 34, 746–758 (2019).PubMed 
Article 

Google Scholar 
Smith, M. D. & Knapp, A. K. Dominant species maintain ecosystem function with non-random species loss. Ecol. Lett. 6, 509–517 (2003).Article 

Google Scholar 
Lennon, J. J., Koleff, P., Greenwood, J. J. & Gaston, K. J. Contribution of rarity and commonness to patterns of species richness. Ecol. Lett. 7, 81–87 (2004).Article 

Google Scholar 
Foster, C. N. et al. Herbivory and fire interact to affect forest understory habitat, but not its use by small vertebrates. Anim. Conserv. 19, 15–25 (2016).Article 

Google Scholar 
Lamont, B. B., Enright, N. J. & He, T. Fitness and evolution of resprouters in relation to fire. Plant Ecol. 212, 1945–1957 (2011).Article 

Google Scholar 
Tolhurst, K. G. & Turvey, N. D. Effects of bracken (Pteridium esculentum (forst. f.) cockayne) on eucalypt regeneration in west-central Victoria. For. Ecol. Manag. 54, 45–67 (1992).Candeias, M. & Warren, R. J. Rareness starts early for disturbance-dependent grassland plant species. Biodivers. Conserv. 25, 2771–2785 (2016).Article 

Google Scholar 
Beadle, N. Soil phosphate and the delimitation of plant communities in eastern Australia. Ecology 35, 370–375 (1954).CAS 
Article 

Google Scholar 
Orians, G. H. & Milewski, A. V. Ecology of Australia: The effects of nutrient-poor soils and intense fires. Biol. Rev. 82, 393–423 (2007).PubMed 
Article 

Google Scholar 
Vesk, P. A. & Westoby, M. Funding the bud bank: A review of the costs of buds. Oikos 106, 200–208 (2004).Article 

Google Scholar 
Wilfahrt, P. et al. Temporal rarity is a better predictor of local extinction risk than spatial rarity. Ecology https://doi.org/10.1002/ecy.3504 (2021).Article 
PubMed 

Google Scholar 
Miller, B. P. et al. Persistence of rare species depends on rare events: Demography, fire response and phenology of two plant species endemic to a semiarid Banded Iron Formation range. Aust. J. Bot. 67, 268–280 (2019).Article 

Google Scholar 
Gillespie, I. G. & Allen, E. B. Fire and competition in a southern California grassland: Impacts on the rare forb Erodium macrophyllum. J. Appl. Ecol. 41, 643–652 (2004).Article 

Google Scholar 
Maire, V. et al. Habitat filtering and niche differentiation jointly explain species relative abundance within grassland communities along fertility and disturbance gradients. New Phytol. 196, 497–509 (2012).PubMed 
Article 

Google Scholar 
Yenni, G., Adler, P. B. & Ernest, S. M. Do persistent rare species experience stronger negative frequency dependence than common species?. Glob. Ecol. Biogeogr. 26, 513–523 (2017).Article 

Google Scholar 
Mayberry, R. J. & Elle, E. Conservation of a rare plant requires different methods in different habitats: Demographic lessons from Actaea elata. Oecologia 164, 1121–1130 (2010).ADS 
PubMed 
Article 

Google Scholar 
Rabinowitz, D. & Rapp, J. K. Dispersal abilities of seven sparse and common grasses froma Missouri prairie. Am. J. Bot. 68, 616–624 (1981).Article 

Google Scholar 
McIntyre, S. Comparison of a common, rare and declining plant species in the Asteraceae: Possible causes of rarity. Pac. Conserv. Biol. 2, 177–190 (1995).Article 

Google Scholar 
Hopfensperger, K. N. A review of similarity between seed bank and standing vegetation across ecosystems. Oikos 116, 1438–1448 (2007).Article 

Google Scholar 
Cross, A. T. et al. Defining the role of fire in alleviating seed dormancy in a rare Mediterranean endemic subshrub. AoB Plants 9, (2017). More

Global vegetation physiological response to increasing atmospheric CO2 and its reduction of mitigation potentialWe calculate different effects of increasing CO2 on mean annual near-surface air temperature change over global vegetated land. We compare the direct atmospheric radiative effect (RAD)-induced climate warming to the temperature reductions caused by the BGC effect. These two global temperature changes are, in turn, then compared against our main focus of aggregated local PHY-induced temperature contributions (PHYall; Fig.1). The main finding is that the spatial aggregation of PHY feedbacks on temperature over global vegetated land (green bars) offsets a substantial amount of the cooling effect through enhanced terrestrial carbon storage because of the BGC effects (blue bars). Terrestrial carbon storage continuously increases with rising atmospheric CO2, and reaches a global total of 621 ± 260 Pg C under 4 × CO2 (Supplementary Fig. 1). This increased land carbon storage is equivalent to 293 ± 122 ppm of CO2 removal from the atmosphere and results in a temperature cooling of −1.24 ± 0.57 °C. The PHYall-induced land temperature increase is 0.83 ± 0.47 °C, from the ensemble of five ESMs we use (here we excluded bcc-csm1-1; see “Methods”), corresponding to a large offset of the cooling effect by terrestrial ecosystems through BGC. However, our estimated temperature cooling induced by the BGC effect is a transient response. This cooling effect will be larger for subsequently stabilised CO2 concentrations, since terrestrial ecosystems continue to fix carbon until reaching equilibrium. We note that BGC-induced cooling may be overestimated in the absence of land-cover change effect in the simulations, as the latter may reduce terrestrial carbon stores (Supplementary Fig. 2). However, inter-model differences, such as different parameterisation or different biogeochemistry module20, may prevent a definitive answer as to how land-cover change influence global temperature change through BGC (Supplementary Fig. 2–5). We also note that we focus our analysis on CMIP5 data mainly due to the fact that the parameters in Eqs. (1) and (2)21,22 (see “Methods”) for calculating BGC-induced cooling are only currently available for CMIP5 ESMs.Fig. 1: Climate warming mitigation potential of terrestrial ecosystems.Global mean land temperature change due to total CO2 physiological forcing (PHYall), increased terrestrial carbon storage (BGC) and CO2 radiative forcing (RAD). Note axes as coloured, and that the vertical blue axis is temperature cooling through BGC. It is straightforward to compare PHYall-, RAD- and BGC-induced temperature change when the same range and directions of bars are used. Three atmospheric CO2 horizons are selected here, for CO2 concentrations of 500 ppm, 800 ppm and 1032 ppm (4 × CO2). Each bar represents the global area-weighted average from the ensemble of five ESMs. The error bars indicate the standard error of these five models. Note, the bcc-csm1-1 model did not provide diagnostics of carbon storage data and so is not included in calculating the temperature change by carbon storage change.Full size imageThe RAD response to rising CO2 projects a global warming of 5.25 ± 0.65 °C, again under 4 × CO2, corresponding to roughly four times the magnitude of BGC-induced cooling. Hence, PHYall-induced temperature increase adds about 16 ± 8% to the RAD-induced warming globally. The relative magnitude of warming/cooling effects is similar for the lower CO2 levels of 500 ppm and 800 ppm (Fig. 1), illustrating the importance of accounting for PHYall-induced warming and how it affects the ability of terrestrial vegetation to mitigate global warming, irrespective of CO2 concentration. Our identified PHYall feedbacks are a combination of different altered land properties. The balance of changed land components will depend on location, and so spatial variations in the overall warming effect may suggest a reappraisal of some climate change adaptation measures. Hence to aid such assessments, we now consider in detail the global contributions of individual drivers of the PHY-based local warming, and then any geographical variations.In Fig. 2a, we show changes in global vegetated land air temperature associated with PHYall with increasing atmospheric CO2 from our ensemble of six ESMs (see “Methods”). PHYdir (see “Methods”; Supplementary Table 1) is based on the direct CO2 effect on vegetation physiology. PHYall represents all the vegetation physiology-related feedbacks, and captures any additional interactions between RAD and PHYdir, due to all effects not being a simple linear addition of RAD and PHY responses. The difference, PHYall minus PHYdir, is termed PHYint. Although inter-model difference exists, warming levels are projected to increase, by all the ESMs, as CO2 concentration rises, and for both PHYdir and PHYall (Fig. 2a). The interactive effects result in PHYall-induced temperature change being higher than PHYdir throughout the period. For the smaller increases in atmospheric CO2 of up to ~450 ppm, interactive effects dominate with PHYdir being almost zero up to that concentration. However, above the CO2 concentration of 450 ppm, PHYdir increases global warming reaching 0.17 ± 0.06 °C for CO2 of 500 ppm, and climbing to 0.62 ± 0.48 °C under quadrupled atmospheric CO2 (Fig. 2a). At that 4 × CO2 level, interaction term PHYint increases warming by approximately one-fifth of that induced by PHYdir. Hence, the overall physiological feedbacks described by PHYall produce a global temperature increase of 0.74 ± 0.47 °C under 4 × CO2 (again from multi-model mean of six ESMs; see “Methods”).Fig. 2: Change in global mean annual land air temperature and individual climate forcing induced by vegetation physiological response to increasing atmospheric CO2.a Global annual area-weighted temperature change of vegetated land induced by total CO2 physiological forcing (PHYall) and the direct CO2 physiological forcing (PHYdir) in response to increasing atmospheric CO2 concentration. Shaded areas are the standard errors of the six Earth System Models (ESMs) used, and the thick curves are their multi-model means. For each atmospheric CO2 concentration in panel (a), values are based on smoothing using a twenty-year running window (to match with the decomposition results in panel (b). The final temperature change induced by PHYall and PHYdir effects under 4 × CO2 are further marked on the righthand side, with the “+” markers indicating multi-model means. b PHYall-induced climate forcing associated with changes in albedo, aerodynamic resistance (ra), evapotranspiration (ET), downwelling shortwave radiation (SW) and near-surface air emissivity (ɛa). Again, the shaded areas are the standard errors of the models, and the mean is the thick continuous lines. The changes in these variables are calculated using a moving average with a 20-year window. The resulting values under 4 × CO2 are plotted on the righthand side.Full size imageTo better understand the factors influencing CO2 physiological drivers of temperature change, we decompose the global PHYall into individual biophysical components6 (see “Methods”). These five aspects are albedo, aerodynamic resistance (ra), evapotranspiration (ET), downwelling shortwave radiation (SW) and near-surface air emissivity (ɛa). These five-component changes result in climate forcings with different signs and magnitudes (Fig. 2b), and thus perturb the surface energy balance, where positive values correspond to an increase in temperature. Specifically, ET, SW, and albedo generates positive climate forcings that increase local temperature, while ra and ɛa produce negative effects and thus offset local temperature increase. The relative role of each biophysical component in influencing PHYall-induced temperature change remains largely invariant in the transition from low to high CO2 concentrations. In addition, the changes in these five quantities affecting PHYdir show similar variations with that of PHYall as atmospheric CO2 rises (Supplementary Fig. 6). Those results are generally valid for CMIP6 results but with a lower magnitude of change (Supplementary Fig. 7).The vegetation physiological response to rising CO2 causes changes in LAI and stomata closure, which adjust in parallel with more detailed attributes of the land surface. The ESM simulations of changes in LAI and transpiration compare moderately well with available field measurements23,24,25,26,27,28,29,30,31,32 (Supplementary Fig. 8). Here we focus on their effects on the near-surface thermal changes, expressed as climate forcings on near-surface energy fluxes (Fig. 2b). The LAI increase due to elevated CO2 (Supplementary Fig. 9a) leads to decreases in albedo11,33 (Supplementary Fig. 10a; Supplementary Table 2) and ra7 (Supplementary Fig. 10c). These changes are continuous with rising atmospheric CO2 and result in positive and negative effects on global land air temperature change, respectively (Fig. 2b). Specifically, the albedo reduction increases solar radiation absorption by the land surface11,33, and imposes a positive forcing of 0.30 ± 0.14 W m−2 for 500 ppm, and of 0.51 ± 0.45 W m−2 for a quadrupling of CO2 (Fig. 2b). The decreased ra favours the turbulent transport of heat from land to atmosphere7, and leads to a persistent surface cooling with increasing atmospheric CO234,35, which is −1.66 ± 0.44 W m−2 for 500 ppm and −3.15 ± 2.05 W m−2 for 4 × CO2. In contrast, ET decreases considerably (Supplementary Fig. 10e; Supplementary Table 2) due to decreased stomatal conductance responding to increasing atmospheric CO2. These reductions in ET reduce evaporative cooling15,17, and therefore result in a strong positive climate forcing on global warming (Fig. 2b), which increases with rising atmospheric CO2 and reaches 0.88 ± 0.25 W m−2 for atmospheric CO2 of 500 ppm and 2.88 ± 1.49 W m−2 with quadrupled atmospheric CO2 concentration (Fig. 2b).Moreover, ET reduction decreases the inflow of evaporative water to the atmosphere, reducing cloud fraction and water vapour content, thereby feeding back to impose indirect effects that influence temperature change. We quantify these indirect effects through their effects on the components of SW and ɛa, (Fig. 2b). The lower ET values reduce cloud fraction (Supplementary Table 2; Supplementary Fig. 9c) thereby increasing the amount of SW to the land surface6,35,36 (Supplementary Fig. 10g), and producing a moderate positive forcing of 0.67 ± 0.16 W m−2 for 500 ppm and 1.76 ± 1.26 W m−2 under 4 × CO2 (Fig. 2b). However, the reduced cloud fraction and water vapour content decreases ɛa (Supplementary Fig. 10i), which weakens the absorption of longwave radiation37 (Supplementary Fig. 9e). The ɛa changes result in a small negative forcing (−0.28 ± 0.11 W m−2 for 500 ppm and −0.57 ± 0.43 W m−2 for 4 × CO2), and thereby lowering global temperatures. Of particular note is that the direct changes in ET produce the strongest positive forcing for warming (Fig. 2b). Hence, ET changes have a dominant role in influencing the magnitude and sign of PHYall feedbacks on global temperature change. The changes in these five factors for PHYdir are very close in relative terms to those of PHYall (Supplementary Fig. 10), driving climate forcings in a similar magnitude and sign to influence PHYdir-induced global temperature change (Supplementary Fig. 6).Overall, as CO2 rises, these five biophysical factors cause vegetation to amplify warming locally, and when aggregated spatially act to raise planetary global warming. This additional warming effect is primarily driven by changes in ET, with smaller warming contributions from changes in SW and albedo, but also compensated with cooling effects from ra and ɛa changes. The substantial role of ET on biophysical climate feedbacks is consistent with a previous study6, which investigates the biophysical feedbacks because of vegetation greening. That study proposed that CO2-driven greening will enhance ET and thereby producing a net cooling effect. We build on that analysis by here additionally including the CO2-induced partial stomatal closure17,18. We find this inclusion overtakes the LAI influence on ET changes, resulting in large reductions in ET that will instead contribute to net warming as CO2 rises. Large inter-model differences in simulating ra dynamics (Supplementary Fig. 10c, d) mostly explains the substantial spread of ra-induced climate forcing (Fig. 2b). The different extents of LAI increase (Supplementary Fig. 9a, b) also contribute to such differences in simulated ra changes and effects on temperature change. Our noted large cooling through ra changes is also indicated in a recent study7, suggesting the importance of ra in affecting vegetation biophysical climate feedbacks, and the importance of constraining this factor in climate models.Spatial patterns and attributions of vegetation physiologically induced warmingWe present in Fig. 3 the area-weighted regional contributions to global temperature change of the physiological responses and BGC under 4 × CO2. PHYall-based warming and BGC-induced cooling both show larger values in East and Central North America (ENA and CNA), North and Central Europe (NEU and CEU), Amazon (AMZ) and North Asia (NAS). Whereas, the smallest changes are in the Sahel region (SAH) (Fig. 3a). PHYall-induced warming reduces large proportions of the temperature cooling through BGC in the northern mid-to-high latitudes (Fig. 3b). Furthermore, this cooling effect is fully offset by warming through vegetation physiological response in Alaska (ALA), Canada/Greenland/Iceland (CGI), NAS, NEU, SAH, Tibetan Plateau (TIB), Central (CAS) and West Asia (WAS), and West North America (WNA), resulting in slight warming in these regions. Presented as a map of the net effects of BGC and PHYall (Fig. 3b), we further illustrate this overall cooling effect from PHY and BGC for the tropical regions and South Hemisphere. Additionally, the balance of PHY and BGC to influence regional temperature change under relatively low atmospheric CO2 level of 500 ppm (Supplementary Fig. 11) is very close to that for 4 × CO2, suggesting consistency of vegetation biophysical and biogeochemical effects on climate irrespective of CO2 levels. In summary, the PHYall feedbacks offset the cooling benefits from ecosystem “draw-down” of CO2 to a large extent, in line with the global average values presented in Fig. 1. Northern mid-to-high latitudes may contribute less than expected to slow global warming (Fig. 3). However, as also noted for the global change values, the estimated temperature cooling by BGC is a transient effect that would keep increasing as the ecosystems approach their equilibrium. However, BGC-induced cooling may be overestimated without consideration of the effect of land-cover change (Supplementary Fig. 2), which is often associated with the deliberate removal of terrestrial carbon.Fig. 3: Regional contributions to temperature change by PHY and BGC under 4 × CO2.a Contribution of regional vegetation physiological responses (PHYall; green bars) and increased carbon storage (BGC; blue bars) to the overall global temperature change for each of the IPCC AR5 SREX regions. Bars represent area-weighted multi-model means, and the error bars indicate the standard errors of the models for each region. b Spatial distribution of the net effects of warming induced by PHYall and cooling through BGC.Full size imageThe regional pattern in Fig. 3b provides a motivation to investigate further the five climate forcings-driven PHYall contributions (Fig. 2a) to near-surface temperature change triggered by rising CO2. Here, we first examine the much finer spatial patterns and component contributions of PHYall, PHYdir and PHYint on warming for 4 × CO2. Focussing on PHYdir and typically for 4 × CO2, there is a larger local warming in the tropical forests (Fig. 4a), where vegetation shows higher potential to stabilise carbon than other ecosystems10 (Fig. 3). Within these tropical regions, the PHYdir-forced temperature increase is the highest in the Amazon forests (Fig. 4a), reaching approximately 30% of that induced by RAD (Supplementary Fig. 12b). Strong warming by PHYdir is also found in the northern mid-to-high latitudes ( >40°N). In contrast, smaller temperature increases due to PHYdir are seen in arid and semi-arid regions, such as Australia and Sahel (Fig. 4a). The spatial variations of PHYall-forced temperature change (Supplementary Fig. 13a) have strong similarities to those of PHYdir (Fig. 4a). These similarities again suggest a relatively small role of interactions on temperature change (Supplementary Fig. 13b) under quadrupled CO2. Additionally, the agreements across the six ESMs are reasonably high, with all the models agreeing that PHYall, PHYdir and PHYint result in local warming for 4 × CO2 across most of the global vegetated land (Supplementary Fig. 14).Fig. 4: Global patterns of local temperature change and climate forcings through vegetation physiological response to 4 × CO2.Spatial distribution of annual mean temperature change from multi-model ensemble induced by (a). direct CO2 physiological forcing (PHYdir) in response to a 4 × CO2 rise since pre-industrial level. The spatial patterns of the individual climate forcing contributing to PHYdir of panel (a) are as follows. In (b). albedo (α), c Aerodynamic resistance (ra), d Evapotranspiration (ET), e Downwelling shortwave radiation (SW) and (f). near-surface air emissivity (ɛa). In all panels, estimates use the mean of the final 20 years of the simulations (atmospheric CO2 at ~1032 ppm).Full size imageWe next analyse spatially the decomposition of PHYall- and PHYdir-induced temperature change into the five biophysical factors shown in Fig. 2b, and with findings shown in Fig. 4b–f. In response to quadrupled CO2, LAI shows increases over most global vegetated land (Supplementary Fig. 15b), leading to a positive climate forcing (and thus warming) through reducing albedo (Fig. 4b), and especially in the south Sahel and Tibetan Plateau. Moreover, a large albedo decrease occurs in the northern mid-to-high latitudes, such as for most of Siberia. This albedo decrease may be due to LAI increase combined with reduced snow cover promoted by PHY-driven warming17. The LAI increase also contributes to a strong negative forcing through ra decrease7. This ra decrease enhances the energy exchange between land and atmosphere, lowering the local warming effect in response to 4 × CO2, particularly in Australia, Sahel, South Africa and South Asia (Fig. 4c). In comparison, the projected forcing from ET reductions (Supplementary Fig. 15d) causes strongly positive warming, and especially in the tropical and boreal forests (Fig. 4d). The larger ET reductions in the tropical and boreal forests also lead to stronger cloud fraction decreases (Supplementary Fig. 15f). These feedbacks induce a positive climate forcing for additional warming by increasing SW reaching the land surface (Fig. 4e). Simultaneously, the decreased cloud cover and water vapour by the ET reductions generate a cooling effect (Fig. 4f) through decreasing net longwave radiation absorption in most locations. The PHYdir-induced warming and the associated climate forcings at an atmospheric CO2 concentration of 500 ppm (Supplementary Fig. 16) show similar spatial distributions, but smaller magnitude, compared with that adjustment for 4 × CO2 increase (Fig. 4). This result manifests that PHY continues to amplify global warming through the combined climate forcings from our identified five biophysical factors, irrespective of CO2 concentration.Substantial differences exist between PHYall- and PHYdir-induced temperature changes in northwestern Eurasia, implying pronounced PHYint-based feedbacks there under 4 × CO2 (Supplementary Fig. 13b). A relatively large change in SW, inducing a warming effect (Supplementary Fig. 13j), mainly contributes to the large PHYint-coupled changes in this region. Elsewhere, for the Amazon forests, the cooling effect in PHYint, through changes in ET (Supplementary Fig. 13h) and SW (Supplementary Fig. 13j), cancels out the warming due to changes in ra (Supplementary Fig. 13f) and ɛa (Supplementary Fig. 13l). This cancellation causes negligible PHYint feedback on temperature there (Supplementary Fig. 13b). In summary, local warming due to biophysical feedback of increasing CO2 is regulated primarily by the positive forcing from ET reductions (with smaller positive forcings from albedo and SW changes), which is partly compensated by reductions in ra (and a small contribution from ɛa).Strong physiologically induced warming in dense ecosystemsWe further investigate the finding presented in Fig. 4 that PHYdir-forced temperature increase is higher in the tropical forests, while lower in the arid and semi-arid ecosystems. This finding suggests that there may be a relationship between the forced temperature change and background baseline LAI, which we confirm in Fig. 5a (for both PHYall and PHYdir). We find widespread warming amplification with increasing LAI for lower background LAI levels. Such rates of increase in temperature per unit of LAI show evidence of flattening at higher LAI values38. That is, the CO2-fertilised LAI increase saturates in ecosystems with dense canopy cover such as tropical forests (Supplementary Fig. 17a). However, the stomatal closure-induced ET decrease varies almost linearly with baseline LAI, although the reduction breaks down when the baseline LAI approaches six (Supplementary Fig. 17b). Hence, we conclude that ET reductions caused by stomatal closure are the primary cause of PHY-induced temperature rises, although changes in ra because of LAI increases remain an important factor (Fig. 4; Supplementary Fig. 13). Taking all the components together, PHYdir-triggered temperature change increases almost linearly with baseline LAI gradients under 4 × CO2 (Fig. 5a). Of particular interest is that when expressing PHYdir-induced temperature increase as a fraction of warming induced by RAD, it also increases with the baseline LAI value (Fig. 5b). Hence, Fig. 5 provides strong evidence that background vegetation functional structure (LAI) influences the feedback of CO2 physiological forcing on warming in response to 4 × CO2. In terms of change, as the trend of global LAI increase slows down as atmospheric CO2 concentration increases to very high levels (Supplementary Fig. 9a, b), this suggests that the magnitude of the cooling effect (via ra) will decrease relative to the ET-based warming.Fig. 5: Effects of vegetation structure on CO2 physiological forcing in response to 4 × CO2.a Variations of local temperature change induced by total (PHYall) and direct (PHYdir) CO2 physiological forcing, for quadrupled CO2 and presented as a function of baseline pre-industrial leaf area index (LAI). b Variations of the ratio of PHYall- and PHYdir-induced temperature change relative to that of RAD warming, again for 4 × CO2 and as a function of baseline LAI. In both panels, the temperature changes along LAI gradients are smoothed using a five-bin running window for the 0.1 increments in LAI bin size. The shaded areas indicate the standard error among different models.Full size imageImplications and conclusionsTerrestrial ecosystems respond to increasing atmospheric CO2 not only through accumulating carbon (the “BGC” effect) but also through their physiological response (the “PHY” feedback), which can have opposing effects on local surface temperature. As such, to fully understand the net climate benefits of terrestrial ecosystems, comprehensive assessments need to account for both PHY and BGC effects39,40. We find that the warming through PHY feedback largely reduces the capacity of vegetation to slow global warming through BGC. In particular, vegetation transiently operates to amplify local warming in the northern mid-to-high latitudes, as PHY-induced warming is larger than BGC-induced cooling (Fig. 3). Tropical forests have net cooling effects, and forests and ecosystem restoration efforts there would have much higher net cooling benefit for mitigation than if they are placed in the northern mid-to-high latitudes. At the global scale, this PHY-induced warming can offset by up to 67% of the cooling gains from the transient BGC effect (Fig. 1). This is likely an upper-bound value, because the steady-state BGC effect is larger than the transient one. Thus, afforestation can still result in a net cooling effect, especially in the tropics, although the cooling achieved may be less than first expected. Of particular interest is that the PHY-based warming is generally stronger for high baseline LAI values (Fig. 5). This result further suggests the PHY produces extra warming through afforestation due to higher LAI for forests than non-forests. However, this does not contradict the finding for the tropics (where LAI is high) as possibly the most beneficial regions for afforestation, due to the compensating much larger BGC potential to lower temperature for such locations.Our analysis highlights the importance of including PHY feedbacks in any assessments of future levels of global warming. Such understanding may be especially important for the northern mid-to-high latitudes. We point out that only one-third of ESMs we use incorporate dynamic global vegetation schemes (DGVMs). Since land-cover change acts to influence PHY and BGC effects (Supplementary Fig. 2–5), future simulations with dynamic vegetation may allow refining these results. However, we note that inter-model structural differences (e.g. different ratio of transpiration to ET15,41) contribute to large uncertainties in these results. In particular, we note the potential implications of our results for climate mitigation policies with substantial reforestation as a mechanism to slow global warming through emissions offsetting. Reforestation may cause a near-surface cooling through the BGC effect. However, the PHY-based warming effect is stronger for higher LAI values as atmospheric CO2 rises (Fig. 5a). This must be accounted for in any major global reforestation plans with consideration of particular location of reforestation measures42. In general, the effects of anthropogenic land use and land cover change on vegetation biophysical and biogeochemical processes are important and complex with substantial spatial heterogeneity42,43,44,45,46,47. Given the lack of exploitable factorial simulations designed to separate land use and land-cover change effects, these are ignored in the idealised simulations, and may introduce a bias in our results. We note, however, that these idealized “4 × CO2” simulations broadly capture the sign and spatial distributions of land-atmosphere coupling effects on temperature change through comparisons to the ESM results under the RCP8.5 scenario (Supplementary Fig. 18) which include anthropogenic land use and land-cover change. We suggest that future climate projections should much more routinely account for the effect of anthropogenic land use and land-cover change together with PHY and BGC effects42,47, and so that this can be considered in reforestation climate mitigation strategies.An additional caveat of our research is that the PHY and RAD factorial ESM simulations that inform our analysis, are at present only available for the illustrative but potentially unrealistic exponential increase in atmospheric CO2 (1% per year). For this scenario, the systems can exhibit a unrealistic linear behaviour48. PHY and RAD simulations under other scenarios, such as abrupt 4 × CO2 or even historical forcing followed by a potential future scenario, would greatly help to improve projections of PHY feedbacks. We suggest this should be made a high priority of climate modelling exercises, especially if additionally including land use and land-cover change representation. We recognise that the newer CMIP6 ESMs contain improved representation of many processes (e.g. atmospheric aerosols, clouds, and land processes)49. We anticipate our findings to remain broadly valid with the CMIP6 models, although the magnitude of PHY-induced warming is lower for CMIP6 than CMIP5 ESMs18, especially in the northern high latitudes. We hope our analysis provides an incentive to others to undertake such analysis as all the calculations (BGC, PHY, and RAD) become available for that newer CMIP6 ensemble. In summary, our results illustrate that vegetation physiological response to increasing atmospheric CO2 has a substantial local warming effect, which requires consideration alongside the cooling effect vegetation offers by “drawing down” atmospheric carbon dioxide. More

Karlsson, O. et al. Pesticide-induced multigenerational effects on amphibian reproduction and metabolism. Sci. Total Environ. 775, 145771 (2021).ADS 
CAS 
PubMed 
Article 

Google Scholar 
IUCN. The IUCN Red List of Threatened Species. Version 2021-3. https://www.iucnredlist.org (2022).Wake, D. B. & Koo, M. S. Amphibians. Curr. Biol. 28, R1237–R1241 (2018).CAS 
PubMed 
Article 

Google Scholar 
Campbell Grant, E. H., Miller, D. A. & Muths, E. A synthesis of evidence of drivers of amphibian declines. Herpetologica 76, 101–107 (2020).Article 

Google Scholar 
Green, D. M., Lannoo, M. J., Lesbarrères, D. & Muths, E. Amphibian population declines: 30 years of progress in confronting a complex problem. Herpetologica 76, 97–100 (2020).Article 

Google Scholar 
Mason, R., Tennekes, H., Sánchez-Bayo, F. & Jepsen, P. U. Immune suppression by neonicotinoid insecticides at the root of global wildlife declines. J. Environ. Immunol. Toxicol. 1, 3–12 (2013).Article 

Google Scholar 
Adams, E., Leeb, C. & Brühl, C. A. Pesticide exposure affects reproductive capacity of common toads (Bufo bufo) in a viticultural landscape. Ecotoxicology 30, 213–223 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Frost, D. R. Amphibian species of the world 6,1, an online reference. Electron. Datab. https://doi.org/10.5531/db.vz.0001 (American Museum of Natural History, 2021).Article 

Google Scholar 
Eterovick, P. C., Souza, A. M. & Sazima, I. Anfíbios da Serra do Cipó [Amphibians from the Serra do Cipó]. http://herpeto.org/wp-content/uploads/2020/11/ANFIBIOS-DA-SERRA-DO-CIPO.pdf (PUCMINAS, 2020).Mijares, A., Rodrigues, M. T. & Baldo, D. Physalaemus cuvieri The IUCN Red List of Threatened Species, version 2014.3. http://www.iucnredlist.org (2010). Accessed 9 Jan 2015.de Sá, F. P., Zina, J. & Haddad, C. F. B. Reproductive dynamics of the Neotropical treefrog Hypsiboas albopunctatus (Anura, Hylidae). J. Herpetol. 48, 181–185 (2014).Article 

Google Scholar 
Herek, J. S. et al. Can environmental concentrations of glyphosate affect survival and cause malformation in amphibians? Effects from a glyphosate-based herbicide on Physalaemus cuvieri and P. gracilis (Anura: Leptodactylidae). Environ. Sci. Pollut. Res. 27, 22619–22630 (2020).CAS 
Article 

Google Scholar 
Silva, F. L. et al. Swimming ability in tadpoles of Physalaemus cf. cuvieri, Scinax x-signatus and Leptodactylus latrans (Amphibia: Anura) exposed to the insecticide chlorpyrifos. Ecotoxicol. Environ. Contam. 16, 13–18 (2021).
Google Scholar 
Pavan, F. A. et al. Morphological, behavioral and genotoxic effects of glyphosate and 2,4-D mixture in tadpoles of two native species of South American amphibians. Environ. Toxicol. Pharmacol. 85, 103637 (2021).CAS 
PubMed 
Article 

Google Scholar 
Simon-Delso, N. et al. Systemic insecticides (Neonicotinoids and fipronil): Trends, uses, mode of action and metabolites. Environ. Sci. Pollut. Res. 22, 5–34 (2015).CAS 
Article 

Google Scholar 
Pietrzak, D., Kania, J., Malina, G., Kmiecik, E. & Wątor, K. Pesticides from the EU first and second watch lists in the water environment. Clean 47, 1–10 (2019).
Google Scholar 
IBAMA: Instituto Brasileiro do Meio Ambiente e dos Recursos Naturais Renováveis. Relatório de comercialização de agrotóxicos 2019 [Brazilian Pesticide Marketing Report 2019] https://www.ibama.gov.br/agrotoxicos/relatorios-de-comercializacao-de-agrotoxicos#boletinsanuais (2021).IBAMA: Instituto Brasileiro do Meio Ambiente e dos Recursos Naturais Renováveis. Vendas de ingredientes ativos por UF [Active ingredient sales by UF in Brazil]. http://ibama.gov.br/phocadownload/qualidadeambiental/relatorios/2019/Vendas_ingredientes_ativos_UF_2019.x (2021).IBAMA – Instituto Brasileiro do Meio Ambiente e dos Recursos Naturais Renováveis. Boletins anuais de produção, importação, exportação e vendas de agrotóxicos no Brasil [Annual bulletins of production, import, export and sales of pesticides in Brazil]. http://ibama.gov.br/index.php?option=com_content&view=article&id=594&Itemid=54 (2021).Pietrzak, D., Kania, J., Kmiecik, E., Malina, G. & Wątor, K. Fate of selected neonicotinoid insecticides in soil–water systems: Current state of the art and knowledge gaps. Chemosphere 255, 126981 (2020).ADS 
CAS 
PubMed 
Article 

Google Scholar 
ANVISA: Agência Nacional de Vigilância Sanitária; Índice Monográfico I13. Imidacloprido. http://portal.anvisa.gov.br/documents/111215/117782/I13+%E2%80%93+Imidacloprido/9d08c7e5-8979-4ee9-b76c-1092899514d7 (2021).Kagabu, S. Discovery of imidacloprid and further developments from strategic molecular designs. J. Agric. Food Chem. 59, 2887–2896 (2011).CAS 
PubMed 
Article 

Google Scholar 
Tomizawa, M. & Casida, J. E. Neonicotinoid insecticide toxicology: Mechanisms of selective action. Annu. Rev. Pharmacol. Toxicol. 45, 247–268 (2005).CAS 
PubMed 
Article 

Google Scholar 
Hashimoto, F. et al. Occurrence of imidacloprid and its transformation product (imidacloprid-nitroguanidine) in rivers during an irrigating and soil puddling duration. Microchem. J. 153, 12 (2020).Article 
CAS 

Google Scholar 
Hladik, M. L. et al. Year-round presence of neonicotinoid insecticides in tributaries to the Great Lakes, USA. Environ. Pollut. 235, 1022–1029 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Jurado, A., Walther, M. & Díaz-Cruz, M. Occurrence, fate and environmental risk assessment of the organic microcontaminants included in the Watch Lists set by EU Decisions 2015/495 and 2018/840 in the groundwater of Spain. Sci. Total Environ. 663, 285–296 (2019).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Montagner, C. C. et al. Ten years-snapshot of the occurrence of emerging contaminants in drinking, surface and ground waters and wastewaters from São Paulo State, Brazil. J. Braz. Chem. Soc. 30, 614–632 (2019).CAS 

Google Scholar 
CCME. Council of Ministers of the Environment. Canadian water quality guidelines for the protection of aquatic life. Imidacloprid. In Canadian water quality guidelines, Council of Ministers of the Environment. Winnipeg. https://ccme.ca/en/res/imidacloprid-en-canadian-water-quality-guidelines-for-the-protection-of-aquatic-life.pdf (2007).RIVM. Water quality standards for imidacloprid: Proposal for an update according to the Water Framework Directive in National Institute for Public Health and the Environment. https://www.rivm.nl/bibliotheek/rapporten/270006001.pdf (2014).PAN. Pesticide Action Network. International Consolidated List of Banned Pesticides. https://pan-international.org/pan-international-consolidated-list-of-banned-pesticides/ (2021).Brazil. Secretaria Estadual da Saúde do Rio Grande do Sul. Portaria SES RS nº 320, de 28 de abril de 2014. https://www.cevs.rs.gov.br/upload/arquivos/201705/11110603-portaria-agrotoxicos-n-320-de-28-de-abril-de-2014.pdf. (2014).Kobashi, K. et al. Comparative ecotoxicity of imidacloprid and dinotefuran to aquatic insects in rice mesocosms. Ecotoxicol. Environ. Saf. 138, 122–129 (2017).CAS 
PubMed 
Article 

Google Scholar 
Islam, M. A., Hossen, M. S., Sumon, K. A. & Rahman, M. M. Acute toxicity of imidacloprid on the developmental stages of common carp Cyprinus carpio. Toxicol. Environ. Health Sci. 11, 244–251 (2019).Article 

Google Scholar 
Pérez-Iglesias, J. M. et al. The genotoxic effects of the imidacloprid-based insecticide formulation Glacoxan Imida on Montevideo tree frog Hypsiboas pulchellus tadpoles (Anura, Hylidae). Ecotoxicol. Environ. Saf. 104, 120–126 (2014).PubMed 
Article 
CAS 

Google Scholar 
Sievers, M., Hale, R., Swearer, S. E. & Parris, K. M. Contaminant mixtures interact to impair predator-avoidance behaviours and survival in a larval amphibian. Ecotoxicol. Environ. Saf. 161, 482–488 (2018).CAS 
PubMed 
Article 

Google Scholar 
USEPA. United States Environmental Protection Agency. Aquatic Life Benchmarks and Ecological Risk Assessments for Registered Pesticides. https://www.epa.gov/pesticide-science-and-assessing-pesticide-risks/aquatic-life-benchmarks-and-ecological-risk. (2021).Feng, S., Kong, Z., Wang, X., Zhao, L. & Peng, P. Acute toxicity and genotoxicity of two novel pesticides on amphibian, Rana N. Hallwell. Chemosphere 56, 457–463 (2004).ADS 
CAS 
PubMed 
Article 

Google Scholar 
De Arcaute, C. R. et al. Genotoxicity evaluation of the insecticide imidacloprid on circulating blood cells of Montevideo tree frog Hypsiboas pulchellus tadpoles (Anura, Hylidae) by comet and micronucleus bioassays. Ecol. Indic. 45, 632–639 (2014).Article 
CAS 

Google Scholar 
Nkontcheu, D. B. K., Tchamadeu, N. N., Ngealekeleoh, F. & Nchase, S. Ecotoxicological effects of imidacloprid and lambda-cyhalothrin (insecticide) on tadpoles of the African common toad, Amietophrynus regularis (Reuss, 1833) (Amphibia: Bufonidae). Emerg. Sci. J. 1, 49–53 (2017).
Google Scholar 
Bortoluzzi, E. C. et al. Contaminação de águas superficiais por agrotóxicos em função do uso do solo numa microbacia hidrográfica de Agudo, RS. Rev. Bras. Eng. Agric. Ambient. 10, 881–887 (2006).Article 

Google Scholar 
Bortoluzzi, E. C. et al. Investigation of the occurrence of pesticide residues in rural wells and surface water following application to tobacco. Quim. Nova 30, 1872–1876 (2007)CAS 
Article 

Google Scholar 
La, N., Lamers, M., Bannwarth, M., Nguyen, V. V. & Streck, T. Imidacloprid concentrations in paddy rice fields in northern Vietnam: measurement and probabilistic modeling. Paddy Water Environ. 13, 191–203 (2015).Article 

Google Scholar 
Sweeney, M. R., Thompson, C. M. & Popescu, V. D. Sublethal, behavioral, and developmental effects of the neonicotinoid pesticide imidacloprid on larval wood frogs (Rana sylvatica). Environ. Toxicol. Chem. 40, 1838–1847 (2021).Article 
CAS 

Google Scholar 
Gibbons, D., Morrissey, C. & Mineau, P. A review of the direct and indirect effects of neonicotinoids and fipronil on vertebrate wildlife. Environ. Sci. Pollut. Res. 22, 103–118 (2015).CAS 
Article 

Google Scholar 
Morrissey, C. A. et al. Neonicotinoid contamination of global surface waters and associated risk to aquatic invertebrates: A review. Environ. Int. 74, 150920 (2015).Article 
CAS 

Google Scholar 
Stinson, S. A. et al. Agricultural surface water, imidacloprid, and chlorantraniliprole result in altered gene expression and receptor activation in Pimephales promelas. Sci. Total Environ. 806, 150920. (2022).ADS 
CAS 
PubMed 
Article 

Google Scholar 
DiGiacopo, D. G. & Hua, J. Evaluating the fitness consequences of plasticity in tolerance to pesticides. Ecol. Evol. 10, 4448–4456 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Carlson, B. E. & Langkilde, T. Body size variation in aquatic consumers causes pervasive community effects, independent of mean body size. Ecol. Evol. 7, 9978–9990 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Phung, T. X., Nascimento, J. C. S., Novarro, A. J. & Wiens, J. J. Correlated and decoupled evolution of adult and larval body size in frogs. Proc. Royal Soc. B 287, 20201474 (2020).Article 

Google Scholar 
Beasley, V. R. Direct and indirect effects of environmental contaminants on amphibians. In Reference Module in Earth Systems and Environmental Sciences https://doi.org/10.1016/b978-0-12-409548-9.11274-6 (Elsevier, 2020).Toledo, L. F., Sazima, I. & Haddad, C. F. B. Behavioural defences of anurans: An overview. Ethol. Ecol. Evol. 23, 1–25 (2011).Article 

Google Scholar 
Hartmann, M. T., Hartmann, P. A. & Haddad, C. F. B. Reproductive modes and fecundity of an assemblage of anuran amphibians in the Atlantic rainforest, Brazil. Inheringia 100, 207–215 (2010).Article 

Google Scholar 
Pupin, N. C., Gasparini, J. L., Bastos, R. P., Haddad, C. F. B. & Prado, C. P. A. Reproductive biology of an endemic Physalaemus of the Brazilian Atlantic forest, and the trade-off between clutch and egg size in terrestrial breeders of the P. signifer group. Herpetol. J. 20, 147–156 (2010).
Google Scholar 
Pereira, G. & Maneyro, R. Size-fecundity relationships and reproductive investment in females of Physalaemus riograndensis Milstead, 1960 (Anura, Leiuperidae) in Uruguay. Herpetol. J. 22, 145–150 (2012).
Google Scholar 
Tolledo, J., Silva, E. T., Nunes-de-Almeida, C. H. L. & Toledo, L. F. Anomalous tadpoles in a Brazilian oceanic archipelago: implications of oral anomalies on foraging behaviour, food intake and metamorphosis. Herpetol. J. 24, 237–243 (2014).
Google Scholar 
Annibale, F. S. et al. Smooth, striated, or rough: how substrate textures affect the feeding performance of tadpoles with different oral morphologies. Zoomorphology 139, 97–110 (2020).Article 

Google Scholar 
Venesky, M. D., Wassersug, R. J. & Parris, M. J. The impact of variation in labial tooth number on the feeding kinematics of tadpoles of southern leopard frog (Lithobates sphenocephalus). Copeia 3, 481–486 (2010).Article 

Google Scholar 
Venesky, M. D. et al. Comparative feeding kinematics of tropical hylid tadpoles. J. Exp. Biol. 216, 1928–1937 (2013).PubMed 

Google Scholar 
Jones, S. K. C., Munn, A. J., Penman, T. D. & Byrne, P. G. Long-term changes in food availability mediate the effects of temperature on growth, development and survival in striped marsh frog larvae: implications for captive breeding programmes. Conserv. Physiol. 3, cov029 (2015).Article 
CAS 

Google Scholar 
Bach, N. C., Natale, G. S., Somoza, G. M. & Ronco, A. E. Effect on the growth and development and induction of abnormalities by a glyphosate commercial formulation and its active ingredient during two developmental stages of the South-American Creole frog, Leptodactylus latrans. Environ. Sci. Pollut. Res. 23, 23959–23971 (2016).CAS 
Article 

Google Scholar 
Capellán, E. & Nicieza, A. G. Non-equivalence of growth arrest induced by predation risk or food limitation: context-dependent compensatory growth in anuran tadpoles. J. Anim. Ecol. 76, 1026–1035 (2007).PubMed 
Article 

Google Scholar 
Chin, A. M., Hill, D. R., Aurora, M. & Spence, J. R. Morphogenesis and maturation of the embryonic and postnatal intestine. Semin. Cell Dev. Biol. 66, 81–93 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Sun, Y., Zhang, J., Song, W. & Shan, A. Vitamin E alleviates phoxim-induced toxic effects on intestinal oxidative stress, barrier function, and morphological changes in rats. Environ. Sci. Pollut. Res. 25, 26682–26692 (2018).
Google Scholar 
Ouellet, M. Amphibian deformities: current state of knowledge. In Ecotoxicology of Amphibians and Reptiles (eds Sparling, D. W. et al.) 617–661 (Society of Environmental Toxicology and Chemistry, 2000).Hussein, M. & Singh, V. Effect on chick embryos development after exposure to neonicotinoid insecticide imidacloprid. J. Anat. Soc. India 65, 83–89 (2016).Article 

Google Scholar 
Crosby, E. B., Bailey, J. M., Oliveri, A. N. & Levin, E. D. Neurobehavioral impairments caused by developmental imidacloprid exposure in zebrafish. Neurotoxicol. Teratol. 49, 81–90 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lonare, M. et al. Evaluation of imidacloprid-induced neurotoxicity in male rats: A protective effect of curcumin. Neurochem. Int. 78, 122–129 (2014).CAS 
PubMed 
Article 

Google Scholar 
Žegura, B., Lah, T. T. & Filipič, M. The role of reactive oxygen species in microcystin-LR-induced DNA damage. Toxicology 200, 59–68 (2004).PubMed 
Article 
CAS 

Google Scholar 
Odetti, L. M., López González, E. C., Romito, M. L., Simoniello, M. F. & Poletta, G. L. Genotoxicity and oxidative stress in Caiman latirostris hatchlings exposed to pesticide formulations and their mixtures during incubation period. Ecotoxicol. Environ. Saf. 193, 110312 (2020).CAS 
PubMed 
Article 

Google Scholar 
Rutkoski, C. F. et al. Morphological and biochemical traits and mortality in Physalaemus gracilis (Anura: Leptodactylidae) tadpoles exposed to the insecticide chlorpyrifos. Chemosphere 250, 126162 (2020).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Herek, J. S. et al. Genotoxic effects of glyphosate on Physalaemus tadpoles. Environ. Toxicol. Pharmacol. 81, 103516 (2021).CAS 
PubMed 
Article 

Google Scholar 
Natale, G. S. et al. Lethal and sublethal effects of the pirimicarb-based formulation Aficida® on Boana pulchella (Duméril and Bibron, 1841) tadpoles (Anura, Hylidae). Ecotoxicol. Environ. Saf. 147, 471–479 (2018)
Google Scholar 
Gilbert, S. F. Developmental Biology, 8th edn. (Sinauer Associates, 2006).Soto, M., García-Santisteban, I., Krenning, L., Medema, R. H. & Raaijmakers, J. A. Chromosomes trapped in micronuclei are liable to segregation errors. J. Cell Sci. 131, 214742 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Crott, J. & Fenech, M. Preliminary study of the genotoxic potential of homocysteine in human lymphocytes in vitro. Mutagenesis 16, 213–217 (2001).CAS 
PubMed 
Article 

Google Scholar 
Benvindo-Souza, M. et al. Micronucleus test in tadpole erythrocytes: Trends in studies and new paths. Chemosphere 240, 124910 (2020).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Fenech, M. The in vitro micronucleus technique. Mutat. Res. 455, 81–95 (2000).CAS 
PubMed 
Article 

Google Scholar 
Podratz, J. L. et al. Drosophila melanogaster: A new model to study cisplatin-induced neurotoxicity. Neurobiol. Dis. 43, 330–337 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Iturburu, F. G. et al. Uptake, distribution in different tissues, and genotoxicity of imidacloprid in the freshwater fish Australoheros facetus. Environ. Toxicol. Chem. 36, 699–708 (2017).CAS 
PubMed 
Article 

Google Scholar 
Vieira, C. E. D., Pérez, M. R., Acayaba, R. D. A., Raimundo, C. C. M. & Martinez, C. B. R. DNA damage and oxidative stress induced by imidacloprid exposure in different tissues of the Neotropical fish Prochilodus lineatus. Chemosphere 195, 125–134 (2018).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Sanchéz-Bayo, F., Goka, K. & Hayasaka, D. Contamination of the aquatic environment with neonicotinoids and its implication for ecosystems. Front. Environ. Sci. 4, 71 (2016).Article 

Google Scholar 
Wood, T. & Goulson, D. The environmental risks of neonicotinoid pesticides: a review of the evidence post-2013. Environ. Sci. Pollut. Res. 24, 17285–17325 (2017).CAS 
Article 

Google Scholar 
Craddock, H. A., Huang, D., Turner, P.C., Quirós-Alcalá, L. & Payne-Sturges, D. C. Trends in neonicotinoid pesticide residues in food and water in the United States, 1999–2015. Environ. Health 18, 7 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Heyer, R. et al. Leptodactylus latrans. IUCN Red List https://doi.org/10.2305/IUCN.UK.2010-2.RLTS.T57151A11592655.en (2010).Ade, C. M., Boone, M. D. & Puglis, H. J. Effects of an insecticide and potential predators on green frogs and northern cricket frogs. J. Herpetol. 44, 591–600 (2010).Article 

Google Scholar 
Sarkar, M. A., Roy, S., Kole, R. K. & Chowdhury, A. Persistence and metabolism of imidacloprid in different soils of West Bengal. Pest Manag. Sci. 57, 598–602 (2001).CAS 
PubMed 
Article 

Google Scholar 
Goulson, D. Review: An overview of the environmental risks posed by neonicotinoid insecticides. J. Appl. Ecol. 50, 977–987 (2013).Article 

Google Scholar 
Mineau, P. Neonic insecticides and invertebrate species endangerment. In Reference Module in Earth Systems and Environmental Sciences https://doi.org/10.1016/B978-0-12-821139-7.00126-4 (2021).Yamamuro, M. et al. Neonicotinoids disrupt aquatic food webs and decrease fishery yields. Science 366, 620–623 (2019).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Gosner. K. L. A simplified table for staging anuran embryos and larvae with notes on identification. Herpetologica 16, 183–190 (1960).
Google Scholar 
Percie-du-Sert, N. et al. The ARRIVE guidelines 2.0: Updated guidelines for reporting animal research. PLoS Biol. 18, e3000410 (2020). CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Herkovits, J. & Pérez-Coll, C. S. AMPHITOX: A customized set of toxicity tests employing amphibian embryos. Symposium on multiple stressor effects in relation to declining amphibian populations. In Multiple Stressor Effects in Relation to Declining Amphibian Populations (eds Linder, G. et al.) 46–60 (ASTM International STP 1443, 2003).Merga, L. B. & Van den Brink, P. J. Ecological effects of imidacloprid on a tropical freshwater ecosystem and subsequent recovery dynamics. Sci. Total Environ. 784, 147167 (2021).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Bonmatin, J.-M. et al. Environmental fate and exposure; neonicotinoids and fipronil. Environ. Sci. Pollut. Res. 22, 35–67 (2015).CAS 
Article 

Google Scholar 
Sumon, K. A. et al. Effects of imidacloprid on the ecology of sub-tropical freshwater microcosms. Environ. Pollut. 236, 432–441 (2018).CONCEA – Conselho Nacional de Controle e Experimentação Animal. Resolução normativa Nº 25, 29 de setembro de 2015. Guia Brasileiro de Produção, Manutenção ou Utilização de Animais para Atividades de Ensino ou Pesquisa Científica do Conselho Nacional de Controle e Experimentação Animal. http://www.mctic.gov.br/mctic/export/sites/institucional/institucional/concea/arquivos/legislacao/resolucoes_normativas/Resolucao-Normativa-CONCEA-n-27-de-23.10.2015-D.O.U.-de-27.10.2015-Secao-I-Pag.-10.pdf. (2015).Rutkoski, C. F. et al. Lethal and sublethal effects of the herbicide atrazine in the early stages of development of Physalaemus gracilis (Anura: Leptodactylidae). Arch. Environ. Contam. Toxicol. 74, 587–593 (2018).CAS 
PubMed 
Article 

Google Scholar 
Pérez-Iglesias, J. M., Soloneski, S., Nikoloff, N., Natale, G. S. & Larramendy, M. L. Toxic and genotoxic effects of the imazethapyr-based herbicide formulation Pivot H® on montevideo tree frog Hypsiboas pulchellus tadpoles (Anura, Hylidae). Ecotoxicol. Environ. Saf. 119, 15–24 (2015).PubMed 
Article 
CAS 

Google Scholar 
Montalvão, M. F. et al. The genotoxicity and cytotoxicity of tannery effluent in bullfrog (Rana catesbeianus). Chemosphere 183, 491–502 (2017).ADS 
PubMed 
Article 
CAS 

Google Scholar  More

Remarks on the residual color patterns in the Kitadani Freshwater BivalvesResidual color patterns in the form of visible pigmentation on fossil molluscan shells are generally uncommon2,3. In the Paleozoic to Mesozoic fossil records, the color patterns were limited to marine species3, which are preserved as black to dark-colored bands running on the shell surface as melanin pigments20,21. The black to dark-colored stripes on the shells of the Kitadani Freshwater Bivalves resemble the color patterns in some extant freshwater bivalves, suggesting that the dark bands are residual color patterns remaining as melanin pigments. Consequently, the Kitadani Freshwater Bivalves represents the oldest and second fossil record of residual color patterns among fossil freshwater bivalves.The residual color patterns of the Kitadani Freshwater Bivalves resemble the color patterns of extant freshwater bivalves in terms of width, number, and distribution of the colored bands. Both the Kitadani Freshwater Bivalves and extant freshwater bivalves examined in this study consist of two types of color patterns: stripes along the growth lines and radial rays tapered toward the umbo. Notably, the former pattern is similar among all the species examined, as it forms in the peripheries of prominent growth lines occurring periodically. In the latter pattern, however, the morphology and distribution of the bands are slightly different between the Kitadani Freshwater Bivalves and the extant species. The Kitadani Freshwater Bivalves exhibits relatively distinct and wide radial rays running roughly parallel to the lengths of the sculpture elements (radial plications and/or wrinkles), while the extant species bear obscure and fine radial rays running diagonally to the lengths of the sculpture elements. Nonetheless, the taxa with V-shaped sculpture elements (wrinkles, ribs or arranged nodules) lack or bear ambiguous radial rays, whether extant (e.g., Triplodon spp., Indochinella spp. and Tritogonia spp.)13,15,22 or extinct (†Trigonioides tetoriensis).Hypothesis I: phylogenetic constraintsThe resemblance of the color patterns between the Kitadani Freshwater Bivalves and the extant unionids possibly resulted from the phylogenetic constrains. Each of the three species of the Kitadani Freshwater Bivalves belongs to a separate family (†Trigonioides tetoriensis: †Trigonioididae, †Plicatounio naktongensis: †Plicatounionidae, and †Matsuomtoina matsumotoi: †Pseudohyriidae) in the order Trigoniida19. Trigoniida in turn, forms the subclass Palaeoheterodonta with Unionida23. This raises a possibility that the color patterns observed in the Kitadani Freshwater Bivalves and the extant unionids is inherited from their most recent common ancestor. In other words, these color patterns, stripes along the growth lines and radial rays tapered toward the umbo, may be the apomorphy for Palaeoheterodonta. In fact, some extant trigoniid species belonging to Neotrigonia exhibit color pattern similar to those in the Kitadani Freshwater Bivalves and extant unionids in this study (e.g. Neotrigonia margaritacea)1.Interestingly, the coloration of color patterns is quite different between unioniids (green to blue colorings) and trigoniids (red to yellow colorings), and the oldest known color patterns of the Palaeoheterodonta (Myophorella nodulosa, a marine species of Trigoniida from the Oxfordian of the Early Jurassic) appears different (concentric rows of patches)10 from those of the Kitadani Freshwater Bivalves or the extant unioniids. These observations suggest that colorations evolved independently, in contrast to the color patterns, between Trigoniida and Unionida, and that Trigoniida more diverse color patterns than Unionida did in the Palaeoheterodont evolutionary history. Although further examination of the fossil record for the residual colors and color patterns in Palaeoheterodonta is essential, it is plausible that the habitat differences may have caused such discrepancy in the colorations and color patterns between Trigoniida (mainly marine) and Unionida (freshwater) in spite of the phylogenetic constrains.Hypothesis II: convergent evolutionThe other possible interpretation of the color pattern similarity between the Kitadani Freshwater Bivalves and extant Unionida, is the convergent evolution. One potential factor that may have caused this convergent evolution of the color patterns is an adaptation to their habitats. In general, much of the convergent evolution in animals occurs through the morphological evolution in response to their habitats24. Similarly in mollusks, shell colors and their patterns are generally influenced by their habitats2,6,25. Considering marine mollusks, the shell colors and their patterns have great diversity due to varying habitat environments, especially in coral reeves that exhibit various colors and complex ecosystem2,6. Conversely, in the freshwater ecosystem, the environmental colors are relatively monotonous with rocks, sand, mud, and green algae8, and such habitat conditions are likely indifferent between the Mesozoic and Cenozoic. As a result, the freshwater bivalves retained simple and monotonous color patterns for adapting to such environments throughout their evolution.Another conceivable factor to explain the convergent evolution in the color patterns of the studied freshwater bivalves is the selection pressure by visual predators. In general, the shell colors and their patterns in bivalves act as camouflages against the predators2,7,8,26,27,28. Previous studies have demonstrated that extant freshwater bivalves are preyed upon by crayfish, fish, birds, reptiles, and mammals29,30. Because shell colors in freshwater bivalves tend to be greenish, such colors may be an adaptation against visual predators for blending into the freshwater sediments on which abundant greenish phytoplanktons occur2,8. Therefore, the evolutionary conservatism in color patterns of freshwater bivalves may result from camouflages into freshwater microenvironments, which has been advantageous against visual predators since the late Early Cretaceous.The above discussion assumes that the visual predators of freshwater bivalves remained similar for at least 120 million years. Which animals could have been potential threads to the Kitadani Freshwater Bivalves, and, in turn, the Early Cretaceous freshwater bivalves? Among the extant visual predators of the freshwater bivalves, those whose lineages were present in the Early Cretaceous include crustaceans (especially brachyuran decapoda31), fish, lizards, turtles, crocodiles, birds, and mammals. Among them, the fossil record of durophagous lizards and mammals can be traced back only to the Late Cretaceous32,33. Conversely, lines of fossil evidence suggest that some fish34,35, turtles36, and crocodiles35 fed on molluscan invertebrates during the Early Cretaceous, and the Kitadani Freshwater Bivalves indeed occurs with abundant lepisosteiform scales, testudinate shells and crocodile teeth. Additionally, at least one Early Cretaceous avian species with crustacean gut contents can be attributed to the durophagous diet37, and the Kitadani Formation has yielded avialan skeletal remains38, and footprints39,40. Therefore, fish, turtles, crocodiles, and birds are likely candidates for visual predators of the Early Cretaceous freshwater bivalves, and have remained so until present. Additionally, while crustaceans have not been identified in the Kitadani Formation, they flourished in the Early Cretaceous and their remains occur with the fossil freshwater bivalves of the time elsewhere31. Thus, crustaceans may have also played a role as visual predators of the freshwater bivalves since the Early Cretaceous.In addition to the crustaceans, fishes, turtles, crocodiles and birds, the visual predators of the Early Cretaceous freshwater bivalves likely include extinct lineages. For example, some pliosauroid plesiosaurs are suggested as being durophagous34, although the freshwater members of the group are considered endemic41 and less likely to be a major thread to the Early Cretaceous freshwater bivalves. Another extinct candidate is non-avian dinosaurs. Ornithischians are suggested to have possessed a dietary flexibility including the durophagy. For instance, well-preserved hadrosaurid coprolites from the Late Cretaceous of Montana, U.S.A. include sizeable crustaceans and mollusks, possibly suggesting that the Cretaceous freshwater mollusks were consumed by these herbivorous dinosaurs42. In addition, some basal ceratopsian psittacosaurids are hypothesized for the durophagy based on the predicted large bite force in the caudal portion of the toothrow43. Among saurischians, some oviraptorosaurian theropods are indicated to consume mollusks with hard shells based on their mandibular features44. While hadrosaurids, psittacosaurids, and oviraptorosaurians have not been identified in the Kitadani Formation, psittacosaurids, and oviraptorosaurians are common elsewhere in the Early Cretaceous of East Asia45,46, and hadrosauroid Koshisaurus is present in the formation47. Because dinosaurs occupied a niche of large terrestrial predators throughout the Mesozoic, they may have acted as one of major mollusk consumers in absence of large lizards and mammals in the Early Cretaceous ecosystem. Thus, the predation pressure by visual predators to the freshwater bivalves in the Early Cretaceous is likely similar to that in the present. Consequently, one of evolutionary adaptations of the freshwater bivalves against such pressure has remained to camouflage in the phytoplankton-rich sediments, leading to the long-term evolutionary conservatism of their color patterns. More

Benchmark regression resultsParallel trend testThe premise of using DID is that the treatment group and control group meet the assumptions of parallel trend, which means that before ETS is officially implemented, the evolution trend of ecologicalization efficiency of industry of the control group and the experimental group is consistent and does not show a systematic difference. This study uses a more rigorous empirical test in parallel trend test: if the interaction coefficient is not significant and is different from zero before the implementation of ETS; and if the interaction coefficient is significant and is different from zero after the implementation of ETS, it indicates that there is no significant difference in ecologicalization efficiency of industry between the control group and the experimental group before the implementation of ETS. Results are shown in Table 4: before ETS was officially implemented, the difference coefficient was not significant; after the official implementation of ETS in 2013, the difference coefficient was significant and not equal to 0, and the ecologicalization efficiency of industry was improved significantly, which met the parallel trend of the DID. Therefore, it is scientific and reasonable to evaluate the effectiveness of ETS with DID.Table 4 Parallel trend test.Full size tableDynamic effect analysisTo compare the conditions of the experimental group and the control group before and after the implementation of ETS, dynamic graphs are drawn in this study, as shown in Fig. 1, which shows the impact of ETS on the regional ecologicalization efficiency of industry. The vertical line represents a 95% confidence interval and the broken line shows the marginal effect of regional ecologicalization efficiency, which means that the confidence interval contains is 0 before ETS’s implementation, and the result is not significant. In contrast, after 2013, the effect of ETS became apparent, the marginal effect gradually increased and the results became significant, perhaps owing to the implementation of ETS.Figure 1Dynamic analysis diagram.Full size imageThe effect of ETS on ecologicalization efficiency of industryControlling time effect and fixed effect, this study collected the data of pilot and non-pilot provinces of ETS from 2007 to 2019 to analyze the impact of ETS on the regional ecologicalization efficiency of industry and regional heterogeneity. The results are shown in Table 5. According to the results in the first column, ETS has significantly promoted the regional ecologicalization efficiency of industry, and the national implementation of ETS has achieved remarkable results. Compared with the regions that are not ETS pilot areas, the ecologicalization efficiency of industry of pilot provinces and cities has increased by 35%. Results also show that ETS has different effects on the ecologicalization efficiency of industry in different regions. Specifically, ETS significantly promoted regional ecologicalization efficiency of industry in the eastern and central regions, and the efficiency in the eastern region was more significant than that of the central region. However, the impact of ETS on the regional ecologicalization efficiency of industry in the western region was negative which may result from the fact that compared to the central and western regions, the east region has better economic development, advanced technology, and lots of talents that can respond to the implementation of ETS, accelerate the upgrade of industries, and improve the utilization level of regional resources. There are many traditional industries in the central and western regions, and the development of scientific and technological levels as well as the resource utilization efficiency there are relatively slow. Besides, it is difficult for the central and western regions to adapt to ETS in a short-term of time leading to the failure of improving the regional ecologicalization efficiency of industry in a short time.Table 5 Influence of ETS on ecologicalization efficiency of industry.Full size tableRobustness testPropensity matching score—double difference method (PSM-DID)The assumption of homogeneity and randomness between the control group and the experimental group is the premise of using the DID model. However, due to the large economic and regional differences among provinces and cities, there may be systematic differences between the experimental group and the control group, which may cause deviations in the results. Therefore, the data after propensity score matching is used in this study, making the matched individuals have no other significant differences unless they have been treated or not. The dual difference is conducted again to avoid self-selection bias, and the robustness of the above results is verified according to the measurement results. Control variables were used to match characteristic variables, the experimental group was matched with the control group, and the Logit model was adopted to delete the samples that fail to meet the matching criteria. After the matching, there are 168 observation values. The regression results of PSM-DID model show that, ETS has positive effects on the regional ecologicalization of industry (0.460***), which again proves that the conclusion that ETS improves regional ecologicalization of industry efficiency is reliable. The results are shown in Table 6.Table 6 The result of the PSM-DID.Full size tableCounterfactual testTo verify the robustness of the results again, six provinces and cities are randomly selected as experimental groups for multiple tests to construct new dummy variables of ETS, and the DID model was used again to verify the credibility of the above results. Four random samples were conducted in this study, and the results are shown in Table 7. It can be seen that the results are not significant, which also reversely proves that ETS improves the regional ecologicalization efficiency of industry.Table 7 Counterfactual test results.Full size tableActing pattern analysis of ETS on the regional ecologicalization efficiency of industryFirst, ETS may improve the regional ecologicalization efficiency of industry through industrial structure optimization and upgrading. Promoting upgrading of the industrial structure is one of the important approaches of social and economic development during the 14th Five-Year Plan formulation and is the only way to promote low-carbon and sustainable development of modern national industries. The upgrading of the industrial structure has been promoted to the national strategic level, contributing to the healthy development of the national economy system. ETS bring costs and benefits to enterprises, forcing them to transform and upgrade, increase investment in environmental protection and use clean energy, and accelerate the pace of energy conservation and emission reduction31. Second, ETS may improve the regional ecologicalization efficiency of industry through the coordinated agglomeration of resources. Marshall’s theory of scale economy, Krugman’s theory of new economic geography, Weber’s theory of agglomeration economy, Coase’s transaction cost theory, and so on reflect the importance of resource aggregation of economic activities through cost-saving, resource sharing, and other ways to improve industrial input–output efficiency, enhance industrial competitiveness, increase regional comprehensive strength and strengthen the competitive advantage of regional industrial clusters32. The benefits generated by resource aggregation far exceed the sum of benefits generated by various industries in the decentralized state. Under the pressure of ETS, enterprises may alleviate the mismatch between labor and capital through the collaborative aggregation of industrial resources, aiming to improve economic benefits and regional resource allocation efficiency and promote regional ecologicalization efficiency of industry. Third, ETS may improve the regional ecologicalization efficiency of industry by supporting ecological optimization. The sustainable development of the ecological environment is closely related to emission reduction policy. To alleviate the bad effects on the ecology, environmental protection is more and more brought to the attention of society and government. Policies for ecological protection have been introduced to reduce pollution20. All regions take effective and targeted measures to control environmental pollution and optimize the investment structure in light of their actual conditions. The purpose of ecological optimization is to improve the regional environment and strengthen pollution control which is one of the important parts of China’s fiscal spending. The government must guide the market to carry out ecological protection and environmental governance according to ETS. Studies have found that a low-carbon pilot policy helps to enhance the level of regional pollution control, promote the harmonious development of regional economy and environment, and then improve the regional ecologicalization efficiency of industry.To explore the transmission mechanism of ETS on the regional ecologicalization of industry efficiency, Baron and Kenny (1986)’s mediating effect model was referred to explore and verify whether there exists a structural optimization upgrade effect, resource synergistic agglomeration effect, ecological optimization support effect when ETC promotes regional ecologicalization efficiency of industry. Table 8 shows the regression results of the influence mechanism of ETS on the regional ecologicalization efficiency of industry. This study refers to the definition and research of industrial optimization and upgrading by Wang Qunwei, Huang Xianglan, and others, and the proportion of tertiary industry added value accounting for industrial added value is selected to measure the effectiveness of industrial optimization and upgrading. For resource synergistic agglomeration effect, this study refers to the calculation methods of Cui Shuhui, Chen Jianjun et al. and adopts the collaborative aggregation index of manufacturing and producer services to measure the collaborative aggregation effect of resources, which effectively avoids the scale difference between different regions. It can be seen from the table that the implementation of ETS has significantly influenced the three effects proposed by this study: the optimization and upgrading effect of industrial structure, the synergistic aggregation effect of resources, and the support effect of ecological optimization. In addition, ETS has a positive and significant impact on the regional ecologicalization efficiency of industry. The results in Columns 3, 5, and 7 of the table show the industrial optimization and upgrading effect, resource synergistic aggregation effect, structural upgrading effect, and resource allocation effect generated in the process of low-carbon pilot policy operation can significantly promote regional ecologicalization efficiency of industry and have an obvious intermediary effect. The mediating effect produced by industrial structure optimization and upgrading is about 0.042, the mediating effect produced by resource synergy agglomeration is about 0.148, and the mediating effect produced by ecological optimization support is about 0.166. According to the Sobal test results, all of them have passed the test, indicating that the above results are reliable.Table 8 Mediating effect test results.Full size table More