Ecology

Wang, Z. & Yuan, J. Typhoon numerical simulation and visualization for disaster risk assessment. Bull. Surv. Mapp. 108–110, 132 (2015).ADS 

Google Scholar 
Tang, J., Xu, L., Li, Y., He, Y. & Cui, S. Assessing the damage caused by typhoon on urban green space ecosystme service based on UAV remote sensing. J. Nat. Disasters 27, 153–161 (2018).
Google Scholar 
Tian, Y., Zhou, W., Qian, Y., Zheng, Z. & Pan, X. The influence of Typhoon Mangkhut on urban green space and biomass in Shenzhen, China. Acta Ecol. Sin. 40, 2589–2598 (2020).
Google Scholar 
Lin, G. Major anti-wind and alkali-resisting landscape plants of south China’s seaside region. For. Invertory Plan. 29, 78–81 (2004).
Google Scholar 
Lin, T. C., Hogan, J. A. & Chang, C. T. Tropical cyclone ecology: a scale-link perspective. Trends Ecol. Evol. 35, 594–604 (2020).Article 

Google Scholar 
Lin, H. et al. Risk assessments of trees and urban spaces under attacks of typhoons: a survey and statistical study in Guangzhou, China. IOP Conf. Ser. Earth Environ. Sci. 588, 032055 (2020).Article 

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

Google Scholar 
Zhu, J., Matsuzaki, T. & Sakioka, K. Wind speeds within a single crown of Japanese black pine (Pinus thunbergii Parl.). For. Ecol. Manag. 135, 19–31 (2000).Article 

Google Scholar 
Ver Planck, N. R. & MacFarlane, D. W. Branch mass allocation increases wind throw risk for Fagus grandifolia. For. Int. J. For. Res. 92, 490–499 (2019).
Google Scholar 
Dunham, R. A. & Cameron, A. D. Crown, stem and wood properties of wind-damaged and undamaged Sitka spruce. For. Ecol. Manag. 135, 73–81 (2000).Article 

Google Scholar 
Burcham, D. C., Autio, W. R., Modarres-Sadeghi, Y. & Kane, B. After pruning, wind-induced bending moments and vibration decrease more on reduced than raised Senegal mahogany (Khaya senegalensis). Urban For. Urban Green. 61, 127100 (2021).Article 

Google Scholar 
Gilman, E. F., Masters, F. & Grabosky, J. C. Pruning affects tree movement in hurricane force wind. Arboric. Urban For. 34, 20–28 (2008).Article 

Google Scholar 
Päätalo, M. L., Peltola, H. & Kellomäki, S. Modelling the risk of snow damage to forests under short-term snow loading. For. Ecol. Manag. 116, 51–70 (1999).Article 

Google Scholar 
North, E., Johnson, G., Murphy, R., Giblin, C. & Rendahl, A. Boulevard tree failures during wind loading events. Arboric. Urban For. 45, 259–269 (2019).
Google Scholar 
Angelou, N., Dellwik, E. & Mann, J. Wind load estimation on an open-grown European oak tree. For. Int. J. For. Res. 92, 381–392 (2019).
Google Scholar 
Kontogianni, A., Tsitsoni, T. & Goudelis, G. An index based on silvicultural knowledge for tree stability assessment and improved ecological function in urban ecosystems. Ecol. Eng. 37, 914–919 (2011).Article 

Google Scholar 
He, D. & Li, Z. Wind tunnel test on wind- induced responses of roadside trees. J. Nat. Disasters 28, 44–53 (2019).
Google Scholar 
Cao, J., Tamura, Y. & Yoshida, A. Wind tunnel study on aerodynamic characteristics of shrubby specimens of three tree species. Urban For. Urban Green. 11, 465–476 (2012).Article 

Google Scholar 
Zhang, W., Kang, L., Zhang, Q., Li, C. & Zou, X. Speed upwind and downwind of a single plant. J. Beijing Norm. Univ. Nat. Sci. 56, 573–581 (2020).
Google Scholar 
Cheng, H. et al. Wind tunnel study of airflow recovery on the lee side of single plants. Agric. For. Meteorol. 263, 362–372 (2018).Article 
ADS 

Google Scholar 
Ma, S. et al. Experimental research of viscous flow around a Nitraria tangutorum boscage. Res. Soil Water Conserv. 13, 147–149 (2006).ADS 

Google Scholar 
Rahman, M. et al. Disentangling the role of competition, light interception, and functional traits in tree growth rate variation in South Asian tropical moist forests. For. Ecol. Manag. 483, 118908 (2021).Article 

Google Scholar 
Forrester, D. I., Benneter, A., Bouriaud, O. & Bauhus, J. Diversity and competition influence tree allometric relationships—Developing functions for mixed-species forests. J. Ecol. 105, 761–774 (2017).Article 

Google Scholar 
Coombes, A., Martin, J. & Slater, D. Defining the allometry of stem and crown diameter of urban trees. Urban For. Urban Green. 44, 1–15 (2019).Article 

Google Scholar 
Stoffel, M. Mechanical stability and growth performance of trees. (2009).Lento, M., Thijs, D., Jonas, A., Dominique, D. & Jan, C. Comparative study of flow field and drag coefficient of model and small natural trees in a wind tunnel. Urban For. Urban Green. 35, 230–239 (2018).Article 

Google Scholar 
West, G. B., Brown, J. H. & Enquist, B. J. A general model for the structure and allometry of plant vascular systems. Nature 400, 664–667 (1999).Article 
ADS 
CAS 

Google Scholar 
Enquist, B. J. Cope’s rule and the evolution of long-distance transport in vascular plants: Allometric scaling, biomass partitioning and optimization. Plant Cell Environ. 26, 151–161 (2003).Article 

Google Scholar 
Bentley, L. P. et al. An empirical assessment of tree branching networks and implications for plant allometric scaling models. Ecol. Lett. 16, 1069–1078 (2013).Article 

Google Scholar 
Eloy, C., Fournier, M., Lacointe, A. & Moulia, B. Wind loads and competition for light sculpt trees into self-similar structures. Nat. Commun. 8, 1–11 (2017).Article 
CAS 

Google Scholar 
Lin, M. Y. & Khlystov, A. Investigation of ultrafine particle deposition to vegetation branches in a wind tunnel. Aerosol Sci. Technol. 46, 465–472 (2012).Article 
ADS 
CAS 

Google Scholar 
Ji, W. & Zhao, B. A wind tunnel study on the effect of trees on PM2.5 distribution around buildings. J. Hazard. Mater. 346, 36–41 (2018).Article 
CAS 

Google Scholar 
Hui, K. K. W. et al. Unveiling falling urban trees before and during Typhoon Higos (2020): Empirical case study of potential structural failure using tilt sensor. Forests 13, 359 (2022).Article 

Google Scholar 
Liao, S., Huang, M., Lou, W., Lin, W. & Xiao, Z. Numerical simulation of a urban wind field under the influence of typhoon “Mangkhut”. Acta Aerodyn. Sin. 39, 107–116 (2021).
Google Scholar 
Gardiner, B., Berry, P. & Moulia, B. Review: Wind impacts on plant growth, mechanics and damage. Plant Sci. 245, 94–118 (2016).Article 
CAS 

Google Scholar 
Lüttge, U. & Buckeridge, M. Trees: structure and function and the challenges of urbanization. Trees Struct. Funct. https://doi.org/10.1007/s00468-020-01964-1 (2020).Article 

Google Scholar 
Hwang, H. M., Fiala, M. J., Park, D. & Wade, T. L. Review of pollutants in urban road dust and stormwater runoff: part 1. Heavy metals released from vehicles. Int. J. Urban Sci. 20, 334–360 (2016).Article 

Google Scholar 
Kuitert, W. The nature of urban Seoul: Potential vegetation derived from the soil map. Int. J. Urban Sci. 17, 95–108 (2013).Article 

Google Scholar 
Stubbs, C. J., Cook, D. D. & Niklas, K. J. A general review of the biomechanics of root anchorage. J. Exp. Bot. 70, 3439–3451 (2019).Article 
CAS 

Google Scholar 
Sterck, F. J. & Bongers, F. Ontogenetic changes in size, allometry, and mechanical design of tropical rain forest trees. Am. J. Bot. 85, 266–272 (1998).Article 
CAS 

Google Scholar 
Sim, V. & Jung, W. Wind fragility for urban street tree in Korea. J. Wetl. Res. 21, 298–304 (2019).
Google Scholar 
Ueda, M. & Shibata, E. Why do trees decline or dieback after a strong wind? Water status of Hinoki cypress standing after a typhoon. Tree Physiol. 24, 701–706 (2004).Article 

Google Scholar 
Jalkanen, A. & Mattila, U. Logistic regression models for wind and snow damage in northern Finland based on the National Forest Inventory data. For. Ecol. Manag. 135, 315–330 (2000).Article 

Google Scholar 
Hale, S. E., Gardiner, B. A., Wellpott, A., Nicoll, B. C. & Achim, A. Wind loading of trees: Influence of tree size and competition. Eur. J. For. Res. 131, 203–217 (2012).Article 

Google Scholar 
Olson, M. E., Aguirre-Hernández, R. & Rosell, J. A. Universal foliage-stem scaling across environments and species in dicot trees: Plasticity, biomechanics and Corner’s Rules. Ecol. Lett. 12, 210–219 (2009).Article 

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

Google Scholar 
King, D. A., Davies, S. J., Tan, S. & Noor, N. S. M. The role of wood density and stem support costs in the growth and mortality of tropical trees. J. Ecol. 94, 670–680 (2006).Article 

Google Scholar 
Petty, J. A. & Worrell, R. Stability of coniferous tree stems in relation to damage by snow. Forestry 54, 115–128 (1981).Article 

Google Scholar 
Gilman, E. F. & Lilly, S. Tree pruning (International Society of Arboriculture, 2002).
Google Scholar 
Zhang, J., Gou, Z., Zhang, F. & Shutter, L. A study of tree crown characteristics and their cooling effects in a subtropical city of Australia. Ecol. Eng. 158, 106027 (2020).Article 

Google Scholar 
Marchi, L. et al. State of the art on the use of trees as supports and anchors in forest operations. Forests 9, 467 (2018).Article 

Google Scholar  More

All statistical analyses were performed through R version 4.1.050. Besides the explicitly mentioned packages, the R packages cowplot51, data.table52, dplyr53, ggplot254, itsadug55, purrr56, raster57, sf58, sfheaders59, tibble60 and tidyr61 were key for data handling, data analysis, and plotting. Posterior distributions were summarised through means and credible intervals (CIs). CIs are the highest density intervals, calculated through the package bayestestR62. To summarise multiple posterior distributions, 5000 Monte Carlo simulations were used.Study regionThe study included data from the whole of Switzerland. As an observation unit for records, we chose 1 × 1 km squares (henceforth squares). Switzerland covers 41,285 km2, spanning a large gradient in elevation, climate and land use. It can be divided into five coarse biogeographic regions based on floristic and faunistic distributions and on institutional borders of municipalities63 (Fig. 1b). The Jura is a mountainous but agricultural landscape in the northwest (~4200 km2, 300–1600 m asl; annual mean temperature: ~9.4 °C, annual precipitation: ~1100 mm (gridded climate data here and in the following from MeteoSwiss (https://www.meteoswiss.admin.ch), average 1980–2020, at sites ~500 m asl.)). The Jura is separated from the Alps by the Plateau, which is the lowland region spanning from the southwest to the northeast (~11,300 km2, 250–1400 m asl, mostly below 1000 m asl; ~9.5 °C, ~1100 mm). It is the most densely populated region with most intensive agricultural use. For the Alps, three regions can be distinguished. The Northern Alps (~10,700 km2, 350–4000 m asl; ~9.2 °C, ~1400 mm), which entail the area from the lower Prealps, which are rather densely populated and largely used agriculturally, up to the northern alpine mountain range. The Central Alps (~11,300 km2, 450–4600 m asl; ~9.5 °C, ~800 mm) comprise of the highest mountain ranges in Switzerland and the inner alpine valleys characterised by more continental climate (i.e., lower precipitation). Intensive agricultural land use is concentrated in the lower elevations and agriculture in higher elevations is mostly restricted to grassland areas used for summer grazing. The Southern Alps (~3800 km2, 200–3800 m asl; ~10.4 °C, 1700 mm) range from the southern alpine mountain range down to the lowest elevations of Switzerland and are clearly distinguished from the other regions climatically, as they are influenced by Mediterranean climate, resulting in, e.g., milder winters. Besides differences between biogeographic regions, climate, land use and changes therein vary greatly between different elevations (Supplementary Fig. S9). To account for these differences, we split the regions in two elevation classes at the level of squares. All squares with a mean elevation of less than 1000 m asl were assigned to the low elevation, whereas squares above 1000 m asl were assigned to the high elevation (no squares in the Plateau fell in the high elevation). This resulted in nine bioclimatic zones (Fig. 1b), for which separate species trends were estimated in the subsequent analyses. The threshold of 1000 m asl enabled a meaningful distinction based on the studied drivers (climate and land-use change) and was also determined by the availability of records data (high coverage in all nine bioclimatic zones).Species detection dataWe extracted records of butterflies (refers here to Papilionoidea as well as Zygaenidae moths), grasshoppers (refers here to all Orthoptera) and dragonflies (refers here to all Odonata) from the database curated by info fauna (The Swiss Faunistic Records Centre; metadata available from the GBIF database at https://doi.org/10.15468/atyl1j, https://doi.org/10.15468/bcthst, https://doi.org/10.15468/fcxtjg). This database unites faunistic records made in Switzerland from various sources including both records by private persons and from projects such as research projects, Red-List inventories or checks of revitalisation measures. Only records with a sufficient precision, both temporally (day of recording) and spatially (place of recording known to the precision of 1 km2 or less), were used for analyses. Besides temporal and spatial information, information on the observer and the project (if any) was obtained for each record. All records made by a person/project on a day in a square were attributed to one visit, which was later used as replication unit to model the observation process (see below).We included records from the focal time range 1980–2020. Additionally, we included records from 1970–1979 for butterflies in occupancy-detection models to increase the robustness of mean occupancy estimates. We excluded the mean occupancy estimates for these additional years from further analyses to cover the same period for all groups. Prior to analyses, following the approach in ref. 26, we excluded observations of non-adult stages and observations from squares that only were visited in 1 year of the studied period, because these would not contain any information on change between years64. This resulted in 18,018 squares (15,248 for butterflies, 9870 for grasshoppers, 5188 for dragonflies) and 1,448,134 records (879,207 butterflies, 272,863 grasshoppers, 296,064 dragonflies) that we included in the analyses (Supplementary Fig. S2). The three datasets for the different groups were treated separately for occupancy-detection modelling, following the same procedures for all three groups. To determine detections and non-detections for each species and visit, which could then be used for occupancy-detection modelling, we only included visits that (a) did not originate from a project, which had a restricted taxonomic focus not including the focal species, (b) were not below the 5% quantile or above the 95% quantile of the day of the year at which the focal species has been recorded26 and (c) were from a bioclimatic zone, from which the focal species was recorded at least once.Occupancy-detection modelsWe used occupancy-detection models65,66 to estimate annual mean occupancy of squares for the whole of Switzerland and for the nine bioclimatic zones for each species (i.e., mean number of squares occupied by a species), mostly following the approach in ref. 26. We fitted a separate model for each species, based on different datasets for the three groups. We included only species that were recorded in any square in at least 25% of all analysed years. Occupancy-detection models are hierarchical models in which two interconnected processes are modelled jointly, one of which describes occurrence probability (ecological process; used to infer mean occupancy), whereas the other describes detection probability (observation process)65. The two processes are modelled through logistic regression models. The occupancy model estimates occurrence probability for all square and year combinations, whereas the observation model estimates the probability that a species has been detected by an observer during a visit. More formally, each square i in the year t has the latent occupancy status zi,t, which may be either 1 (present) or 0 (absent). zi,t depends on the occurrence probability ψi,t as follows$${z}_{i,t}sim {{{mbox{Bern}}}}left({psi }_{i,t}right)$$
(1)
The occupancy status is linked to the detection/non-detection data yi,t,j at square i in year t at visit j as$${y}_{i,t,, j}{{|}}{z}_{i,t}sim {mathrm {Bern}}({z}_{i,t}{p}_{i,t,j})$$
(2)
where pi,t,j is the detection probability.The regression model for occurrence probability (occupancy model) looked as follows$${{mbox{logit}}}({psi }_{i,t})={mu }_{o}+{beta }_{o1}{{{{{rm{elevatio}}}}}}{{{{{{rm{n}}}}}}}_{i}+{beta }_{o2}{{{{{rm{elevatio}}}}}}{{{{{{rm{n}}}}}}}_{i}^{2}+{alpha }_{o1,i}+{alpha }_{o2,i}+{gamma }_{r(i),t}$$
(3)
with μo being the global intercept, elevationi being the scaled elevation above sea level and αo1,i, αo2,i and γr(i),t being the random effects for fine biogeographic region (12 levels, Supplementary Fig. S10; these were again defined based on floristic and faunistic distributions and followed institutional borders63), square and year. The random effects for fine biogeographic region and square were modelled as follows:$${alpha }_{o1}sim {{{{{rm{Normal}}}}}}left(0,{sigma }_{o1}right)$$
(4)
and$${alpha }_{o2}sim {{{{{rm{Normal}}}}}}left(0,{sigma }_{o2}right)$$
(5)
The random effect of the year was implemented with separate random walks per zone following ref. 67, which allowed the effect to vary between the nine bioclimatic zones, while accounting for dependencies among consecutive years. Conceptually, in random walks, the effect of 1 year is dependent on the previous year’s effect, resulting in trajectories with less sudden changes between consecutive years. This was implemented as follows:$${gamma }_{r,t}sim left{begin{array}{c}{{{{{rm{Normal}}}}}}left(0,{1.5}^{2}right){{{{rm{for}}}}},t=1\ {{{{{rm{Normal}}}}}}left({gamma }_{r,t-1},{sigma }_{gamma r}^{2}right){{{{rm{for}}}}},t , > ,1end{array}right.$$
(6)
with$${sigma }_{gamma r}sim {{mbox{Cauchy}}}left(0,1right)$$
(7)
The regression model for detection probability (observation model) looked as follows$${{{{rm{logit}}}}}({p}_{i,t,j}) =, {mu }_{d}+{beta }_{d1}{{{{{rm{yda}}}}}}{{{{{{rm{y}}}}}}}_{j}+{beta }_{d2}{{{{{rm{yda}}}}}}{{{{{{rm{y}}}}}}}_{j}^{2}+{beta }_{d3}{{{{{rm{shortlis}}}}}}{{{{{{rm{t}}}}}}}_{j}+{beta }_{d4}{{{{{rm{longlis}}}}}}{{{{{{rm{t}}}}}}}_{j} \ quad+ {beta }_{d5}{{{{{rm{exper}}}}}}{{{{{{rm{t}}}}}}}_{j}+{beta }_{d6}{{{{{rm{projec}}}}}}{{{{{{rm{t}}}}}}}_{j}+{beta }_{d7}{{{{{rm{targeted}}}}}}_{{{{{rm{projec}}}}}}{{{{{{rm{t}}}}}}}_{j} \ quad+ {beta }_{d8}{{{{{rm{redlis}}}}}}{{{{{{rm{t}}}}}}}_{j}+{alpha }_{d1,t}$$
(8)
where μd is the global intercept, ydayj is the scaled day of the year of visit j, shortlistj and longlistj are dummies of a three-level factor denoting the number of species recorded during the visit (1; 2–3; >3), and expertj, projectj, targeted_projectj and redlistj are dummies of a five-level factor denoting the source of the data. The source might either be a common naturalist observation (reference level), an observation by an expert naturalist, an observation made during a not further specified project, an observation made in a project targeted at the focal species or an observation made during a Red-List inventory. An expert naturalist was defined as an observer that contributed a significant number of records, which was defined as the upper 2.5% quantile of all observers arranged by their total number of records, and that made at least one visit with an exceptionally long species list, which was defined as a visit in the upper 2.5% quantile of all visits arranged by the number of records. The proportions of records originating from these different sources changed across years, but change was not unidirectional and differed among the investigated groups (Supplementary Fig. S11), such that accounting for data source in the model should suffice to yield reliable estimates of occupancy trends. αd1,t is a random effect for year, which was modelled as$${alpha }_{d1}sim {{{{{rm{Normal}}}}}}left(0,{sigma }_{d1}right)$$
(9)
The occupancy and observation models were fitted jointly in Stan through the interface CmdStanR68. Four Markov chain Monte Carlo chains with 2000 iterations each, including 1000 warm-up iterations, were used. Priors of the model parameters are specified in Supplementary Table S5. For the prior distribution of global intercepts, a standard deviation of 1.5 was chosen to not overweight the extreme values on the probability scale. To ensure that chains mixed well, Rhat statistics for annual mean occupancy estimates were calculated through the package rstan69. For Switzerland-wide annual estimates (n = 18,140), 98.0% of values met the standard threshold of 1.1 (99.9% of values More

Tularemia affects animal welfare, human health, and the environment and is thus better approached from a one-health perspective27. Several studies in the Northern hemisphere28, and more recently in Australia15,16, have provided a vital research track in the epidemiology of this disease. In contrast, studies in Africa are too limited and scarce. The aim of this study was to investigate the presence of tularemia in wild leporids collected in Northern Algeria. These animals are highly susceptible to F. tularensis infection and considered sentinel hosts for surveillance of tularemia. The strategy we used to detect F. tularensis in leporids mainly used molecular, histological and immunohistochemical analyzes of tissues taken from animals found dead or hunted. To the best of our knowledge, detection of F. tularensis by PCR or culture has not been previously reported in wild leporidae in Algeria or other African countries.Animal tissue samples were tested using three qPCR assays of variable sensitivity and specificity. The Type B-qPCR test targets a specific junction between ISFtu2 and a flanking 3′ region, which is considered specific for F. tularensis subsp. holarctica26, the only tularemia agent found in Europe and Asia. The Tul4-qPCR assay targets a simple copy gene encoding a surface protein, which can be found in the genome of all F. tularensis subspecies causing tularemia and that of the aquatic bacterium F. novicida. Because F. novicida has never been isolated from lagomorphs or other animal species, and very rarely from human29, a positive Tul4 qPCR for the studied tissue samples likely indicated the presence of F. tularensis DNA. The ISFtu2 qPCR is considered highly sensitive because multiple copies of this insertion sequence are found in the F. tularensis genome. However, it lacks specificity because ISFtu2 is also found in many other Francisella species25.Two animals were considered “probable” tularemia cases because some of their samples were positive for the three qPCR tests. Ten animals were considered “possible” tularemia cases because their samples were positive for the ISFtu2 and Tul4 qPCRs but not the Type B qPCR. Finally 19 leporids were “uncertain” cases because only samples positive for the ISFtu2 qPCR were found. For the remaining 43 animals, all the tested samples were negative for the three qPCRs. Overall, we detected F. tularensis DNA-positive samples in 12/74 (16.21%) leporids, which strongly suggest that tularemia is present in the lagomorph population of the study area. The positive Type B qPCR tests in two animals suggested that F. tularensis subsp. holarctica could be the involved subspecies. We did not confirm these data by isolating F. tularensis from the studied leporids. However, the isolation of this pathogen from human or animal samples is tedious and has a low sensitivity13. Moreover, most of our samples were not appropriate for F. tularensis culture because of their long-term preservation in ethanol 70° or 10% formalin. Further study using fresh (non-fixed) tissue samples from dead leporids collected in the same study area is needed to definitively confirm the presence of tularemia in these animals and characterize the F. tularensis subspecies and genotypes involved.Although PCR is usually more sensitive than culture for detecting F. tularensis, it also has some limitations. Firstly, the DNA extraction from organs preserved in ethanol for several months was difficult although easier for spleen than for liver samples. Some tissue samples could be lysed only after overnight incubation with proteinase K. Secondly, tissue samples contained PCR inhibitors as demonstrated by better DNA amplification from some samples after their dilution in PCR grade water. To reduce the effect of PCR inhibitors, organ samples with negative qPCR were retested using Bovine Serum Albumin (BSA) and the Real-time PCR system TaqMan (Applied Biosystems, Munich, Germany)30. Finally, DNA regions to be amplified were optimized to obtain high sensitivity and specificity of qPCR tests.IHC detection of F. tularensis in formalin-fixed tissue can be helpful for tularemia diagnosis31,32. For one possible tularemia case, F. tularensis could be detected on immunohistochemical (IHC) examination of a liver sample using a specific anti-F. tularensis antibody. The intensity and localization of positive staining were comparable to those previously recorded for other animals32,33. IHC did not provide interpretable findings for four other tested specimens. Such negative results might be explained by an inhomogeneous distribution of infectious foci in the involved organs as well as a low bacterial inoculum in infected tissues. This has been previously demonstrated in tularemia granulomatous lesions in cell types like epithelial cells of the kidney, testis, and epididymis, hepatocytes, and bronchiolar epithelial cells31. Besides, IHC is a delicate technology whose results are highly dependent on the quality and fixation time of the organ tissues34. IHC analysis of dead animal tissues remains challenging, especially in case of tissue necrosis34.In our limited case series we found a F. tularensis infection prevalence in leporids of 2.7% (2/74) for probable tularemia cases and 16.2% (12/74) when considering both probable and posible cases. We cannot make a guess about the prevalence of tularemia because our series is not representative of the general lagomorph population in the study area. In Germany, F. tularensis DNA was detected in 1.1% of European Brown hares and 2.4% of wild rabbits collected between 2009 and 201435. Higher infection rates were reported in the same country, including 11.8% (100/848 animals) in hares collcted in the North Rhine-Westphalia region36 and 30% (55/179) in brown hares collected between 2010 and 2016 in Baden-Wuerttemberg37. In Hungary, the prevalence of tularemia in hares was evaluated at 4.9–5.3%38. In Portugal, prevalences of 4.3% and 6.3% were reported in brown hares and wild rabbits, respectively39. However, the comparison of the reported tularemia prevalences in leporids is irrelevant because studies involved different animal species and geographic areas, and used different methods for F. tularensis detection.Two possibilities could explain the lack of detection of tularemia in Algeria before this study. The first hypothesis is that this disease was not searched for in previous years, while it could have been present in this country for decades. The second hypothesis is that tularemia was recently imported in Algeria. Migratory birds may have been involved in the long-distance spread of F. tularensis40. These hosts can be infested by ectoparasites such as ticks which are the primary vectors of tularemia41,42. They can also spread the bacteria in the hydro-telluric environment through their secretions and feces18,43,44. An alternative possibility is that F. tularensis-infected animals (especially game animals) have been imported in Algeria from endemic countries. Whatever the mode of introduction of tularemia in Algeria, the dissemination of this disease over time might have been facilitated by the ability of F. tularensis to infect multiple hosts and its better survival in a cool environment45, which characterizes Northern Algeria climate. The emergence or re-emergence of tularemia in other countries has been related to climate change, human-mediated movement of infected animals, and wartime resulting in a significant rise of F. tularensis infections in the rodent populations39,46.In our study, infected animals were collected throughout 4 years, although more frequently in autumn. Probable and possible tularemia cases were mainly collected during the hunting season (i.e., September, October, November, and December). Animals could not be collected in February because of heavy rains and in May and June because it corresponds to female leporids’ lactation period. In most endemic countries, tularemia cases are typically more frequent in late spring, the summer months, and early autumn37,47,48,49,50. Occasionally, fatal tularemia cases in hares have been predominantly reported during the cold season11,51. The climatic conditions can affect tularemia outbreaks in animals, depending on the reservoir involved and the predominant modes of infection52.We detected tularemia more frequently in female than in male hares, and the reverse was true for wild rabbits. The prevalence of tularemia in male or female lagomorphs varies between studies. In Sweden, Morener et al.50 reported a tularemia case series only involving male hares. In the same country, Borg et al.50 observed an overrepresentation of females in the epizootic of 1967. They suggested that, compared to males, females had a higher risk of exposure to infected mosquitoes or were more vulnerable to tularemia because they were pregnant or had just given birth to a litter50. Tularemia was found in a few juveline leporids, which might be explained by a shorter exposure time to F. tularensis, a higher death rates due to higher susceptibility to F. tularensis infection or easier predation by their natural enemies, or more frequent hunting of adults compared to the juveniles53.Tularemia is usually more frequently detected in leporids found dead than in hunted animals. As an example, a German study reported a higher prevalence of tularemia in hares found dead (2.9%) than in hunted ones (0.7%)35. In our study, most qPCR-positive animals were hunted. Our study might not be representative of the prevalence of tularemia in either population because most collected animals had been hunted.The incubation period and clinical presentation of tularemia in leporids vary according to the species considered. Tularemia is typically an acute disease in mountain hares (Lepus timidus) in Scandinavia and has a chronic pattern in European brown hares (Lepus europaeus) in Central Europe50. The incubation time and clinical presentation of tularemia can be different in Cape hares (Lepus capensis). Wild rabbits are less sensitive to F. tularensis infection than hares31,39,54. An extended incubation period and chronic evolution of tularemia would facilitate the detection of F. tularensis in infected animals. In our study, a similar tularemia prevalence was found in the Cape hares and wild rabbits, which might reflect exposure to a same biotope area and environmental reservoirs of F. tularensis.The pathological lesions of tularemiia in leporids can vary according to the F. tularensis strain involved, the mode and route of infection, and the susceptibility and immune status of the host32,50. In the European brown hares, granulomas with central necrosis have been reported in the lungs and kidneys and occasionally in the liver, spleen, bone marrow, and lymph nodes50. In contrast, only acute necrosis in the liver, spleen, bone marrow, and lymph nodes have been found in Lepus timudus hares in Sweden50. The lesions in the Japanese hare (Lepus brachyurus angustidens) are comparable to those of Lepus timidus, except for cutaneous, lung, brain, and adrenal gland lesions32. In the European rabbit, Oryctolagus cuniculus, tularemia is not associated with identifiable macroscopic tissue lesions39,55. To our knowledge, no reports describing post-mortem lesions in Cape hares with tularemia are available. In this study, similar lesions were found in hares and wild rabbits except necrotic foci only observed in some wild rabbit organs (such as liver, lungs, kidney, ovary). Most animals had pathological lesions of pneumonia, gastritis and enteritis. Kidney lesions and adrenal glands enlargment were oberved. Necrotic lesions were occasionally found in the lungs, liver, spleen and ovary and hemorrhages in the lungs, liver, and intestines.Tularemia is an arthropod-born disease in most endemic areas14,22,28. In our study, 50% of positive leporids were infested by known tularemia vectors such as ticks (Ixodes ricinus56,57, Rhipicephalus sanguineus39), fleas (Spillopsylus cuniculi58), and lice of lagomorphs (Haemodipsus lepori and Haemodipsus setoni59,60). Ticks are the most significant arthropod vectors of tularemia61. Ticks are frequently involved in the transmission of tularemia in North America, including Dermacentor andersoni, D. variabilis, and Amblyomma americanum57,62,63. In Europe, tick-borne tularemia represents 13% to 26% of human cases57,64. The involved species include D. marginatus, D. reticulatus, I. ricinus, R. sanguineus, and Haemaphysalis concinna65,66. Further research on wild leporid sucking arthropods is needed to confirm the presence and clarify the ecology of F. tularensis in Algeria.Our study reports for the first time the detection of F. tularensis DNA in leporids from Northern Algeria. The markers most in favor of tularemia in the animals studied are the positivity of qPCR tests, in particular, the “type B” qPCR test which amplifies a specific DNA sequence of F. tularensis subsp. holarctica, and a positive immunohistological examination in one animal. Further investigation is needed to confirm our results by the isolation of this pathogen from animal samples and determine the F. tularensis subspecies and genotypes involved. This would allow the characterization of the F. tularensis subspecies and genotypes present in Algeria. Furthermore, our findings push us in future studies to seek tularemia in the Algerian human population. To achieve this, interdisciplinary or trans-disciplinary collaborative efforts underpinned by the One Health concept will be necessary. More

According to the information that exists so far regarding reproduction in echinoderms, this is the first work in which the occurrence of trioecy in sea urchins is reported. This is also the first report of trioecy among members of the phylum Echinodermata, one of the most widespread taxa, both latitudinally and bathymetrically. Our results show that trioecy in this population of T. roseus is temporally stable, since the three sexes were observed together throughout the year in each month of sampling. Hermaphroditic individuals also presented the same gametogenic developmental pattern as females and males. Finally, during the spawning period of the population they contributed to the reproductive process by releasing mature gametes, which evidenced their full functionality within the studied population.We were unable to obtain evidence of self-fertilization in the studied hermaphrodites; but self- fertilization in the gonads and gonadal ducts of a hermaphrodite individual of Echinocardium cordatum was recorded in 193543. However, the embryos produced did not complete development successfully, probably due to the premature fertilization within the gonad43. Also, the cases of fully functional hermaphrodites of Arbacia punctulata have been reported44,45. The gametes of the hermaphrodites were fertilized as soon as they were released into seawater and the development of self-fertilized eggs was absolutely normal in time and morphology. After nine days, typical pluteus larvae were obtained and both the eggs and sperm of the hermaphrodites functioned ordinarily with gametes from other males and females.Therefore, we consider that there are no reasons to think that in the case of Toxopneustes roseus hermaphrodites cannot carry out self-fertilization. According to the analysis of the gonad developmental stages, their gametes were released into seawater. Theoretically, those gametes would be able to follow the normal course of fertilization, interacting among them and with gametes of females and males.The trioecic condition has been recorded so far only in some animals, such as a few nematode species and a hydra9,10,14,46,47,48. In marine invertebrates, it has been reported in one anemone under laboratory conditions and in one bivalve mollusk15,16. The coexistence of males, females and hermaphrodites has been considered an evolutionarily transitory state; for example, androdioecy (male / hermaphrodite) in nematodes such as Caenorhabditis elegans is believed to have evolved from dioecy (male / female) through a trioic intermediate. Consequently, it is very difficult to find the ecological or evolutionary causes that lead a species or population to present three sexes simultaneously49.In the species in which trioecy has been studied and monitored, it is noticeable that their populations are subjected to strong environmental stress in situ or under laboratory manipulation50,51,52. For example, some nematodes of the genus Tokorhabditis are extremophilic species that live in the Californian Mono Lake, which is characterized by being hypersaline and exhibiting high levels of arsenic10,50. In the case of Auanema freiburgensis the flexible sex determination and mating system and, consequently, its trioecy can be critical for resilience at the population level in patchy, resource-limited environments49. These results thus demonstrate that life-history, ecology and environment can play defining roles in the development of sexual systems and determine the continued presence of trioecy in the nematode. In the case of Hydra viridissima, it unlike most European species, is a “warm crisis” hydra, since it usually reproduces asexually, but when the temperatures rise to, or are maintained at high levels (≥ 20 °C), it reproduces sexually14,53. In experimental conditions, the population studied essentially behaved as androdioecic and only at the end of the research period, when the temperature was the highest (~ 25 °C), a few females appeared and joined the other existing sexes, thus generating the condition of trioecy14. Trioecy has been identified in another non-described species (e.g., Rhabditis sp. JU1783) isolated from star fruit, although it is closely related to A. rhodensis and A. freiburgensis and likely to belong to the same genus11,12. Little is known about the ecology of Auanema, as A. rhodensis has been isolated from a tick and a beetle, and A. freiburgensis from dung and a rotting plant of the genus Petasites12,47,51.Regarding the sea anemone Aiptasia diaphana, it is mainly found in isolated fouling communities, and no hermaphrodites exist in natural populations that could reproduce asexually or sexually54. However, under laboratory conditions, a single founder individual (asexual clone) produced not only males and females, but also hermaphroditic individuals. In addition, A. diaphana can fertilize within and between cloning lines, producing larval-swimming planules, which could explain the success of the species as an invader of artificial marine substrates. The condition of trioecy was also identified in individuals of this anemone manipulated in the laboratory, to create age-homogeneous populations of asexual propagules (pedal lacerations) and ontogenetic patterns of sexual differentiation were documented15.In the case of the marine bivalve Semimytilus algosus, there was not an obvious explanation for the occurrence of its trioecy, despite the intense analyses of factors such as motility versus a sessile way of life or reproductive density within a population, which could have relevance for gamete interactions16. In many respects, S. algosus is a “typical” marine intertidal mussel, since it is sessile in adulthood, occurs at high densities in wild populations, and has a very large population. S. algosus also co-occurs with other species that are close relatives within the Mytilidae family and have evolved and conserved their dioecy16.Toxopneustes roseus is another typical species of sea urchin, which has a wide latitudinal distribution throughout the tropical eastern Pacific and co-inhabits with other species of sea urchins and echinoderms that have a similar distribution and in which hermaphroditism has not been reported40,55,56,57. Regarding its population density, T. roseus is not considered among the most abundant species in the study area and its densities are relatively low (between 0.04 and 1.2 ind.m2). However, it cannot be considered a rare species in terms of abundance58,59.All of the above makes it difficult to clearly explain the reasons for the occurrence of trioecy in this species; however, certain aspects of its early development are known that could indicate the factors behind the development of this reproductive mating system in the pink sea urchin. In recent experiments carried out with gametes, larvae, and embryos of a population of T. roseus from the same area as our study, it was found that the increase in temperature above the normal values of its habitat has a deleterious effect on the success of early development60. There exists experimental evidence that at an increase of temperature to 32 °C, which is 2 °C above the maximum values registered in the study area, fertilization occurred at a very low percentage. There was also a deleterious effect on embryos, resulting in abnormal development and the lowest percentage of larval survival also occurred at 32 °C60. The same kind of experiments has been performed on other species from the study area, such as the irregular sea urchin Ryncholampas pacificus and the intertidal Echinometra vanbrunti. The deleterious effects on these species were observed only at 34 °C, which was the highest temperature tested (unpublished data). At 32 °C, however, there was no evidence of negative effects in the case on E. vanbrunti, and there was just arrested development, but no abnormalities in the case of R. pacificus. These results indicate that T. roseus is much more sensitive to the rise in temperature than other cohabiting sea urchins, and probably lives near its upper thermal limit. In that context, the continuous ocean warming could threaten the permanence of the species in the study area, since the early stages of development constitute a bottleneck for successful recruitment and later population maintenance in populations that carry out reproduction by means of external fertilization.Within the phylum Echinodermata, when stressful conditions appear in the habitat or the environment becomes hostile, the species can generally resort to asexual reproduction by fission (ophiuroids) or fission and autotomy (holothuroids and asteroids) to increase the abundance of populations in a relatively short time or counteract a threat with numbers61. This does not apply to sea urchins since they are unable to reproduce asexually. The only way for sea urchins to reproduce asexually would be by cloning larvae, but this process would also require that sexual reproduction occurs first62. Therefore, any reproductive strategy that a sea urchin population could develop to respond to drastic changes in their area must involve sexual reproduction. In this regard, in an experimental evolution study with the nematode Caenorhabditis elegans, in which partial selfing, exclusive selfing, and predominant outcrossing were compared, it was evidenced that monoecious populations only have hermaphrodites and, therefore, reproduction is carried out exclusively by self-fertilization. However, in trioic populations that have males, females, and a small number of hermaphrodites, reproduction is predominantly carried out by external crossing49. Also populations that underwent some degree of interbreeding during the evolutionary experiments (trioic and androdioic populations), maintained more genetic diversity than expected solely under genetic drift or under genetic drift and directional selection49. In this sense, it is possible that high levels of interbreeding, such as that which occurs in trioic populations, develop with populations that have sufficient deleterious recessive alleles to avoid extinction, since selection is less efficient to purge them. Trioecy, therefore, becomes an efficient system to select characteristics of the genome that allows a population that only reproduces sexually to adequately cope with significant changes in the environment that could threaten the permanence of the species in that habitat. Interbreeding (gonochorism, self-incompatible hermaphroditism) also favors genetic diversity and offers greater potential to adapt to changing environments63. The costs and advantages of crossing over selfing depend on environmental factors and, therefore, selection may favor transitions between mating systems. Androdioecy, gynodioecy, and trioecy are evolutionarily unstable intermediate strategies, but they offer important systems for testing models of the causes and consequences of the mating system in the evolution of populations63.However, the question remains why T. roseus has developed trioecy, when in the same habitat there are other sea urchins with very similar life-histories that only maintain dioecy. In the case of the bivalve Semimytilus algosus; which presents the same situation as we have with T. roseus, it was proposed that the trioecy of the species may be related to the sex determination mechanism, considering what it is known about the nematodes of the genus Auanema10,16,46. In Auanema, the male versus non-male (hermaphrodite or female) decision is determined genetically (XO for males, and XX for females and hermaphrodites)9,64. The hermaphrodite versus female decision, however, is determined by the environment of the mother. For A. freiburgensis the maternal social environment is determinant, whereas for A. rhodensis it is the age of the mother9,12,51,65. Therefore, in Auanema, environmental sex determination and genetic sex determination interact to produce trioecy.Although there is apparently no clear cause of strong, stressful conditions in the habitat of T. roseus that could threaten the survival of this species, according to the United States Environmental Protection Agency (EPA, 2021), sea surface temperature increased during the twentieth century and continues to rise. From 1901 to 2020, the global temperature rose at an average rate of 0.004 °C per decade, resulting in a total increase of 0.5 °C to date. Additionally, regional studies based on continuous monitoring, which have not yet been published, have shown that between 2002 and 2020 there has been an increase of approximately 1 °C above the historical average of the sea surface temperature in the study area.The foregoing discussion leads us to speculate that the studied population of T. roseus lives at the limit of its thermal tolerance, and the constant increase in ocean temperature due to global warming constitutes a threat to its survival and a constant source of stress for the population. This is because its early-development stages are more vulnerable to high temperature than other sea urchins that live in the same area and its population density is also significantly lower58.Phylogenetically T. roseus belongs to Family Toxopneustidae and although no other species within the genus Toxopneustes has shown hermaphroditism, this condition was reported in Tripneustes gratilla, which belongs to the same family36. Toxopneustids belong to the Order Camarodonta, and almost all the species of sea urchins in which hermaphroditism has been reported belong to this Order except for a couple that belong to the Arbacioida. At the same time, this order is contained in the Superorder Echinacea along with Camarodonta, according to the last exhaustive analysis resolving the position of the clades within Echinoidea66. In this context, theoretically T. roseus at some point underwent the environmental pressure of its early stage living under constantly rising temperatures, along with its low population densities in the study area. Consequently, it was able to develop hermaphroditism and, therefore, trioecy, similarly to what occurred to Hydra viridissima under conditions of extreme high temperature14. We hypothesize that these permanent conditions generate a constant source of strong environmental stress, which is the determining factor that keeps trioecy stable in the species in which it has been studied, and, thus, trioecy remains stable in this population of T. roseus.The mechanism of sex determination in echinoids, as well as in other echinoderms, is still unknown, although the sex ratio, which is generally close to 1:1, suggests that it occurs through sex chromosomes67. It is known that in mammals, sex determination is dictated by the presence or absence of the Y-chromosomal gene SRY. SRY functions as the primary sex-determining gene by activating testis formation, and in its absence, the embryo will form ovaries. SRY only exists in mammals; however it evolved as a duplication of the Sox gene family, which exists in all metazoans68.In vertebrates, Sox genes are involved in sex determination, neurogenesis, skeletonogenesis, eye development, pituitary development, pancreas formation, and neural crest and notochord formation69. In invertebrates, they are involved in processes such as metamorphosis, eye development, neural crest formation, and ectoderm formation70. In the sea urchin Strongylocentrotus purpuratus, SoxB1 was determined to be expressed in the primordial gut during development and is closely related in sequence to Sox genes of the mouse embryo71. An investigation of sex determination was carried out in the sea urchin Strongylocentrotus purpuratus using RNA-seq and quantitative mRNA measurements, but the mechanisms that govern sexual determination of the species could not be clearly established72. However; the results show that the male fate factors Dmrt and SoxH are expressed early and meiosis initiates early. Also, gonad-specific transcripts involved in egg and sperm biology, are first activated before rudiment formation in the larvae of this sea urchin. The study provided additional evidence for the hypothesis that in sea urchins, sex determination occurs genetically72. Another research with the sea cucumber Apostichopus japonicus, which integrated genome-wide association study and analyzes of sex-specific variations evidenced that the species exhibits genetic sexual determination73. Furthermore, analysis of homozygous and heterozygous genotypes of abundant sex-specific SNPs in females and males, confirmed that A.japonicus might have a XX/XY sex determination system73.On the other hand, it has been proposed that a deviation from the 1:1 sex ratio in echinoids could reflect environmental conditions that influence sex determination67. For example, a relatively large proportion of Lytechinus variegatus and Tripneustes ventricosus (as Tripneustes esculentus) hermaphrodites was recorded in southern Florida during an unusually cold winter, suggesting that adverse winter conditions in some way affected sex determination in juveniles74,75. Also relatively large number of Strongylocentrotus purpuratus hermaphrodites was reported in Bahía de Todos los Santos, Mexico, where extreme seasonal fluctuations in temperature (from about 12–24 °C) are recorded76. However, posterior studies did not find a single hermaphrodite of Strongylocentrotus purpuratus in more than 500 individuals analyzed77,78.Considering that sex determination in sea urchins is highly probable to occur genetically and the possibility that the environment may also influence sex determination, we think that in the case of Toxopneustes roseus, genetic sex determination and environmental sex determination are interacting to maintain the condition of trioecy stable. We propose that, especially because the cases in which environmental conditions have assumed to influence sex determination, extreme temperatures are invoked as the main affecting factor. However, more detailed studies are needed in terms of sexual determination and experimental evolution to be able to verify our assumption.In general, the efforts that have been made to explain the evolution of the sexes and the origin of hermaphroditism and trioecy are still scarce, and critical questions remain to be answered. The case of trioecy detected in T. roseus may constitute an important model to seek these answers about the evolution of sexual systems and the environmental mechanisms that trigger trioecy in marine macroinvertebrates and, in particular, in echinoderms. More

Savenije, H. H. G. The importance of interception and why we should delete the term evapotranspiration from our vocabulary. Hydrol. Process. 18, 1507–1511 (2004).Article 
ADS 

Google Scholar 
Gerrits, A. M. J., Pfister, L. & Savenije, H. H. G. Spatial and temporal variability of canopy and forest floor interception in a beech forest. Hydrol. Process. 24, 3011–3025 (2010).Article 
ADS 

Google Scholar 
Porada, P., Van Stan, J. T. & Kleidon, A. Significant contribution of non-vascular vegetation to global rainfall interception. Nat. Geosci. 11, 563–567 (2018).Article 
ADS 
CAS 

Google Scholar 
van der Ent, R. J., Wang-Erlandsson, L., Keys, P. W. & Savenije, H. H. G. Contrasting roles of interception and transpiration in the hydrological cycle – Part 2: moisture recycling. Earth Syst. Dyn. 5, 471–489 (2014).Article 
ADS 

Google Scholar 
Lian, X. et al. Partitioning global land evapotranspiration using CMIP5 models constrained by observations. Nat. Clim. Change 8, 640–646 (2018).Article 
ADS 

Google Scholar 
Coenders-Gerrits, A. M. et al. Uncertainties in transpiration estimates. Nature 506, E1–E2 (2014).Article 
CAS 

Google Scholar 
Chang, L.-L. et al. Why do large-scale land surface models produce a low ratio of transpiration to evapotranspiration? J. Geophys. Res. Atmos. 123, 9109–9130 (2018).Article 

Google Scholar 
Zwieback, S., Chang, Q., Marsh, P. & Berg, A. Shrub tundra ecohydrology: rainfall interception is a major component of the water balance. Environ. Res. Lett. 14, 055005 (2019).Article 
ADS 

Google Scholar 
Cuartas, L. A. et al. Interception water-partitioning dynamics for a pristine rainforest in Central Amazonia: Marked differences between normal and dry years. Agric. For. Meteorol. 145, 69–83 (2007).Article 
ADS 

Google Scholar 
Yue, K. et al. Global patterns and drivers of rainfall partitioning by trees and shrubs. Glob. Change Biol. 27, 3350–3357 (2021).Article 

Google Scholar 
Pastorello, G. et al. The FLUXNET2015 dataset and the ONEFlux processing pipeline for eddy covariance data. Sci. Data 7, 225 (2020).Article 

Google Scholar 
Tramontana, G. et al. Predicting carbon dioxide and energy fluxes across global FLUXNET sites with regression algorithms. Biogeosciences 13, 4291–4313 (2016).Article 
ADS 
CAS 

Google Scholar 
Jung, M., Reichstein, M. & Bondeau, A. Towards global empirical upscaling of FLUXNET eddy covariance observations: validation of a model tree ensemble approach using a biosphere model. Biogeosciences 6, 2001–2013 (2009).Article 
ADS 
CAS 

Google Scholar 
Li, X. et al. Spatiotemporal pattern of terrestrial evapotranspiration in China during the past thirty years. Agric. For. Meteorol. 259, 131–140 (2018).Article 
ADS 

Google Scholar 
Koppa, A., Rains, D., Hulsman, P., Poyatos, R. & Miralles, D. G. A deep learning-based hybrid model of global terrestrial evaporation. Nat. Commun. 13, 1912 (2022).Article 
ADS 
CAS 

Google Scholar 
Zheng, C. & Jia, L. Global canopy rainfall interception loss derived from satellite Earth observations. Ecohydrology 13, e2186 (2019).
Google Scholar 
Muzylo, A. et al. A review of rainfall interception modelling. J. Hydrol. 370, 191–206 (2009).Article 
ADS 

Google Scholar 
Miralles, D. G., Gash, J. H., Holmes, T. R. H., de Jeu, R. A. M., & Dolman, A. J. Global canopy interception from satellite observations. J. Geophys. Res. 115, D16122 (2010).Article 
ADS 

Google Scholar 
Martens, B. et al. GLEAM v3: satellite-based land evaporation and root-zone soil moisture. Geosci. Model Dev. 10, 1903–1925 (2017).Article 
ADS 

Google Scholar 
Oleson, K. et al. Technical Description of Version 4.5 of the Community Land Model (CLM) Report NCAR/TN-503+STR, https://doi.org/10.5065/D6RR1W7M (2013).Gash, J. An analytical model of rainfall interception by forests. Q. J. Roy. Meteor. Soc. 105, 43–55 (1979).Article 
ADS 

Google Scholar 
Fan, Y. et al. Reconciling canopy interception parameterization and rainfall forcing frequency in the Community Land Model for simulating evapotranspiration of rainforests and oil palm plantations in Indonesia. J. Adv. Model. Earth Syst. 11, 732–751 (2019).Article 
ADS 

Google Scholar 
Návar, J. Modeling rainfall interception loss components of forests. J. Hydrol. 584, 124449 (2019).Article 

Google Scholar 
Kang, M., Kwon, H., Cheon, J. H. & Kim, J. On estimating wet canopy evaporation from deciduous and coniferous forests in the Asian monsoon climate. J. Hydrometeorol. 13, 950–965 (2012).Article 
ADS 

Google Scholar 
Llorens, P., Domingo, F., Garcia-Estringana, P., Muzylo, A. & Gallart, F. Canopy wetness patterns in a Mediterranean deciduous stand. J. Hydrol. 512, 254–262 (2014).Article 
ADS 

Google Scholar 
Czikowsky, M. J. & Fitzjarrald, D. R. Detecting rainfall interception in an Amazonian rain forest with eddy flux measurements. J. Hydrol. 377, 92–105 (2009).Article 
ADS 

Google Scholar 
Renninger, H. J., Phillips, N. & Salvucci, G. D. Wet- vs. dry-season transpiration in an Amazonian rain forest palm iriartea deltoidea. Biotropica 42, 470–478 (2010).Article 

Google Scholar 
Zhao, W. et al. Physics-constrained machine learning of evapotranspiration. Geophys. Res. Lett. 46, 14496–14507 (2019).Article 
ADS 

Google Scholar 
Zabret, K. & Šraj, M. How characteristics of a rainfall event and the meteorological conditions determine the development of stemflow: A case study of a birch tree. Front. Glob. Change 4, 663100 (2022).Article 

Google Scholar 
Calder, I. R. Dependence of rainfall interception on drop size: 1. Development of the two-layer stochastic model. J. Hydrol. 185, 363–378 (1996).Article 
ADS 

Google Scholar 
Niinemets, Ü. A review of light interception in plant stands from leaf to canopy in different plant functional types and in species with varying shade tolerance. Ecol. Res. 25, 693–714 (2010).Article 

Google Scholar 
Gordon, D. A. R., Coenders-Gerrits, M., Sellers, B. A., Sadeghi, S., & Van Stan II, J. T. Rainfall interception and redistribution by a common North American understory and pasture forb, Eupatorium capillifolium (Lam. dogfennel). Hydrol. Earth Syst. Sci. 24, 4587–4599 (2020).Article 
ADS 

Google Scholar 
Ciruzzi, D. M. & Loheide, S. P. II Monitoring tree sway as an indicator of interception dynamics before, during, and following a storm. Geophys. Res. Lett. 48, e2021GL094980 (2021).Article 
ADS 

Google Scholar 
Karimi, P., Bastiaanssen, W. G. & Molden, D. Water Accounting Plus (WA+)–a water accounting procedure for complex river basins based on satellite measurements. Hydrol. Earth Syst. Sci. 17, 2459–2472 (2013).Article 
ADS 

Google Scholar 
del Campo, A. D., González-Sanchis, M., Lidón, A., Ceacero, C. J. & García-Prats, A. Rainfall partitioning after thinning in two low-biomass semiarid forests: Impact of meteorological variables and forest structure on the effectiveness of water-oriented treatments. J. Hydrol. 565, 74–86 (2018).Article 

Google Scholar 
Lian, X. et al. Multifaceted characteristics of dryland aridity changes in a warming world. Nat. Rev. Earth Environ. 2, 232–250 (2021).Article 
ADS 

Google Scholar 
Piao, S. et al. Characteristics, drivers and feedbacks of global greening. Nat. Rev. Earth Environ. 1, 14–27 (2020).Article 
ADS 

Google Scholar 
Feng, X. et al. Revegetation in China’s Loess Plateau is approaching sustainable water resource limits. Nat. Clim. Change 6, 1019–1022 (2016).Article 
ADS 

Google Scholar 
Dawson, T. E. & Goldsmith, G. R. The value of wet leaves. N. Phytol. 219, 1156–1169 (2018).Article 

Google Scholar 
Aparecido, L. M. T., Miller, G. R., Cahill, A. T. & Moore, G. W. Comparison of tree transpiration under wet and dry canopy conditions in a Costa Rican premontane tropical forest. Hydrol. Process. 30, 5000–5011 (2016).Article 
ADS 

Google Scholar 
Huang, L. & Zhang, Z. Effect of rainfall pulses on plant growth and transpiration of two xerophytic shrubs in a revegetated desert area: Tengger Desert, China. CATENA 137, 269–276 (2016).Article 

Google Scholar 
Fathizadeh, O., Hosseini, S., Zimmermann, A., Keim, R. & Boloorani, A. D. Estimating linkages between forest structural variables and rainfall interception parameters in semi-arid deciduous oak forest stands. Sci. Total Environ. 601, 1824–1837 (2017).Article 
ADS 

Google Scholar 
Zhang, Z.-S., Zhao, Y., Li, X.-R., Huang, L. & Tan, H.-J. Gross rainfall amount and maximum rainfall intensity in 60-minute influence on interception loss of shrubs: a 10-year observation in the Tengger Desert. Sci. Rep. 6, 26030 (2016).Article 
ADS 
CAS 

Google Scholar 
de Groen, M. M. & Savenije, H. H. G. A monthly interception equation based on the statistical characteristics of daily rainfall. Water Resour. Res. 42, W12417 (2006).Article 
ADS 

Google Scholar 
Chinita, M. J., Richardson, M., Teixeira, J. & Miranda, P. M. A. Global mean frequency increases of daily and sub-daily heavy precipitation in ERA5. Environ. Res. Lett. 16, 074035 (2021).Article 
ADS 

Google Scholar 
Donat, M. G., Lowry, A. L., Alexander, L. V., O’Gorman, P. A. & Maher, N. More extreme precipitation in the world’s dry and wet regions. Nat. Clim. Change 6, 508–513 (2016).Article 
ADS 

Google Scholar 
IPCC. The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (eds Masson-Delmotte, V. et al) (Cambridge Univ. Press, 2021).Ficklin, D. L., Null, S. E., Abatzoglou, J. T., Novick, K. A. & Myers, D. T. Hydrological intensification will increase the complexity of water resource management. Earth’s Futur. 10, e2021EF002487 (2022).Article 
ADS 

Google Scholar 
Haslwanter, A., Hammerle, A. & Wohlfahrt, G. Open-path vs. closed-path eddy covariance measurements of the net ecosystem carbon dioxide and water vapour exchange: a long-term perspective. Agric. For. Meteorol. 149, 291–302 (2009).Article 
ADS 

Google Scholar 
Migliavacca, M. et al. The three major axes of terrestrial ecosystem function. Nature 598, 468–472 (2021).Article 
ADS 
CAS 

Google Scholar 
Zhang, W. et al. The effect of relative humidity on eddy covariance latent heat flux measurements and its implication for partitioning into transpiration and evaporation. Preprint at https://doi.org/10.2139/ssrn.4106267 (2022).van Dijk, A. I. J. M. et al. Rainfall interception and the coupled surface water and energy balance. Agric. For. Meteorol. 214–215, 402–415 (2015).Article 

Google Scholar 
Barr, A. G., Morgenstern, K., Black, T. A., McCaughey, J. H. & Nesic, Z. Surface energy balance closure by the eddy-covariance method above three boreal forest stands and implications for the measurement of the CO2 flux. Agric. Meteorol. 140, 322–337 (2006).Article 

Google Scholar 
Reichstein, M. et al. Deep learning and process understanding for data-driven Earth system science. Nature 566, 195–204 (2019).Article 
ADS 
CAS 

Google Scholar 
Zhi, W. et al. From hydrometeorology to river water quality: can a deep learning model predict dissolved oxygen at the continental scale? Environ. Sci. Technol. 55, 2357–2368 (2021).Article 
ADS 
CAS 

Google Scholar 
Kraft, B., Jung, M., Körner, M. & Reichstein, M. Hybrid modeling: fusion of a deep approach and physics-based model for global hydrological modeling. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 43, 1537–1544 (2020).Article 

Google Scholar 
Hoffmann, L. et al. From ERA-Interim to ERA5: the considerable impact of ECMWF’s next-generation reanalysis on Lagrangian transport simulations. Atmos. Chem. Phys. 19, 3097–3124 (2019).Article 
ADS 
CAS 

Google Scholar 
Wang, D., Wang, G. & Anagnostou, E. N. Evaluation of canopy interception schemes in land surface models. J. Hydrol. 347, 308–318 (2007).Article 
ADS 

Google Scholar 
Wang, G. & Eltahir, E. A. Modeling the biosphere–atmosphere system: The impact of the subgrid variability in rainfall interception. J. Clim. 13, 2887–2899 (2000).Article 
ADS 

Google Scholar 
Sitch, S. et al. Recent trends and drivers of regional sources and sinks of carbon dioxide. Biogeosciences 12, 653–679 (2015).Article 
ADS 

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

Google Scholar  More

Ravindran, S. Coral reefs at a tipping point. Proc. Natl Acad. Sci. 113, 5140–5141 (2016).CAS 

Google Scholar 
Lenton, T. M. et al. Climate tipping points—too risky to bet against. Nature 575, 592–595 (2019).CAS 

Google Scholar 
Veron, J. E. N. et al. Delineating the Coral Triangle. Galaxea J. Coral Reef. Stud. 11, 91–100 (2009).
Google Scholar 
Hoegh-Guldberg, O. et al. Coral Reefs Under Rapid Climate Change and Ocean Acidification. Science 318, 1737–1742 (2007).CAS 

Google Scholar 
Brown, C., Corcoran, E. & Herkenrath, P. Marine and coastal ecosystems and human well-being: a synthesis report based on the findings of the Millennium Ecosystem Assessment. (2006).Heinze, C. et al. The quiet crossing of ocean tipping points. Proc. Natl Acad. Sci. 118, e2008478118 (2021).CAS 

Google Scholar 
Barber, P. H. The challenge of understanding the Coral Triangle biodiversity hotspot. J. Biogeogr. 36, 1845–1846 (2009).
Google Scholar 
Ekman, S. Zoogeography of the Sea. (Sidgwick & Jackson, 1953).Ladd, H. S. Origin of the Pacific island molluscan fauna. Am. J. Sci. 256, 137–150 (1960).
Google Scholar 
Woodland, D. J. Zoogeography of the Siganidae (Pisces): an interpretation of distribution and richness patterns. Bull. Mar. Sci. 33, 713–717 (1983).
Google Scholar 
Loveland, T. R. & Merchant, J. M. Ecoregions and ecoregionalization: geographical and ecological perspectives. Environ. Manag. 34, S1–S13 (2004).
Google Scholar 
Levins, R. Some Demographic and Genetic Consequences of Environmental Heterogeneity for Biological Control. Bull. Entomol. Soc. Am. 15, 237–240 (1969).
Google Scholar 
Obura, D. The Diversity and Biogeography of Western Indian Ocean Reef-Building Corals. PLoS One. 7, e45013 (2012).CAS 

Google Scholar 
Fontoura, L. et al. Protecting connectivity promotes successful biodiversity and fisheries conservation. Science 375, 336–340 (2022).CAS 

Google Scholar 
Roberts, C. M. Connectivity and Management of Caribbean Coral Reefs. Science 278, 1454–1457 (1997).CAS 

Google Scholar 
Ayre, D. J. & Hughes, T. P. Climate change, genotypic diversity and gene flow in reef-building corals: Gene flow in reef building corals. Ecol. Lett. 7, 273–278 (2004).
Google Scholar 
Graham, N. A. et al. Dynamic fragility of oceanic coral reef ecosystems. Proc. Natl Acad. Sci. 103, 8425–8429 (2006).CAS 

Google Scholar 
McClanahan, T. R. et al. Prioritizing Key Resilience Indicators to Support Coral Reef Management in a Changing Climate. PLoS One. 7, e42884 (2012).CAS 

Google Scholar 
Gilmour, J. P., Smith, L. D., Heyward, A. J., Baird, A. H. & Pratchett, M. S. Recovery of an Isolated Coral Reef System Following Severe Disturbance. Science 340, 69–71 (2013).
Google Scholar 
Grayson, N., Clements, C. S., Towner, A. A., Beatty, D. S. & Hay, M. E. Did the historic overharvesting of sea cucumbers make coral more susceptible to pathogens? Coral Reefs. 41, 447–453 (2022).
Google Scholar 
Spalding, M. D. et al. Marine Ecoregions of the World: A Bioregionalization of Coastal and Shelf Areas. BioScience 57, 573–583 (2007).
Google Scholar 
Berline, L., Rammou, A.-M., Doglioli, A., Molcard, A. & Petrenko, A. A Connectivity-Based Eco-Regionalization Method of the Mediterranean Sea. PLoS ONE. 9, e111978 (2014).
Google Scholar 
Ser-Giacomi, E., Rossi, V., López, C. & Hernández-García, E. Flow networks: A characterization of geophysical fluid transport. Chaos Interdiscip. J. Nonlinear Sci. 25, 036404 (2015).
Google Scholar 
Thompson, D. M. et al. Variability in oceanographic barriers to coral larval dispersal: Do currents shape biodiversity? Prog. Oceanogr. 165, 110–122 (2018).
Google Scholar 
Treml, E. A., Halpin, P. N., Urban, D. L. & Pratson, L. F. Modeling population connectivity by ocean currents, a graph-theoretic approach for marine conservation. Landsc. Ecol. 23, 19–36 (2008).
Google Scholar 
Liu, G., Bracco, A., Quattrini, A. M. & Herrera, S. Kilometer-Scale Larval Dispersal Processes Predict Metapopulation Connectivity Pathways for Paramuricea biscaya in the Northern Gulf of Mexico. Front. Mar. Sci. 8, 790927 (2021).
Google Scholar 
Fountalis, I., Dovrolis, C., Bracco, A., Dilkina, B. & Keilholz, S. δ-MAPS: from spatio-temporal data to a weighted and lagged network between functional domains. Appl. Netw. Sci. 3, 21 (2018).
Google Scholar 
Falasca, F., Bracco, A., Nenes, A. & Fountalis, I. Dimensionality Reduction and Network Inference for Climate Data Using δ‐MAPS: Application to the CESM Large Ensemble Sea Surface Temperature. J. Adv. Model. Earth Syst. 11, 1479–1515 (2019).
Google Scholar 
Novi, L., Bracco, A. & Falasca, F. Uncovering marine connectivity through sea surface temperature. Sci. Rep. 11, 8839 (2021).CAS 

Google Scholar 
Kleypas, J. A., Castruccio, F. S., Curchitser, E. N. & Mcleod, F. The impact of ENSO on coral heat stress in the western equatorial Pacific. Glob. Change Biol. 21, 2525–2539 (2015).
Google Scholar 
GLOBAL_REANALYSIS_001_030. Global Ocean Physics Reanalysis GLORYS12V1 1/12° product. MERCATOR GLORYS12V1 (global-reanalysis-001-030-monthly). E.U. Copernicus Marine Service Information (CMEMS). https://doi.org/10.48670/moi-00021.Lellouche, J.-M. et al. The Copernicus Global 1/12° Oceanic and Sea Ice GLORYS12 Reanalysis. Front. Earth Sci. 9, 698876 (2021).
Google Scholar 
Treml, E. A. & Halpin, P. N. Marine population connectivity identifies ecological neighbors for conservation planning in the Coral Triangle: Ecological neighbors in conservation. Conserv. Lett. 5, 441–449 (2012).
Google Scholar 
Meyers, G. Variation of Indonesian throughflow and the El Niño-Southern Oscillation. J. Geophys. Res. Oceans 101, 12255–12263 (1996).
Google Scholar 
Wolfram Research (2012), FindGraphCommunities, Wolfram Language function. https://reference.wolfram.com/language/ref/FindGraphCommunities.html (updated 2015).MacArthur, R. H. & Wilson, E. O. The theory of island biogeography. In The Theory of Island Biogeography (Princeton university press, 2016).Brin, S. & Page, L. The anatomy of a large-scale hypertextual Web search engine. Comput. Netw. ISDN Syst. 30, 107–117 (1998).
Google Scholar 
Wolfram Research (2010), PageRankCentrality, Wolfram Language function. https://reference.wolfram.com/language/ref/PageRankCentrality.html (Updated 2015).NOAA Coral Reef Watch program, 20180813, NOAA Coral Reef Watch Version 3.1 Daily Global 5km Satellite Coral Bleaching Heat Stress Monitoring Product Suite: NOAA Coral Reef Watch program, College Park, Maryland, USA. https://coralreefwatch.noaa.gov/product/5km/.Liu, G. et al. Reef-Scale Thermal Stress Monitoring of Coral Ecosystems: New 5-km Global Products from NOAA Coral Reef Watch. Remote Sens. 6, 11579–11606 (2014).
Google Scholar 
Liu, G. et al. NOAA Coral Reef Watch’s 5km Satellite Coral Bleaching Heat Stress Monitoring Product Suite Version 3 and Four-Month Outlook Version 4. 32, 7 (2017).Claar, D. C., Szostek, L., McDevitt-Irwin, J. M., Schanze, J. J. & Baum, J. K. Global patterns and impacts of El Niño events on coral reefs: A meta-analysis. PLOS ONE 13, e0190957 (2018).
Google Scholar 
Sully, S., Burkepile, D. E., Donovan, M. K., Hodgson, G. & van Woesik, R. A global analysis of coral bleaching over the past two decades. Nat. Commun. 10, 1264 (2019).CAS 

Google Scholar 
Darling, E. S. et al. Social–environmental drivers inform strategic management of coral reefs in the Anthropocene. Nat. Ecol. Evol. 3, 1341–1350 (2019).
Google Scholar 
Dance, A. These corals could survive climate change—and help save the world’s reefs. Nature 575, 580–582 (2019).CAS 

Google Scholar 
Renema, W. et al. Hopping Hotspots: Global Shifts in Marine Biodiversity. Science 321, 654–657 (2008).CAS 

Google Scholar 
Weiss, T. L., Denniston, R. F., Wanamaker, A. D., Villarini, G. & von der Heydt, A. S. El Niño–Southern Oscillation–like variability in a late Miocene Caribbean coral. Geology 45, 643–646 (2017).
Google Scholar 
Watanabe, T. et al. Permanent El Niño during the Pliocene warm period not supported by coral evidence. Nature 471, 209–211 (2011).CAS 

Google Scholar 
Von Der Heydt, A. S. & Dijkstra, H. A. The impact of ocean gateways on ENSO variability in the Miocene. Geol. Soc. Lond. Spec. Publ. 355, 305–318 (2011).
Google Scholar 
Yasuhara, M. et al. Past and future decline of tropical pelagic biodiversity. Proc. Natl Acad. Sci. 117, 12891–12896 (2020).CAS 

Google Scholar 
Falasca, F., Crétat, J., Bracco, A., Braconnot, P. & Marti, O. Climate change in the Indo-Pacific basin from mid- to late Holocene. Clim. Dyn. 59, 753–766 (2022).
Google Scholar 
Treml, E. A., Ford, J. R., Black, K. P. & Swearer, S. E. Identifying the key biophysical drivers, connectivity outcomes, and metapopulation consequences of larval dispersal in the sea. Mov. Ecol. 3, 17 (2015).
Google Scholar 
Hackerott, S., Martell, H. A. & Eirin-Lopez, J. M. Coral environmental memory: causes, mechanisms, and consequences for future reefs. Trends Ecol. Evol. 36, 1011–1023 (2021).
Google Scholar 
Ogle, K. et al. Quantifying ecological memory in plant and ecosystem processes. Ecol. Lett. 18, 221–235 (2015).
Google Scholar 
Peterson, G. D. Contagious Disturbance, Ecological Memory, and the Emergence of Landscape Pattern. Ecosystems 5, 329–338 (2002).
Google Scholar 
Thomas, L., López, E. H., Morikawa, M. K. & Palumbi, S. R. Transcriptomic resilience, symbiont shuffling, and vulnerability to recurrent bleaching in reef‐building corals. Mol. Ecol. 28, 3371–3382 (2019).
Google Scholar 
Dziedzic, K. E., Elder, H., Tavalire, H. & Meyer, E. Heritable variation in bleaching responses and its functional genomic basis in reef‐building corals (Orbicella faveolata). Mol. Ecol. 28, 2238–2253 (2019).
Google Scholar 
Ainsworth, T. D. et al. Climate change disables coral bleaching protection on the Great Barrier Reef. Science 352, 338–342 (2016).CAS 

Google Scholar 
Harrison, H. B., Bode, M., Williamson, D. H., Berumen, M. L. & Jones, G. P. A connectivity portfolio effect stabilizes marine reserve performance. Proc. Natl Acad. Sci. 117, 25595–25600 (2020).CAS 

Google Scholar 
Leeuwenburgh, O. & Stammer, D. The Effect of Ocean Currents on Sea Surface Temperature Anomalies. J. Phys. Oceanogr. 31, 2340–2358 (2001).
Google Scholar 
Box, G. E., Jenkins, G. M. & Reinsel, G. C. Time series analysis: forecasting and control. (Wiley, 2011).Falasca, F. & Bracco, A. Exploring the tropical Pacific manifold in models and observations. Phys. Rev. X 12, 021054 (2022).CAS 

Google Scholar 
NOAA (National Oceanic and Atmospheric Administration), (2019a). Nino regions. https://www.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/nino_regions.shtml.NOAA (National Oceanic and Atmospheric Administration), (2019b). Cold and warm episodes by season. https://origin.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ONI_v5.php.Baird, A. et al. Coral Spawning Database. 10552719 Bytes https://doi.org/10.25405/DATA.NCL.13082333.V1 (2020).UNEP-WCMC, WorldFish Centre, WRI, TNC (2021). Global distribution of warm-water coral reefs, compiled from multiple sources including the Millennium Coral Reef Mapping Project. Version 4.1. Includes contributions from IMaRS-USF and IRD (2005), IMaRS-USF (2005) and Spalding et al. (2001). Cambridge (UK): UN Environment World Conservation Monitoring Centre. Data https://doi.org/10.34892/t2wk-5t34.IMaRS-USF, IRD (Institut de Recherche pour le Developpement) (2005). Millennium Coral Reef Mapping Project. Validated maps. Cambridge (UK): UNEP World Conservation Monitoring Centre.IMaRS-USF (Institute for Marine Remote Sensing-University of South Florida) (2005). Millennium Coral Reef Mapping Project. Unvalidated maps. These maps are unendorsed by IRD, but were further interpreted by UNEP World Conservation Monitoring Centre. Cambridge (UK): UNEP World Conservation Monitoring Centre.Spalding, M., Ravilious, C. & Green, E. World atlas of coral reefs. Choice Rev. Online. 39, 39-2540–39–2540 (2002).
Google Scholar  More

Yu, Z., Loisel, J., Brosseau, D. P., Beilman, D. W. & Hunt, S. J. Global peatland dynamics since the Last Glacial Maximum. Geophys. Res. Lett. 37, L13402 (2010).
Google Scholar 
Joosten, H., Tapio-BiströmM, L. & Susanna, T. Peatlands: guidance for climate change mitigation through conservation, rehabilitation and sustainable use. Food and Agriculture Organization of the United Nations and Wetlands International. FAO (2012).IUCN. Issues brief: peatlands and climate change. www.icun.org (2017).Joosten, H. Peatlands, Climate Change Mitigation and Biodiversity Conservation. An Issue Brief on the Importance of Peatlands for Carbon and Biodiversity Conservation and the Role of Drained Peatlands as Greenhouse Gas Emission Hotspots (Nordic Council of Ministers, 2015).Moore, T. R. & Knowles, R. The influence of water table levels on methane and carbon dioxide emissions from peatland soils. Can. J. Soil Sci. 69, 33–38 (1989).CAS 

Google Scholar 
Tfaily, M. M. et al. Organic matter transformation in the peat column at Marcell Experimental Forest: humification and vertical stratification. J. Geophys. Res. Biogeosci. 119, 661–675 (2014).CAS 

Google Scholar 
Clymo, R. S. & Bryant, C. L. Diffusion and mass flow of dissolved carbon dioxide, methane, and dissolved organic carbon in a 7-m deep raised peat bog. Geochim. Cosmochim. Acta 72, 2048–2066 (2008).CAS 

Google Scholar 
Clymo, R. S. The limits to peat bog growth. Philos. Trans. R. Soc. B 303, 605–654 (1984).
Google Scholar 
Qin, S. et al. Temperature sensitivity of SOM decomposition governed by aggregate protection and microbial communities. Sci. Adv. 5, eaau1218. 1211–1219 (2019).
Google Scholar 
Dorrepaal, E. et al. Carbon respiration from subsurface peat accelerated by climate warming in the subarctic. Nature 460, 616–619 (2009).CAS 

Google Scholar 
Luo, Z. K., Wang, G. C. & Wang, E. L. Global subsoil organic carbon turnover times dominantly controlled by soil properties rather than climate. Nat. Commun. 10, 3688 (2019).
Google Scholar 
Wilson, R. M. et al. Stability of peatland carbon to rising temperatures. Nat. Commun. 7, 13723 (2016).CAS 

Google Scholar 
Sihi, D., Inglett, P. W. & Inglett, K. S. Carbon quality and nutrient status drive the temperature sensitivity of organic matter decomposition in subtropical peat soils. Biogeochemistry 131, 103–119 (2016).CAS 

Google Scholar 
Wang, Q., Liu, S. & Tian, P. Carbon quality and soil microbial property control the latitudinal pattern in temperature sensitivity of soil microbial respiration across Chinese forest ecosystems. Glob. Chang. Biol. 24, 2841–2849 (2018).
Google Scholar 
Cheng, L. et al. Warming enhances old organic carbon decomposition through altering functional microbial communities. ISME J. 11, 1825–1835 (2017).
Google Scholar 
Luan, J., Wu, J., Liu, S., Roulet, N. & Wang, M. Soil nitrogen determines greenhouse gas emissions from northern peatlands under concurrent warming and vegetation shifting. Commun. Biol. 2, 132 (2019).
Google Scholar 
Meyer, N. et al. Nitrogen and phosphorus supply controls soil organic carbon mineralization in tropical topsoil and subsoil. Soil Biol. Biochem. 119, 152–161 (2018).CAS 

Google Scholar 
Fontaine, S. et al. Stability of organic carbon in deep soil layers controlled by fresh carbon supply. Nature 450, 277–280 (2007).CAS 

Google Scholar 
Moni, C. et al. Temperature response of soil organic matter mineralisation in arctic soil profiles. Soil Biol. Biochem. 88, 236–246 (2015).CAS 

Google Scholar 
Xu, X., Sherry, R. A., Niu, S., Zhou, J. & Luo, Y. Long-term experimental warming decreased labile soil organic carbon in a tallgrass prairie. Plant Soil 361, 307–315 (2012).CAS 

Google Scholar 
Broder, T., Blodau, C., Biester, H. & Knorr, K. H. Peat decomposition records in three pristine ombrotrophic bogs in southern Patagonia. Biogeosciences 9, 1479–1491 (2012).CAS 

Google Scholar 
Adamczyk, M., Perez-Mon, C., Gunz, S. & Frey, B. Strong shifts in microbial community structure are associated with increased litter input rather than temperature in High Arctic soils. Soil Biol. Biochem. 151, 108054 (2020).CAS 

Google Scholar 
Hug, L. A. et al. Community genomic analyses constrain the distribution of metabolic traits across the Chloroflexi phylum and indicate roles in sediment carbon cycling. Microbiome 1, 22 (2013).
Google Scholar 
Yun, J. L., Ju, Y. W., Deng, Y. C. & Zhang, H. X. Bacterial community structure in two permafrost wetlands on the Tibetan Plateau and Sanjiang Plain, China. Microb. Ecol. 68, 360–369 (2014).
Google Scholar 
Zhong, Q. et al. Water table drawdown shapes the depth-dependent variations in prokaryotic diversity and structure in Zoige peatlands. FEMS Microbiol. Ecol. https://doi.org/10.1093/femsec/fix049 (2017).Article 

Google Scholar 
Karhu, K. et al. Temperature sensitivity of soil respiration rates enhanced by microbial community response. Nature 513, 81–84 (2014).CAS 

Google Scholar 
Thiessen, S., Gleixner, G., Wutzler, T. & Reichstein, M. Both priming and temperature sensitivity of soil organic matter decomposition depend on microbial biomass – An incubation study. Soil Biol. Biochem. 57, 739–748 (2013).CAS 

Google Scholar 
Walker, T. W. N. et al. Microbial temperature sensitivity and biomass change explain soil carbon loss with warming. Nat. Clim. Chang. 8, 885–899 (2018).CAS 

Google Scholar 
Dungait, J. A. J., Hopkins, D. W., Gregory, A. S. & Whitmore, A. P. Soil organic matter turnover is governed by accessibility not recalcitrance. Glob. Chang. Biol. 18, 1781–1796 (2012).
Google Scholar 
Conant, R. T. et al. Temperature and soil organic matter decomposition rates – synthesis of current knowledge and a way forward. Global Chang. Biol. 17, 3392–3404 (2011).
Google Scholar 
Hietz, P. et al. Long-term change in the nitrogen cycle of tropical forests. Science 4, 334 (2011).
Google Scholar 
Manzoni, S., Taylor, P., Richter, A., Porporato, A. & Agren, G. I. Environmental and stoichiometric controls on microbial carbon-use efficiency in soils. New Phytol. 196, 79–91 (2012).CAS 

Google Scholar 
Sistla, S. A., Asao, S. & Schimel, J. P. Detecting microbial N-limitation in tussock tundra soil: Implications for Arctic soil organic carbon cycling. Soil Biol. Biochem. 55, 78–84 (2012).CAS 

Google Scholar 
Chen, L. et al. Nitrogen availability regulates topsoil carbon dynamics after permafrost thaw by altering microbial metabolic efficiency. Nat. Commun. 9, 3951 (2018).
Google Scholar 
Soong, J. L. et al. Five years of whole-soil warming led to loss of subsoil carbon stocks and increased CO2 efflux. Sci. Adv. 7, eabd1343 (2021).Chen, L. et al. Determinants of carbon release from the active layer and permafrost deposits on the Tibetan Plateau. Nat. Commun. 7, 13046 (2016).CAS 

Google Scholar 
Girkin, N. T. et al. Interactions between labile carbon, temperature and land use regulate carbon dioxide and methane production in tropical peat. Biogeochemistry 147, 87–97 (2019).
Google Scholar 
Swails, E. et al. Will CO2 emissions from drained tropical peatlands decline over time? Links between soil organic matter quality, nutrients, and C mineralization rates. Ecosystems 21, 868–885 (2017).
Google Scholar 
Ismawi, S., Gandaseca, S. & Ahmed, O. Effects of deforestation on soil major macro-nutrient and other selected chemical properties of secondary tropical peat swamp forest. Int. J. Phys. Sci. 7, 2225–2228 (2012).CAS 

Google Scholar 
Kimura, S., Melling, L. & Goh, K. Influence of soil aggregate size on greenhouse gas emission and uptake rate from tropical peat soil in forest and different oil palm development years. Geoderma 185, 1–5 (2012).
Google Scholar 
Takakai, F. et al. Effects of agricultural land-use change and forest fire on N2O emission from tropical peatlands, Central Kalimantan, Indonesia. Soil Sci. Plant Nutr. 52, 662–674 (2006).CAS 

Google Scholar 
Knoblauch, C., Beer, C., Sosnin, A., Wagner, D. & Pfeiffer, E. M. Predicting long-term carbon mineralization and trace gas production from thawing permafrost of Northeast Siberia. Glob. Chang. Biol. 19, 1160–1172 (2013).
Google Scholar 
Treat, C. C. et al. Temperature and peat type control CO2 and CH4 production in Alaskan permafrost peats. Glob. Chang. Biol. 20, 2674–2686 (2014).CAS 

Google Scholar 
Hobbie, S. E., Schimel, J. P., Trumbore, S. E. & Randerson, J. Controls over carbon storage and tureover in high-latitude soils. Glob. Chang. Biol. 6, 196–210 (2000).
Google Scholar 
Keller, J. K., Bauers, A. K., Bridgham, S. D., Kellogg, L. E. & Iversen, C. M. Nutrient control of microbial carbon cycling along an ombrotrophic-minerotrophic peatland gradient. J. Geophys. Res. https://doi.org/10.1029/2005jg000152 (2006).Chen, H. et al. A historical overview about basic issues and studies of mires (in Chinese). Sci. Sin. 51, 15–26 (2020).
Google Scholar 
Ridl, J. et al. Plants rather than mineral fertilization shape microbial community structure and functional potential in legacy contaminated soil. Front. Microbiol. 7, 1–10 (2016).
Google Scholar 
Kane, E. S. et al. Response of anaerobic carbon cycling to water table manipulation in an Alaskan rich fen. Soil Biol. Biochem. 58, 50–60 (2013).CAS 

Google Scholar 
Carrell, A. A. et al. Experimental warming alters the community composition, diversity, and N2 fixation activity of peat moss (Sphagnum fallax) microbiomes. Glob. Chang. Biol. 25, 2993–3004 (2019).
Google Scholar 
Lamit, L. J. et al. Patterns and drivers of fungal community depth stratification in Sphagnum peat. FEMS Microbiol. Ecol. 93, fix082 (2017).
Google Scholar 
Harrison, R. B., Footen, P. W. & Strahm, B. D. Deep soil horizons: contribution and importance to soil carbon pools and in assessing whole-ecosystem response to management and global change. Forest Sci. 57, 67–76 (2011).
Google Scholar 
Krüger, J. P., Leifeld, J., Glatzel, S., Szidat, S. & Alewell, C. Biogeochemical indicators of peatland degradation – a case study of a temperate bog in northern Germany. Biogeosciences 12, 2861–2871 (2015).
Google Scholar 
Franzén, L. G. Increased decomposition of subsurface peat in Swedish raised bogs: are temperate peatlands still net sinks of carbon? Mires Peat 1, 3 (2006).
Google Scholar 
Eilers, K. G., Lauber, C. L., Knight, R. & Fierer, N. Shifts in bacterial community structure associated with inputs of low molecular weight carbon compounds to soil. Soil Biol. Biochem. 42, 896–903 (2010).CAS 

Google Scholar 
de Graaff, M. A., Jastrow, J. D., Gillette, S., Johns, A. & Wullschleger, S. D. Differential priming of soil carbon driven by soil depth and root impacts on carbon availability. Soil Biol. Biochem. 69, 147–156 (2014).
Google Scholar 
Peay, K. G., Kennedy, P. G. & Brun, T. D. Fungal community ecology: a hybrid beast with a molecular master. BioScience 58, 799–810 (2008).
Google Scholar 
Gillabel, J., Cebrian, B., Six, J. & Merckx, R. Experimental evidence for the attenuating effect of SOM protection on temperature sensitivity of SOM decomposition. Glob. Chang. Biol. 16, 2789–2798 (2010).
Google Scholar 
Pries, C. E. H., Castanha, C., Porras, R. C. & Torn, M. S. The whole-soil carbon flux in response to warming. Science 355, 1420–1423 (2017).
Google Scholar 
Hicks Pries, C. E., Schuur, E. A. G. & Crummer, K. G. Thawing permafrost increases old soil and autotrophic respiration in tundra: partitioning ecosystem respiration using δ13C and ∆14C. Global Chang. Biol. 19, 649–661 (2013).
Google Scholar 
Tian, J. et al. Aerobic environments in combination with substrate additions to soil significantly reshape depth-dependent microbial distribution patterns in Zoige peatlands, China. Appl.Soil Ecol. 170, 104252 (2022).
Google Scholar 
Feng, W. et al. Enhanced decomposition of stable soil organic carbon and microbial catabolic potentials by long-term field warming. Glob. Chang. Biol. 00, 1–12 (2017).
Google Scholar 
Feng, W. et al. Methodological uncertainty in estimating carbon turnover times of soil fractions. Soil Biol. Biochem. 100, 118–124 (2016).CAS 

Google Scholar 
Liang, J. et al. Methods for estimating temperature sensitivity of soil organic matter based on incubation data: A comparative evaluation. Soil Biol. Biochem. 80, 127–135 (2015).CAS 

Google Scholar 
Cai, A., Feng, W., Zhang, W. & Xu, M. Climate, soil texture, and soil types affect the contributions of fine-fraction-stabilized carbon to total soil organic carbon in different land uses across China. J. Environ. Manag. 172, 2–9 (2016).CAS 

Google Scholar 
Liu, L. et al. Response of anaerobic mineralization of different depths peat carbon to warming on Zoige plateau. Geoderma 337, 1218–1226 (2019).CAS 

Google Scholar 
Waldrop, M. et al. Molecular investigations into a globally important carbon pool: permafrost protected carbon in Alaskan soils. Glob. Chang. Biol. 16, 2543–2554 (2014).
Google Scholar 
Mooshammer, M., Wanek, W., Zechmeister-Boltenstern, S. & Richter, A. Stoichiometric imbalances between terrestrial decomposer communities and their resources: mechanisms and implications of microbial adaptations to their resources. Front. Microbiol. 5, 22 (2014).
Google Scholar 
Blagodatskaya, E. & Kuzyakov, Y. Mechanisms of real and apparent priming effects and their dependence on soil microbial biomass and community structure: critical review. Biol. Fertil. Soils 45, 115–131 (2008).
Google Scholar 
Chen, H. et al. The carbon stock of alpine peatlands on the Qinghai–Tibetan Plateau during the Holocene and their future fate. Quat. Sci. Rev. 95, 151–158 (2014).
Google Scholar 
Sun, G. A study on the mineral formation law, classifictation and reserves of the peat in the Rouergai Plateau. J. Nat. Res. 7, 334–345 (1992).
Google Scholar 
Liu, L. et al. Responses of peat carbon at different depths to simulated warming and oxidizing. Sci. Total Environ. 548-549, 429–440 (2016).CAS 

Google Scholar 
Liu, L. et al. Water table drawdown reshapes soil physicochemical characteristics in Zoige peatlands. Catena 170, 119–128 (2018).CAS 

Google Scholar 
Liu, L. et al. Carbon stock stability in drained peatland after simulated plant carbon addition: Strong dependence on deeper soil. Sci. Total Environ. 848, 157539 (2022).CAS 

Google Scholar 
Yang, Z. et al. Soil properties and species composition under different grazing intensity in an alpine meadow on the eastern Tibetan Plateau, China. Environ. Monit. Assess 188, 678 (2016).
Google Scholar 
Simpson, M. J. & Simpson, A. J. The chemical ecology of soil organic matter molecular constituents. J. Chem. Ecol. 38, 768–784 (2012).CAS 

Google Scholar 
Lalonde, K., Mucci, A., Ouellet, A. & Gelinas, Y. Preservation of organic matter in sediments promoted by iron. Nature 483, 198–200 (2012).CAS 

Google Scholar 
Deforest, J. L., zak, D. R., Pregitzer, K. S. & Burtonf, A. J. Atomspheric nitrate deposition and enhanced dissolved organic carbon leaching: test of a potential mechanism. Soil Sci. Soc. Am. J. 69, 1233–1237 (2005).CAS 

Google Scholar 
Schadel, C. et al. Circumpolar assessment of permafrost C quality and its vulnerability over time using long-term incubation data. Glob. Chang. Biol. 20, 641–652 (2014).
Google Scholar 
Bell, M. & Lawrence, D. Soil carbon sequestration – myths and mysteries. Department of Primary Industries and Fisheries, Queensland Government (2009).Schadel, C., Luo, Y., David Evans, R., Fei, S. & Schaeffer, S. M. Separating soil CO2 efflux into C-pool-specific decay rates via inverse analysis of soil incubation data. Oecologia 171, 721–732 (2013).
Google Scholar 
Caporaso, J. G. et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl Acad. Sci. USA 108, 4516–4522 (2011).CAS 

Google Scholar 
Gardes, M. & Bruns, T. D. ITS primers with enhanced specificity for basidiomycetes – application to the identification of mycorrhizae and rusts. Mol. Ecol. 2, 113–118 (1993).CAS 

Google Scholar 
White, T. J. in PCR-Protocols: A Guide to Methods and Applications (Academic Press, 1990).Bell, C. et al. High-throughput fluorometric measurement of potential soil extracellular enzyme activities. J. Vis. Exp. 81, e50961 (2013).
Google Scholar 
DeForest, J. L. The influence of time, storage temperature, and substrate age on potential soil enzyme activity in acidic forest soils using MUB-linked substrates and L-DOPA. Soil Biol. Biochem. 41, 1180–1186 (2009).CAS 

Google Scholar 
Amundson, R. The carbon budget in soils. Annu. Rev. Earth Planet. Sci. 29, 535–562 (2001).CAS 

Google Scholar 
Trumbore, S. E. Potential responses of soil organic carbon to global environmental change. Proc. Natl Acad. Sci. USA 94, 8284–8291 (1997).CAS 

Google Scholar 
R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing. https://www.r-project.org (2017).Oksanen, J. et al. vegan: community ecology package. R Packag version 24-1 (2016).Asshauer, K. P., Wemheuer, B., Daniel, R. & Meinicke, P. Tax4Fun: predicting functional profiles from metagenomic 16S rRNA data. Bioinformatics 31, 2882–2884 (2015).CAS 

Google Scholar  More

Wood, E. M. & Esaian, S. The importance of street trees to urban avifauna. Ecol. Appl. 30(7), e02149 (2020).Article 

Google Scholar 
Li, Z. & Ma, J. Discussing street tree planning based on pedestrian volume using machine learning and computer vision. Build. Environ. 219, 109178 (2022).Article 

Google Scholar 
Tan, X. & Shibata, S. Factors influencing street tree health in constrained planting spaces: Evidence from Kyoto City, Japan. Urban For. Urban Green. 67, 127416 (2022).Article 

Google Scholar 
Plant, L. & Sipe, N. Adapting and applying evidence gathering techniques for planning and investment in street trees: A case study from Brisbane. Australia. Urban For. Urban Green. 19, 79–87 (2016).Article 

Google Scholar 
Dümpelmann, S. Urban trees in times of crisis: Palliatives, mitigators, and resources. One Earth 2, 402–404 (2020).Article 
ADS 

Google Scholar 
Liu, J. & Slik, F. Are street trees friendly to biodiversity?. Landsc. Urban Plan. 218, 104304 (2022).Article 

Google Scholar 
Suchocka, M. et al. Old trees are perceived as a valuable element of the municipal forest landscape. PeerJ 10, 12700 (2022).Article 

Google Scholar 
Marselle, M. R. et al. Urban Street tree biodiversity and antidepressant prescriptions. Sci. Rep. 10, 22445 (2020).Article 
ADS 
CAS 

Google Scholar 
Radu, S. The ecological role of deadwood in natural forests. In Nature Conservation. Environmental Science and Engineering (eds Gafta, D. & Akeroyd, J.) (Springer, 2006).
Google Scholar 
Piovesan, G. & Biondi, F. On tree longevity. New Phytol. 231, 1318–1337 (2021).Article 

Google Scholar 
Ferenc, M., Sedláček, O. & Fuchs, R. How to improve urban greenspace for woodland birds: Site and local-scale determinants of bird species richness. Urban Ecosyst. 17, 625–640 (2014).Article 

Google Scholar 
Birch, J. D., Lutz, J. A., Turner, B. L. & Karst, J. Divergent, age-associated fungal communities of Pinus flexilis and Pinus longaeva. For. Ecol. Manage. 494, 119277 (2021).Article 

Google Scholar 
Siitonen, J., Ranius, T. The importance of veteran trees for saproxylic insects. In Europe’s Changing Woods and Forests: From Wildwood to Managed Landscapes (2015).Polyakov, A. Y., Weller, T. J. & Tietje, W. D. Remnant trees increase bat activity and facilitate the use of vineyards by edge-space bats. Agr. Ecosyst. Environ. 281, 56–63 (2019).Article 

Google Scholar 
Hall, S. J. G. & Bunce, R. G. H. Mature trees as keystone structures in Holarctic ecosystems – a quantitative species comparison in a northern English park. Plant Ecol. Divers. 4, 243–250 (2011).Article 

Google Scholar 
Suchocka, M. et al. Transit versus Nature. Depreciation of environmental values of the road alleys. Case study: Gamerki-Jonkowo, Poland. Sustain. 11(6), 1816 (2019).Article 

Google Scholar 
What Are Ancient & Veteran Trees. Ancient Tree Forum | Championing the Biological, Cultural And Heritage Value Of The UK’s Ancient Trees. URL https://www.ancienttreeforum.org.uk/ancient-trees/what-are-ancient-veteran-trees/ (2022).Fay, N. Environmental arboriculture, tree ecology and veteran tree management. Arbor. J. 26, 213–236 (2002).Article 

Google Scholar 
Dujesiefken, D., Fay, N., De Groot, J. W. & De Berker, N. Trees—a lifespan approach. Contributions to arboriculture from European practitioners (eds. Witkoś-Gnach, K., Tyszko-Chmielowiec, P.) (Fundacja EkoRozwoju, 2016).Roman, L. How many trees are enough? Tree death and the urban canopy. Scenar. J. 04, 8 (2014).
Google Scholar 
Roman, L. A. & Scatena, F. N. Street tree survival rates: Meta-analysis of previous studies and application to a field survey in Philadelphia, PA, USA. Urban For. Urban Green. 10(4), 269–274 (2011).Article 

Google Scholar 
Czaja, M., Kołton, A. & Muras, P. The complex issue of urban trees—stress factor accumulation and ecological service possibilities. Forests 11, 932 (2020).Article 

Google Scholar 
Olchowik, J., Suchocka, M., Jankowski, P., Malewski, T. & Hilszczańska, D. The ectomycorrhizal community of urban linden trees in Gdańsk, Poland. PlosOne. 16(4), e0237551 (2021).Article 
CAS 

Google Scholar 
Nilsson, K., Konijnendijk, C. C. & Nielsen, A. B. Urban forest function, design and management. In Encyclopedia of Sustainability Science and Technology (ed. Meyers, R. A.) https://doi.org/10.1007/978-1-4419-0851-3_218 (Springer, New York, NY, 2013).Chapter 

Google Scholar 
Pokorny, J.D. Urban tree risk management, a Community Guide to Program Design and Implementation. USDA Forest Service Northeastern Area State and Private Forestry (2003).James, K. R., Haritos, N. & Ades, P. K. Mechanical stability of trees under dynamic loads. Am. J. Bot. 93(10), 1361–1369 (2006).Article 

Google Scholar 
Hickman, G. W., Perry, E. & Evans, R. Validation of a tree failure evaluation system. J. Arboric. 21(5), 233–234 (1995).
Google Scholar 
Klein, R., Koeser, A., Hauer, R., Hansen, G. & Escobedo, F. Risk assessment and risk perception of trees: A review of literature relating to arboriculture and urban forestry. Arboric. Urban For. 45(1), 26–38 (2019).
Google Scholar 
Smiley, E. T. Root pruning and stability of young willow oak. Arboric. Urban For. 34(2), 123–128 (2008).Article 

Google Scholar 
Terho, M. & Hallaksela, A.-M. Decay characteristics of hazardous Tilia, Betula, and Acer trees felled by municipal urban tree managers in the Helsinki city area. Forestry 81(2), 151–159. https://doi.org/10.1093/forestry/cpn002 (2008).Article 

Google Scholar 
Terho, M. An assessment of decay among urban Tilia, Betula, and Acer trees felled as hazardous. Urban For. Urban Green. 8, 77–85 (2009).Article 

Google Scholar 
Koeser, A. K., Klein, R. W., Hasing, G. & Northrop, R. J. Factors driving professional and public urban tree risk perception. Urban For. Urban Green. 14(4), 968–974 (2015).Article 

Google Scholar 
Johnson, G. R. Storms over Minnesota. Minn. Shade Tree Advocate 2(1), 1–12 (1999).ADS 

Google Scholar 
Zhang, Y., Hussain, A., Deng, J. & Letson, L. Public attitudes toward urban trees and supporting urban tree programs. Environ. Behav. 39(6), 797–814 (2007).Article 

Google Scholar 
Suchocka, M., Swoczyna, T., Kosno-Jończy, J. & Kalaji, H. M. Impact of heavy pruning on development and photosynthesis of Tilia cordata Mill Trees. PLoS ONE 16(8), e0256465. https://doi.org/10.1371/journal.pone.0256465 (2021).Article 
CAS 

Google Scholar 
Gilman, E. F. & Knox, G. Pruning type affects ecay and structure of crape myrtle. J. Arboric. 31, 38–47 (2005).
Google Scholar 
Gilman, E. F. & Lilly, S. J. Best Management Practices: Tree Pruning (International Society of Arboriculture, 2008).
Google Scholar 
Perrette, G., Delagrange, S., Ramirez, J. A. & Messier, C. Optimisingreduction pruning under electrical lines: The influence of tree vitality before pruning on traumatic responses. Urban For. Urban Green. 63, 127139 (2021).Article 

Google Scholar 
von Döhren, P. & Haase, D. Risk assessment concerning urban ecosystem disservices: The example of street trees in Berlin. Germany. Ecosyst. Serv. 40, 101031 (2019).Article 

Google Scholar 
Papandrea, S. F., Cataldo, M. F., Zimbalatti, G. & Proto, A. R. Comparative evaluation of inspection techniques for decay detection in urban trees. Environ. Sci. Proc. 3, 14 (2021).
Google Scholar 
McPherson, G. & Peper, P. P. Costs of street tree damage to infrastructure. Arbor. J. 20, 143–160 (1996).Article 

Google Scholar 
Mullaney, J., Lucke, T. & Trueman, S. J. A review of benefits and challenges in growing street trees in paved urban environments. Landsc. Urban Plan. 134, 157–166 (2015).Article 

Google Scholar 
Vogt, J., Hauer, R. J. & Fischer, B. C. The costs of maintaining and not maintaining the urban forest: A review of the urban forestry and arboriculture literature. Arboric. Urban For. 41(6), 293–323 (2015).
Google Scholar 
Mattheck, C. & Breloer, H. Field guide for visual tree assessment (VTA). Arboric. J. 18(1), 1–23 (1994).Article 

Google Scholar 
Smiley E.T., Matheny N., & Lilly S. Best management practices: Tree risk assessment. In International Society of Arboriculture, 86 (Champaign, Illinois, 2011).Dunster J.A., Smiley E.T., Matheny N., Lilly S. Tree risk assessment manual. International Society of Arboriculture 194 (Champaign, Illinois, 2013).Li, H., Zhang, X., Li, Z., Wen, J. & Tan, X. A review of research on tree risk assessment methods. Forests 13, 1556 (2022).Article 

Google Scholar 
Koeser, A. K., Hauer, R. J., Klein, R. W. & Miesbauer, J. W. Assessment of likelihood of failure using limited visual, basic, and advanced assessment techniques. Urban For. 24, 71–79 (2017).
Google Scholar 
TRAQ [URL TRAQandOtherTreeRiskAssessmentMethodsforEvaluationandPrioritizingTreeRiskConditions(forestmetrix.com) (2021).TRAQ Tree Risk Assessment Qualification Application Guide https://www.isa-arbor.com/Portals/0/Assets/PDF/Certification-Applications/TRAQ-App-Guide.pdf (2021).Matheny N. P., Clark J. R. A photographic guide to the evaluation of hazard trees in urban areas. In International Society of Arboriculture 85 (Champaign, 1994).Linhares, C. S. F., Gonçalves, R., Martins, L. M. & Knapic, S. Structural stability of urban trees using visual and instrumental techniques: A review. Forests 12, 1752. https://doi.org/10.3390/f12121752 (2021).Article 

Google Scholar 
Ellison, M. Quantified tree risk assessment: Nota De procedimiento V5.2.3 (ES)2018-01 Quantified Tree Risk Assessment Limited (2018).Forbes-Laird, J. THREATS – tree hazard risk evaluation and treatment system – Guidance note for users Retrieved March 27th, 2020 from Forbes-Laird Arboricultural Consultancy http://www.flac.uk.com/wp-content/uploads/2010/07/THREATS-GN-June-2010.pdf, (2010).Guyon C. Cleaver M. Jackson A. Saavedra P. Zambino A. Guide to Identifying, Assessing, and Managing Hazard Trees in Developed Recreational Sites of the Northern Rocky Mountains and the Intermountain West Retrieved March 31st, 2020 from USDA Forest Service, Northern and Intermountain Regions (2017). https://www.fs.usda.gov/Internet/FSE_DOCUMENTS/fseprd571021.pdfBlodgett, J. T., Burns, K. S., Worrall J. J.Guide to hazard tree management Retrieved March 31st, 2020 from USDA Forest Service, Rocky Mountain Region (2017) https://www.fs.usda.gov/Internet/FSE_DOCUMENTS/fseprd572690.pdf (2017).Norris M. A review of methods used to undertake risk assessments of urban trees. MSc. Thesis (2007).Smiley, E. T., Matheny, N., Lilly, S. Best management practices: Tree risk assessment. International Society of Arboriculture 86 (Champaign, Illinois, 2011).ALARP – Hart, A, 2013, ALARP – Recent Developments, ALARP: Learning from the Experiences of Others, London: IMechE, 4th June 2013 (2013).HSE, 2001 Reducing risks, protecting people, HSE’s decision making process, Liverpool: Health and Safety Executive. (2001).Rinn, F. Holzanatomische Grundlagen mechanischer impuls – Tomographie an Baumen [Wood anatomy background through mechanical pulses – tomografy of trees]. Allg. Forstwirtsch. 8, 450–456 (2003).
Google Scholar 
Gilbert, E. A. & Smiley, E. T. Picus sonic tomography for the quantification of decay in white oak (Quercus alba) and hickory (Carya spp.). J. Arboric 30, 277–281 (2004).
Google Scholar 
Wang, X. & Allison, R. B. Decay detection in red oak trees using a combination of visual inspection, acoustic testing, and resistance microdrilling. Arboric. Urban For. 34(1), 1–4 (2008).Article 

Google Scholar 
Wu, Y. & Shao, Z. Measurement and mechanical analysis of the strains–stresses induced by tree-pulling experiments in tree stems. Trees 30, 675–684 (2016).Article 

Google Scholar 
Schindler, D. & Kolbe, S. Assessment of the response of a scots pine tree to effective wind loading. Forests 11(2), 145 (2020).Article 

Google Scholar 
Koeser, A. K. & Smiley, E. T. Impact of assessor on tree risk assessment ratings and prescribed mitigation measures. Urban For. 24, 109–115 (2017).
Google Scholar 
Klein, R. W. et al. Assessing the consequences of tree failure. Urban Forestry & Urban Greening 65, 127307 (2021).Article 

Google Scholar 
Renn, O. Perception of risks. Toxicol. Lett. 149(1), 405–413 (2004).Article 
CAS 

Google Scholar 
Hasan, R., Othman, N. & Ismail, F. Roadside tree management in urban area for public safety and properties. Asian J. Quality Life 3, 10–21834 (2018).Article 

Google Scholar 
Williams, V. How do You Decide When to Remove a Tree? (University Of Maryland extension, 2018).Rhoades, H. Filling holes in tree trunks: how to patch a hole in a tree trunk or a hollow tree. https://www.gardeningknowhow.com/ornamental/trees/tgen/patching-tree-hole.htm (2020).Terho, M. & Hallaksela, A. M. Potential hazard characteristics of Tilia, Betula, and Acer trees removed in the Helsinki City Area during 2001–2003. Urban For. Urban Green. 3, 113–120 (2005).Article 

Google Scholar 
Nagendra, H. & Gopal, D. Tree diversity, distribution, history and change in urban parks: Studies in Bangalore India. Urban Ecosyst. 14, 211–223 (2011).Article 

Google Scholar 
Lindenmayer, D. B., Blanchard, W., Blair, D. & McBurney, L. The road to oblivion – Quantifying pathways in the decline of large old trees. For. Ecol. Manage. 430, 259–264 (2018).Article 

Google Scholar 
Lusk, A. C., da Silva Filho, D. F. & Dobbert, L. Pedestrian and cyclist preferences for tree locations by sidewalks and cycle tracks and associated benefits: Worldwide implications from a study in Boston, MA. Cities 106, 102111 (2020).Article 

Google Scholar 
Galenieks, A. Importance of urban street tree policies: A comparison of neighboring southern California Cities. Urban For. Urban Green. 22, 105–110 (2017).Article 

Google Scholar 
Wessolly, L. Material and structural features of trees Contribution to the Stargardt strength catalogue. In Proceedings of the 15th Bad Goteborg Tree Seminar (1992).Schwarze, F. Diagnosis and prognosis of the development of wood decay in urban trees. Agrios GN 1997 Plant Patology. (Academic Press, San Diego, 2008).Footway. Cycling Embassy Of Great Britain [https://www.cycling-embassy.org.uk/dictionary/footway] (2022).Roloff, A. Handbuch Baumdiagnostik Baum-Korpersprache und Baum-Beurtailung (Ulmer Verlag, 2015).
Google Scholar 
Koeser, A. K., Hasing, G., McLean, D., Northrop R. Tree risk assessment methods: A comparison of three common evaluation forms Retrieved March 24th, 2020 from https://edis.ifas.ufl.edu/ep487 (2016).Smiley, E. T. & Kumamoto, H. Qualitative Tree Risk Assessment. 12–18 (2012).Mattheck, C. Trees: The Mechanical Design (Springer, 1991).Book 

Google Scholar 
R Core Team (2020) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.r-project.org/ (2020).Olchowik, J. et al. The ectomycorrhizal community of crimean linden trees in Warsaw, Poland. Forests 11(9), 926 (2020).Article 

Google Scholar 
Dupre, S., Thiebaut, B. & Tessier du Cros, E. Morphologie architecture des jeunes hfitres (Fagus sylvatica L.). Influence du milieu variability genetique. Ann. Sci. For. 43, 85–102 (1986).Article 

Google Scholar 
Power, S. A., Ashmore, M. R. & Ling, K. A. Recent trends in beech tree health in southern Britain and the influence of soil type. Water Air Soil Pollut. 85, 1293–1298 (1995).Article 
ADS 
CAS 

Google Scholar 
Masarovičova, E. & Štefančik, L. Some ecophysiological features in sun and shade leaves of tall beech trees. Biol. Plant 32, 374–387 (1990).Article 

Google Scholar 
Nicolini, E. & Caraglio, Y. L’influence de divers caracteres architecturaux sur l’apparition de la fourche chez le Fagus sylvatica, en fonction de l’absence ou de la presence d’un couvert. Botany 72, 1723–1734 (1994).
Google Scholar 
van Wassenaer, P. V. & Richardson, M. A review of tree risk assessment using minimally invasive technologies and two case studies. Arboric. J. 32, 275–292 (2009).Article 

Google Scholar 
dos Reis, M. N., Gonçalves, R., Brazolin, S. & de Assis Palma, S. S. Strength loss inference due to decay or cavities in tree trunks using tomographic imaging data applied to equations proposed in the literature. Forests 13, 596 (2022).Article 

Google Scholar 
Kanea, B., Warrena, P. S. & Lermanab, S. B. A broad scale analysis of tree risk, mitigation and potential habitat for cavity-nesting birds. Urban For. 14, 1137–1146 (2015).
Google Scholar 
Wolf, K. L. Roadside urban trees—balancing safety and community values. Arborist News 15, 25–27 (2006).
Google Scholar 
Hightshoe, G. L. Native Trees, Shrubs and Vines for Urban and Rural America (Wiley and Sons, 1988).
Google Scholar 
Costello, L. R. & Jones, K. S. Western chapter of the international society of arboriculture. In Reducing Infrastructure Damage by The Tree Roots: A Compendium of Strategies. 64–65 (2003).Kjaer, E. D. Introduction part 2. Consequences of ash dieback: Damage level, resistance and resilience of European Ash Forests. Balt. For. 23, 141–143 (2017).
Google Scholar 
Timmermann, V., Nagy, N., Hietala, A., Børja, I. & Solheim, H. Progression of ash dieback in Norway related to tree age, disease history and regional aspects. Balt. For. 23, 150–158 (2017).
Google Scholar 
Zajączkowska, U., Kaczmarczyk, K. & Liana, J. Birch sap exudation: influence of tree position in a forest stand on birch sap production, trunk wood anatomy and radial bending strength. Silva Fennica 53(2), 10048. https://doi.org/10.14214/sf.10048 (2019).Article 

Google Scholar 
Reed, H. J. Veteran Trees: A Guide to Good Management (England Nature, 2000).
Google Scholar  More