More stories

  • in

    Potential negative effects of the installation of video surveillance cameras in raptors’ nests

    Ribic, C. A., Thompson, F. R. & Pietz, P. J. Video Surveillance of Nesting Birds (University of California Press, 2012).Book 

    Google Scholar 
    O’Brien, T. G. & Kinnaird, M. F. A picture is worth a thousand words: The application of camera trapping to the study of birds. Bird Conserv. Int. 18, S144–S162 (2008).Article 

    Google Scholar 
    Kristan, D. M., Golightly Jr, R. T. & Tomkiewicz Jr, S. M. A solar-powered transmitting video camera for monitoring raptor nests. Wildl. Soc. Bull. 24, 284–290 (1996).
    Google Scholar 
    Grubb, T. An infrared video camera system for monitoring diurnal and nocturnal raptors. J. Raptor Res. 32, 290–296 (1998).
    Google Scholar 
    Margalida, A. et al. A solar-powered transmitting video camera for monitoring cliff-nesting raptors. J. Field Ornithol. 77, 7–12 (2006).Article 

    Google Scholar 
    Bolton, M., Butcher, N., Sharpe, F., Stevens, D. & Fisher, G. Remote monitoring of nests using digital camera technology. J. Field Ornithol. 78, 213–220 (2007).Article 

    Google Scholar 
    Pierce, A. J. & Pobprasert, K. A portable system for continuous monitoring of bird nests using digital video recorders. J. Field Ornithol. 78, 322–328 (2007).Article 

    Google Scholar 
    Benson, T. J., Brown, J. D. & Bednarz, J. C. Identifying predators clarifies predictors of nest success in a temperate passerine. J. Anim. Ecol. 79, 225–234 (2010).Article 

    Google Scholar 
    Lewis, S. B., Fuller, M. R. & Titus, K. A comparison of 3 methods for assessing raptor diet during the breeding season. Wildl. Soc. Bull. 32, 373–385 (2004).Article 

    Google Scholar 
    Rogers, A. S. Quantifying Northern Goshawk diets using remote cameras and observations from blinds. J. Raptor Res. 39, 303–309 (2005).ADS 

    Google Scholar 
    Tornberg, R. & Reif, V. Assessing the diet of birds of prey: A comparison of prey items found in nests and images. Ornis Fenn. 84, 21 (2007).
    Google Scholar 
    López-López, P. & Urios, V. Use of digital trail cameras to study Bonelli’s eagle diet during the nestling period. Ital. J. Zool. 77, 289–295 (2010).Article 

    Google Scholar 
    Harrison, J. T., Kochert, M. N., Pauli, B. P. & Heath, J. A. Using motion-activated trail cameras to study diet and productivity of cliff-nesting Golden Eagles. J. Raptor Res. 53, 26–37 (2019).Article 

    Google Scholar 
    McRae, S. B., Weatherhead, P. J. & Montgomerie, R. American robin nestlings compete by jockeying for position. Behav. Ecol. Sociobiol. 33, 101–106 (1993).Article 

    Google Scholar 
    Nathan, A., Legge, S. & Cockburn, A. Nestling aggression in broods of a siblicidal kingfisher, the laughing kookaburra. Behav. Ecol. 12, 716–725 (2001).Article 

    Google Scholar 
    Grivas, C. et al. An audio–visual nest monitoring system for the study and manipulation of siblicide in bearded vultures Gypaetus barbatus on the island of Crete (Greece). J. Ethol. 27, 105–116 (2009).Article 

    Google Scholar 
    Gula, R., Theuerkauf, J., Rouys, S. & Legault, A. An audio/video surveillance system for wildlife. Eur. J. Wildl. Res. 56, 803–807 (2010).Article 

    Google Scholar 
    Sanaiotti, T. M., Seixas, G. H., Duleba, S. & Martins, F. D. Camera trapping at harpy eagle nests: Interspecies interactions under predation risk. J. Raptor Res. 51, 72–78 (2017).Article 

    Google Scholar 
    Allen, M. L., Inagaki, A. & Ward, M. P. Cannibalism in raptors: A review. J. Raptor Res. 54, 424–430 (2020).Article 

    Google Scholar 
    Academia, M. H. & Dalgleish, H. J. Use of nest web cameras and citizen science to quantify osprey prey delivery rate and nest success. J. Raptor Res. 56, 212–219 (2022).Article 

    Google Scholar 
    Gysel, L. W. & Davis, E. M. A simple automatic photographic unit for wildlife research. J. Wildl. Manag. 20, 451–453 (1956).Article 

    Google Scholar 
    Royama, T. A device of an auto-cinematic food-recorder. Jpn. J. Ornithol. 15, 172–176 (1959).Article 

    Google Scholar 
    Cox, W. A. et al. Development of camera technology for monitoring nests. In Chapter 15. Video Surveill. Nesting Birds Stud. Avian Biol. (eds Ribic, C. A. et al.) 185–210 (Univ. Calif. Press, 2012).Sanders, M. D. & Maloney, R. F. Causes of mortality at nests of ground-nesting birds in the Upper Waitaki Basin, South Island, New Zealand: A 5-year video study. Biol. Conserv. 106, 225–236 (2002).Article 

    Google Scholar 
    Reif, V. & Tornberg, R. Using time-lapse digital video recording for a nesting study of birds of prey. Eur. J. Wildl. Res. 52, 251–258 (2006).Article 

    Google Scholar 
    McKinnon, L. & Bêty, J. Effect of camera monitoring on survival rates of High-Arctic shorebird nests. J. Field Ornithol. 80, 280–288 (2009).Article 

    Google Scholar 
    Powell, L. A., Giovanni, M. D., Groepper, S. R., Reineke, M. & Schacht, W. H. Attendance Patterns and Survival of Western Meadowlark Nests (University of California Press, 2012).Book 

    Google Scholar 
    Herranz, J., Yanes, M. & Suárez, F. Does Photo-Monitoring Affect Nest Predation? J. Field Ornithol. 73, 97–101 (2002).Article 

    Google Scholar 
    Richardson, T. W., Gardali, T. & Jenkins, S. H. Review and meta-analysis of camera effects on avian nest success. J. Wildl. Manag. 73, 287–293 (2009).Article 

    Google Scholar 
    Cain, S. L. Nesting activity time budgets of Bald Eagles in southeast Alaska. (1985).García-Salgado, G. et al. Evaluation of trail-cameras for analyzing the diet of nesting raptors using the Northern Goshawk as a model. PLoS ONE 10, e0127585 (2015).Article 

    Google Scholar 
    Swann, D. E., Kawanishi, K. & Palmer, J. Evaluating types and features of camera traps in ecological studies: a guide for researchers. In Camera Traps in Animal Ecology 27–43 (Springer, 2011).Dykstra, C., Meyer, M. & Warnke, D. Bald Eagle reproductive performance following video camera placement. J. Raptor Res. 36, 136–139 (2002).
    Google Scholar 
    Del Moral, J. C. & Molina, B. El águila perdicera en España, población reproductora en 2018 y método de censo (SEO/BirdLife, 2018).
    Google Scholar 
    Generalitat Valenciana. Orden 2/2022, de 16 de febrero, de la Conselleria de Agricultura, Desarrollo Rural, Emergencia Climática y Transición Ecológica, por la que se actualizan los listados valencianos de especies protegidas de flora y fauna (2022).Real, J. & Mañosa, S. Demography and conservation of western European Bonelli’s eagle Hieraaetus fasciatus populations. Biol. Conserv. 79, 59–66 (1997).Article 

    Google Scholar 
    Hernández-Matías, A. et al. From local monitoring to a broad-scale viability assessment: A case study for the Bonelli’s Eagle in western Europe. Ecol. Monogr. 83, 239–261 (2013).Article 

    Google Scholar 
    Rollan, A. et al. Guiding local-scale management to improve the conservation of endangered populations: The example of Bonelli’s Eagle Aquila fasciata. Bird Conserv. Int. 31, 395–409 (2021).Article 

    Google Scholar 
    López-López, P., García-Ripollés, C. & Urios, V. Population size, breeding performance and territory quality of Bonelli’s Eagle Hieraaetus fasciatus in eastern Spain. Bird Study 54, 335–342 (2007).Article 

    Google Scholar 
    López-López, P. Informe científico valoración de la inclusión del águila perdicera como especie en peligro de extinción en el Catálogo Valenciano de Especies de Fauna Amenazadahttps://doi.org/10.13140/RG.2.2.32806.04166 (2021).López-López, P., Perona, A., Egea-Casas, O., Morant, J. & Urios, V. Tri-axial accelerometry shows differences in energy expenditure and parental effort throughout the breeding season in long-lived raptors. Curr. Zool. 68, 57–67 (2022).Article 

    Google Scholar 
    Perona, A. M., Urios, V. & López-López, P. Holidays? Not for all Eagles have larger home ranges on holidays as a consequence of human disturbance. Biol. Conserv. 231, 59–66 (2019).Article 

    Google Scholar 
    Morollón, S., Urios, V. & López-López, P. Fifteen days are enough to estimate home-range size in some long-lived resident eagles. J. Ornithol. 163, 849–854 (2022).Article 

    Google Scholar 
    Stewart-Oaten, A., Murdoch, W. W. & Parker, K. R. Environmental impact assessment:” Pseudoreplication” in time?. Ecology 67, 929–940 (1986).Article 

    Google Scholar 
    Underwood, A. Beyond BACI: The detection of environmental impacts on populations in the real, but variable, world. J. Exp. Mar. Biol. Ecol. 161, 145–178 (1992).Article 

    Google Scholar 
    López-López, P., García-Ripollés, C., García-López, F., Aguilar, J. M. & Verdejo, J. Patrón de distribución del águila real Aquila chrysaetos y del águila-azor perdicera Hieraaetus fasciatus en la provincia de Castellón. Ardeola 51, 275–283 (2004).
    Google Scholar 
    López-López, P., García-Ripollés, C., Aguilar, J. M., Garcia-López, F. & Verdejo, J. Modelling breeding habitat preferences of Bonelli’s eagle (Hieraaetus fasciatus) in relation to topography, disturbance, climate and land use at different spatial scales. J. Ornithol. 147, 97–106 (2006).Article 

    Google Scholar 
    Gil-Sánchez, J. Effects of altitude and prey availability on the laying date of Bonelli’s Eagles (Hieraaetus fasciatus) in Granada (SE Spain). Ardeola 47, 1–8 (2000).
    Google Scholar 
    Forsman, D. Flight Identification of Raptors of Europe, North Africa and the Middle East (Bloomsbury Publishing, 2016).
    Google Scholar 
    Zuur, A. F., Ieno, E. N., Walker, N. J., Saveliev, A. A. & Smith, G. M. Mixed Effects Models and Extensions in Ecology with R Vol. 574 (Springer, 2009).Book 
    MATH 

    Google Scholar 
    Harrison, X. A. et al. A brief introduction to mixed effects modelling and multi-model inference in ecology. PeerJ 6, e4794 (2018).Article 

    Google Scholar 
    Brooks, M. E. et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 9, 378–400 (2017).Article 

    Google Scholar 
    Nakagawa, S. & Schielzeth, H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol. Evol. 4, 133–142 (2013).Article 

    Google Scholar 
    Barton, K. MuMIn: Multi-Model Inference. R package version 1.46.0. (2022).Cutler, T. L. & Swann, D. E. Using remote photography in wildlife ecology: A review. Wildl. Soc. Bull. 27, 571–581 (1999).
    Google Scholar 
    Richardson, C. T. & Miller, C. K. Recommendations for protecting raptors from human disturbance: A review. Wildl. Soc. Bull. 25, 634–638 (1997).
    Google Scholar 
    Balbontin, J., Penteriani, V. & Ferrer, M. Variations in the age of mates as an early warning signal of changes in population trends? The case of Bonelli’s eagle in Andalusia. Biol. Conserv. 109, 417–423 (2003).Article 

    Google Scholar 
    Martinez, J. A. et al. Breeding performance, age effects and territory occupancy in a Bonelli’s eagle Hieraaetus fasciatus population. Ibis 150, 223–233 (2008).Article 

    Google Scholar 
    Sánchez-Zapata, J., Calvo, J., Carrete, M. & Martínez, J. Age and breeding success of a Golden Eagle Aquila chrysaetos population in southeastern Spain. Bird Study 47, 235–237 (2000).Article 

    Google Scholar 
    Ferrer, M., Penteriani, V., Balbontin, J. & Pandolfi, M. The proportion of immature breeders as a reliable early warning signal of population decline: Evidence from the Spanish imperial eagle in Donana. Biol. Conserv. 114, 463–466 (2003).Article 

    Google Scholar 
    Cano, A. & Parrinder, E. Studies of less familiar birds, Bonelli’s Eagle. Br. Birds 54, 422–427 (1961).
    Google Scholar 
    Blondel, J., Coulon, L., Girerd, B. & Hortigue, M. Deux cents heures d’observation aupre‘s de l’aire de l’Aigle de Bonelli Hieraaetus fasciatus. Nos Oiseaux 30, 37–60 (1969).
    Google Scholar 
    Vaucher, C. Notes sur 1’ethologie de I’Aigle de Bonelli. Nos Oiseaux 31, 101–111 (1971).
    Google Scholar 
    Elósegui, J. Informe preliminar sobre alimentación de aves rapaces en Navarra y provincias limítrofes. Ardeola 19, 249–256 (1974).
    Google Scholar 
    Cheylan, G. L. place trophique de l’Aigle de Bonelli Hieraaetus fasciatus dans les biocénoses méditerranéennes. Alauda 45, 1–15 (1977).
    Google Scholar 
    Palma, L., Cancela da Fonseca, L. & Oliveira, L. L’alimentation de l’aigle de Bonelli Hieraaetus fasciatus dans la coˆte portugaise. Rapinyaires Mediterranis 2, 87–96 (1984).
    Google Scholar 
    Real, J. Biases in diet study methods in the Bonelli’s eagle. J. Wildl. Manag. 60, 632–638 (1996).Article 

    Google Scholar 
    Gil-Sánchez, J. M., Molino, F., Valenzuela, G. & Moleón, M. Demografía y alimentación del Águila-azor Perdicera (Hieraaetus fasciatus) en la provincia de Granada. Ardeola 47, 69–75 (2000).
    Google Scholar 
    Ontiveros, D., Pleguezuelos, J. M. & Caro, J. Prey density, prey detectability and food habits: The case of Bonelli’s eagle and the conservation measures. Biol. Conserv. 123, 19–25 (2005).Article 

    Google Scholar 
    Moleón, M. et al. Large-scale spatio-temporal shifts in the diet of a predator mediated by an emerging infectious disease of its main prey. J. Biogeogr. 36, 1502–1515 (2009).Article 

    Google Scholar 
    Resano-Mayor, J. et al. Diet–demography relationships in a long-lived predator: From territories to populations. Oikos 125, 262–270 (2016).Article 

    Google Scholar 
    Di Vittorio, M. et al. Long-term changes in the breeding period diet of Bonelli’s eagle (Aquila fasciata) in Sicily, Italy. Wildl. Res. 46, 409–414 (2019).Article 

    Google Scholar  More

  • in

    Global predictions for the risk of establishment of Pierce’s disease of grapevines

    Thermal requirements to develop PDWe examined the response of a wide spectrum of European grapevine varieties to XfPD infection in three independent experiments conducted in 2018, 2019, and 2020. Overall, 86.1% (n = 764) of 886 inoculated plants, comprising 36 varieties and 57 unique scion/rootstock combinations, developed PD symptoms 16 weeks after inoculation. European V. vinifera varieties exhibited significant differences in their susceptibility to XfPD (Supplementary Table S1). All varieties, however, showed PD symptoms to some extent, confirming previous field observations of general susceptibility to XfPD9,12,37. We also found significant differences in virulence (χ2 = 68.73, df = 1, P = 2.2 × 10−16) between two XfPD strains isolated from grapevines in Majorca across grapevine varieties (Supplementary Fig. S1). Full details on the results of the inoculation tests are available in “Methods”, Supplementary Note 1, Supplementary Table S1 and Supplementary Data 1.Growing degree days (GDD) have traditionally been used to describe and predict phenological events of plants and insect pests, but rarely in plant diseases58. We took advantage of data collated in the inoculation trials together with temperature to relate symptom development to the accumulated heat units at weeks eight, 10, 12, 14, and 16 after inoculation (Supplementary Data 1). Rather than counting GDDs linearly above a threshold temperature, we consider Xf ’s specific growth rate in vitro within its cardinal temperatures. The empirical growth rates come from the seminal work by Feil & Purcell38 shown in the inset of Fig. 1a. This Arrhenius plot was transformed, as explained in Supplementary Note 2A, to obtain a a piece-wise function f(T) Eq. (1). Our model and risk maps are based on f(T) (red line in Fig. 1a) because it provides the best fit to the experimental data when compared with the commonly used Beta function (blue line) for representing the thermal response in biological processes59,60. This Modified Growing Degree Day (MGDD) profile Eq. (1) enables to measure the thermal integral from hourly average temperatures, improving the prediction scale of the biological process61. MGDD also provides an excellent metric to link XfPD growth in culture with PD development as, once the pathogen is injected into the healthy vine, symptoms progression mainly depends upon the bacterial load (i.e., multiplication) and the movement through the xylem vessel network, which are fundamentally temperature-dependent processes38,62. Moreover, MGDD can be mathematically related to the exponential or logistic growth of the pathogen within the plant (Supplementary Note 2B).Fig. 1: Climatic and transmission layers composing the epidemiological model.a MGDD profile fitted to the in vitro data of Xf growth rate in Feil & Purcell 200138. The original Arrhenius plot in Kelvin degrees (inset) was converted to Celsius, as explained in (Supplementary Note 2A), to obtain the fit shown in the main plot red line; the blue line represents the fit with a Beta function. b Correlation between CDD and the average ({T}_{min }) of the coldest month between 1981 and 2019. Plotted black dots (worldwide) and yellow dots (main wine-producing zones) depict climatic data from 6,487,200 cells at 0.1∘ × 0.1∘ resolution, spread globally and retrieved from ERA5-Land dataset. The red solid line depicts the fitted exponential function for worldwide data and the blue solid line for main vineyard zones. c Nonlinear relationship between MGDD (red line) and CDD (blue line) and the likelihood of developing chronic infections. Black dots depict the cumulative proportion of grapevine plants in the population of 36 inoculated varieties showing five or more symptomatic leaves at each of the 15 MGDD levels (see Supplementary Information). Vertical bars are the 95% CI. d Combined ranges of MGDD and CDD on the likelihood of developing chronic infection. e Transmission layer in the dynamic equation (1) of the SIR compartmental model. f Relationship between the exponential growth of the number of infected plants with the risk index and their ranks.Full size imageInterannual infection survival in grapevines plays a relevant role when modelling PD epidemiology. In our model, we assumed a threshold of five or more symptomatic leaves for these chronic infections based on the relationship between the timing and severity of the infection during the growing season and the likelihood of winter recovery38,39,42. This five-leaf cut-off was grounded on: (i) the bimodal distribution in the frequency of the number of symptomatic leaves among the population of inoculated grapevines (Supplementary Fig. S1), whereby vines that generally show less than five symptomatic leaves at 12 weeks after inoculation remain so in the following weeks, while those that pass that threshold continue to produce symptomatic leaves, and (ii) the observed correlation between the acropetal and basipetal movement of Xf along the cane (Supplementary Fig. S1). The likelihood of developing chronic infections as a function of accumulated MGDD among the population of grapevine varieties was modelled using survival analysis with data fitted to a logistic distribution ({{{{{{{mathcal{F}}}}}}}}({{{{{rm{MGDD}}}}}})). A minimum window of MGDD = 528 was needed to develop chronic infections (var. Tempranillo), about 975 for a median estimate, while a cumulative MGDD  > 1159 indicate over 90% probability within a growing season (red curve in Fig. 1c and “Methods”).Next, we intended to model the probability of disease recovery by exposure to cold temperatures. Previous works had specifically modelled cold curing on Pinot Noir and Cabernet Sauvignon varieties in California as the effect of temperature and duration39 by assuming a progressive elimination of the bacterial load with cold temperatures42. In the absence of appropriate empirical data to formulate a general average pattern of winter curing among grapevine varieties, we combined the approach of Lieth et al.39 and the empirical observations of Purcell on the distribution of PD in the US related to the average minimum temperature of the coldest month, Tmin, isolines41. To consider the accumulation of cold units in an analogy of the MGDD, we searched for a general correlation between Tmin and the cold degree days (CDDs) with base temperature = 6 ∘C (see “Methods”). We found an exponential relation, ({{{{{rm{CDD}}}}}} sim 230exp (-0.26cdot {T}_{min })), where specifically, CDD ≳ 306 correspond to ({T}_{min } < -1.{1},^{circ }{{{{{rm{C}}}}}}) (Fig. 1b). To transform this exponential relationship to a probabilistic function analogous to ({{{{{{{mathcal{F}}}}}}}}({{{{{rm{MGDD}}}}}})), hereafter denoted ({{{{{{{mathcal{G}}}}}}}}({{{{{rm{CDD}}}}}})), ranging between 0 and 1, we considered the sigmoidal family of functions (f(x)=frac{A}{B+{x}^{C}}) with A = 9 × 106, B = A and C = 3 (Fig. 1c), fulfilling the limit ({{{{{{{mathcal{G}}}}}}}}({{{{{rm{CDD}}}}}}=0)=1), i.e., no winter curing when no cold accumulated, and a conservative 75% of the infected plants recovered at ({T}_{min }=-1.{1},^{circ }{{{{{rm{C}}}}}}) instead of 100% to reflect uncertainties on the effect of winter curing.MGDD/CDD distribution mapsMGDD were used to compute annual risk maps of developing PD during summer for the period 1981–2019 (see “Methods”). The resulting averaged map identifies all known areas with a recent history of severe PD in the US corresponding to ({{{{{{{mathcal{F}}}}}}}}({{{{{rm{MGDD}}}}}}) , > , 90 %) (i.e., high-risk), such as the Gulf Coast states (Texas, Alabama, Mississippi, Louisiana, Florida), Georgia and Southern California sites (e.g., Temecula Valley) (Fig. 2a), while captures areas with a steep gradation of disease endemicity in the north coast of California (({{{{{{{mathcal{F}}}}}}}}({{{{{rm{MGDD}}}}}} , > , 50 % )). Overall, more than 95% of confirmed PD sites (n = 155) in the US (Supplementary Data 2) fall in grid cells with ({{{{{{{mathcal{F}}}}}}}}({{{{{rm{MGDD}}}}}}) , > , 50 %).Fig. 2: Average thermal-dependent maps for Pierce’s disease (PD) development and recovery in North America and Europe.PD development during the growing season based on average ({{{{{{{mathcal{F}}}}}}}}({{{{{rm{MGDD}}}}}})) estimations between 1981 and 2019 in North America (a) and Europe (b) derived from the results of the inoculation experiments on 36 grapevine varieties. Large differences in the areal extension with favourable MGDDs can be observed between the US and Europe. The winter curing effect is reflected in the distribution of the average ({{{{{{{mathcal{G}}}}}}}}({{{{{rm{CDD}}}}}})) for the 1981–2019 period in the United States (c) and Europe (d). A snapshot of the temperature-driven probability of chronic infection averaged for the 1981–2019 period is obtained from the joint effect of MGDD and CDD in North America (e) and Europe (f). Warmer colours indicate more favourable conditions for chronic PD and the dashed line highlights the threshold of chronic infection probability being 0.5.Full size imageThe average MGDD-projected map for Europe during 1981–2019 spots a high risk for the coast, islands and major river valleys of the Mediterranean Basin, southern Spain, the Atlantic coast from Gibraltar to Oporto, and continental areas of central and southeast Europe (Fig. 2b). Of these, however, only some Mediterranean islands, such as Cyprus and Crete, show ({{{{{{{mathcal{F}}}}}}}}({{{{{rm{MGDD}}}}}}) , > , 99 %) comparable to areas with high disease incidence in the Gulf Coast states of the US and California. Almost all the Atlantic coast from Oporto (Portugal) to Denmark are below suitable MGDD, with an important exception in the Garonne river basin in France (Bordeaux Area) with low to moderate MGDD (Fig. 2b).Figure 2a shows how the area with high-risk MGDD values extends further north of the current known PD distribution in the southeastern US, suggesting that winter temperatures limit the expansion of PD northwards9. A comparison between MGDD and CDD maps (Fig. 2a vs. Fig. 2c, Fig. 2e) further supports the idea that winter curing is restricting PD northward migration from the southeastern US. However, consistent with growing concern among Midwest states winegrowers on PD northward migration led by climate change63, we found a mean increase of 0.12% y−1 in the areal extent with CDD  0.075) in 22.3% of the vineyards in Europe. However, no vineyard is in epidemic-risk zones with a high-risk index and only 2.9% of the vineyard surface is at moderate risk (Supplementary Table S8). The areas with the highest risk index (r(t) between 0.70 and 0.88) are mainly located in the Mediterranean islands of Crete, Cyprus and the Balearic Islands or at pronounced peninsulas like Apulia (Italy) and Peloponnese (Greece) in the continent (Fig. 6a and Supplementary Table S8). Most vineyards are in non-risk zones (42.1%), whereas 35.6% are located in transition zones with presently non-risk but where XfPD could become established in the next decades causing some sporadic outbreaks. In Supplementary Data 4 and Supplementary Table S8, we provide full details of the total vineyard areas currently at risk for each country and region.Fig. 6: Intersection between Corine-land-cover vineyard distribution map and PD-risk maps for 2020 and 2050.Data were obtained from Corine-land-cover (2018) and the layer of climatic suitability forP. spumarius in Europe from35. The surface of the vineyard contour has been enlarged to improve the visualisation of the risk zones and disease-incidence growth-rate ranks. a PD risk map for 2019 and its projection for 2050 (b). Blue colours represent non-risk zones and transient risk zones for chronic PD (R0  More

  • in

    Multistressor global change drivers reduce hatch and viability of Lingcod embryos, a benthic egg layer in the California Current System

    IPCC Climate Change The physical science basis. In Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (eds Masson-Delmotte, V. et al.) (Cambridge University Press, 2021).
    Google Scholar 
    Doney, S. C. et al. Climate change impacts on marine ecosystems. Annu. Rev. Mar. Sci. 4, 11–37 (2012).Article 
    ADS 

    Google Scholar 
    Song, H. et al. Thresholds of temperature change for mass extinctions. Nat. Commun. 12, 4694 (2021).Article 
    ADS 
    CAS 

    Google Scholar 
    Free, C. M. et al. Impacts of historical warming on marine fisheries production. Science 363, 979–983 (2019).Article 
    ADS 
    CAS 

    Google Scholar 
    Cheung, W. W. L. et al. Large-scale redistribution of maximum fisheries catch potential in the global ocean under climate change: Climate change impacts on catch potential. Glob. Change Biol. 16, 24–35 (2010).Article 
    ADS 

    Google Scholar 
    Harley, C. D. G. et al. The impacts of climate change in coastal marine systems: Climate change in coastal marine systems. Ecol. Lett. 9, 228–241 (2006).Article 
    ADS 

    Google Scholar 
    Dahlke, F. T., Wohlrab, S., Butzin, M. & Pörtner, H.-O. Thermal bottlenecks in the life cycle define climate vulnerability of fish. Science 369, 65–70 (2020).Article 
    ADS 
    CAS 

    Google Scholar 
    Hodgson, E. E., Essington, T. E. & Kaplan, I. C. Extending vulnerability assessment to include life stages considerations. PLoS ONE 11, e0158917 (2016).Article 

    Google Scholar 
    Peck, M. A., Reglero, P., Takahashi, M. & Catalán, I. A. Life cycle ecophysiology of small pelagic fish and climate-driven changes in populations. Prog. Oceanogr. 116, 220–245 (2013).Article 
    ADS 

    Google Scholar 
    Tsoukali, S., Visser, A. W. & MacKenzie, B. R. Functional responses of North Atlantic fish eggs to increasing temperature. Mar. Ecol. Prog. Ser. 555, 151–165 (2016).Article 
    ADS 

    Google Scholar 
    Pörtner, H. O. & Peck, M. A. Climate change effects on fishes and fisheries: Towards a cause-and-effect understanding. J. Fish Biol. 77, 1745–1779 (2010).Article 

    Google Scholar 
    Pankhurst, N. W. & Munday, P. L. Effects of climate change on fish reproduction and early life history stages. Mar. Freshw. Res. 62, 1015–1026 (2011).Article 
    CAS 

    Google Scholar 
    Brauner, C. J. Acid-base balance. In Fish Larval physiology (eds Finn, R. N. & Kapoor, B. G.) 185–198 (Science Publishers, 2008).
    Google Scholar 
    Dahlke, F. T. et al. Effects of ocean acidification increase embryonic sensitivity to thermal extremes in Atlantic cod, Gadus morhua. Glob. Chang. Biol. 23, 1499–1510 (2017).Article 
    ADS 

    Google Scholar 
    Shelbourne, J. E. Significance of the subdermal space in pelagic fish embryos and larvae. Nature 176, 743–744 (1955).Article 
    ADS 

    Google Scholar 
    Sundby, S. & Kristiansen, T. The principles of buoyancy in marine fish eggs and their vertical distributions across the world oceans. PLoS ONE 10, e0138821 (2015).Article 

    Google Scholar 
    Shei, M., Mies, M. & Olivotto, I. Other demersal spawners and mouthbrooders. Marine ornamental species aquaculture, 223–250 (2017).Beaudreau, A. H. The predatory role of lingcod (Ophiodon elongatus) in the San Juan Archipelago, Washington. (University of Washington, 2009).Love, M. Certainly More Than You Want to Know About the Fishes of the Pacific Coast: A Postmodern Experience. (Really Big Press, 2011).Pauly, D. & Christensen, V. Primary production required to sustain global fisheries. Nature 374, 255–257 (1995).Article 
    ADS 
    CAS 

    Google Scholar 
    Reum, J. C. et al. Interpretation and design of ocean acidification experiments in upwelling systems in the context of carbonate chemistry co-variation with temperature and oxygen. ICES J. Mar. Sci. 73, 582–595 (2016).Article 

    Google Scholar 
    Cheresh, J. & Fiechter, J. Physical and biogeochemical drivers of alongshore pH and oxygen variability in the California Current System. Geophys. Res. Lett. 47, e2020089553 (2020).Article 
    ADS 

    Google Scholar 
    Gruber, N. et al. Rapid progression of ocean acidification in the California Current System. Science 337, 220–223 (2012).Article 
    ADS 
    CAS 

    Google Scholar 
    Hauri, C. et al. Spatiotemporal variability and long-term trends of ocean acidification in the California Current System. Biogeosciences 10, 193–216 (2013).Article 
    ADS 

    Google Scholar 
    Pepin, P. Effect of temperature and size on development, mortality, and survival rates of the pelagic early life history stages of marine fish. Can. J. Fish. Aquat. Sci. 48, 503–518 (1991).Article 

    Google Scholar 
    Lauel, B. J. & Blood, D. M. The Effects of Temperature on Hatching and Survival of Northern Rock Sole Larvae (Lepidopsetta polyxystra) (Springer, 2011).
    Google Scholar 
    Guevara-Fletcher, C., Alvarez, P., Sanchez, J. & Iglesias, J. Effect of temperature on the development and mortality of European hake (Merluccius merluccius L.) eggs from southern stock under laboratory conditions. J. Exp. Mar. Biol. Ecol. 476, 50–57 (2016).Article 

    Google Scholar 
    Collins, L. A. & Nelson, S. G. Effects of temperature on oxygen consumption, growth, and development of embryos and yolk-sac larvae of Siganus randalli (Pisces: Siganidae). Mar. Biol. 117, 195–204 (1993).Article 

    Google Scholar 
    Cook, M. A., Guthrie, K. M., Rust, M. B. & Plesha, P. D. Effects of salinity and temperature during incubation on hatching and development of lingcod Ophiodon elongatus Girard, embryos. Aquac. Res. 36, 1298–1303 (2005).Article 

    Google Scholar 
    Pörtner, H. Integrating climate-related stressor effects on marine organisms: Unifying principles linking molecule to ecosystem-level changes. Mar. Ecol. Prog. Ser. 470, 273–290 (2012).Article 
    ADS 

    Google Scholar 
    Laurel, B. J., Copeman, L. A., Spencer, M. & Iseri, P. Comparative effects of temperature on rates of development and survival of eggs and yolk-sac larvae of Arctic cod (Boreogadus saida) and walleye pollock (Gadus chalcogrammus). ICES J. Mar. Sci. 75, 2403–2412 (2018).Article 

    Google Scholar 
    Jordaan, A., Hayhurst, S. E. & Kling, L. J. The influence of temperature on the stage at hatch of laboratory reared Gadus morhua and implications for comparisons of length and morphology. J. Fish Biol. 68, 7–24 (2006).Article 

    Google Scholar 
    Peña, R., Dumas, S., Zavala-Leal, I. & Contreras-Olguín, M. Effect of incubation temperature on the embryonic development and yolk-sac larvae of the Pacific red snapper Lutjanus peru (Nichols & Murphy, 1922). Aquac Res 45, 519–527 (2014).Article 

    Google Scholar 
    Breitburg, D. Effects of hypoxia, and the balance between hypoxia and enrichment, on coastal fishes and fisheries. Estuaries 25, 767–781 (2002).Article 

    Google Scholar 
    Hassell, K. L., Coutin, P. C. & Nugegoda, D. Hypoxia impairs embryo development and survival in black bream (Acanthopagrus butcheri). Mar. Pollut. Bull. 57, 302–306 (2008).Article 
    CAS 

    Google Scholar 
    Giorgi, A. E. The Environmental Biology of the Embryos, Egg Masses and Nesting Sites of the Lingcod, Ophiodon elongatus. (University of Washington, 1981).Oseid, D. M. & Smith, L. L. Survival and hatching of walleye eggs at various dissolved oxygen levels. Progress. Fish-Cult. 33, 81–85 (1971).Article 
    CAS 

    Google Scholar 
    Shumway, D. L., Warren, C. E. & Doudoroff, P. Influence of oxygen concentration and water movement on the growth of steelhead trout and coho salmon embryos. Trans. Am. Fish. Soc. 93, 342–356 (1964).Article 

    Google Scholar 
    Baumann, H., Talmage, S. C. & Gobler, C. J. Reduced early life growth and survival in a fish in direct response to increased carbon dioxide. Nat. Clim Change 2, 38–41 (2012).Article 
    ADS 
    CAS 

    Google Scholar 
    Faria, A. M. et al. Effects of high pCO2 on early life development of pelagic spawning marine fish. Mar. Freshw. Res. 68, 2106–2114 (2017).Article 
    CAS 

    Google Scholar 
    Frommel, A. Y. et al. Severe tissue damage in Atlantic cod larvae under increasing ocean acidification. Nat. Clim. Change 2, 42–46 (2012).Article 
    ADS 
    CAS 

    Google Scholar 
    Munday, P. L. et al. Effects of elevated CO2 on early life history development of the yellowtail kingfish, Seriola lalandi, a large pelagic fish. ICES J. Mar. Sci. 73, 641–649 (2016).Article 

    Google Scholar 
    Hurst, T. P., Fernandez, E. R. & Mathis, J. T. Effects of ocean acidification on hatch size and larval growth of walleye pollock (Theragra chalcogramma). ICES J. Mar. Sci. 70, 812–822 (2013).Article 

    Google Scholar 
    Wang, X., Song, L., Chen, Y., Ran, H. & Song, J. Impact of ocean acidification on the early development and escape behavior of marine medaka (Oryzias melastigma). Mar. Environ. Res. 131, 10–18 (2017).Article 

    Google Scholar 
    Franke, A. & Clemmesen, C. Effect of ocean acidification on early life stages of Atlantic herring (Clupea harengus L.). Biogeosciences 8, 3697–3707 (2011).Article 
    ADS 
    CAS 

    Google Scholar 
    Forsgren, E., Dupont, S., Jutfelt, F. & Amundsen, T. Elevated CO 2 affects embryonic development and larval phototaxis in a temperate marine fish. Ecol. Evol. 3, 3637–3646 (2013).Article 

    Google Scholar 
    Bromhead, D. et al. The potential impact of ocean acidification upon eggs and larvae of yellowfin tuna (Thunnus albacares). Deep Sea Res. II 113, 268–279 (2015).Article 
    CAS 

    Google Scholar 
    Garrido, S. et al. Born small, die young: Intrinsic, size-selective mortality in marine larval fish. Sci. Rep. 5, 17065 (2015).Article 
    ADS 
    CAS 

    Google Scholar 
    Sampaio, E. et al. Impacts of hypoxic events surpass those of future ocean warming and acidification. Nat. Ecol. Evol. 5, 311–321 (2021).Article 

    Google Scholar 
    Crain, C. M., Kroeker, K. & Halpern, B. S. Interactive and cumulative effects of multiple human stressors in marine systems. Ecol. Lett. 11, 1304–1315 (2008).Article 

    Google Scholar 
    Pörtner, H. O. Synergistic effects of temperature extremes, hypoxia, and increases in CO2 on marine animals: From Earth history to global change. J. Geophys. Res. 110, 0910 (2005).Article 

    Google Scholar 
    Piggott, J. J., Townsend, C. R. & Matthaei, C. D. Reconceptualizing synergism and antagonism among multiple stressors. Ecol. Evol. 5, 1538–1547 (2015).Article 

    Google Scholar 
    Boyd, P. W. et al. Experimental strategies to assess the biological ramifications of multiple drivers of global ocean change: A review. Glob. Change Biol 24, 2239–2261 (2018).Article 
    ADS 

    Google Scholar 
    Giorgi, A. E. & Congleton, J. L. Effects of current velocity on development and survival of lingcod, Ophiodon elongatus, embryos. Environ. Biol. Fish 10, 15–27 (1984).Article 

    Google Scholar 
    Liu, G., Zhu, S., Liu, D. & Ye, Z. Effect of the C/N ratio on inorganic nitrogen control and the growth and physiological parameters of tilapia s fingerlings, Oreochromis niloticu reared in biofloc systems. Aquac. Res. 49, 2429–2439 (2018).Article 
    CAS 

    Google Scholar 
    Houde, E. D. Fish early life dynamics and recruitment variability. Am. Fish. Soc. Symp. 2, 17–29 (1987).ADS 

    Google Scholar 
    Miller, T. J., Crowder, L. B., Rice, J. A. & Marschall, E. A. Larval size and recruitment mechanisms in fishes: Toward a conceptual framework. Can. J. Fish. Aquat. Sci. 45, 1657–1670 (1988).Article 

    Google Scholar 
    Doi, H., Akamatsu, F. & González, A. L. Starvation effects on nitrogen and carbon stable isotopes of animals: An insight from meta-analysis of fasting experiments. R. Soc. open sci. 4, 170633 (2017).Article 
    ADS 

    Google Scholar 
    Pimentel, M. S. et al. Defective skeletogenesis and oversized otoliths in fish early stages in a changing ocean. J. Exp. Biol. 1, 092635. https://doi.org/10.1242/jeb.092635 (2014).Article 

    Google Scholar 
    Politis, S. N., Dahlke, F. T., Butts, I. A., Peck, M. A. & Trippel, E. A. Temperature, paternity and asynchronous hatching influence early developmental characteristics of larval Atlantic cod, Gadus morhua. J. Exp. Mar. Biol. Ecol. 459, 70–79 (2014).Article 

    Google Scholar 
    Appelbaum, S. et al. Studies on rearing of lingcod Ophiodon elongatus. Aquaculture 135, 219–227 (1995).Article 

    Google Scholar 
    Hempel, G. Early life history of marine fish: The egg stage. Washington Sea Grant. (University of Washington Press, 1979)Gadomski, D. M. & Caddell, S. M. Effects of temperature on the development and survival of eggs of four coastal California fishes. Fish. Bull. 94, 41–48 (1996).
    Google Scholar 
    Parker, L. M. et al. Adult exposure influences offspring response to ocean acidification in oysters. Glob. Change Biol. 18, 82–92 (2012).Article 
    ADS 

    Google Scholar 
    Rombough, P. The effects of temperature on embryonic and larval development. In Global Warming: Implications for Freshwater and Marine Fish (Society for Experimental Biology Seminar Series) (eds Wood, C. & McDonald, D.) 177–224 (Cambridge University Press, 1997).Chapter 

    Google Scholar 
    Bownds, C., Wilson, R. & Marshall, D. J. Why do colder mothers produce larger eggs? An optimality approach. J. Exp. Biol. 213, 3796–3801 (2010).Article 

    Google Scholar 
    Longo, G. C. et al. Strong population differentiation in lingcod ( Ophiodon elongatus ) is driven by a small portion of the genome. Evol. Appl. 13, 2536–2554 (2020).Article 
    CAS 

    Google Scholar 
    Silberberg, K. R., Laidig, T. E., Adams, P. B. & Albin, D. Analysis of maturity in lingcod, Ophiodon elongatus. California Fish Game 87, 139–152 (2001).
    Google Scholar 
    Palumbi, S. R. Why mothers matter. Nature 430, 621–622 (2004).Article 
    ADS 
    CAS 

    Google Scholar 
    Berkeley, S. A., Chapman, C. & Sogard, S. M. Maternal age as a determinant of larval growth and survival in a marine fish, Sebastes melanops. Ecology 85, 1258–1264 (2004).Article 

    Google Scholar 
    Miller, D. J., & Geibel, J. J. Summary of Blue Rockfish and Lingcod Life Histories, a Reef Ecology Study, and Giant Kelp, Macrocystis Pyrifera, Experiments in Monterey Bay, California. (State of California, Resources Agency, Department of Fish and Game, 1973).Low, C. J. & Beamish, R. J. A study of the nesting behavior of lingcod (Ophiodon elongatus) in the strait of Georgia, British Columbia. Can. Fish. Mar. Serv. Tech. Rep. 843, 1–10 (1978).
    Google Scholar 
    King, J. R. & Withler, R. E. Male nest site fidelity and female serial polyandry in lingcod (Ophiodon elongatus, Hexagrammidae): Lingcod nest site fidelity. Mol. Ecol. 14, 653–660 (2005).Article 

    Google Scholar 
    Withler, R. E. et al. Polygamous mating and high levels of genetic variation in lingcod, Ophiodon elongatus of the Strait of Georgia, British Columbia. In Genetics of Subpolar Fish and Invertebrates 345–357 (Springer, 2004).
    Google Scholar 
    Perkins, M. J. et al. Application of nitrogen and carbon stable isotopes (δ15N and δ13C) to quantify food chain length and trophic structure. PLoS ONE 9, e93281 (2014).Article 
    ADS 

    Google Scholar 
    Earth Systems Research Laboratory (ESRL). NOAA’s Ocean Climate Change Web Portal. http://www.esrl.noaa.gov/psd/ipcc/ocn/ (2019).Feely, R., Doney, S. & Cooley, S. Ocean acidification: Present conditions and future changes in a high-CO2 world. Oceanography 22, 36–47 (2009).Article 

    Google Scholar 
    Frieder, C. A., Nam, S. H., Martz, T. R. & Levin, L. A. High temporal and spatial variability of dissolved oxygen and pH in a nearshore California kelp forest. Biogeosciences 9, 3917–3930 (2012).Article 
    ADS 
    CAS 

    Google Scholar 
    Olito, C., White, C. R., Marshall, D. J. & Barneche, D. R. Estimating monotonic rates from biological data using local linear regression. J. Exp. Biol. 1, 148775. https://doi.org/10.1242/jeb.148775 (2017).Article 

    Google Scholar 
    Brooks, M. E. et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 9(2), 378–400 (2017).Article 

    Google Scholar  More

  • in

    Nations forge historic deal to save species: what’s in it and what’s missing

    National negotiators inked a deal to protect nature in the early hours of 19 December in Montreal.Credit: Julian Haber/UN Biodiversity (CC BY 2.0)

    Despite earlier signals of possible failure, countries around the world have cemented a deal to safeguard nature — and for the first time, the agreement sets quantitative biodiversity targets akin to the one that nations set seven years ago to limit global warming to 1.5–2 ºC above pre-industrial levels.In the early hours of 19 December, more than 190 countries eked out the deal, known as the Kunming-Montreal Global Biodiversity Framework, during the COP15 international biodiversity summit in Montreal, Canada. A key target it sets is for nations to protect and restore 30% of the world’s land and seas globally by 2030, while also respecting the rights of Indigenous peoples who depend on and steward much of Earth’s remaining biodiversity. Another target is for nations to reduce the extinction rate by 10-fold for all species by 2050.
    10 startling images of nature in crisis — and the struggle to save it
    Steven Guilbeault, the Canadian environment minister, described COP15 as the most significant biodiversity conference ever held. “We have taken a great step forward in history,” he said at a plenary session where the framework was adopted.At several points during the United Nations summit, which ran from 7–19 December, arguments over details threatened to derail a deal. In the final hours of negotiations, the Democratic Republic of the Congo (DRC) objected to how the framework would be funded. Nonetheless, Huang Runqiu, China’s environment minister and president of COP15, brought the gavel down on the agreement.Negotiators from several African countries, which are home to biodiversity hotspots but say they need funding to preserve those areas, thought that China’s presidency strong-armed the deal. Uganda called it “fraud”. A source who spoke to Nature from the African delegation, and who asked not to be named to maintain diplomacy, said the negotiating process was not equitable towards developing countries and that the deal will not enable significant progress towards stemming biodiversity loss. “It was a coup d’état,” they say. However, a legal expert for the Convention on Biological Diversity — the treaty within which the framework now sits — told COP15 attendees that the adoption of the framework is legitimate.Concerns and disappointmentsScientists and conservation groups have welcomed the deal, emphasizing that there has never been an international agreement to protect nature on this scale. Kina Murphy, an ecologist and chief scientist at the Campaign for Nature, a conservation group, says, “It’s a historic moment for biodiversity.”

    Huang Runqiu, China’s environment minister and president of COP15, brought the gavel down on the biodiversity deal, despite objections from representatives of the Democratic Republic of the Congo.Credit: Julian Haber/UN Biodiversity (CC BY 2.0)

    But some concerns and disappointments remain. For one, the deal lacks a mandatory requirement for companies to track and disclose their impact on biodiversity. “Voluntary action is not enough,” says Eva Zabey, executive director of Business for Nature, a global coalition of 330 businesses seeking such a requirement so that firms can compete on a level playing field. Nevertheless, it sends a powerful signal to industry that it will need to reduce negative impacts over time, says Andrew Deutz, an environmental law and finance specialist at the Nature Conservancy, a conservation group in Arlington, Virginia.In addition, the deal is weak on tackling the drivers of biodiversity loss, because it does not specifically call out the most ecologically damaging industries, such as commercial fishing and agriculture, or set precise targets for them to put biodiversity conservation at the centre of their operations, researchers say.
    Can the world save a million species from extinction?
    “I would have liked more ambition and precision in the targets” to address those drivers, says Sandra Diaz, an ecologist at the National University of Córdoba, in Argentina.The deal is not legally binding, but countries will have to demonstrate progress towards achieving the framework’s goals through national and global reviews. Countries failed to meet the previous Aichi Biodiversity Targets, which were set in 2010 and expired in 2020; scientists have suggested that this failure occurred because of a lack of an accountability mechanism.With the reviews included, the framework “is a very good start, with clear quantitative targets” that will allow us to understand progress and the reasons for success and failure, says Stuart Pimm, an ecologist at Duke University in Durham, North Carolina, and head of Saving Nature, a non-profit conservation organization.A long time comingScientists have estimated that one million species are under threat because of habitat loss, mainly through converting land for agriculture. And they have warned that this biodiversity loss could threaten the health of ecosystems on which humans depend for clean water and disease prevention, and called for a new international conservation effort.
    Crucial biodiversity summit will go ahead in Canada, not China: what scientists think
    The new agreement took 4 years to resolve, in part because of delays caused by the COVID-19 pandemic (the summit was supposed to take place in Kunming, China, in 2020), but also because of arguments over how to finance conservation efforts. Nations finally agreed that by 2030, funding for biodiversity from all public and private sources must rise to at least US$200 billion per year. This includes at least $30 billion per year, contributed from wealthy to low-income nations. These figures fall short of the approximately $700 billion that researchers say is needed to fully safeguard and restore nature, but represents a tripling of existing donations.Low- and middle-income countries (LMICs), including the DRC, had called for a brand-new, independent fund for biodiversity financing. Lee White, environment minister from Gabon, told Nature that biodiversity-rich LMICs have difficulty accessing the Global Environment Facility (GEF), the current fund held by the World Bank in Washington DC, and that it is slow to distribute funds.But France and the European Union strongly objected to a new fund, arguing it would take too long to set up. The framework instead compromises by establishing a trust fund by next year under the GEF. The final agreement also calls on the GEF to reform its process to address the concerns of LMICs.Progress without drastic changeAnother sticking point during negotiations was how to fairly and equitably share the benefits of ‘digital sequence information’ — genetic data collected from plants, animals and other organisms. Communities in biodiversity-rich regions where genetic material is collected have little control over the commercialization of the data, and no way to recoup financial or other benefits from them. But countries came to an agreement to set up a mechanism to share profits, the details of which will be worked out by the next international biodiversity summit, COP16, in 2024.Overall, the deal marks progress toward tackling biodiversity loss, but it is not the drastic change scientists say they were hoping for. “I am not so sure that it has enough teeth to curb the activities that do most of the harm,” Diaz says. More

  • in

    Number of simultaneously acting global change factors affects composition, diversity and productivity of grassland plant communities

    Study species and pre-cultivationTo create the mesocosm communities, we selected nine herbaceous grassland species that are native to and widespread in Central Europe (Supplementary Table 9), where they can also co-occur. The species were Alopecurus pratensis L., Diplotaxis tenuifolia (L.) DC., Lolium perenne L., Poa pratensis L., Prunella vulgaris L., Sinapis arvensis L., Sonchus oleraceus L., Vicia cracca L., Vicia sativa L. To increase generalizability54, the species were selected from three functional groups (grasses, annual forbs, perennial forbs), and they represent five families.Seeds were obtained from different sources (Supplementary Table 9). For the transplanted-seedling community (see section ‘Experimental lay-out), seedlings were pre-cultivated in a greenhouse of the Botanical Garden of the University of Konstanz. As the species require different times for germination, they were sown on different dates (Supplementary Table 10) to ensure that seedlings of all species were at a similar developmental stage at transplantation. Seeds were sown separately per species in plastic trays filled with potting soil (Einheitserde®, Pikiererde CL P). The greenhouse had a regular day-night rhythm of c. 16:8 hours, and its ventilation windows automatically opened at 21 °C during the day and at 18 °C during the night. Two days before transplanting, the seedlings were placed outdoors to acclimatize. For the sown community, we sowed a seed mixture of the nine study species directly into the outdoor mesocosm pots.Experimental setupGlobal change treatmentsWe imposed six global change treatments: climate warming, light pollution, microplastic pollution, soil salinization, eutrophication, and fungicide accumulation, all of which frequently occur in the environment. These GCFs were chosen because they differ in their nature (i.e., physical, chemical), are likely to differ in their mode of action and effect direction21, and can be easily implemented. Each of the six GCFs have been shown to impact plants and their environment when applied on their own10,13,17,19,20,55,56,57,58,59,60. Furthermore, all of the chosen GCFs are likely to continue to increase in magnitude or extent in the near future61,62,63,64,65. For the climate-warming treatment, we used infrared-heater lamps (HS-2420; 240 V, 2000 W; Kalglo Electronics Co., Bethlehem, USA) set to 70% of their maximum capacity to achieve an average temperature increase of 2.0 °C (±SD = 0.2 °C) at plant level. This is within the range of temperature increases predicted by the RCP 4.5 scenario for the year 2100 [+1.1 − 2.6 °C; 63]. For the light-pollution treatment, we used LED spotlights (LED-Strahler Flare 10 W, IP 65, 900 lm, cool white 6500 K; REV Ritter GmbH, Mömbris, Germany), which were switched on daily from 9 pm to 5 am, corresponding to the times of sunset and -rise. The average light intensity was 24.5 lx at ground level, which is within the range of light intensities found below street lights, and matches the light intensities used in other light-pollution experiments14,56. For the microplastic pollution treatment, we used granules (1.0–2.5 mm diameter) of the synthetic rubber ethylene propylene diene monomer (EPDM Granulat, Gummi Appel GmbH + Co. KG, Kahl am Main, Germany) at a concentration of 1% (w/w, granules/dry soil, approximately corresponding to 1.5% v/v). EPDM granules are, for example, used in artificial sport turfs, from where they easily spread into the surroundings, and have been used previously to investigate the effects of microplastics on plants18. The chosen concentration is well within the range of concentrations used in previous studies18,66,67, and is at the low to intermediate range of concentrations found in sites polluted with plastics68. For the soil-salinization treatment, dissolved NaCl was added to the soil. Soil salinity is commonly measured as electrical conductivity, with a conductivity between 4 and 8 dS m−1 considered to be moderately saline69. For the experiment, we used a salinity of 6 dS m−1. To maintain a more or less constant salinity level, electrical conductivity was measured weekly, and, if required, adjusted by adding dissolved NaCl. For the eutrophication treatment, 3 g of a dissolved NPK fertilizer (Universol® blue oxide, ICL SF Germany & Austria, Nordhorn, Germany) was added per pot. For N, this corresponds to an input of 100 kg N ha−1, comparable to the yearly amounts of atmospheric N deposition in large parts of Europe52 and the yearly nitrogen input on agricultural field in the European Union70. To ensure a more or less constant nutrient availability during the experiments, we split total fertilizer input into three applications (directly after, 3 weeks after, and 6 weeks after starting the experiments) of 1 g fertilizer per pot per application. In addition, to avoid severe nutrient limitation in the other pots, all pots (irrespective of the eutrophication treatment) received basic fertilization. This was applied four times to the transplanted-seedling-community pots and five times to the sown community pots, with 0.2 g fertilizer per pot per application. For the pesticide treatment, we used the fungicide Landor® CT (Syngenta Agro GmbH, Maintal, Germany). This fungicide was chosen because it contains three azoles as active agents, which belong to the most widely used class of antifungal agents71. To each pot in this treatment, we added 1.5 μl fungicide dissolved in water (1‰). This corresponds to 60% of the maximum amount that should be used per hectare of cropland. A summary of the levels of the individual GCFs used in our experiment is provided in Supplementary Table 8.Combinations of simultaneously acting GCFsTo examine the potential effects of the numbers of simultaneously acting GCFs, we created five levels of increasing GCF numbers. These levels were: zero (i.e., the control without any GCF application), one (single), two, four and six GCFs. For the one-, two- and four-GCF levels, there were six different combinations, so that each of these levels included either six different GCFs in case of the one-factor, or six different GCF combinations in case of the two- and four-GCF levels. In the six-GCF level, all six factors were combined, so that there was only one combination. To avoid potential biases due to unequal representation of the different GCFs in each GCF-number level, we created the GCF combinations randomly but with the restriction that each GCF was present in an equal number of combinations for each GCF-number level (i.e., each GCF was included once in GCF-number levels 1 and 6, respectively, twice in GCF-number level 2, and four times in GCF-number level 4; Supplementary Table 11).Experimental lay-outThe experiment was conducted outdoors in the climate-warming-simulation facility of the Botanical Garden of the University of Konstanz, Germany (N: 47°69’19.56”, E: 9°17’78.42”). Twenty of the 2 m × 2 m plots of this facility were used for our experiment. As the climate-warming and light-pollution treatments could not be applied to each individual pot separately, we applied those treatments at the plot level. Therefore, we assigned four of the 20 plots to the climate-warming treatment, four plots to the light-pollution treatment and four plots to both climate-warming and light-pollution treatment combination. Each plot had a 145 cm high metal frame. The eight plots assigned to the climate-warming treatment were equipped with a 1.80 m long, horizontally hanging infrared-heating lamp at the top of the metal frame (i.e., at 145 cm above soil level). The heating lamp slowly oscillated along its longitudinal axis to ensure uniform heating of the whole 2 m × 2 m plot. The eight plots assigned to the light-pollution treatment, each had a LED spotlight attached to one of the sides of the metal frame at a height of 120 cm. To reduce illumination of the neighboring plots, light-pollution was only applied to the outer plots of the climate-warming-simulation facility (Supplementary Fig. 5), and LEDs were pointing away from the inner plots and were equipped with lamp shades made of black plastic pots (18 cm × 18 cm × 25.5 cm). Furthermore, to reduce the light intensity to a realistic light-pollution level (24.5 lx) as found below street lights, we covered the spotlight with a layer of white cloth (Supplementary Fig. 6). For further details on the artificial light treatment, see Supplementary Fig. 7.To create mesocosms with the transplanted-seedling and sown communities, we filled 10-L pots (CEP- Container, 10.0 F, Burger GmbH, Renningen-Malmsheim, Germany) with a mixture of 40% potting soil (see above), 40% quartz sand (0.5–0.8 mm), and—to inoculate the substrate with a natural soil community—20% top soil excavated from a seminatural grassland patch in the botanical garden. In total, the experiments with the transplanted-seedling and sown communities, each included 120 pots (i.e., 20 treatment combinations × six replicates × 2 experiments = 240 pots in total; see Supplementary Table 11), which were distributed across the 20 plots. To prevent leakage of fertilizer or salt solutions, each pot was placed onto a plastic dish. To reduce differences due to environmental variation within plots, the positions of pots within each plot were re-randomized every 14 days. Plants were watered regularly to avoid drought stress and to avoid differences in soil moisture due to application of fertilizer- and salt-solutions.For the sown community, we randomly distributed five seeds of each of the nine species on the substrate in each pot on 3 July 2020. For the transplanted-seedling community, two seedlings of each of the nine species were transplanted into each pot (i.e., 18 seedlings per pot) according to a fixed pattern (Supplementary Fig. 8) on 6 July 2020. Since there were a few seedlings missing for S. arvensis (six seedlings) and V. cracca (four seedlings), we re-sowed these species in germination trays on 6 July 2020. On 13 July 2020, dead seedlings, and the missing seedlings for S. arvensis were replaced. Since V. cracca took longer to germinate, the missing seedlings were transplanted on 17 July 2020.MeasurementsTo investigate the effects of single-GCFs and their number on the sown and transplanted-seedling communities, we used plant biomass as an indicator for plant performance72. As it was impossible to disentangle the roots, we only used aboveground biomass. On 14 and 15 September 2020, i.e., 10 weeks after transplanting, we harvested the transplanted-seedling communities. On 28 and 29 September, i.e., twelve weeks after sowing, we harvested the sown communities. For both community types, we harvested the plants separately by species. The harvested plants were stored in paper bags, dried at 70 °C for at least 72 hours and weighed.Statistical analysisAll analyses were done in R 3.6.273. As the transplanted-seedling and sown communities were harvested at different times, we treated them as separate experiments, and therefore analyzed them separately (but see the subsection “Community type specific responses” below).Community aboveground biomassTo analyze the effects GCF number on plant-community productivity, we fitted linear mixed-effects models separately for the transplanted-seedling and sown communities, using the lmer function in the “lme4” package74. Total aboveground biomass per pot was the response variable. To improve normality of the residuals, biomass of the transplanted-seedling and sown communities was square-root- and natural-log-transformed, respectively. We included GCF number as a continuous fixed variable. To account for non-independence of pots in the same GCF combination and of pots in the same plot, GCF combination and plot were included as random effects. The effects of the individual GCFs on biomass production were also assessed by fitting linear mixed-effects models, using only the data of the control and single-GCF treatments, and including GCF identity as fixed effect.Community compositionTo assess potential effects of single-GCFs and GCF number on the final composition of the transplanted-seedling and sown communities, we first assessed variation in species composition, based on biomass proportions, among pots using principal component analysis (PCA) [rda function of the “vegan” package75,]. For each PCA (Supplementary Fig. 1), we extracted the PC1 and PC2 values, which together explained more than 65% of the variation in community composition and included them as response variables in separate linear mixed models, as described above for community biomass.To evaluate whether GCF number affects the diversity and evenness of plant communities, we calculated the Shannon index (H)76, using the diversity function in the “vegan” package, and evenness index (J)77 based on species biomass proportions. Subsequently, the single-GCF and GCF-number effects on diversity and evenness of the sown and transplanted-seedling communities were analyzed using linear mixed-effects models, or—if adding random effects did not improve the model—more parsimonious linear models78,79. For all models, we used type II analysis of variance (ANOVA) tests (Anova function in the “car” package) to assess the significance of fixed effects.Hierarchical diversity-interaction modelingWhen there is a significant GCF-number effect, this could reflect that with increasing numbers of co-acting GCFs, there is a higher chance that it will include a GCF with a strong and dominant effect (i.e., sampling or selection effects). However, it could also be that the GCF-number effect is driven by interactions among the GCFs, and the effects of these interactions could be GCF-specific or general. As our experiment does not include all possible combinations of GCFs, it does not allow to test the contributions of each possible multi-way GCF interaction. Therefore, to gain insights into whether the GCF identities and specific or general GCF interactions underlie the significant GCF-number effects, we applied the hierarchical diversity-interaction modeling framework of Kirwan et al.80. This framework was originally developed for estimating contributions of species identities and their interactions to ecosystem functions, but we here applied it to GCF identities and interactions. For each of the response variables showing a significant GCF-number effect, we ran five hierarchical models specifying different assumptions about the potential contributions of individual GCFs and their interactions to the GCF-number effect, and compared them using likelihood ratio tests (Fig. 4). For these analyses, the data of the control treatment (i.e., GCF number zero) was excluded. Each of the five models specified different assumptions about the potential contributions of individual GCFs and their interactions to the GCF-number effect. The first model is the null model, which assumed that there were no GCF-specific contributions (i.e., all GCFs contributed equally) and that there were no contributions of GCF interactions. Therefore, the null model only included the centered sum of the GCFs of each treatment (M) as fixed effect. M accounts for differences in ‘initial abundances’ of GCFs—meaning that the other model terms are interpreted based on the average initial abundance—and was also included in the four other models80. This way, we could include the GCFs’ relative proportions in each GCF combination, instead of just considering GCF presence, while taking into account that, with increasing GCF number, the relative proportion of each individual GCF is automatically reduced. In the second model, the GCF identities (i.e., their proportions in the respective GCF combination) were added, assuming that individual GCFs contribute differently to the effect of GCF number. In the third model, separate-pairwise interactions between the GCFs were added, considering that, in addition to contributions of individual GCFs, specific pairwise interactions contributed to the GCF-number effect. In the fourth model, the average GCF-interaction model (which is also called the evenness model in Kirwan et al. 2009), the separate-pairwise GCF interactions were replaced by an average interaction effect. Thus, the average GCF-interaction model assumed equivalent contributions of all pairwise GCF interactions. In the fifth model, the additive GCF-specific interaction contributions model, the average interaction effect of the fourth model was replaced by average GCF-specific interaction effects. This model assumed that each GCF’s contribution to a pairwise interaction remains constant. For the calculation of the average GCF-specific and average interaction effect, we used the equations provided by Kirwan et al.80. For each of the response variables, we generally included the same random terms as in the main analyses of the GCF-number effect. However, as this resulted in singularity warnings for some of the hierarchical diversity-interaction models, e.g., those for species diversity and evenness measures, we used for these cases linear models instead of linear mixed models.Fig. 4: Hierarchical diversity-interaction-modeling framework to assess contributions of GCF identities and GCF interactions to GCF-number effects.The framework was adapted from Kirwan et al.80. The null model assumes equivalent contributions of all GCFs and no interactions between them. The subsequent models assume more complex effects of how the individual GCFs and their interactions determine the GCF-number effects. The questions that can be answered by comparing specific models are depicted next to the arrows connecting the two models.Full size imageCommunity type-specific responsesAs the transplanted-seedling and sown communities were harvested at different times, we treated them as separate experiments, and therefore analyzed them separately. However, to test explicitly whether both community types differed in their responses to single-GCFs and GCF number, we also analyzed them jointly. To this end, we fitted linear mixed-effects models for each response variable including GCF number (or single-factor treatments), community type and their interaction as fixed effects (Supplementary Table 5).Final number of plants per speciesTo test for effects of individual GCFs and GCF number on species presence, i.e., the number of individuals per species present at harvest, we fitted generalized linear mixed-effects models for the transplanted-seedling and sown communities separately. We included the survival rate (number of individuals present at harvest divided by the number of planted/sown individuals) as response variables. For the models testing the effects of GCF number, we included GCF combination, species, pot, and plot as random effects. For the models testing the effects of single-GCFs, the same random effects were included, except for GCF combination. Specific random effects were removed from the model if their incorporation resulted in singular fit warnings due to low variation. We assessed the effects of individual GCFs or GCF number using type III ANOVA tests (Anova function in the “car” package, Supplementary Table 7).Eutrophication effectsIn addition to the general assessment of individual GCF effects in the hierarchical diversity-interaction models, we specifically assessed the effects of eutrophication. This was done because eutrophication had the strongest effect on productivity as individual GCF, and this might also have dominated the GCF-number effect, indicating a sampling effect. To this end, we added a binary-coded variable to include information on whether eutrophication was included in the different GCF combinations. Subsequently, we fitted linear mixed-effects models for all response traits that were affected by GCF number. In these models, we included GCF number, community type, eutrophication, and the respective two-way interactions as fixed effects, and plot and GCF combination as random effects. Effects of fixed factors were assessed using type III ANOVA tests (Anova function in the “car” package; Supplementary Table 6).Reporting summaryFurther information on research design is available in the Nature Portfolio Reporting Summary linked to this article. More

  • in

    Wood structure explained by complex spatial source-sink interactions

    Overall frameworkCells in our model are arranged along independent radial files, with each cell in one of either the proliferation, enlargement-only, wall thickening, or mature zones, depending on the distance of the cell’s centre from the inside edge of the phloem and the time of year. Only cells that contribute to the formation of xylem tracheids are treated explicitly. A daily timestep is used, on which cells in the proliferation and enlargement-only zones can enlarge in the radial direction if these zones are non-dormant, and on which secondary-wall thickening can occur in the wall thickening zone. Cells in the proliferation zone divide periclinally if they reach a threshold radial length. Cell-size control at division is intermediate between a critical size and a critical increment22. Mother cells divide asymmetrically, with the subsequent relative growth rates of the daughters inversely proportional to their relative sizes. Size at division and asymmetry of division are computed with added statistical noise22, and therefore the model is run for an ensemble of independent radial files with perturbed initial conditions.Equations and parametersCell enlargement and divisionCells in the proliferation and enlargement-only zones, when not dormant, enlarge in the radial direction at a rate dependent on temperature and relative sibling birth size. A Boltzmann-Arrhenius approach is used for the temperature dependence30:$$mu={mu }_{0}{e}^{frac{{E}_{a}}{k}left(frac{1}{{T}_{0}}-frac{1}{T}right)}$$
    (1)
    where μ is the relative rate of radial cell growth at temperature T (μm μm−1 day−1), μ0 is μ at temperature T0 (=283.15 K), Ea is the effective activation energy for cell enlargement, k is the Boltzmann constant (i.e. 8.617 x 10−5 eV K−1), and T is temperature (K). μ0 was calibrated to an observed mean radial file length at the end of the elongation period dataset23 (Table 1; see “Observations”), and Ea was calibrated to an observed temperature dependence of annual ring width dataset31 (Table 1; Supplementary Fig. 4; see “Observations”).Table 1 Model parameters calibrated to observationsFull size tableRadial length of an individual cell then increases according to:$${{Delta }}{L}_{r}={L}_{r}({e}^{epsilon mu }-1)$$
    (2)
    where ΔLr is the radial increment of the cell (μm day−1), Lr is the radial length of the cell (μm), and ϵ is the cell’s growth dependence on relative birth size, given by22:$$epsilon=1-{g}_{asym}{alpha }_{b}$$
    (3)
    where gasym is the strength of the dependence of relative growth rate on asymmetric division (Table 2; unitless), and αb is the degree of asymmetry relative to the cell’s sister22 (scalar):$${alpha }_{b}=frac{{L}_{r}{,}_{b}-{L}_{r}{,}_{b}^{sis}}{{L}_{r}{,}_{b}+{L}_{r}{,}_{b}^{sis}}$$
    (4)
    where Lr,b is the radial length of the cell at birth (μm) and ({L}_{r}{,}_{b}^{sis}) is the radial length of its sister at birth (μm), which are calculated stochastically22:$${L}_{r}{,}_{b}={L}_{r}{,}_{d}(0.5-{Z}_{a})$$
    (5)
    $${L}_{r}{,}_{b}^{sis}={L}_{r}{,}_{d}(0.5+{Z}_{a})$$
    (6)
    where Lr,d is the length of the mother cell when it divides (μm) and Za is Gaussian noise with mean zero and standard deviation σa (Table 2; −0.49 ≤Za≤ 0.49 for numerical stability).Table 2 Parameters used in the model that are taken directly from literatureFull size tableLength at division is derived as22:$${L}_{r}{,}_{d}=f{L}_{r}{,}_{b}+{chi }_{b}(2-f+Z)$$
    (7)
    where f is the mode of cell-size regulation (Table 2; unitless), χb is the mean cell birth size (Table 3; μm), and Z is Gaussian noise with mean zero and standard deviation σ (Table 2).Table 3 Parameters used in the model that are calculated from observationsFull size tableThe first cell in each radial file is an initial, which produces phloem mother cells outwards and xylem mother cells inwards. It grows and divides as other cells in the proliferation zone, but on division one of the daughters is stochastically assigned to phloem or xylem, the other remaining as the initial. The probability of the daughter being on the phloem side is fphloem (Table 3).Cell-wall growthBoth primary and secondary cell-wall growth are influenced by temperature, carbohydrate concentration, and lumen volume. A Michaelis-Menten equation is used to relate the rate of wall growth to the concentration of carbohydrates in the cytoplasm:$${{Delta }}M=frac{{{Delta }}{M}_{max}theta }{theta+{K}_{m}}$$
    (8)
    where ΔM is the rate of cell-wall growth (mg cell−1 day−1), ΔMmax is the carbohydrate-saturated rate of wall growth (mg cell−1 day−1), θ is the concentration of carbohydrates in the cell’s cytoplasm (mg ml−1), and Km is the effective Michaelis constant (mg ml−1; Table 1).The maximum rate of cell-wall growth, ΔMmax, is assumed to depend linearly on lumen volume (a proxy for the amount of machinery for wall growth), and on temperature as in Eq. (1):$${{Delta }}{M}_{max}=omega {V}_{l}{e}^{frac{{E}_{aw}}{k}left(frac{1}{{T}_{0}}-frac{1}{T}right)}$$
    (9)
    where ω is the normalised rate of cell-wall mass growth (i.e. the rate at T0; Table 1; mg ml−1 day−1), Vl is the cell lumen volume (ml cell−1), and Eaw is the effective activation energy for wall building (eV; Table 1). ω and Km were calibrated to an observed distribution of carbohydrates23 (see next section). Eaw was calibrated to an observed temperature dependence of maximum density31 (Table 1; see “Observations”).Lumen volume is given by:$${V}_{l}={V}_{c}-{V}_{w}$$
    (10)
    where Vc is total cell volume (ml cell−1) and Vw is total wall volume (ml cell−1). Cells are assumed cuboid and therefore Vc is given by:$${V}_{c}={L}_{a}{L}_{t}{L}_{r}/1{0}^{12}$$
    (11)
    where La is axial length (μm; Table 2) and Lt is tangential length (μm; Table 3). Vw is given by:$${V}_{w}=M/rho$$
    (12)
    where M is wall mass (mg cell−1) and ρ is wall-mass density (Table 2; mg[DM] ml−1).Cells in the proliferation and enlargement-only zones only have primary cell walls. ΔMmax (Eq. (9)) is therefore given the following limit:$${{Delta }}{M}_{max}=min ({{Delta }}{M}_{max},rho {V}_{{w}_{p}}-M)$$
    (13)
    where ({V}_{{w}_{p}}) is the required primary wall volume:$${V}_{{w}_{p}}={V}_{c}-({L}_{a}-2{W}_{p})({L}_{t}-2{W}_{p})({L}_{r}-2{W}_{p})/1{0}^{12}$$
    (14)
    where Wp is primary cell-wall thickness (Table 3; μm).Carbohydrate distributionThe distribution of carbohydrates across each radial file is calculated independently from the balance of diffusion from the phloem and the uptake into primary and secondary cell walls. The carbohydrate concentration in the phloem is prescribed at the mean value observed across the three observational dates in23, as described below in “Observations”, and the resulting concentration in the cytoplasm of the furthest living cell from the phloem is solved numerically. The inside wall of this cell is assumed to be impermeable to carbohydrates and therefore provides the inner boundary to the solution. It is assumed that the rate of diffusion across each file is rapid relative to the rate of cell-wall building, and therefore concentrations are assumed to be in equilibrium on each day. Carbohydrate diffusion between living cells is assumed to be proportional to the concentration gradient:$${q}_{i}=({theta }_{i-1}-{theta }_{i})/eta$$
    (15)
    where qi is the rate of carbohydrate diffusion from cell i − 1 to cell i (mg day−1) and η is the resistance to flow between cells (calibrated to the observed distribution of carbohydrates23, see next section; Table 1; day ml−1).As it is assumed that carbohydrates cannot diffuse between radial files, at equilibrium the flux into a given cell must equal the rate of wall growth in that cell plus the wall growth in all cells further along the radial file away from the phloem. From this it can be shown that the equilibrium carbohydrate concentration in the furthest living cell from the phloem in each radial file is given by:$${theta }_{n}={theta }_{p}-eta mathop{sum }limits_{i=1}^{n}mathop{sum }limits_{j=i}^{n}{{Delta }}{M}_{j}$$
    (16)
    where θp is the concentration of carbohydrates in the phloem (Table 3; mg ml−1) and n is the number of living cells in the file (phloem mother cells are ignored for simplicity). The rate of wall growth in each cell depends on the concentration of carbohydrates (Eq. (8)), and therefore θn must be found that results in an equilibrium flow across the radial file. This is done using Brent’s method41 as implemented in the “ZBRENT” function42.Zone widthsThe widths of the proliferation, enlargement-only, and secondary wall thickening zones vary through the year, and are fitted to observations on three dates23 (see Supplementary Fig. 2 and “Observations”). Linear responses to daylength were found, which are therefore used to determine widths for the observational period and later days:$${z}_{k}={a}_{k}+{b}_{k}{{{{{{{rm{dl}}}}}}}};{{mathrm{DOY}}}ge 185$$
    (17)
    where zk is the distance of the inner edge of the zone from the inner edge of the phloem (μm), k is proliferation (p), secondary wall thickening (t), or enlargement-only (e), ak and bk are constants (Table 3), dl is daylength (s), and DOY is day-of-year. The proliferation zone width on earlier days when non-dormant was fixed at its DOY 185 width (assuming this to be its maximum, and that it would reach its maximum very soon after cambial dormancy is broken in the spring). During dormancy, the proliferation zone width is fixed at its value on DOY 231 (the first day of dormancy23). The enlargement-only zone width prior to DOY 185, the first observational day, is assumed to be a linear extension of the rate of change after DOY 185. The wall-thickening zone width plays little role prior to DOY 185 at the focal site, and so was set to its Eq. (17) value each earlier day. On all days the condition zt≥ze≥zp is imposed, and zone widths cannot exceed their values at 24 h daylength (necessary for sites north of the Arctic circle). Supplementary Figure 2 shows the resulting progression of zone widths through the year, together with the observed values.DormancyProliferation was observed to be finished by DOY 23123, and so the proliferation and enlargement-only zones are assumed to enter dormancy then. Release from dormancy in the spring is calculated using an empirical thermal time/chilling model33. It was necessary to adjust the asymptote and temperature threshold of the published model because the heat sum on the day of release calculated from observations in Sweden (see “Observations”) was much lower than reported for Sitka spruce buds in Britain in the original work:$${{{{{{{{rm{dd}}}}}}}}}_{req}=15+4401.8{e}^{-0.042{{{{{{{rm{cd}}}}}}}}}$$
    (18)
    where ddreq is the required sum of degree-days (°C) from DOY 32 for dormancy release and cd is the chill-day sum from DOY 306. The degree-day sum is the sum of daily mean temperatures above 0 °C, and the chill-day sum is the number of days with mean temperatures below 0 °C. Dormancy can only be released during the first half of the year.Simulation protocolsEach simulation consisted of an ensemble of 100 independent radial files. Each radial file was initialised by producing a file of 100 cells with radial lengths χb(1+Za), allowing these to divide once, ignoring the second daughter from each division, and then limiting the remaining daughters to only those falling inside the proliferation zone on DOY 1. Values for ϵ (the relative growth of daughter cells) and Lr,d (the radial length at division) were derived for each cell. The main simulations were conducted at the observation site in boreal Sweden (64.35°N, 19.77°E) over 1951–1995 to capture the growth period of the observed trees23. Results are mostly presented for 1995 when the observations were made. Simulations for calibration of the effective activation energies (i.e. Ea and Eaw) were performed at 68.26°N, 19.63°E in Arctic Sweden over 1901–200431. Daily mean temperatures for both sites were derived from the appropriate gridbox in a 6 h 1/2 degree global-gridded dataset43.ObservationsObservations of cellular characteristics and carbohydrate concentrations23 were used to derive a number of model parameters, and to test model output (model calibration and testing were performed using different outputs). According to the published study we used, samples were cut from three 44 yr old Scots pine trees growing in Sweden (64°21’ N; 19°46’ E) at different times during the growing season. 30 μm thick longitudinal tangential sections of the cambial region were made, and the radial distributions of soluble carbohydrates measured using microanalytical techniques23. Cell sizes, wall thicknesses, and positions in their Fig. 123, an image of transverse sections on three sampling dates, were digitised using “WebPlotDigitizer”44. These, together with the numbers of cells in each zone and their sizes given in the text of that paper, were used to estimate zone widths, which were then regressed against daylength to give the parameters for Eq. (17) (Table 3), mean cell size in the proliferation zone on the first sampling date (used to derive χb; Table 3), mean cell tangential length (Table 3), and final ring width (used to calibrate μ0; Table 1). The thickness of the primary cell wall (Table 3) was derived by plotting cell-wall thickness against time and taking the low asymptote.The distributions of carbohydrates along the radial files on the last sampling date for “Tree 1” and “Tree 3” (results for “Tree 2” were not shown for this date) shown in Fig. 2 of the observational paper23 were calculated. The masses for each of sucrose, glucose, and fructose in each 30 μm section were digitised using the same method as for cell properties and then summed and converted to concentrations, with the results shown in Supplementary Figure 5. Mean observed carbohydrate concentrations and cell masses at four points were used to calibrate values for the η, ω, and Km parameters in Table 1. Calibration was performed by minimising the summed relative error across the observations.The calibration target for the effective activation energy for wall deposition (i.e. Eaw) was the observed relationship between maximum density and mean June-July-August temperature over 1901-2004 in northern Sweden31 (Supplementary Fig. 3), and for the effective activation energy for cell enlargement the relationship between ring width and temperature (i.e. Ea) target was the same study (Supplementary Fig. 4). These observations were made on living and subfossil Scots pine sample material from the Lake Tornesträsk area (68.21–68.31°N; 19.45–19.80°E; 350–450 m a.s.l.) using X-ray densitometry for maximum density, and standardised to remove non-climatic information31.Reporting summaryFurther information on research design is available in the Nature Portfolio Reporting Summary linked to this article. More

  • in

    Out-of-date datasets hamper conservation of species close to extinction

    Scheffers, B. R., Yong, D. L., Harris, J. B. C., Giam, X. & Sodhi, N. S. The world’s rediscovered species: back from the brink? PloS ONE 6, e22531 (2011).Article 
    CAS 

    Google Scholar 
    Abeli, T., Albani Rocchetti, G., Barina, Z., Bazos, I. & Draper, D. et al. Seventeen ‘extinct’ plant species back to conservation attention in Europe. Nat. Plants 7, 282–286 (2021).Article 

    Google Scholar 
    Guidelines for Using the IUCN Red List Categories and Criteria Version 14 (IUCN Standards and Petitions Committee, 2019); http://www.iucnredlist.org/documents/RedListGuidelines.pdfDalrymple, S. E. & Abeli, T. Ex situ seed banks and the IUCN Red List. Nat. Plants 5, 122–123 (2019).Article 

    Google Scholar 
    Albani Rocchetti, G. et al. Selecting the best candidates for resurrecting extinct-in-the-wild plants from herbaria. Nat. Plants. https://doi.org/10.1038/s41477-022-01296-7 (2022).The IUCN Red List of Threatened Species Version 2022-1 (IUCN, accessed 264 October 2022); https://www.iucnredlist.orgHumphreys, A. M., Govaerts, R., Ficinski, S. Z., Lughadha, E. N. & Vorontsova, M. S. Global dataset shows geography and life form predict modern plant extinction and rediscovery. Nat. Ecol. Evol. 3, 1043–1047 (2019).Article 

    Google Scholar 
    Knapp, W. M., Frances, A., Noss, R., Naczi, R. F. & Weakley, A. et al. Vascular plant extinction in the continental United States and Canada. Conserv. Biol. 35, 360–368 (2021).Article 

    Google Scholar 
    Sasidharan, N. Cynometra beddomei. The IUCN Red List of Threatened Species 2020 (IUCN, accessed 27 October 2021); https://www.iucnredlist.org/species/31184/115932185Cronk, Q. C. B. A new species and hybrid in the St Helena endemic genus Trochetiopsis. Edinb. J. Bot. 52, 205–213 (1995).Article 

    Google Scholar 
    Loizeau, P. A. & Jackson, P. W. World Flora Online mid-term update. Ann. Missouri Bot. Gard. 102, 341–346 (2017).Article 

    Google Scholar 
    Edwards, C., Bassüner, B., Birkinshaw, C., Camara, C. & Lehavana, A. et al. A botanical mystery solved by phylogenetic analysis of botanical garden collections: the rediscovery of the presumed-extinct Dracaena umbraculifera. Oryx 52, 427–436 (2018).Article 

    Google Scholar 
    MosaChristas, K., Karthick, R., Kowsalya, E. & Jaquline, C. R. I. Musa kattuvazhana (Musaceae): rediscovery and additional notes on a critically endangered species from Western Ghats of Tamil Nadu, India. Feddes Repert. 132, 263–268 (2021).Article 

    Google Scholar 
    Van Hoi, Q. U. A. C. H., Doudkin, R. V., Cuong, T. Q., Le Van, S. O. N. & Van Dung, L. U. O. N. G. et al. Rediscovery of Camellia langbianensis (Theaceae) in Vietnam. Phytotaxa 480, 85–90 (2021).Article 

    Google Scholar 
    Wilkinson, M. D., Dumontier, M., Aalbersberg, I. J., Appleton, G. & Axton, M. et al. The FAIR guiding principles for scientific data management and stewardship. Sci. Data 3, 160018 (2016).Costello, M. J. & Wieczorek, J. Best practice for biodiversity data management and publication. Biol. Conserv. 173, 68–73 (2014).Article 

    Google Scholar 
    Wieczorek, J., Bloom, D., Guralnick, R., Blum, S., Döring, M., Giovanni, R., Robertson, T. & Vieglais, D. Darwin Core: an evolving community-developed biodiversity data standard. PloS ONE 7, e29715 (2012).Article 
    CAS 

    Google Scholar 
    de Lange, P.J. Lepidium obtusatum Fact Sheet (content continuously updated) (New Zealand Plant Conservation Network, accessed 16 December 2021); https://www.nzpcn.org.nz/flora/species/lepidium-obtusatum/Knapp, W.M., Poindexter, D.B. & Weakley, A.S. The true identity of Marshallia grandiflora an extinct species and the description of Marshallia pulchra (Asteraceae Helenieae Marshalliinae). Phytotaxa 447, 1–15 (2020).Article 

    Google Scholar  More

  • in

    Environmentally driven phenotypic convergence and niche conservatism accompany speciation in hoary bats

    Orr, M. R. & Smith, T. B. Ecology and speciation. Trends Ecol. Evol. 13, 502–506 (1998).Article 
    CAS 

    Google Scholar 
    Coyne, J. A. & Orr, H. A. Speciation (Sinauer Associates, 2004).
    Google Scholar 
    Gillespie, R. G. Adaptive radiation: Convergence and non-equilibrium. Curr. Biol. 23, R71–R74 (2013).Article 
    CAS 

    Google Scholar 
    Price, T. Speciation in Birds (Roberts and Company Publishers, 2008).
    Google Scholar 
    Schluter, D. Evidence for ecological speciation and its alternative. Science 323, 737–741 (2009).Article 
    ADS 
    CAS 

    Google Scholar 
    Stroud, J. T. & Losos, J. B. Ecological opportunity and adaptive radiation. Annu. Rev. Ecol. Evol. Syst. 47, 507–532 (2016).Article 

    Google Scholar 
    Jønsson, K. A. et al. Ecological and evolutionary determinants for the adaptive radiation of the Madagascan vangas. Proc. Natl. Acad. Sci. 109, 6620–6625 (2012).Article 
    ADS 

    Google Scholar 
    Wiens, J. J. Speciation and ecology revisited: Phylogenetic niche conservatism and the origin of species. Evolution 58, 193–197 (2004).
    Google Scholar 
    Barve, N. et al. The crucial role of the accessible area in ecological niche modeling and species distribution modeling. Ecol. Model. 222, 1810–1819 (2011).Article 

    Google Scholar 
    Wiens, J. J. & Graham, C. H. Niche Conservatism: Integrating evolution, ecology, and conservation biology. Annu. Rev. Ecol. Evol. Syst. 36, 519–539 (2005).Article 

    Google Scholar 
    Petitpierre, B. et al. Climatic niche shifts are rare among terrestrial plant invaders. Science 335, 1344–1348 (2012).Article 
    ADS 
    CAS 

    Google Scholar 
    Winger, B. M., Barker, F. K. & Ree, R. H. Temperate origins of long-distance seasonal migration in New World songbirds. Proc. Natl. Acad. Sci. 111, 12115–12120 (2014).Article 
    ADS 
    CAS 

    Google Scholar 
    Alerstam, T., Hedenström, A. & Åkesson, S. Long-distance migration: Evolution and determinants. Oikos 103, 247–260 (2003).Article 

    Google Scholar 
    Gómez, C., Tenorio, E. A., Montoya, P. & Cadena, C. D. Niche-tracking migrants and niche-switching residents: Evolution of climatic niches in New World warblers (Parulidae). Proc. R. Soc. B Biol. Sci. 283, 20152458 (2016).Article 

    Google Scholar 
    Menchaca, A., Arteaga, M. C., Medellin, R. A. & Jones, G. Conservation units and historical matrilineal structure in the tequila bat (Leptonycteris yerbabuenae). Glob. Ecol. Conserv. 23, e01164 (2020).Article 

    Google Scholar 
    Medellín, R. A. et al. Follow me: Foraging distances of Leptonycteris yerbabuenae (Chiroptera: Phyllostomidae) in Sonora determined by fluorescent powder. J. Mammal. 99, 306–311 (2018).Article 

    Google Scholar 
    Broennimann, O. et al. Evidence of climatic niche shift during biological invasion. Ecol. Lett. 10, 701–709 (2007).Article 
    CAS 

    Google Scholar 
    Martínez-Meyer, E., Peterson, A. T. & Hargrove, W. W. Ecological niches as stable distributional constraints on mammal species, with implications for Pleistocene extinctions and climate change projections for biodiversity. Glob. Ecol. Biogeogr. 13, 305–314 (2004).Article 

    Google Scholar 
    Soto-Centeno, J. A. & Steadman, D. W. Fossils reject climate change as the cause of extinction of Caribbean bats. Sci. Rep. 5, 7971 (2015).Article 
    ADS 
    CAS 

    Google Scholar 
    Avise, J. C. Phylogeography: The History and Formation of Species (Harvard University Press, 2000).Book 

    Google Scholar 
    Hickerson, M. J. et al. Phylogeography’s past, present, and future: 10 years after Avise, 2000. Mol. Phylogenet. Evol. 54, 291–301 (2010).Article 
    CAS 

    Google Scholar 
    Pahad, G., Montgelard, C. & Jansen van Vuuren, B. Phylogeography and niche modelling: Reciprocal enlightenment. Mammalia 84, 10–25 (2019).Article 

    Google Scholar 
    Flanders, J. et al. Phylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum: Contrasting results from mitochondrial and microsatellite data. Mol. Ecol. 18, 306–318 (2009).Article 
    CAS 

    Google Scholar 
    Machado, A. F. et al. Integrating phylogeography and ecological niche modelling to test diversification hypotheses using a Neotropical rodent. Evol. Ecol. 33, 111–148 (2019).Article 

    Google Scholar 
    Kalkvik, H. M., Stout, I. J., Doonan, T. J. & Parkinson, C. L. Investigating niche and lineage diversification in widely distributed taxa: Phylogeography and ecological niche modeling of the Peromyscus maniculatus species group. Ecography 35, 54–64 (2012).Article 

    Google Scholar 
    Wang, Y. et al. Ring distribution patterns—diversification or speciation? Comparative phylogeography of two small mammals in the mountains surrounding the Sichuan Basin. Mol. Ecol. 30, 2641–2658 (2021).Article 

    Google Scholar 
    Soto-Centeno, J. A., Barrow, L. N., Allen, J. M. & Reed, D. L. Reevaluation of a classic phylogeographic barrier: New techniques reveal the influence of microgeographic climate variation on population divergence. Ecol. Evol. 3, 1603–1613 (2013).Article 

    Google Scholar 
    Amador, L. I., Moyers Arévalo, R. L., Almeida, F. C., Catalano, S. A. & Giannini, N. P. Bat systematics in the light of unconstrained analyses of a comprehensive molecular supermatrix. J. Mamm. Evol. 25, 37–70 (2018).Article 

    Google Scholar 
    Rojas, D., Warsi, O. M. & Dávalos, L. M. Bats (Chiroptera: Noctilionoidea) challenge a recent origin of extant neotropical diversity. Syst. Biol. 65, 432–448 (2016).Article 

    Google Scholar 
    Shi, J. J. & Rabosky, D. L. Speciation dynamics during the global radiation of extant bats. Evolution 69, 1528–1545 (2015).Article 

    Google Scholar 
    Dumont, E. R. et al. Morphological innovation, diversification and invasion of a new adaptive zone. Proc. Biol. Sci. 279, 1797–1805 (2012).
    Google Scholar 
    Leiser-Miller, L. B. & Santana, S. E. Morphological diversity in the sensory system of phyllostomid bats: Implications for acoustic and dietary ecology. Funct. Ecol. 34, 1416–1427 (2020).Article 

    Google Scholar 
    Hedrick, B. P. & Dumont, E. R. Putting the leaf-nosed bats in context: A geometric morphometric analysis of three of the largest families of bats. J. Mammal. 99, 1042–1054 (2018).Article 

    Google Scholar 
    Clare, E. L. Cryptic species? Patterns of maternal and paternal gene flow in eight neotropical bats. PLoS One 6, e21460 (2011).Article 
    ADS 
    CAS 

    Google Scholar 
    Chaverri, G. et al. Unveiling the hidden bat diversity of a neotropical montane forest. PLoS One 11, e0162712 (2016).Article 

    Google Scholar 
    Calahorra-Oliart, A., Ospina-Garcés, S. M. & León-Paniagua, L. Cryptic species in Glossophaga soricina (Chiroptera: Phyllostomidae): Do morphological data support molecular evidence?. J. Mammal. 102, 54–68 (2021).Article 

    Google Scholar 
    Lim, B. K., Loureiro, L. O. & Garbino, G. S. T. Cryptic diversity and range extension in the big-eyed bat genus Chiroderma (Chiroptera, Phyllostomidae). Zookeys 918, 41–63 (2020).Article 

    Google Scholar 
    Loureiro, L. O., Engstrom, M., Lim, B., González, C. L. & Juste, J. Not all Molossus are created equal: Genetic variation in the mastiff bat reveals diversity masked by conservative morphology. Acta Chiropterologica 21, 51 (2019).Article 

    Google Scholar 
    Morales, A., Villalobos, F., Velazco, P. M., Simmons, N. B. & Piñero, D. Environmental niche drives genetic and morphometric structure in a widespread bat. J. Biogeogr. 43, 1057–1068 (2016).Article 

    Google Scholar 
    Hedrick, B. P. et al. Morphological diversification under high integration in a hyper diverse mammal clade. J. Mamm. Evol. 27, 563–575 (2020).Article 

    Google Scholar 
    Morales, A. E. & Carstens, B. C. Evidence that myotis lucifugus “subspecies” are five nonsister species, despite gene flow. Syst. Biol. 67, 756–769 (2018).Article 

    Google Scholar 
    Simmons, N. B. & Cirranello, A. L. Bat species of the world: A taxonomic and geographic database. https://batnames.org.Russell, A. L., Pinzari, C. A., Vonhof, M. J., Olival, K. J. & Bonaccorso, F. J. Two tickets to paradise: Multiple dispersal events in the founding of hoary bat populations in Hawai’i. PLoS One 10, 1–13 (2015).
    Google Scholar 
    Shump, K. A. & Shump, A. U. Lasiurus cinereus. Mamm. Species 185, 1–5 (1982).
    Google Scholar 
    Ziegler, A. C., Howarth, F. G. & Simmons, N. B. A second endemic land mammal for the Hawaiian Islands: A new genus and species of fossil bat (Chiroptera: Vespertilionidae). Am. Museum Novit. 1–52 (2016).Bonaccorso, F. J. & McGuire, L. P. Modeling the colonization of Hawaii by hoary bats (Lasiurus cinereus). In Bat Evolution, Ecology, and Conservation (eds Adams, R. A. & Pedersen, S. C.) 187–205 (Springer, 2013).Chapter 

    Google Scholar 
    Baird, A. B. et al. Molecular systematic revision of tree bats (Lasiurini): Doubling the native mammals of the Hawaiian Islands. J. Mammal. 96, 1255–1274 (2015).Article 

    Google Scholar 
    Jacobs, D. S. Morphological divergence in an insular bat, Lasiurus cinereus semotus. Funct. Ecol. 10, 622–630 (1996).Article 

    Google Scholar 
    Baird, A. B. et al. Nuclear and mtDNA phylogenetic analyses clarify the evolutionary history of two species of native Hawaiian bats and the taxonomy of Lasiurini (Mammalia: Chiroptera). PLoS One 12, e0186085 (2017).Article 

    Google Scholar 
    Kumar, S. & Subramanian, S. Mutation rates in mammalian genomes. Proc. Natl. Acad. Sci. U.S.A. 99, 803–808 (2002).Article 
    ADS 
    CAS 

    Google Scholar 
    Gillespie, R. G. et al. Comparing adaptive radiations across space, time, and taxa. J. Hered. 111, 1–20 (2020).Article 

    Google Scholar 
    Fišer, C., Robinson, C. T. & Malard, F. Cryptic species as a window into the paradigm shift of the species concept. Mol. Ecol. 27, 613–635 (2018).Article 

    Google Scholar 
    Espíndola, A. et al. Identifying cryptic diversity with predictive phylogeography. Proc. R. Soc. B Biol. Sci. 283, 20161529 (2016).Article 

    Google Scholar 
    Padial, J. M., Miralles, A., De la Riva, I. & Vences, M. The integrative future of taxonomy. Front. Zool. 7, 1–14 (2010).Article 

    Google Scholar 
    Fujita, M. K., Leaché, A. D., Burbrink, F. T., McGuire, J. A. & Moritz, C. Coalescent-based species delimitation in an integrative taxonomy. Trends Ecol. Evol. 27, 480–488 (2012).Article 

    Google Scholar 
    Solari, S., Sotero-Caio, C. G. & Baker, R. J. Advances in systematics of bats: Towards a consensus on species delimitation and classifications through integrative taxonomy. J. Mammal. 100, 838–851 (2018).Article 

    Google Scholar 
    Mayr, E. Geographical character gradients and climatic adaptation. Evolution 10, 105–108 (1956).
    Google Scholar 
    Morales, A. E., De-la-Mora, M. & Piñero, D. Spatial and environmental factors predict skull variation and genetic structure in the cosmopolitan bat Tadarida brasiliensis. J. Biogeogr. 45, 1529–1540 (2018).Article 

    Google Scholar 
    Pavan, A. C. & Marroig, G. Integrating multiple evidences in taxonomy: Species diversity and phylogeny of mustached bats (Mormoopidae: Pteronotus). Mol. Phylogenet. Evol. 103, 184–198 (2016).Article 

    Google Scholar 
    Kozlov, A. M., Darriba, D., Flouri, T., Morel, B. & Stamatakis, A. RAxML-NG: A fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 35, 4453–4455 (2019).Article 
    CAS 

    Google Scholar 
    Robinson, D. & Foulds, L. Comparison of phylogenetic trees. Math. Biosci. 53, 131–147 (1981).Article 
    MathSciNet 
    MATH 

    Google Scholar 
    Pattengale, N. D., Alipour, M., Bininda-Emonds, O. R., Moret, B. M. & Stamatakis, A. How many bootstrap replicates are necessary?. J. Comput. Biol. 17, 337–354 (2010).Article 
    MathSciNet 
    CAS 

    Google Scholar 
    Lemoine, F. et al. Renewing Felsenstein’s phylogenetic bootstrap in the era of big data. Nature 556, 452–456 (2018).Article 
    ADS 
    CAS 

    Google Scholar 
    Ronquist, F. et al. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 61, 539–542 (2012).Article 

    Google Scholar 
    Rambaut, A., Drummond, A. J., Xie, D., Baele, G. & Suchard, M. A. Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst. Biol. 67, 901–904 (2018).Article 
    CAS 

    Google Scholar 
    Kapli, P. et al. Multi-rate Poisson Tree Processes for single-locus species delimitation under Maximum Likelihood and Markov Chain Monte Carlo. Bioinformatics 33, 1630–1638 (2017).CAS 

    Google Scholar 
    Yang, Z. & Rannala, B. Unguided species delimitation using DNA sequence data from multiple loci. Mol. Biol. Evol. 31, 3125–3135 (2014).Article 
    CAS 

    Google Scholar 
    Flouri, T., Jiao, X., Rannala, B. & Yang, Z. Species tree inference with BPP using genomic sequences and the multispecies coalescent. Mol. Biol. Evol. 35, 2585–2593 (2018).Article 
    CAS 

    Google Scholar 
    Van Buuren, S. & Groothuis-Oudshoorn, K. Multivariate imputation by chained equations. J. Stat. Softw. 45, 1–67 (2011).Article 

    Google Scholar 
    Penone, C. et al. Imputation of missing data in life-history trait datasets: Which approach performs the best?. Methods Ecol. Evol. 5, 961–970 (2014).Article 

    Google Scholar 
    Berner, D. Size correction in biology: How reliable are approaches based on (common) principal component analysis?. Oecologia 166, 961–971 (2011).Article 
    ADS 

    Google Scholar 
    Simmons, N. B. Order Chiroptera. In Mammal Species of the World: A Taxonomic and Geographic Reference (eds Wilson, D. E. & Reeder, D. M.) 312–529 (The John Hopkins University Press, 2005).
    Google Scholar 
    Wilson, D. E. & Mittermeier, R. A. Handbook of the Mammals of the World. Vol. 9. Bats (Lynx Editions, 2019).
    Google Scholar 
    R Core Team. R: A language and environment for statistical computing (2022).Kuhn, M. caret: Classification and Regression Training. R package version 6.0-86. https://CRAN.R-project.org/package=caret (2020).Venables, W. N. & Ripley, B. D. Modern Applied Statistics with S (Springer, 2002).Book 
    MATH 

    Google Scholar 
    Kuhn, M. & Johnson, K. Applied Predictive Modeling (Springer, 2013).Book 
    MATH 

    Google Scholar 
    Fick, S. E. & Hijmans, R. J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315 (2017).Article 

    Google Scholar 
    Hijmans, R. J. raster: Geographic Data Analysis and Modeling (2022).Barker, B. S., Rodríguez-Robles, J. A. & Cook, J. A. Climate as a driver of tropical insular diversity: Comparative phylogeography of two ecologically distinctive frogs in Puerto Rico. Ecography 38, 769–781 (2015).Article 

    Google Scholar 
    Petitpierre, B., Broennimann, O., Kueffer, C., Daehler, C. & Guisan, A. Selecting predictors to maximize the transferability of species distribution models: Lessons from cross-continental plant invasions. Glob. Ecol. Biogeogr. 26, 275–287 (2017).Article 

    Google Scholar 
    Akinwande, M. O., Dikko, H. G. & Samson, A. Variance inflation factor: As a condition for the inclusion of suppressor variable(s) in regression analysis. Open J. Stat. 05, 754–767 (2015).Article 

    Google Scholar 
    Izenman, A. J. Linear discriminant analysis. in Modern Multivariate Statistical Techniques 237–280 (2013).Lever, J., Krzywinski, M. & Altman, N. Points of significance: Principal component analysis. Nat. Methods 14, 641–642 (2017).Article 
    CAS 

    Google Scholar 
    Guisan, A., Petitpierre, B., Broennimann, O., Daehler, C. & Kueffer, C. Unifying niche shift studies: Insights from biological invasions. Trends Ecol. Evol. 29, 260–269 (2014).Article 

    Google Scholar 
    Di Cola, V. et al. ecospat: An R package to support spatial analyses and modeling of species niches and distributions. Ecography 40, 774–787 (2017).Article 

    Google Scholar 
    Broennimann, O. et al. Measuring ecological niche overlap from occurrence and spatial environmental data. Glob. Ecol. Biogeogr. 21, 481–497 (2012).Article 

    Google Scholar 
    Liu, C., Wolter, C., Xian, W. & Jeschke, J. M. Most invasive species largely conserve their climatic niche. Proc. Natl. Acad. Sci. 117, 23643–23651 (2020).Article 
    ADS 
    CAS 

    Google Scholar 
    Warren, D. L., Glor, R. E. & Turelli, M. Environmental niche equivalency versus conservatism: Quantitative approaches to niche evolution. Evolution 62, 2868–2883 (2008).Article 

    Google Scholar 
    Warren, D. L., Glor, R. E. & Turelli, M. ENMTools: A toolbox for comparative studies of environmental niche models. Ecography 33, 607–611 (2010).
    Google Scholar  More