Philippa Kaur

Tsunekawa, A., Liu, G., Yamanaka, N. & Du, S. Restoration and Development of the Degraded Loess Plateau China 3–21 (Springer Press, 2017).
Google Scholar 
Kimura, R., Kamichika, M., Takayama, N., Matsuoka, N. & Zhang, X. C. Heat balance and soil moisture in the Loess Plateau, China. J. Agric Meteorol. 60(2), 103–113 (2004).
Google Scholar 
Deng, X. P., Lun, S., Zhang, H. P. & Turner, N. C. Improving agricultural water use efficiency in arid and semiarid areas of China. Agric. Water Manage. 80(1), 23–40 (2006).
Google Scholar 
Liu, C. A. et al. Maize yield and water balance is affected by nitrogen application in a film-mulching ridge-furrow system in a semiarid region of China. Eur. J. Agron. 52, 103–111 (2014).CAS 

Google Scholar 
Bu, L. D. et al. The effects of mulching on maize growth, yield and water use in a semi-arid region. Agric. Water Manage. 123, 71–78 (2013).
Google Scholar 
Hou, F. Y. et al. Effect of plastic mulching on the photosynthetic capacity, endogenous hormones and root yield of summer-sown sweet potato (Ipomoea batatas (L.) Lam.) in Northern China. Acta Physiol. Plant. 37, 164 (2015).
Google Scholar 
Jensen, K., Kimball, E. R. & Ricketson, C. L. Effect of perforated plastic row covers on residues of the herbicide DCPA in soil and broccoli. Environ Contam. Toxicol. B 35(6), 716–722 (1985).CAS 

Google Scholar 
Li, F. M., Guo, A. H. & Wei, H. Effects of clear plastic film mulch on yield of spring wheat. Field Crop. Res. 63(1), 79–86 (1999).
Google Scholar 
Liu, J. L. et al. Response of nitrogen use efficiency and soil nitrate dynamics to soil mulching in dryland maize (Zea mays L.) fields. Nutr. Cycl. Agroecosyst. 101(2), 271–283 (2015).CAS 

Google Scholar 
Li, R. et al. Effects on soil temperature, moisture, and maize yield of cultivation with ridge and furrow mulching in the rained area of the Loess Plateau, China. Agric. Water Manage. 116, 101–109 (2013).
Google Scholar 
Anzalone, A., Cirujeda, A., Aibar, J., Pardo, G. & Zaragoza, C. Effect of biodegradable mulch materials on weed control in processing tomatoes. Weed Technol. 24(3), 369–377 (2010).
Google Scholar 
Summers, C. G. & Stapleton, J. J. Use of UV reflective mulch to delay the colonization and reduce the severity of Bemisia argentifolii (Homoptera: Aleyrodidae) infestations in cucurbits. Crop Prot. 21(10), 921–928 (2002).
Google Scholar 
Chen, Y. S. et al. Empirical estimation of pollution load and contamination levels of phthalate esters in agricultural soils from plastic film mulching in China. Environ. Earth Sci. 70(1), 239–247 (2013).CAS 

Google Scholar 
Wang, S. Y. et al. Occurrence of macroplastic debris in the long-term plastic film-mulched agricultural soil: A case study of Northwest China. Sci. Total Environ. https://doi.org/10.1016/j.scitotenv.2022.154881 (2003).Article 
PubMed 

Google Scholar 
Hu, X. Y., Wen, B. & Shan, X. Q. Survey of phthalate pollution in arable soils in China. J. Environ Monit. 5(4), 649–653 (2003).CAS 
PubMed 

Google Scholar 
Zhou, X. Y. et al. Effects of residual mulch film on the growth and fruit quality of tomato (Lycopersicon esculentum Mill.). Water Air Soil. Pollut. 228(2), 1–18 (2017).ADS 

Google Scholar 
Hu, Q. et al. Effects of residual plastic-film mulch on field corn growth and productivity. Sci. Total Environ. 729, 1–10 (2020).
Google Scholar 
Wang, J. Z. et al. Crop yield and soil organic matter after long-term straw return to soil in China. Nutr. Cycl Agroecosyst. 102(3), 371–381 (2015).
Google Scholar 
Huang, Y. L., Chen, L. D., Fu, B. J., Huang, Z. L. & Gong, J. The wheat yields and water-use efficiency in the Loess Plateau: Straw mulch and irrigation effects. Agric. Water Manage. 72(3), 209–222 (2005).
Google Scholar 
Su, Z. Y. et al. Effects of conservation tillage practices on winter wheat water-use efficiency and crop yield on the Loess Plateau, China. Agric. Water Manage. 87(3), 307–314 (2007).
Google Scholar 
Ibrahim, A., Abaidoo, R. C., Fatondji, D. & Opoku, A. Integrated use of fertilizer micro-dosing and Acacia tumida mulching increases millet yield and water use efficiency in Sahelian semi-arid environment. Nutr. Cycl Agroecosys. 103(3), 375–388 (2015).CAS 

Google Scholar 
Zhang, D. K. et al. Suitable furrow mulching material for maize and sorghum production with ridge-furrow rainwater harvesting in semiarid regions of China. Agric. Water Manage. 228, 105928 (2020).
Google Scholar 
Myint, T. et al. Mulching improved soil water, root distribution and yield of maize in the Loess Plateau of Northwest China. Agric. Water Manage. 241, 106340 (2020).
Google Scholar 
Bai, Y. L. et al. Effects of long-term full straw return on yield and potassium response in wheat-maize rotation. J. Integr. Agric. 14(012), 2467–2476 (2015).CAS 

Google Scholar 
Liu, Z. J., Meng, Y., Cai, M. & Zhou, J. B. Coupled effects of mulching and nitrogen fertilization on crop yield, residual soil nitrate, and water use efficiency of summer maize in the Chinese Loess Plateau. Environ Sci. Pollut R. 24(33), 25849–25860 (2017).CAS 

Google Scholar 
Thomas, F. D., Michael, B., Jürgen, H., Maria, R. F. & Helmut, S. Effects of straw mulch on soil nitrate dynamics, weeds, and yield and soil erosion in organically grown potatoes. Field Crop. Res. 94(2–3), 238–249 (2005).
Google Scholar 
Tu, C., Ristaino, J. B. & Hu, S. J. Soil microbial biomass and activity in organic tomato farming systems: Effects of organic inputs and straw mulching. Soil Biol. Biochem. 38(2), 247–255 (2006).CAS 

Google Scholar 
Rao, Z. X. et al. Effect of rice straw mulching on migration and transportation of Cd, Cu, Zn, and Ni in surface runoff under simulated rainfall. J. Soils Sediment. 16(8), 2021–2029 (2016).CAS 

Google Scholar 
Ma, J., Xu, H., Yagi, K. & Cai, Z. C. Methane emission from paddy soils as affected by wheat straw returning mode. Plant Soil. 313, 167–174 (2008).CAS 

Google Scholar 
Xue, L. L. et al. Influence of straw mulch on yield, chlorophyll contents, lipid peroxidation and antioxidant enzymes activities of soybean under drought stress. J. Food Agric. Environ. 9(2), 699–704 (2011).
Google Scholar 
Wu, Y., Huang, F. Y., Jia, Z. K., Ren, X. R. & Cai, T. Response of soil water, temperature, and maize (Zea mays L.) production to different plastic film mulching patterns in semi-arid areas of Northwest China. Soil Tillage Res. 166, 113–121 (2017).
Google Scholar 
Blake, G. R. & Hartge, K. H. Bulk density. In Methods of Soil Analysis Part 1: Physical and Mineralogical Methods (ed. Klute, A.) 363–375 (American Society of Agronomy, Soil Science Society of America, 1986).
Google Scholar 
Li, F. M., Song, Q. H., Jjemba, P. & Shi, Y. Dynamics of soil microbial biomass and soil fertility in cropland mulched with plastic film in a semiarid agro-ecosystem. Soil Biol. Biochem. 36(11), 1893–1902 (2004).CAS 

Google Scholar 
Zhang, P. et al. Plastic-film mulching for enhanced water-use efficiency and economic returns from maize fields in semiarid China. Front. Plant Sci. 8, 512 (2017).PubMed 
PubMed Central 

Google Scholar 
Custodio, R. P. et al. The effects of increased temperature on crop growth and yield of soybean grown in a temperature gradient chamber. Field Crop. Res. 154, 74–81 (2013).
Google Scholar 
He, G. et al. Wheat yield affected by soil temperature and water under mulching in dryland. Agron. J. 109(6), 2998–3006 (2017).CAS 

Google Scholar 
Zhang, S. L. et al. Effects of mulching and catch cropping on soil temperature, soil moisture and wheat yield on the Loess Plateau of China. Soil Tillage Res. 102(1), 78–86 (2008).
Google Scholar 
Li, R., Hou, X. Q., Jia, Z. K. & Han, Q. F. Soil environment and maize productivity in semi-humid regions prone to drought of Weibei Highland are improved by ridge-and-furrow tillage with mulching. Soil Tillage Res. 196, 104476 (2020).
Google Scholar 
Wang, S. H. et al. Change in the bio-uptake of soil phthalates with increasing mulching years: Underlying mechanism and response to temperature rise. J. Clean Prod. 287(2021), 125049 (2020).
Google Scholar 
Li, W. W., Xiong, L., Wang, C. J., Liao, Y. C. & Wu, W. Optimized ridge–furrow with plastic film mulching system to use precipitation efficiently for winter wheat production in dry semi-humid areas. Agric. Water Manage. 218, 211–221 (2019).
Google Scholar 
Kader, M. A., Nakamura, K., Senge, M., Mojid, M. A. & Kawashima, S. Effects of colored plastic mulch on soil hydrothermal characteristics, growth and water productivity of rain-fed soybean. Irrig. Drain. 69(3), 483–494 (2020).
Google Scholar 
Luo, C. L. et al. Dual plastic film and straw mulching boosts wheat productivity and soil quality under the El Nino in semiarid Kenya. Sci. Total Environ. 738, 139808 (2020).ADS 
CAS 
PubMed 

Google Scholar 
Gouranga, K. & Ashwani, K. Effects of irrigation and straw mulch on water use and tuber yield of potato in Eastern India. Agric. Water Manage. 94(1), 109–116 (2007).
Google Scholar 
Lu, X. J., Li, Z. Z., Sun, Z. G. & Bu, Q. G. Straw mulching reduces maize yield, water, and nitrogen use in Northeastern China. Agron. J. 107(1), 406–414 (2015).
Google Scholar 
Zhou, L. F., Zhao, W. Z., He, J. Q., Flerchinger, G. N. & Feng, H. Simulating soil surface temperature under plastic film mulching during seedling emergence of spring maize with the RZ–SHAW and DNDC models. Soil Tillage Res. 197, 104517 (2020).
Google Scholar 
Chang, L. et al. Straw strip mulching affects soil moisture and temperature for potato yield in semiarid regions. Agron. J. 112(2), 1126–1139 (2020).CAS 

Google Scholar 
Zhang, P. et al. Effects of straw mulch on soil water and winter wheat production in dryland farming. Sci. Rep. 5(1), 209–222 (2015).
Google Scholar 
Ren, X. L., Zhang, P., Chen, X. L., Guo, J. J. & Jia, Z. K. Effect of different mulches under rainfall concentration system on corn production in the semi-arid areas of the Loess Plateau. Sci. Rep. 6(1), 47–50 (2016).
Google Scholar 
Akhtar, K. et al. Integrated use of straw mulch with nitrogen fertilizer improves soil functionality and soybean production. Environ. Int. 132, 105092 (2019).CAS 
PubMed 

Google Scholar 
Eden, G. R. S. M. The impact of organic amendments, mulching and tillage on plant nutrition, Pythium root rot, root-knot nematode and other pests and diseases of capsicum in a subtropical environment, and implications for the development of more sustainable vegetable farming. Australas. Plant Path. 37(2), 123–131 (2008).
Google Scholar 
Kader, M. A., Senge, M., Mojid, M. A., Takeo, O. & Kengo, I. Effects of plastic-hole mulching on effective rainfall and readily available soil moisture under soybean (Glycine max) cultivation. Paddy Water Environ. 15(3), 659–668 (2017).
Google Scholar 
Zhang, Z. et al. Plastic film cover during the fallow season preceding sowing increases yield and water use efficiency of rain-fed spring maize in a semi-arid climate. Agric. Water Manage. 212, 203–210 (2019).
Google Scholar 
Kader, M. A., Nakamura, K., Senge, M., Mojid, M. A. & Kawashima, S. Numerical simulation of water- and heat-flow regimes of mulched soil in rain-fed soybean field in central Japan. Soil Tillage Res. 191, 142–155 (2019).
Google Scholar 
Ryu, J. H. et al. Effects of straw mulching on soil physicochemical properties in Saemangeum reclaimed land. Korean J. Soil Sci. Fert. 49(1), 12–16 (2016).CAS 

Google Scholar 
Yin, W. et al. Growth trajectories of wheat–maize intercropping with straw and plastic management in arid conditions. Agron. J. 112(4), 2777–2790 (2020).
Google Scholar 
Wang, J. et al. Responses of runoff and soil erosion to planting pattern, row direction, and straw mulching on sloped farmland in the corn belt of northeast China. Agric. Water Manage. 25, 106935 (2021).
Google Scholar 
Cao, B. et al. Future landscape of renewable fuel resources: Current and future conservation and utilization of main biofuel crops in China. Sci. Total Environ. 806, 150946 (2022).ADS 
CAS 
PubMed 

Google Scholar 
Khawar, J. et al. Economic assessment of different mulches in conventional and water-saving rice production systems. Environ. Sci. Pollut. Res. 23, 9156–9163 (2016).
Google Scholar  More

Characteristics of LVD dataset and assembled vOTUsThe land use virome dataset LVD was derived from 2.6 billion paired clean reads of sequences across 50 viromes of 25 samples with five types of land uses (Supplementary Data 2). A total of 6,442,065 contigs ( >1500 bp) were yielded, of which 764,466 (11.8%) contigs were identified as putative viral genomes through VIBRANT. Subsequently, putative false positive viral genomes were removed (see Methods section), and 27,951 and 48,936 bona fide viral genomes were retained from the 25 intracellular VLPs (iVLPs) and 25 extracellular VLPs (eVLPs) viromes, respectively. These genomes were clustered into 25,941 and 45,152 vOTUs for iVLPs and eVLPs viromes, respectively, in which the iVLPs and eVLPs viromes shared 11,467 (19.2%) vOTUs. Subsequently, they were merged and dereplicated, resulting in 59,626 vOTUs (Supplementary Data 3) for the following analysis. A total of 8112 (13.6%) vOTUs genomes were classified as complete, in which the median length of all and circular vOTUs were 25,183 bp and 45,511 bp, respectively (Supplementary Fig. 4).To explore the taxonomic affiliation of vOTUs in family and genus-level, a gene-sharing network consist of 59,626 vOTUs genomes from this study and 3502 reference phage genomes (from NCBI Viral RefSeq version 201) revealed 6009 VCs comprising of 37,224 vOTUs, of which 34,417 vOTUs were from LVD, besides 2794 singletons (2653 from LVD dataset), 16,056 outliers (15,833 from LVD) and 8492 overlaps (8061 from LVD) were detected (Supplementary Data 4). Of these, only 157 VCs contained genomes from both the RefSeq and LVD dataset (1864 viral genomes) (Supplementary Data 4). Most of VCs (1837, 30.4%) included only two members.At the family level, most of vOTUs were classified into Siphoviridae (712 by vConTACT2 and 29,671 (50.9%) by Demovir, tailed dsDNA), Podoviridae (610 by vConTACT2 and 9923 (17.6 %) by Demovir, tailed dsDNA), Myoviridae (485 by vConTACT2 and 5445 (9.9%) by Demovir, tailed dsDNA), Tectiviridae (50 by vConTACT2 and 10 (0.10%) by Demovir, non-tailed dsDNA) (Fig. 1). Besides, the Eukaryotic viruses Herpesviridae (159 by Demovir, 0.26%, dsDNA), Phycodnaviridae (120 (0.20%) by Demovir, dsDNA); the Virophage Family Lavidaviridae (15 (0.03%) by Demovir) were detected as well, but a majority of vOTUs were unclassified in genus-level.Fig. 1: The taxonomic assignment of LVD.Pie charts showing the affiliation of 56,870 vOTUs at family level assigned by script Demovir (a). and the affiliation of 1864 vOTUs at family level assigned by package vConTACT2 (b). Source data are provided in the Source Data file.Full size imageViral community structures differ across land use typesBray–Curtis dissimilarity of viral communities (median 0.9951) showed strong heterogeneity of viral communities among different sites (Fig. 2a). While, the Bray–Curtis dissimilarity (median: 0.5109) between paired viral communities of iVLPs and eVLPs from each site have a significant lower heterogeneity than inter-sites (Wilcox.test, p  0.05; Fig. 2b). Therefore, the paired iVLPs and eVLPs viromes from each site were merged for subsequently viral community analysis.Fig. 2: The macrodiversity of soil viral communities.a Boxplot showing Bray–Curtis dissimilarity of viral communities of intra-sites (between the corresponding community of iVLPs and eVLPs, n = 25) and inter-sites (between different sample sites, n = 300). The minima, maxima, center, bounds of box and whiskers in boxplots from bottom to top represented percentile 0, 10, 25, 50, 75, 90, and 100, respectively, the difference between different zones was tested using the two-sided Wilcox.test, ****p  More

Walas, Ł, Ganatsas, P., Iszkuło, G., Thomas, P. A. & Dering, M. Spatial genetic structure and diversity of natural populations of Aesculus hippocastanum L. in Greece. PLoS ONE 14, e0226225 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Song, Y. G. et al. Past, present and future suitable areas for the relict tree Pterocarya fraxinifolia (Juglandaceae): Integrating fossil records, niche modeling, and phylogeography for conservation. Eur. J. For. Res. 140, 1323–1339 (2021).Article 

Google Scholar 
Dyderski, M. K., Paź, S., Frelich, L. E. & Jagodziński, A. M. How much does climate change threaten European forest tree species distributions?. Glob. Change Biol. 24, 1150–1163 (2018).ADS 
Article 

Google Scholar 
Chakraborty, D., Móricz, N., Rasztovits, E., Dobor, L. & Schueler, S. Provisioning forest and conservation science with high-resolution maps of potential distribution of major European tree species under climate change. Ann. For. Sci. 78, 1–18 (2021).Article 

Google Scholar 
Williams, J. N. et al. Using species distribution models to predict new occurrences for rare plants. Divers. Distrib. 15, 565–576 (2009).Article 

Google Scholar 
Watling, J. I. et al. Performance metrics and variance partitioning reveal sources of uncertainty in species distribution models. Ecol. Modell. 309, 48–59 (2015).ADS 
Article 

Google Scholar 
Phillips, S. J., Anderson, R. P. & Schapire, R. E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 190, 231–259 (2006).Article 

Google Scholar 
Phillips, S. J., Dudík, M. & Schapire, R. E. [Internet] Maxent software for modeling species niches and distributions. url: http://biodiversityinformatics.amnh.org/open_source/maxent/. Accessed 13 July 2022.Elith, J. et al. A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 17, 43–57 (2011).Article 

Google Scholar 
Marcer, A., Sáez, L., Molowny-Horas, R., Pons, X. & Pino, J. Using species distribution modelling to disentangle realised versus potential distributions for rare species conservation. Biol. Conserv. 166, 221–230 (2013).Article 

Google Scholar 
Rigling, A., Eilmann, B., Koechli, R. & Dobbertin, M. Mistletoe-induced crown degradation in Scots pine in a xeric environment. Tree Physiol. 30, 845–852 (2010).PubMed 
Article 

Google Scholar 
Sangüesa-Barreda, G., Linares, J. C. & Camarero, J. J. Mistletoe effects on Scots pine decline following drought events: Insights from within-tree spatial patterns, growth and carbohydrates. Tree Physiol. 32, 585–598 (2012).PubMed 
Article 

Google Scholar 
Kollas, C., Gutsch, M., Hommel, R., Lasch-Born, P. & Suckow, F. Mistletoe-induced growth reductions at the forest stand scale. Tree Physiol. 38, 735–744 (2018).PubMed 
Article 

Google Scholar 
Schulze, E. D. & Ehleringer, J. R. The effect of nitrogen supply on growth and water-use efficiency of xylem-tapping mistletoes. Planta 162, 268–275 (1984).PubMed 
Article 

Google Scholar 
Escher, P. et al. Transpiration, CO2 assimilation, WUE, and stomatal aperture in leaves of Viscum album L: Effect of abscisic acid (ABA) in the xylem sap of its host (Populus x euamericana). Plant Physiol. Biochem. 46, 64–70 (2008).PubMed 
Article 

Google Scholar 
Zweifel, R., Bangerter, S., Rigling, A. & Sterck, F. J. Pine and mistletoes: How to live with a leak in the water flow and storage system?. J. Exp. Bot. 63, 2565–2578 (2012).PubMed 
Article 

Google Scholar 
Mutlu, S., Osma, E., Ilhan, V., Turkoglu, H. I. & Atici, O. Mistletoe (Viscum album) reduces the growth of the Scots pine by accumulating essential nutrient elements in its structure as a trap. Trees 30, 815–824 (2016).Article 

Google Scholar 
Tsopelas, P., Angelopoulos, A., Economou, A. & Soulioti, N. Mistletoe (Viscum album) in the fir forest of Mount Parnis Greece. For. Ecol. Manag. 202, 59–65 (2004).Article 

Google Scholar 
Dobbertin, M. & Rigling, A. Pine mistletoe (Viscum album ssp. austriacum) contributes to Scots pine (Pinus sylvestris) mortality in the Rhone valley of Switzerland. For. Pathol. 36, 309–322 (2006).Article 

Google Scholar 
Lech, P., Żółciak, A. & Hildebrand, R. Occurrence of European mistletoe (Viscum album L.) on forest trees in Poland and its dynamics of spread in the period 2008–2018. Forests 11, 83 (2020).Article 

Google Scholar 
Iszkuło, G. et al. Jemioła jako zagrożenie dla zdrowotności drzewostanów iglastych. Sylwan 164, 226–236 (2020) ([In Polish]).
Google Scholar 
Mellado, A., Morillas, L., Gallardo, A. & Zamora, R. Temporal dynamic of parasite-mediated linkages between the forest canopy and soil processes and the microbial community. New Phytol. 211, 1382–1392 (2016).PubMed 
Article 

Google Scholar 
Mellado, A. & Zamora, R. Generalist birds govern the seed dispersal of a parasitic plant with strong recruitment constraints. Oecologia 176, 139–147 (2014).ADS 
PubMed 
Article 

Google Scholar 
Hódar, J. A., Lázaro-González, A. & Zamora, R. Beneath the mistletoe: parasitized trees host a more diverse herbaceous vegetation and are more visited by rabbits. Ann. For. Sci. 75, 1–8 (2018).Article 

Google Scholar 
Zuber, D. Biological flora of Central Europe: Viscum album L. Flora Morphol. Distrib Funct. Ecol. Plants 199, 181–203 (2004).Article 

Google Scholar 
Urech, K. & Baumgartner, S. Chemical constituents of Viscum album L.: Implications for the pharmaceutical preparation of mistletoe. In: Mistletoe: From mythology to evidence-based medicine. (eds. Zänker, K.S. & Kaveri, S. V.), 11–23. (S. Karger AG, Basel, Switzerland, 2015).Singh, B. N. et al. European Viscum album: a potent phytotherapeutic agent with multifarious phytochemicals, pharmacological properties and clinical evidence. RSC Adv. 6, 23837–23857 (2016).ADS 
Article 

Google Scholar 
Jeffree, C. E. & Jeffree, E. P. Redistribution of the potential geographical ranges of mistletoe and colorado beetle in Europe in response to the temperature component of climate change. Funct. Ecol. 10, 562–577 (1996).Article 

Google Scholar 
Troels-Smith, J. Ivy, mistletoe and elm climate indicators-fodder plants. A contribution to the interpretation of the pollen zone border VII-VIII. Dan. Geol. Undersøg. IV Række 4, 1–32 (1960).
Google Scholar 
Dobbertin, M. et al. The upward shift in altitude of pine mistletoe (Viscum album ssp. austriacum) in Switzerland—the result of climate warming?. Int. J. Biometeorol. 50, 40–47 (2005).ADS 
PubMed 
Article 

Google Scholar 
Zamora, R. & Mellado, A. Identifying the abiotic and biotic drivers behind the elevational distribution shift of a parasitic plant. Plant Biol. 21, 307–317 (2019).PubMed 
Article 

Google Scholar 
Barney, C. W., Hawksworth, F. G. & Geils, B. W. Hosts of Viscum album. Eur. J. Plant Pathol. 28, 187–208 (1998).
Google Scholar 
Böhling, N. et al. Notes on the Cretan mistletoe, Viscum album subsp. creticum subsp. nova (Loranthaceae/Viscaceae). Isr. J. Plant Sci. 50, 77–84 (2002).
Google Scholar 
Plants of the World Online [Internet] url: https://powo.science.kew.org/taxon/urn:lsid:ipni.org:names:921668-1. Accessed 13 July 2022.Zuber, D. & Widmer, A. Phylogeography and host race differentiation in the European mistletoe (Viscum album L.). Mol. Ecol. 18, 1946–1962 (2009).PubMed 
Article 

Google Scholar 
Schaller, G., Urech, K., Grazi, G. & Giannattasio, M. Viscotoxin composition of the three European subspecies of Viscum album. Planta Med 64, 677–678 (1998).PubMed 
Article 

Google Scholar 
Kahle-Zuber, D. Biology and evolution of the European mistletoe (Viscum album). Doctoral Thesis. ETH Zurich. (2008).Zuber, D. & Widmer, A. Genetic evidence for host specificity in the hemi-parasitic Viscum album L. (Viscaceae). Mol. Ecol. 9, 1069–1073 (2000).PubMed 
Article 

Google Scholar 
Mejnartowicz, L. Relationship and genetic diversity of mistletoe [Viscum album L.] subspecies. Acta Soc. Bot. Pol. Pol. 75, 39–49 (2006).Article 

Google Scholar 
Xie, W., Adolf, J. & Melzig, M. F. Identification of Viscum album L. miRNAs and prediction of their medicinal values. PLoS ONE 12, e0187776 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Valle, A. C. V., de Carvalho, A. C. & Andrade, R. V. Viscum album-literature review. Int. J. Sci. Res 10, 63–71 (2021).
Google Scholar 
Schröder, L. et al. The gene space of European mistletoe (Viscum album). Plant J. 109, 278–294 (2022).PubMed 
Article 

Google Scholar 
Sangüesa-Barreda, G. et al. Delineating limits: Confronting predicted climatic suitability to field performance in mistletoe populations. J. Ecol. 106, 2218–2229 (2018).Article 

Google Scholar 
GBIF.org [Internet] GBIF Occurrence Download Doi: https://doi.org/10.15468/dl.zw6f5q. Accessed 27 July 2021.GBIF.org [Internet] GBIF Occurrence Download Doi: https://doi.org/10.15468/dl.6wmc9d. Accessed 6 August 2021.FloraWeb [Internet] url: https://www.floraweb.de. Accessed 10 December 2021.Pladias – Database of the Czech Flora and Vegetation. [Internet] url: www.pladias.cz. Accessed 14 July 2022.Zając, A., Zając, M., Tertil, R. & Harman, I. Atlas rozmieszczenia roślin naczyniowych w Polsce. 593 (Instytut Botaniki Uniwersytetu Jagiellońskiego, Kraków, 2001) [In Polish].Idžojtić, M., Kogelnik, M., Franjić, J. & Škvorc, Ž. Hosts and distribution of Viscum album L. ssp. album in Croatia and Slovenia. Plant Biosyst. 140, 50–55 (2006).Article 

Google Scholar 
Varga, I. et al. Changes in the Distribution of European Mistletoe (Viscum album) in Hungary During the Last Hundred Years. Folia Geobot 49, 559–577 (2014).Article 

Google Scholar 
Wild, J. et al. Plant distribution data for the Czech Republic integrated in the Pladias database. Preslia 91, 1–24 (2019).Article 

Google Scholar 
Krasylenko, Y. et al. The European mistletoe (Viscum album L.): Distribution, host range, biotic interactions, and management worldwide with special emphasis on Ukraine. Botany 98, 499–516 (2020).Article 

Google Scholar 
Karger, D. N. et al. Climatologies at high resolution for the Earth land surface areas. Sci. Data 4, 170122 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Karger D. N., et al. Data from: Climatologies at high resolution for the earth’s land surface areas. Dryad Digital Repository (2018).Gutjahr, O. et al. Max planck institute earth system model (MPI-ESM1. 2) for the high-resolution model intercomparison project (HighResMIP). Geosci. Model Dev. 12, 3241–3281 (2019).ADS 
Article 

Google Scholar 
Hijmans, R. J., & van Etten, J. raster: Geographic analysis and modeling with raster data. R package version 2.0-12. (2012).R Core Team. The Comprehensive R Archive Network. [Internet] url: https://cran.r-project.org/ Accessed 14 July 2022.Chakraborty, D., Móricz, N., Rasztovits, E., Dobor, L. & Schueler, S. Provisioning forest and conservation science with European tree species distribution models under climate change (Version v1). Zenodo https://doi.org/10.5281/zenodo.3686918 (2020).Wang, Z., Chang, Y. I., Ying, Z., Zhu, L. & Yang, Y. A parsimonious threshold-independent protein feature selection method through the area under receiver operating characteristic curve. Bioinformatics 23, 2788–2794 (2007).PubMed 
Article 

Google Scholar 
Lobo, J. M., Jiménez-Valverde, A. & Hortal, J. The uncertain nature of absences and their importance in species distribution modelling. Ecography 33, 103–114 (2010).Article 

Google Scholar 
QGIS Development Team. QGIS Geographic Information Sys-tem. Open Source Geospatial Foundation Project. [Internet]. url: https://www.qgis.org/en/site/. Accessed 14 July 2022.Fischer, J. T. Water relations of mistletoes and their hosts. In: The biology of mistletoes. (eds. Calder, M., & Bernhard, T.), 163–184 (Academic Press, Sydney, 1983).Skre, O. The regional distribution of vascular plants in Scandinavia with requirements for high summer temperatures. Norweg. J. Bot. 26, 295–318 (1979).
Google Scholar 
Wangerin, B. Loranthaceae. In: Lebensgeschichte der Blütenpflanzen Mitteleuropas (eds. Kirchner, O. V., Loew, E., & Schroeter, C.) 2, 953–1146 (E. Ulmer, Stuttgart, 1937).Rybalka, I. A. Relationship between density of the white mistletoe (Viscum album L.) and some landscape and environmental characteristics of urban areas in the case of Kharkiv. Ekologicheskiy Vestnik 1, 87–97 (2017).
Google Scholar 
Patykowski, J. & Kołodziejek, J. Comparative analysis of antioxidant activity in leaves of different hosts infected by mistletoe (Viscum album L. subsp. album). Arch. Biol. Sci. 65, 851–861 (2013).Article 

Google Scholar 
Skrypnik, L., Maslennikov, P., Feduraev, P., Pungin, A. & Belov, N. Ecological and landscape factors affecting the spread of European mistletoe (Viscum album L.) in urban areas (A Case Study of the Kaliningrad City, Russia). Plants 9, 394 (2020).PubMed Central 
Article 

Google Scholar 
Kunick, W. Veränderungen von Flora und Vegetation einer Grosstadt dargestellt am Beispiel von Berlin (West). PhD Thesis, Technische Universität (1974). [In German].Kołodziejek, J., Patykowski, J. & Kołodziejek, R. Distribution, frequency and host patterns of European mistletoe (Viscum album subsp. album) in the major city of Lodz Poland. Biol. 68, 55–64 (2013).
Google Scholar 
Caudullo, G., Welk, E. & San-Miguel-Ayanz, J. Chorological maps for the main European woody species. Data Brief 12, 662–666 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
O’Donnell, M. S. & Ignizio, D. A. Bioclimatic predictors for supporting ecological applications in the conterminous United States. US Geol. Surv. Data Ser. 691, 4–9 (2012).
Google Scholar 
Luther, P., Becker, H. & Leroi, R. Die Mistel: Botanik, Lektine, medizinische Anwendung. Springer (1987).Gazol, A. et al. Distinct effects of climate warming on populations of silver fir (Abies alba) across Europe. J. Biogeogr. 42, 1150–1162 (2015).Article 

Google Scholar 
Tikkanen, O. P. et al. Freezing tolerance of seeds can explain differences in the distribution of two widespread mistletoe subspecies in Europe. For. Ecol. Manag. 482, 118806 (2021).Article 

Google Scholar 
Pilichowski, S. et al. Wpływ Viscum album ssp. austriacum (Wiesb.) Vollm. na przyrost radialny Pinus sylvestris L. Sylwan 162, 452–459 (2018) ([In Polish]).
Google Scholar 
Szmidla, H., Tkaczyk, M., Plewa, R., Tarwacki, G. & Sierota, Z. Impact of common mistletoe (Viscum album L.) on scots pine forests—A call for action. Forests 10, 847 (2019).Article 

Google Scholar 
Wójcik, R. & Kędziora, W. Abundance of Viscum in central Poland: Results from a large-scale mistletoe inventory. Environ. Sci. Proc. 3, 98 (2020).
Google Scholar 
Sangüesa-Barreda, G., Linares, J. C. & Camarero, J. J. Drought and mistletoe reduce growth and water-use efficiency of Scots pine. For. Ecol. Manag. 296, 64–73 (2013).Article 

Google Scholar 
Mathiasen, R. L., Nickrent, D. L., Shaw, D. C. & Watson, D. M. Mistletoes: Pathology, systematics, ecology, and management. Plant Dis. 92, 988–1006 (2008).PubMed 
Article 

Google Scholar 
Catal, Y. & Carus, S. Effect of pine mistletoe on radial growth of crimean pine (Pinus nigra) in Turkey. J. Environ. Biol. 32, 263 (2011).PubMed 

Google Scholar 
Skre, O. High temperature demands for growth and development in Norway Spruce [Picea abies (L.) Karst.] in Scandinavia. Meld Nor Landbrukshøgsk 51, 1–29 (1971).
Google Scholar 
Utaaker, K. A temperature-growth index—the respiration equivalent—used in climatic studies on the meso-scale in Norway. Agric. Meteorol. 5, 351–359 (1968).Article 

Google Scholar 
Iversen, J. Viscum, Hedera and Ilex as climate indicators: A contribution to the study of the post-glacial temperature climate. Geol. fören. Stockh. förh. 66, 463–483 (1944).Article 

Google Scholar 
Briggs, J. Mistletoe, Viscum album (Santalaceae), in Britain and Ireland; a discussion and review of current status and trends. Brit. Ir. Bot. 3, 419–454 (2021).
Google Scholar  More

Bartz, K. K. & Naiman, R. J. Effects of Salmon-Borne nutrients on riparian soils and vegetation in Southwest Alaska. Ecosystems 8, 529–545 (2005).Article 

Google Scholar 
Erskine, P. D. et al. Subantarctic Macquarie Island—a model ecosystem for studying animal-derived nitrogen sources using 15N natural abundance. Oecologia 117, 187–193 (1998).ADS 
PubMed 
Article 

Google Scholar 
Hocking, M. D. & Reimchen, T. E. Salmon species, density and watershed size predict magnitude of marine enrichment in riparian food webs. Oikos 118(9), 1307–1318 (2009).Article 

Google Scholar 
Hocking, M. D. & Reynolds, J. D. Impacts of salmon on riparian plant diversity. Science 331, 1609–1612 (2011).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Hocking, M. D., & Reimchen, T. E. Salmon-derived nitrogen in terrestrial invertebrates from coniferous forests of the Pacific Northwest. BMC Ecol. 2, 4. https://doi.org/10.1186/1472-6785-2-4 (2002).Bilby, R. E., Fransen, B. R. & Bisson, P. A. Incorporation of nitrogen and carbon from spawning coho salmon into the trophic system of small streams: Evidence from stable isotopes. Can. J. Fish Aquat. Sci. 53, 164–173 (1996).Article 

Google Scholar 
Talley, D. M. et al. Research challenges at the land–sea interface. Estuar. Coast. Shelf Sci. 58, 699–702 (2003).ADS 
Article 

Google Scholar 
Mizutani, H. & Wada, E. Nitrogen and carbon isotope ratios in seabird rookeries and their ecological implications. Ecology 69(2), 340–349 (1988).Article 

Google Scholar 
Rowe, J. A., Litton, C. M., Lepczyk, C. A. & Popp, B. N. Impacts of endangered seabirds on nutrient cycling in montane forest ecosystems of Hawai’i. Pac. Sci. 71(4), 495–509 (2017).Article 

Google Scholar 
Sanchez-Pinero, F. & Polis, G. A. Bottom-up dynamics of allochthonous input: Direct and indirect effects of seabirds on islands. Ecology 81(11), 3117–3132 (2000).Article 

Google Scholar 
Wait, D. A., Aubrey, D. P. & Anderson, W. B. Seabird guano influences on desert islands: Soil chemistry and herbaceous species richness and productivity. J. Arid Environ. 60, 681–695 (2005).ADS 
Article 

Google Scholar 
Stapp, P., Polis, G. A. & Pinero, F. S. Stable isotopes reveal strong marine and El Nino effects on island food webs. Nature 401, 467–469 (1999).ADS 
CAS 
Article 

Google Scholar 
Anderson, W. B., Wait, D. A. & Stapp, P. Resources from another place and time: Responses to pulses in a spatially subsidized system. Ecology 89(3), 660–670 (2008).PubMed 
Article 

Google Scholar 
Ellis, J. C. Marine birds on land: A review of plant biomass, species richness, and community composition in seabird colonies. Plant Ecol. 181(2), 227–241 (2005).Article 

Google Scholar 
Fukami, T. et al. Above- and below-ground impacts of introduced predators in seabird-dominated island ecosystems. Ecol. Lett. 9, 1299–1307 (2006).PubMed 
Article 

Google Scholar 
Wootton, J. T. Direct and indirect effects of nutrients on intertidal community structure: Variable consequences of seabird guano. J. Exp. Mar. Biol. Ecol. 151, 139–153 (1991).Article 

Google Scholar 
McCauley, D. J., et al., From wing to wing: the persistence of long ecological interaction chains in less-disturbed ecosystems. Sci. Rep. 2, 409. https://doi.org/10.1038/srep00409 (2012).Young, H. S., McCauley, D. J., Dunbar, R. B. & Dirzo, R. Plants cause ecosystem nutrient depletion via the interruption of bird-derived spatial subsidies. Proc. Natl. Acad. Sci. U.S.A. 107(5), 2072–2077 (2010).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lindeboom, H. J. The nitrogen pathway in a Penguin rookery. Ecology 65(1), 269–277 (1984).CAS 
Article 

Google Scholar 
Mizutani, H., Kabaya, Y. & Wada, E. Ammonia volatilization and high 15N/14N ratio in a penguin rookery in Antarctica. Geochem. J. 19(6), 323–327 (1985).ADS 
CAS 
Article 

Google Scholar 
Anderson, W. B. & Polis, G. A. Nutrient fluxes from water to land: seabirds affect plant nutrient status on Gulf of California islands. Oecologia 118, 324–332 (1999).ADS 
PubMed 
Article 

Google Scholar 
Polis, G. A. & Hurd, S. D. Linking marine and terrestrial food webs: Allochthonous input from the ocean supports high secondary productivity on small islands and coastal land communities. Am. Nat. 147, 396–423 (1996).Article 

Google Scholar 
Goss, N. S. New and rare birds found breeding on the San Pedro Martir Isle. University of California Press 5, 240–244 (1888).
Google Scholar 
Velarde, E., et al., Nesting seabirds of the Gulf of California’s Offshore islands: Diversity, ecology and conservation. in Biodiversity, Ecosystems, and Conservation in Northern Mexico, Carton, J.-L. E., Ceballos, G., Felger, R. S. Eds. (Oxford University Press, 2005) pp. 452–470.Wilder, B. T., Felger, R. S. & Ezcurra, E. Controls of plant diversity and composition on a desert archipelago. PeerJ 7, e7286. https://doi.org/10.7717/peerj.7286 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Ellis, J., Fariña, J. & Witman, J. Nutrient transfer from sea to land: the case of gulls and cormorants in the Gulf of Maine. J. Anim. Ecol. 75, 565–574 (2006).PubMed 
Article 

Google Scholar 
Wilder, B. T., Felger, R. S. & Morales, H. R. Succulent plant diversity of the Sonoran Islands, Gulf of California Mexico. Haseltonia 2008(14), 127–160 (2008).Article 

Google Scholar 
Lucassen, F. et al. The stable isotope composition of nitrogen and carbon and elemental contents in modern and fossil seabird guano from Northern Chile—Marine sources and diagenetic effects. PLoS ONE 12(6), e0179440. https://doi.org/10.1371/journal.pone.0179440 (2017).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Robinson, D. δ15N as an integrator of the nitrogen cycle. Trends Ecol. Evol. 16(3), 153–162 (2001).CAS 
PubMed 
Article 

Google Scholar 
Szpak, P., Longstaffe, F. J., Millaire, J.-F. & White, C. D. Stable isotope biogeochemistry of seabird guano fertilization: Results from growth chamber studies with maize (Zea mays). PLoS ONE 7(3), e33741. https://doi.org/10.1371/journal.pone.0033741 (2012).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Ezcurra, E., et al. Natural History and Evolution of the World’s Deserts. Global Deserts Outlook. United Nations Environment Programme (UNEP), 1–26 (2006).Yetman, D. The Great Cacti: Ethnobotany and biogeography (University of Arizona Press, 2007).
Google Scholar 
Álvarez-Borrego, S. Physical oceanography. in A New Island Biogeography of the Sea of Cortés, Case, T. J., Cody, M. L., Ezcurra, E. Eds. (Oxford University Press, 2002), pp. 41–59.Douglas, R., Gonzalez-Yajimovich, O., Ledesma-Vazquez, J. & Staines-Urias, F. Climate forcing, primary production and the distribution of Holocene biogenic sediments in the Gulf of California. Quatern. Sci. Rev. 26, 115–129 (2007).ADS 
Article 

Google Scholar 
Urbán, J. Marine mammals of the Gulf of California: An overview of diversity and conservation status. in The Gulf of California: Biodiversity and conservation, R. C. Brusca, Ed. (The University of Arizona Press and the Arizona-Sonora Desert Museum, 2010), pp. 188–209.Hastings, P. A., Findley, L. T., & Van der Heiden, A. M. Fishes of the Gulf of California. in: Brusca, R. C., (eds) The Gulf of California: Biodiversity and conservation 96–118, The University of Arizona Press and the Arizona-Sonora Desert Museum (2010).
Google Scholar 
Polis, G. A., Hurd, S. D., Jackson, C. T. & Sanchez Piñero, F. El Niño effects on the dynamics and control of an Island ecosystem in the Gulf of California. Ecology 78, 1884–1897 (1997).
Google Scholar 
Wilder, B. T. & Felger, R. S. Dwarf giants, guano, and isolation: The flora and vegetation of San Pedro Mártir Island, Gulf of California, Mexico. Proc. San Diego Soc. Nat. Hist. 42, 1–24 (2010).
Google Scholar 
Medel-Narvaez, A., Leon Luz, J. L., Freaner-Martinez, F. & Molina-Freaner, F. Patterns of abundance and population structure of Pachycereus pringlei (Cactaceae), a columnar cactus of the Sonoran Desert. Plant Ecol. 187, 1–14 (2006).Article 

Google Scholar 
Felger, R.S., Wilder, B.T. in collaboration with Romero-Morales, H. Plant Life of a Desert Archipelago: Flora of the Sonoran Islands in the Gulf of California. Tucson, University of Arizona Press (2012).Wilkinson, C. E., Hocking, M. D. & Reimchen, T. E. Uptake of salmon-derived nitrogen by mosses and liverworts in Coastal British Columbia. Oikos 108, 85–98 (2005).CAS 
Article 

Google Scholar 
Barrett, K., Wait, D. A. & Anderson, W. B. Small island biogeography in the Gulf of California: Lizards, the subsidized island biogeography hypothesis, and the small island effect. J. Biogeogr. 30, 1575–1581 (2003).Article 

Google Scholar 
Young, H. S., McCauley, D. J. & Dirzo, R. Differential responses to guano fertilization among tropical tree species with varying functional traits. Am. J. Bot. 98, 207–214 (2011).PubMed 
Article 

Google Scholar 
Nobel, P. S. Environmental Biology of Agaves and Cacti. Cambridge University Press (2003).Ramirez, K. S., Craine, J. M. & Fierer, N. Consistent effects of nitrogen amendments on soil microbial communities and processes across biomes. Glob. Change Biol. 18(6), 1918–1927 (2012).ADS 
Article 

Google Scholar 
Craine, J. M. et al. Ecological interpretations of nitrogen isotope ratios of terrestrial plants and soils. Plant Soil 396, 1–26 (2015).CAS 
Article 

Google Scholar 
Schoeninger, M. J. & DeNiro, M. J. Nitrogen and carbon isotope composition of bone collagen from marine and terrestrial animals. Geochim. Cosmochim. Acta 48(4), 625–639 (1984).ADS 
CAS 
Article 

Google Scholar 
Amundson, R. et al. Global patterns of the isotopic composition of soil and plant nitrogen. Global Biogeochem. Cycles 17(1), 1031. https://doi.org/10.1029/2002GB001903 (2003).ADS 
CAS 
Article 

Google Scholar 
Kahmen, A., Wanek, W. & Buchmann, N. Foliar δ15N values characterize soil N cycling and reflect nitrate or ammonium preference of plants along a temperate grassland gradient. Oecologia 156, 861–870 (2008).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bowen, T. Unknown Island: Seri Indians, Europeans, and San Esteban Island in the Gulf of California (University of New Mexico Press, 2000).
Google Scholar 
Evans, R. D. Physiological mechanisms influencing plant nitrogen isotope composition. Trends Plant Sci. 6(3), 121–126 (2001).CAS 
PubMed 
Article 

Google Scholar 
Dolby, G., Bennett, S. E. K., Lira-Noriega, A., Wilder, B. T. & Munguia-Vega, A. Assessing the geological and climatic forcing of biodiversity and evolution surrounding the Gulf of California. J. Southw. 57, 391–455 (2015).Article 

Google Scholar 
Case, T. J., Cody, M. L., & Ezcurra, E. A New Island Biogeography of the Sea of Cortés (Oxford University Press, 2002).Book 

Google Scholar 
Tershy, B. R. & Breese, D. The birds of San Pedro Mártir Island, Gulf of California Mexico. West. Birds 28, 96–107 (1997).
Google Scholar 
Tershy, B. R., Breese, D. & Croll, D. A. Human perturbations and conservation strategies for San Pedro Mártir Island, Islas de Golfo de California Reserve México. Environ. Conserv. 24, 261–270 (1997).Article 

Google Scholar 
Wilder, B. T. Historical biogeography of the Midriff Islands in the Gulf of California, Mexico. Dissertation. Riverside: UC, Riverside (2014).Post, D. M. et al. Getting to the fat of the matter: Models, methods and assumptions for dealing with lipids in stable isotope analyses. Oecologia 152, 179–189 (2007).ADS 
PubMed 
Article 

Google Scholar 
Kiljunen, M. et al. A revised model for lipid-normalizing δ13C values from aquatic organisms, with implications for isotope mixing models. J. Appl. Ecol. 43, 1213–1222 (2006).CAS 
Article 

Google Scholar 
Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. B. lmerTest Package: Tests in linear mixed effects models. J. Stat. Softw. 82(13), 1–26. https://doi.org/10.18637/jss.v082.i13 (2017).Article 

Google Scholar 
R Core Team, R: A language and environment for statistical computing. https://www.R-project.org/ (R Foundation for Statistical Computing, Vienna, Austria, 2022). More

Study site and field data collectionPermanent 1 m2 photoquadrats were established on Heron Reef in 1962/63, using 9 mm diameter mild steel (rebar) pegs, which were replaced over time. From the 1990’s, replacement pegs were stainless steel for greater longevity. Four sites were established, the protected (south) crest, inner flat, exposed (north) crest and exposed pools. Co-ordinates for each site are presented in Table 1, the layout shown in Fig. 2, and sites have been well described previously5,6. At each census, a 1 m2 frame divided into a 5 × 5 grid using string was placed over the pegs, and the quadrat photographed from directly above at low tide. From 1963 until 2003, a 35 mm camera and colour slide film were used. The camera was attached to a tripod affixed to the 1 m2 frame, and captured around 2/3 of the quadrat. The frame (and camera) were then rotated 180 degrees to capture the remainder of the quadrat. After 2003, a hand-held digital camera was used, with the entire quadrat being captured in a single image. Concurrent with each census, mud maps of each quadrat were hand drawn in the field, and all colonies identified in situ by someone with expertise in coral taxonomy.Table 1 Coordinates of the study sites on Heron Island Reef (WGS84).Full size tableFig. 2Quadrat layouts for each of the four sites respectively, noting that the north crest and north ridge have been treated as a single north crest site in previous publications. Underlining indicates original 1962/63 quadrats. Other quadrats were added in or after 2008, as indicated in the text. Contiguous quadrats are pictured bordering each other. Spacing between separate quadrats or groups of quadrats is not shown to scale. Note that up until 2005, NRNW was known as NR. The acronyms in each quadrat represent its name.Full size imageAt the protected (south) crest, a set of six contiguous quadrats were established in 1963 in a 2 × 3 arrangement parallel to the waterline, and about 420 m southeast of the island. This site is exposed at low tide, and was photographed once all water had drained off it. Images of quadrats A, C & E (the shoreward row) from 1963 to 2012 have been fully processed, and the data have been through QA/QC. Data after 2012 exist as images only. These quadrats form the basis of previous analyses1,4,5,6 for this site. Photographs are available for quadrats B, D & F, but apart from 2003–2010, have not been processed. In 2010, an additional two quadrats were established either side of the original six, leading to a 2 × 5 arrangement. Again, only imagery is available for these additional quadrats.At the inner flat, two pairs of contiguous quadrats were established in 1962, 44 m apart, about 70 m south of the island. This site is covered by ~10 cm of water at low tide, so could only be photographed on a still day. Imagery for this site is only available to 2012, after which the marker stakes appear to have been removed in a cleanup of the area. Images for one quadrat in each pair have been processed, but have not been subject to full QA/QC.At the exposed (north) crest main site, a set of four contiguous quadrats was established about 1100 m northeast of the island in 1963. An additional single quadrat (north ridge) was established 326 m to the east. Images from 1963 to 2012 have been fully processed, and the data have been through QA/QC. Data after 2012 exist as images only. In 2005, the single north ridge quadrat was expanded to 4 m2, and in 2008, both subsites were expanded to six quadrats in a 2 × 3 arrangement. These additional quadrats have been digitised up to 2012, but have not been through full QA/QC.The exposed pools are two individual quadrats about 5 m apart about 30 m north of the eastern (north ridge) exposed crest site. These are on the edge of a natural pool, and range from ~5–50 cm deep at low tide, and so could only be photographed on a calm day. Imagery for this site is only available until 2005, after which the marker stakes could not be relocated. Images from 1963 to 1998 have been processed, but have not been through full QA/QC.Retrieval of coral composition data from the photoquadratsProcessing of the images involved scanning the colour slides to produce digital images, and then orthorectifying each image to a 1 m2 basemap in ArcGIS (ESRI Ltd). The corners of the frame, and the holes for the string grid, were used as control points for the orthorectification. For images that originated as colour slides, each half of the quadrat was individually orthorectified to the same basemap, producing a single image of the entire quadrat (see Fig. 3). While contiguous quadrats were orthorectified individually, they were done so against a basemap containing all quadrats in the group, meaning that the resulting images can be easily merged to create a single image of the group. The outlines of all visible coral colonies ( >~1 cm2), and other benthic organisms such as algae and clams, were then digitised in ArcGIS to create a single shapefile for each quadrat for each year. Each colony was represented as an individual feature within the shapefile, and was assigned a unique colony number and species based on the mud maps drawn in the field. Colony numbers were consistent across years, allowing individual colonies to be tracked over time. If a colony underwent fission, the original colony number was retained for each, with the addition of a unique identifier after a decimal point. For example, if colony 35 split in two, the resultant colonies were identified as 35.1 and 35.2. If 35.2 later split again, the resultant colonies were identified as 35.2.1 and 35.2.2. If the colony overlapped the edge of the quadrat, only the area within the quadrat was digitised, and a flag was applied to indicate that only part of the colony was included (edgestatus = 1 in the data). Upon completion of digitisation, ArcGIS was used to calculate the area and perimeter of all colonies. While multiple census were conducted in 1963, 1971 and 1983, only a single census in each year has been processed. There are currently no plans to undertake further digitisation or QA/QC of this data set.Fig. 3Example orthorectified and stitched (prior to 2001) images from the NCNE quadrat, showing the effects of a cyclone that removed all colonies in 1972, and slow recovery over subsequent decades.Full size image More

Pecl, G. T. et al. Biodiversity redistribution under climate change: impacts on ecosystems and human well-being. Science 355, (2017).IPBES. Global assessment report on biodiversity and ecosystem services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services. (IPBES secretariat, 2019). https://doi.org/10.5281/zenodo.3831673.Elton, C. S. The Ecology of Invasions by Animals and Plants. (University of Chicago Press, 1958).Lockwood, J. L., Hoopes, M. F. & Marchetti, M. P. Invasion Ecology. (Wiley-Blackwell, 2013).O’Dowd, D. J., Green, P. T. & Lake, P. S. Invasional “meltdown” on an oceanic island. Ecol. Lett. 6, 812–817 (2003).
Google Scholar 
Doherty, T. S., Glen, A. S., Nimmo, D. G., Ritchie, E. G. & Dickman, C. R. Invasive predators and global biodiversity loss. Proc. Natl Acad. Sci. 113, 11261–11265 (2016).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Spatz, D. R. et al. Globally threatened vertebrates on islands with invasive species. Sci. Adv. 3, (2017).Pimentel, D. et al. Economic and environmental threats of alien plant, animal, and microbe invasions. Agriculture, Ecosyst. Environ. 84, 1–20 (2001).
Google Scholar 
Hoffmann, B. D. & Broadhurst, L. M. The economic cost of managing invasive species in Australia. NeoBiota 31, 1–18 (2016).
Google Scholar 
Kolar, C. S. & Lodge, D. M. Progress in invasion biology: predicting invaders. Trends Ecol. Evolution 16, 199–204 (2001).
Google Scholar 
Jeschke, J. M. & Strayer, D. L. Invasion success of vertebrates in Europe and North America. Proc. Natl Acad. Sci. 102, 7198–7202 (2005).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Lovell, R. S. L., Blackburn, T. M., Dyer, E. E. & Pigot, A. L. Environmental resistance predicts the spread of alien species. Nat. Ecol. Evolution 5, 322–329 (2021).
Google Scholar 
Blackburn, T. M. et al. A proposed unified framework for biological invasions. Trends Ecol. Evolution 26, 333–339 (2011).
Google Scholar 
Chapple, D. G., Simmonds, S. M. & Wong, B. B. M. Can behavioral and personality traits influence the success of unintentional species introductions? Trends Ecol. Evolution 27, 57–64 (2012).
Google Scholar 
Chapple, D. G. & Wong, B. B. M. The role of behavioural variation across different stages of the introduction process. in Biological Invasions and Animal Behaviour (eds. Weis, Judith, S. & Sol, Daniel.) 7–25 (Cambridge University Press, 2016).Holway, D. & Suarez, A. Animal behavior: an essential component of invasion biology. Trends Ecol. Evolution 14, 328–330 (1999).CAS 

Google Scholar 
Felden, A. et al. Behavioural variation and plasticity along an invasive ant introduction pathway. J. Anim. Ecol. 87, 1653–1666 (2018).PubMed 

Google Scholar 
D’Amore, D. M., Popescu, V. D. & Morris, M. R. The influence of the invasive process on behaviours in an intentionally introduced hybrid, Xiphophorus helleri-maculatus. Anim. Behav. 156, 79–85 (2019).
Google Scholar 
Perkins, T. A., Boettiger, C. & Phillips, B. L. After the games are over: life‐history trade‐offs drive dispersal attenuation following range expansion. Ecol. Evolution 6, 6425–6434 (2016).
Google Scholar 
Phillips, B. L., Brown, G. P., Travis, J. M. J. & Shine, R. Reid’s Paradox revisited: the evolution of dispersal kernels during range expansion. Am. Naturalist 172, S34–S48 (2008).
Google Scholar 
Shine, R., Brown, G. P. & Phillips, B. L. An evolutionary process that assembles phenotypes through space rather than through time. Proc. Natl Acad. Sci. 108, 5708–5711 (2011).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Lindström, T., Brown, G. P., Sisson, S. A., Phillips, B. L. & Shine, R. Rapid shifts in dispersal behavior on an expanding range edge. Proc. Natl Acad. Sci. 110, 13452–13456 (2013).ADS 
PubMed 
PubMed Central 

Google Scholar 
Heger, T. & Jeschke, J. M. The enemy release hypothesis as a hierarchy of hypotheses. Oikos 123, 741–750 (2014).
Google Scholar 
Colautti, R. I., Ricciardi, A., Grigorovich, I. A. & MacIsaac, H. J. Is invasion success explained by the enemy release hypothesis? Ecol. Lett. 7, 721–733 (2004).
Google Scholar 
Wilson, J. R. U., Dormontt, E. E., Prentis, P. J., Lowe, A. J. & Richardson, D. M. Something in the way you move: dispersal pathways affect invasion success. Trends Ecol. Evolution 24, 136–144 (2009).
Google Scholar 
Wilson, S. & Swan, G. A complete guide to reptiles of Australia. (New Holland Publishers, 2021).Chapple, D. G., Miller, K. A., Kraus, F. & Thompson, M. B. Divergent introduction histories among invasive populations of the delicate skink (Lampropholis delicata): has the importance of genetic admixture in the success of biological invasions been overemphasized? Diversity Distrib. 19, 134–146 (2013).
Google Scholar 
Chapple, D., Knegtmans, J., Kikillus, H. & van Winkel, D. Biosecurity of exotic reptiles and amphibians in New Zealand: building upon Tony Whitaker’s legacy. J. R. Soc. N.Z. 46, 66–84 (2016).
Google Scholar 
Chapple, D. G., Whitaker, A. H., Chapple, S. N. J., Miller, K. A. & Thompson, M. B. Biosecurity interceptions of an invasive lizard: Origin of stowaways and human-assisted spread within New Zealand. Evolut. Appl. 6, 324–339 (2013).
Google Scholar 
Tingley, R., Thompson, M. B., Hartley, S. & Chapple, D. G. Patterns of niche filling and expansion across the invaded ranges of an Australian lizard. Ecography 39, 270–280 (2016).
Google Scholar 
Chapple, D. G. et al. Biology of the invasive delicate skink (Lampropholis delicata) on Lord Howe Island. Aust. J. Zool. 62, 498–506 (2014).
Google Scholar 
Moule, H. et al. A matter of time: temporal variation in the introduction history and population genetic structuring of an invasive lizard. Curr. Zool. 61, 456–464 (2015).CAS 

Google Scholar 
Chapple, D. G., Simmonds, S. M. & Wong, B. B. M. Know when to run, know when to hide: can behavioral differences explain the divergent invasion success of two sympatric lizards? Ecol. Evolution 1, 278–289 (2011).
Google Scholar 
Cromie, G. L. & Chapple, D. G. Impact of tail loss on the behaviour and locomotor performance of two sympatric Lampropholis skink species. PLoS ONE 7, e34732 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Brand, J. A. et al. Rapid shifts in behavioural traits during a recent fish invasion. Behav. Ecol. Sociobiol. 75, 134 (2021).
Google Scholar 
Myles-Gonzalez, E., Burness, G., Yavno, S., Rooke, A. & Fox, M. G. To boldly go where no goby has gone before: boldness, dispersal tendency, and metabolism at the invasion front. Behav. Ecol. 26, 1083–1090 (2015).
Google Scholar 
Pintor, L. M., Sih, A. & Bauer, M. L. Differences in aggression, activity and boldness between native and introduced populations of an invasive crayfish. Oikos 117, 1629–1636 (2008).
Google Scholar 
Mueller, J. C. et al. Selection on a behaviour-related gene during the first stages of the biological invasion pathway. Mol. Ecol. 26, 6110–6121 (2017).MathSciNet 
CAS 
PubMed 

Google Scholar 
Snell-Rood, E. C. An overview of the evolutionary causes and consequences of behavioural plasticity. Anim. Behav. 85, 1004–1011 (2013).
Google Scholar 
Niemelä, P. T., Niehoff, P. P., Gasparini, C., Dingemanse, N. J. & Tuni, C. Crickets become behaviourally more stable when raised under higher temperatures. Behav. Ecol. Sociobiol. 73, 81 (2019).
Google Scholar 
Polverino, G. et al. Psychoactive pollution suppresses individual differences in fish behaviour. Proc. R. Soc. B: Biol. Sci. 288, 20202294 (2021).
Google Scholar 
Royauté, R., Garrison, C., Dalos, J., Berdal, M. A. & Dochtermann, N. A. Current energy state interacts with the developmental environment to influence behavioural plasticity. Anim. Behav. 148, 39–51 (2019).
Google Scholar 
Michelangeli, M., Chapple, D. G., Goulet, C. T., Bertram, M. G. & Wong, B. B. M. Behavioral syndromes vary among geographically distinct populations in a reptile. Behav. Ecol. 30, 393–401 (2019).
Google Scholar 
Nicolaus, M., Tinbergen, J. M., Ubels, R., Both, C. & Dingemanse, N. J. Density fluctuations represent a key process maintaining personality variation in a wild passerine bird. Ecol. Lett. 19, 478–486 (2016).PubMed 

Google Scholar 
Lapiedra, O., Schoener, T. W., Leal, M., Losos, J. B. & Kolbe, J. J. Predator-driven natural selection on risk-taking behavior in anole lizards. Science 360, 1017–1020 (2018).ADS 
CAS 
PubMed 

Google Scholar 
Gruber, J., Brown, G., Whiting, M. J. & Shine, R. Geographic divergence in dispersal-related behaviour in cane toads from range-front versus range-core populations in Australia. Behav. Ecol. Sociobiol. 71, 38 (2017).
Google Scholar 
Gruber, J., Brown, G., Whiting, M. J. & Shine, R. Is the behavioural divergence between range-core and range-edge populations of cane toads (Rhinella marina) due to evolutionary change or developmental plasticity? R. Soc. Open Sci. 4, 170789 (2017).ADS 
PubMed 
PubMed Central 

Google Scholar 
Morgan, D., Waas, J. R. & Innes, J. Do territorial and non-breeding Australian Magpies Gymnorhina tibicen influence the local movements of rural birds in New Zealand? Ibis 148, 330–342 (2006).
Google Scholar 
O’leary, R. A. & Jones, D. N. Foraging by suburban Australian magpies during dry conditions. Corella 26, 53–54 (2002).
Google Scholar 
Wright, T. F., Eberhard, J. R., Hobson, E. A., Avery, M. L. & Russello, M. A. Behavioral flexibility and species invasions: the adaptive flexibility hypothesis. Ethol. Ecol. Evolution 22, 393–404 (2010).
Google Scholar 
Dingemanse, N. J. & Wolf, M. Between-individual differences in behavioural plasticity within populations: causes and consequences. Anim. Behav. 85, 1031–1039 (2013).
Google Scholar 
Ducatez, S., Sol, D., Sayol, F. & Lefebvre, L. Behavioural plasticity is associated with reduced extinction risk in birds. Nat. Ecol. Evolution 4, 788–793 (2020).
Google Scholar 
Cole, E. F. & Quinn, J. L. Personality and problem-solving performance explain competitive ability in the wild. Proc. R. Soc. B: Biol. Sci. 279, 1168–1175 (2012).
Google Scholar 
Webster, M. M., Ward, A. J. W. & Hart, P. J. B. Individual boldness affects interspecific interactions in sticklebacks. Behav. Ecol. Sociobiol. 63, 511–520 (2009).
Google Scholar 
McGhee, K. E., Pintor, L. M. & Bell, A. M. Reciprocal behavioral plasticity and behavioral types during predator-prey interactions. Am. Naturalist 182, 704–717 (2013).
Google Scholar 
Ioannou, C. C., Payne, M. & Krause, J. Ecological consequences of the bold–shy continuum: the effect of predator boldness on prey risk. Oecologia 157, 177–182 (2008).ADS 
CAS 
PubMed 

Google Scholar 
Moran, N. P., Wong, B. B. M. & Thompson, R. M. Weaving animal temperament into food webs: implications for biodiversity. Oikos 126, 917–930 (2017).
Google Scholar 
Bellard, C., Cassey, P. & Blackburn, T. M. Alien species as a driver of recent extinctions. Biol. Lett. 12, 20150623 (2016).PubMed 
PubMed Central 

Google Scholar 
Moule, H., Michelangeli, M., Thompson, M. B. & Chapple, D. G. The influence of urbanization on the behaviour of an Australian lizard and the presence of an activity–exploratory behavioural syndrome. J. Zool. 298, 103–111 (2016).
Google Scholar 
Michelangeli, M., Wong, B. B. M. & Chapple, D. G. It’s a trap: sampling bias due to animal personality is not always inevitable. Behav. Ecol. 27, 62–67 (2016).
Google Scholar 
Michelangeli, M., Melki-Wegner, B., Laskowski, K., Wong, B. B. M. & Chapple, D. G. Impacts of caudal autotomy on personality. Anim. Behav. 162, 67–78 (2020).
Google Scholar 
Shine, R. Locomotor speeds of gravid lizards: Placing “costs of reproduction” within an ecological context. Funct. Ecol. 17, 526–533 (2003).
Google Scholar 
Naimo, A. C., Jones, C., Chapple, D. G. & Wong, B. B. M. Has an invasive lizard lost its antipredator behaviours following 40 generations of isolation from snake predators? Behav. Ecol. Sociobiol. 75, 131 (2021).
Google Scholar 
Brand, J. A. et al. Population differences in the effect of context on personality in an invasive lizard. Behav. Ecol. 32, 1363–1371 (2021).
Google Scholar 
Goulet, C. T., Thompson, M. B., Michelangeli, M., Wong, B. B. M. & Chapple, D. G. Thermal physiology: a new dimension of the pace‐of‐life syndrome. J. Anim. Ecol. 86, 1269–1280 (2017).PubMed 

Google Scholar 
Michelangeli, M., Goulet, C. T., Kang, H. S., Wong, B. B. M. & Chapple, D. G. Integrating thermal physiology within a syndrome: locomotion, personality and habitat selection in an ectotherm. Funct. Ecol. 32, 970–981 (2018).
Google Scholar 
Bell, A. M. Randomized or fixed order for studies of behavioral syndromes? Behav. Ecol. 24, 16–20 (2013).PubMed 

Google Scholar 
Friard, O. & Gamba, M. BORIS: a free, versatile open-source event-logging software for video/audio coding and live observations. Methods Ecol. Evolution 7, 1325–1330 (2016).
Google Scholar 
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.r-project.org/. (2019).Bürkner, P. C. brms: an R package for Bayesian multilevel models using Stan. J. Stat. Softw. 80, 1–28 (2017).
Google Scholar 
Munson, A. A., Michelangeli, M. & Sih, A. Stable social groups foster conformity and among-group differences. Anim. Behav. 174, 197–206 (2021).
Google Scholar 
Royauté, R. & Dochtermann, N. A. Comparing ecological and evolutionary variability within datasets. Behav. Ecol. Sociobiol. 75, 127 (2021).
Google Scholar 
Dalos, J., Royauté, R., Hedrick, A. V. & Dochtermann, N. A. Phylogenetic conservation of behavioural variation and behavioural syndromes. J. Evolut. Biol. 35, 311–321 (2022).
Google Scholar 
Miller, K. A., Duran, A., Melville, J., Thompson, M. B. & Chapple, D. G. Sex-specific shifts in morphology and colour pattern polymorphism during range expansion of an invasive lizard. J. Biogeogr. 44, 2778–2788 (2017).
Google Scholar 
Michelangeli, M., Chapple, D. G. & Wong, B. B. M. Are behavioural syndromes sex specific? Personality in a widespread lizard species. Behav. Ecol. Sociobiol. 70, 1911–1919 (2016).
Google Scholar 
Vehtari, A., Gelman, A. & Gabry, J. Practical Bayesian model evaluation using leave-one-out cross-validation and WAIC. Stat. Comput. 27, 1413–1432 (2017).MathSciNet 
MATH 

Google Scholar 
Nakagawa, S. & Schielzeth, H. Repeatability for Gaussian and non-Gaussian data: a practical guide for biologists. Biol. Rev. 85, 935–956 (2010).PubMed 

Google Scholar 
Chapple, D. G. et al. Data from Chapple et al. “Biological invasions as a selective filter driving behavioral divergence”. Monash University. Dataset. https://doi.org/10.26180/18851036.v2 (2022). More

DefinitionWe present here a new statistical approach to measure and characterize light pollution. The objective is to define an indicator which is not limited to clear sky measurements and does not require a precise calibration of a photometer. The key attributes of the indicator are the following:

It requires the automated acquisition of a large number of zenithal NSB measures when the Sun is below (-18^{circ }) and the Moon below (-5^{circ });

The acquisitions must at least cover a period of 6 months in order to record a wide range of possible weather conditions from perfectly clear to totally overcast skies. The objective is to obtain a significant sample of every type of cloud conditions (e.g. cloud density and ceiling altitude) as well as a good characterization of the average clear sky ;

It is based on the analysis of the zenithal NSB measure dispersion which is directly linked to the level of light pollution a site experiences.

As presented above in Fig. 2, the NSB density histograms, which are assembled from a large number of NSB measures, display a higher density zone which denotes a characteristic clear sky level that we name nominal NSB in the scope of the indicator calculation. On both sides of the clear sky level (above and below), NSB measures are distributed in a way that reflect the zenithal night sky luminance in cloudy conditions: NSB measures above the clear sky level mean that the light pollution is amplified by clouds while those below the clear sky level indicate a darker environment where clouds mask light pollution from distant sources as well as natural light sources. The calculation of the indicator is based on the evaluation of the NSB measure dispersion on both sides of the nominal NSB (i.e. characteristic clear sky level). Since there can have strong variations of artificial light emitted into the environment at the beginning and end of the night (decrease then increase of human activity, extinction of public lighting, etc.), the range of NSB measures retained for calculating the indicator is restricted to a portion in the middle of the night, typically 2 h.Figure 8 shows a typical NSB density histogram for a site which is quite severely impacted by light pollution. It covers a 2 h time range between 23:00 UTC and 01:00 UTC and one can easily see that the zone above the nominal NSB is much higher and denser than the one below, i.e. cloud conditions create more often a brighter environment than a darker one and with a greater amplitude.Figure 8NSB density histogram where the nominal NSB that represents the most common clear sky conditions is identified. It delimits two areas, the NSB bright dispersion above the nominal NSB and the NSB dark dispersion below.Full size imageBased on the determination of the nominal NSB, a quantitative indicator, called NDR for NSB Dispersion Ratio, is calculated in the following way:$$begin{aligned} NDR = (N_b cdot MAD_b) / (N_d cdot MAD_d) end{aligned}$$where (N_b) is the number of measures above the nominal NSB (brighter sky), (N_d) is the number of measures below the nominal NSB (darker sky), (MAD_b) is the median absolute deviation of the measures in the bright dispersion zone (above the nominal NSB) and (MAD_d) is the median absolute deviation of the measures in the dark dispersion zone (below the nominal NSB). The median absolute deviation is a statistical tool used to measure the variability of a data set, which is exactly what we try to achieve with the two NSB extensions above and below the nominal NSB. It is formally defined as (MAD = median(|X_i – tilde{{mathbf {X}}}|)) where (X_i) in our case represents an NSB value and (tilde{{mathbf {X}}}) is (median(X_i)). The median absolute deviation is a better choice than the usual standard deviation to measure the spread of NSB measures since the data does not follow a normal distribution.In order to make the determination of the NSB Dispersion Ratio stronger from a statistical standpoint, we use a bootstrapping with replacement resampling method on the set of night portions used to compute the indicator. Assuming we have N night portions at our disposal, we randomly select a sample of N items in this set of night portions knowing that a given item can appear multiple times in the sample (hence the bootstrapping with replacement). The NDR value is then computed for the considered sample. This process is repeated 1000 times and the average NDR value if eventually computed. This average value represents the actual NDR indicator of the considered site.The NDR indicator takes into account both the number of NSB values on each side of the nominal NSB and the dispersion of these values. This is what makes it relevant as an indicator of light pollution which encompasses all kinds of meteorological conditions experienced at a particular site. On that aspect, it is therefore not an astronomical light pollution indicator since it is not focused on clear sky conditions. On the opposite, it requires to have a significant number of NSB measures in all sorts of cloudy conditions so that a valid NDR indicator can be derived.A key aspect of the NDR calculation methodology is to determine the level of the nominal NSB, i.e. the typical clear sky level, since it will be used to differentiate the NSB measures that go in each of the two sets to calculate the bright and dark dispersions. As we have seen earlier in the article, such a determination can be biased by natural light sources that raise or lower clear sky NSB at different times of the night. This can result into a “blurry” high density zone which makes the determination of the nominal NSB difficult or even impossible depending on the observation period. Based on the quantitative estimate of the different natural light sources presented above, the most important bias to address is the contribution of the Galactic plane. This contribution must be eliminated for all the NSB measures which are used to calculate the NDR indicator. In order to do that, Noxi, the Ninox processing software developed by DarkSkyLab, calculates for each NSB measure the corresponding Galactic plane and star fields contribution using the galactic coordinates of the zenith and integrating the combined flux of all stars in the field of view using the UCAC4 astrometry and photometry star catalogue. It is not possible to cancel the contribution of the airglow due to its unpredictable nature, but since it only appears in rare occasions, it is not seen as a problem and is ignored. Regarding the contribution of the zodiacal light, it is considered as minimal at the zenith and it is also ignored.As an example, Fig. 9 shows on the left an NSB density histogram where the Galactic plane bias has not been corrected in the data, and on the right the same data but with the Galactic plane bias corrected. It is easy to see that in the latter the nominal NSB is much easier to determine, providing a more accurate reference level to calculate the NDR indicator. Once the Galactic plane bias has been corrected, the nominal NSB is determined as the highest density zone of the NSB histogram. It must be noted that, as of today, all the NSB measures are corrected from the Galactic plane bias without regards to the presence of clouds or high levels of light pollution. This results into an additional source of inaccuracy that will be addressed in the future through the implementation of two heuristics within the Noxi software:

1.

A first heuristic will determine if a night portion is considered as having a clear sky or not so that the Galactic bias correction is applied only if the sky is clear. In order to do that, we have developed an indicator called the NSS (for Night Sky Stability). To determine the NSS for a full night of measures or just a night portion, we fit the NSB curve with a degree 10 polynomial and we then compute the difference between each NSB measure and it polynomial counterpart. As a result, we obtain a set of residuals. The variance of all the residuals defines the NSS for the considered NSB dataset. Below a given value, the sky is considered as clear knowing that the NSS indicator has been calibrated on several NSB data sets for which the corresponding weather conditions are known;

2.

A second heuristic will allow us to weight the Galactic bias correction to be applied to NSB measures according to their value. For non-polluted skies with high values of NSB, the full Galactic bias correction will be applied while below a certain NSB threshold (for instance 21 mag(_{mathrm{SQM}})/arcsec(^{2}) which corresponds roughly to the brightest parts of the Milky Way) no correction will be applied.

Figure 9NSB density histograms of the same data set with no correction of the Galactic plane bias applied on the left and a full correction applied on the right.Full size imageThe NDR indicator is unitless since it is the ratio of two quantities with the same unit (mag(_{mathrm{SQM}})/arcsec(^{2})). For the data set presented in Fig. 9, the NDR value which is obtained is 25 (which is justified by the fact that the bright extension in the density histogram is much higher and denser than the dark extension). This denotes a quite high level of light pollution despite the fact that the nominal NSB is at a level of 21.6 mag(_{mathrm{SQM}})/arcsec(^{2}). This highlights the fact that there is not always a strict correlation between the typical clear sky NSB obtained for a given site and its NDR indicator, i.e. the presence of clouds decreases the NSB more than we could have expected just by knowing the clear sky NSB. On that respect, the NDR ratio brings more information that the clear sky NSB alone.In addition to provide an indicator which is representative of light pollution in all possible atmospheric conditions, the NDR provides a tool to compare locations in a more meaningful way than just using a set of standalone NSB evaluations. First it is not dependent of an inter-calibration between different systems and second its statistical nature makes it more robust when it comes to perform comparisons.NDR into practiceThe NDR indicator has been calculated for several different sites by DarkSkyLab during various projects in France that involved NSB measuring sessions in the field. To demonstrate some of the results that have been obtained, Fig. 10 provides the density histograms of 4 different sites which have quite different light pollution profiles.Figure 10NSB density histograms of 4 different sites used to compute the NDR indicator. The nominal NSB (which corresponds to the most common clear sky conditions) is noted with a white tick mark next to the vertical axis. Relative levels of the bright and dark dispersion terms (((N_b cdot MAD_b)) and ((N_d cdot MAD_d))) are noted respectively with an orange tick mark and a green tick mark. The computed values of the NDR indicator and nominal NSB are provided in the top-left corner of each figure.Full size imageTo build these diagrams, only the measures acquired during a few hours in the middle of the nights have been used to ensure the maximum stability of the NSB curves and avoid lighting extinctions that create large gaps in NSB profiles. The Galactic plane bias is corrected on all plots and the same NSB scale is used in order to perform comparisons between the 4 sites. One can notice that the number of measures and nights for the 4 sites are quite different. However, they are all sufficient to derive a meaningful value of the NDR indicator using the bootstrapping with replacement resampling method described above, but it is clear that the more NSB measures used, the more accurate the NDR indicator.Table 2 summarizes the NDR indicators as well as the nominal NSB for the 4 sites which are sorted in the order of decreasing NDR indicator values.Table 2 Summary of the nominal NSB and NDR indicators of the 4 different sites.Full size tableOne can see that the NDR indicator values are not strictly correlated to the nominal NSB values, e.g. despite the fact that the nominal NSB of site (a) is slightly better than the one of site (b), the NDR indicator value is much larger for site (a) than for site (b). This can be explained if we consider the specificities of each site:

Cervières (a) is a small village in the Haut-Forez area, France, which is surrounded by large cities (Lyon, Saint-Etienne and Clermont-Ferrand at a distance between 50 to 80 km) and a closer mid-size city (Roanne at 30 km). At the top of that, the town of Noirétable and a large highway rest area are just 2 km away without any nocturnal extinction applied (as opposed to the village of Cervières itself for which public lighting is turned off from 23:00 to 05:00 local time). These conditions are favourable to the presence of a constant light pollution background which has a negative impact on the zenithal NSB measures in most cloudy conditions (distant large cities for high elevation clouds and Noirétable and the highway rest area for lower elevation clouds). Only rare cloud conditions actually protect the site from the effect of mid-distance light sources. In clear sky conditions, however, the fact that there is no close light sources provides reasonably good NSB levels;

The Copernic Association Observatory (b) is located 6 km from the large town of Gap in the mountain area of Hautes-Alpes in the south of France. There is no significant short distant light sources but in many cloud conditions the contribution of Gap has a very negative impact on the zenithal luminance. However, due to the fact that the observatory is at a higher altitude on the hills surrounding the city of Gap, there are cloud conditions that make the site darker. In clear sky conditions, the proximity of Gap does not permit a quality better than that of a rural sky;

The Astrièves Observatory (c) is located near the center of the small town of Gresse-en-Vercors in the Parc Naturel Régional du Vercors. There is a full nocturnal extinction of the village for a large part of the night resulting in a good sky quality in clear sky conditions. The large city of Grenoble is at a distance of 30 km in a valley at the north-east, and the two locations are separated by a few mountains which efficiently help masking the light pollution as soon as the cloud ceiling is below a certain altitude, resulting into a dark environment. On the opposite, high elevation clouds reflect the light from Grenoble and increase the zenithal luminance;

Eourres (d) is a small and isolated village located 20 km west of Sisteron in the department of Hautes-Alpes, France, which is surrounded by mountains. There is no significant light sources closer than those of Sisteron and this results into a very good night sky quality with, most of the times, a very dark environment in cloudy conditions.

Figure 11 provides a graphical representation of the NDR indicator values for the 4 sites. On the NDR scale, the value 1 indicates that the bright and dark dispersion terms (respectively ((N_b cdot MAD_b)) and ((N_d cdot MAD_d))) are equal, which means there is a balance between dark and bright conditions at the zenith on the considered site with reference to the most common clear sky level.Figure 11Summary of the NDR indicators obtained for the 4 sites. The diagram uses 1 as the pivotal value to delineate sites according to the two bright and dark dispersion terms ((N_b cdot MAD_b)) and ((N_d cdot MAD_d)).Full size imageThe NDR can theoretically vary between 0 (totally dark site) and several hundreds (extremely bright site) but in practice the best sites can reach NDR indicator values down to 0.3 in the best preserved locations and up to 200 for very large and polluted cities.Robustness of the NDR indicatorIt is important to evaluate how the NDR indicator is dependant on the number of measures used to compute it and to figure out what would be the minimum number of night sessions required to obtain a meaningful NDR indicator value at a given site. To achieve that, we have used the data from two of the four sites presented above (the two which have the largest number or recorded nights: Cervières with 424 nights and the Astrièves Observatory with 373 nights). The 1000-step bootstrapping procedure has been repeatedly executed on each data set with a regularly decreasing sample of nights: starting from the full number of nights, a decrement of 10 nights is applied at each step until only 20 nights are remaining. At every bootstrap step, each sample is composed of n nights randomly chosen among the N available ones knowing that any night can be selected several times.Figure 12 shows the NDR indicator values that have been obtained for each of the two sites as a function of the night sample considered. The 95% confidence interval is plotted against each NDR indicator value (it is preferred to the standard deviation since the NSB distribution in the data sets is not normal). In the right plot of Fig. 12, the last confidence interval for the 24 night sample is too wide to fit in the y-axis NDR range (the top value is 195).Figure 12Results of the NDR resampling on the two data sets of Cervières and Astrièves Observatory. The horizontal axis is the number of nights considered into the night sample and the vertical axis provides the NDR indicator obtained for each sampling set through a 1000-iteration bootstrapping with replacement calculation.Full size imageDiscussion on the required number of nightsWe can see in Fig. 12 that the NDR indicator and the confidence interval remain stable down to 200 nights. Below this threshold, the NDR starts to become unstable with growing confidence intervals. Based on this data, we can estimate that the minimum number of nights required to compute a robust NDR indicator is 200 (therefore between 7 and 8 months since there are periods around the full moons where there is no night portions recorded).However, depending on the measuring session objectives, the NDR indicator can be considered as accurate enough even when using a smaller number of nights. If the goal is simply to get a first estimate of the light pollution level at a given site, we can consider that 90 nights (a little more than 3 months of measures) are enough. On the opposite, if we want to perform a comparison between several sites for evaluating the impact of light pollution on a particular species, we might want to perform at least 200 nights of measurement to get a better accuracy for the NDR indicator. The experience from DarkSkyLab through many NSB measuring sessions is that 3 to 4 months of measures are required to get a meaningful density histogram, hence an accurate enough NDR indicator, so that a site can be sufficiently characterized from a light pollution perspective. Such a measuring period usually guarantees that the clear sky nominal NSB is well defined and that various cloud conditions have been observed. This estimate is sustained by the results obtained in Fig. 12.Value of the NDR indicator for ecological researchThe study of the impact of light pollution on biodiversity is currently in full expansion, amplifying a political and citizen demand for the reclamation of the night2,32,33.We identify three main contributions of the NDR indicator for ecological research. First, it overcomes the limits of an old problem of communication in terms of measurement units between disciplines and potentially limits the use of units without real meaning from a biodiversity point of view34,35. Secondly, the use of the NDR indicator limits the common biases linked to a characterization of the effects of anthropogenic light which is too limited in time and space35. Indeed, the life history traits of species are not only shaped by the intensity of light emitted into the nocturnal environment but also by its variation over time34,35,36. Currently, the characterization of light pollution is too often limited in time and space, which can lead to misinterpretation37. Thirdly, the NDR indicator provides ecological researchers with a unit of measurement that integrates a sufficiently long time step to study the impact of light pollution on the evolutionary processes at work in the life of species and particularly on population dynamics and animal behavior36,38,39.Limitations and future improvements of the NDR indicatorThe main limitation of the NDR indicator resides in the possible difficulty to identify a well defined value for the nominal NSB, i.e. the NSB value that represents the most common clear sky conditions of a given site. For the most part, this is due to the contribution of the Galactic plane to the zenithal sky brightness and, to a lower extent, to the contribution of other natural light sources (dense star fields, airglow and zodiacal light). The residual spread of NSB measures is due to changing atmospheric conditions at various time scales, but, for this particular contribution, we can expect a statistical compensation to eliminate a systematic associated bias.At the moment, the contribution of the Galactic plane and star fields is canceled into the NSB measures by calculating in the Noxi software the integrated flux of all the stars that belong to the field of view (using the UCAC4 star catalogue). However, this approach has proven some limitations, especially in the southern hemisphere where the Galactic center goes through the zenith and is particularly bright. A probable explanation for that lack of predictability is the fact that the Galactic plane contains diffuse sources such as nebulae which are not accounted for into the star catalogues and which actually cannot be ignored. To address this issue, DarkSkyLab has the project to create a brightness map of the Galactic plane with a square degree resolution or better so that the contribution of all sources can be correctly accounted for.In addition to better correcting the Galactic plane bias, an improvement must be made with regards to the NSB measures that need to be corrected. At the moment, all NSB measures are corrected from the Galactic plane bias without regards to the presence of clouds or high levels of light pollution. So a first heuristic must be implemented to only apply the bias correction to clear sky NSB measures. An other heuristic must also be developed to reduce the correction applied as a function of the NSB level.A third limitation of the NDR indicator is related to a possible lack of cloudy conditions at some sites (e.g. in the Atacama desert in Chile with more than 320 clear nights per year), the reason simply being that the NDR indicator requires the presence of clouds to differentiate the bright and dark extensions into the NSB density histograms. This means that the NDR indicator can hardly be used for such astronomy-oriented sites which experience rare cloudy conditions. More

Garibaldi, L. A. et al. Stability of pollination services decreases with isolation from natural areas despite honey bee visits. Ecol. Lett. 14(10), 1062–1072 (2011).PubMed 
Article 

Google Scholar 
Kremen, C. et al. Pollination and other ecosystem services produced by mobile organisms: A conceptual framework for the effects of land-use change. Ecol. Lett. 10(4), 299–314 (2007).PubMed 
Article 

Google Scholar 
Kluser, S. & Peduzzi, P. Global pollinator decline: A literature review. Preprint at http://archive-ouverte.unige.ch/unige 32258 (2007).Potts, S. G. et al. Global pollinator declines: Trends, impacts and drivers. Trends Ecol. Evol. 25(6), 345–353 (2010).PubMed 
Article 

Google Scholar 
Rhodes, C. J. Pollinator decline—an ecological calamity in the making?. Sci. Prog. 101(2), 121–160 (2018).PubMed 
Article 

Google Scholar 
Huang, H. & D’Odorico, P. Critical transitions in plant-pollinator systems induced by positive inbreeding-reward-pollinator feedbacks. Iscience 23(2), 100819 (2020).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Krishnan, N. et al. Assessing field-scale risks of foliar insecticide applications to monarch butterfly (Danaus plexippus) larvae. Environ. Toxicol. Chem. 39(4), 923–941 (2020).CAS 
PubMed 
Article 

Google Scholar 
Bargar, T. A., Hladik, M. L. & Daniels, J. C. Uptake and toxicity of clothianidin to monarch butterflies from milkweed consumption. PeerJ 8, e8669 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Emmel, T. C. & Tucker, J. C. In Mosquito Control Pesticides: Ecological Impacts and Management Alternatives (eds Emmel, T. C. & Tucker, J. C.) 105 (Scientific Publishers, 1991).Johnson, R. M., Ellis, M. D., Mullin, C. A. & Frazier, M. Pesticides and honey bee toxicity–USA. Apidologie 41(3), 312–331 (2010).CAS 
Article 

Google Scholar 
Olaya-Arenas, P., Scharf, M. E. & Kaplan, I. Do pollinators prefer pesticide-free plants? An experimental test with monarchs and milkweeds. J. Appl. Ecol. 57(10), 2019–2030 (2020).CAS 
Article 

Google Scholar 
Berryman, A. A. What causes population cycles of forest Lepidoptera?. Trends Ecol. Evol. 11(1), 28–32 (1996).CAS 
PubMed 
Article 

Google Scholar 
Elkinton, J. & Boettner, G. Benefits and harm caused by the introduced generalist tachinid, Compsilura concinnata North America. Biol. Control 57(2), 277–288 (2012).
Google Scholar 
Beschta, R. L. & Ripple, W. J. Riparian vegetation recovery in Yellowstone: The first two decades after wolf reintroduction. Biol. Conserv. 198, 93–103 (2016).Article 

Google Scholar 
Oberhauser, K. et al. Lacewings wasps and fliesoh my insect enemies take a bite out of monarchs. In Monarchs in a Changing World: Biology and Conservation of an iconic insect (eds Oberhauser, K. S. et al.) 71–82 (Cornell University Press, 2015).Chapter 

Google Scholar 
Zalucki, M. P., Clarke, A. R. & Malcolm, S. B. Ecology and behavior of first instar larval Lepidoptera. Annu. Rev. Entomol. 47(1), 361–393 (2002).CAS 
PubMed 
Article 

Google Scholar 
Hermann, S. L., Blackledge, C., Haan, N. L., Myers, A. T. & Landis, D. A. Predators of monarch butterfly eggs and neonate larvae are more diverse than previously recognised. Sci. Rep. 9(1), 1–9 (2019).CAS 
Article 

Google Scholar 
McCoshum, S. M., Andreoli, S. L., Stenoien, C. M., Oberhauser, K. S. & Baum, K. A. Species distribution models for natural enemies of monarch butterfly (Danaus plexippus) larvae and pupae: Distribution patterns and implications for conservation. J. Insect Conserv. 20(2), 223–237 (2016).Article 

Google Scholar 
Geest, E. A., Wolfenbarger, L. L. & McCarty, J. P. Recruitment, survival and parasitism of monarch butterflies (Danaus plexippus) in milkweed gardens and conservation areas. J. Insect Conserv. 23(2), 211–224 (2019).Article 

Google Scholar 
Stenoien, C. et al. Monarchs in decline: A collateral landscape-level effect of modern agriculture. Insect Sci. 25(4), 528–541 (2018).PubMed 
Article 

Google Scholar 
Crone, E. E., Pelton, E. M., Brown, L. M., Thomas, C. C. & Schultz, C. B. Why are monarch butterflies declining in the west? Understanding the importance of multiple correlated drivers. Ecol. Appl. 29(7), e01975 (2019).PubMed 
Article 

Google Scholar 
Brower, L. P. et al. Effect of the 2010–2011 drought on the lipid content of monarchs migrating through Texas to overwintering sites in Mexico. In The Monarchs in a Changing World: Biology and Conservation of an Iconic Butterfly (eds Oberhauser, K. S. et al.) 117–129 (Cornell University Press, 2015).
Google Scholar 
Thogmartin, W. E. et al. Monarch butterfly population decline in North America: Identifying the threatening processes. R. Soc. Open Sci. 4(9), 170760 (2017).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Olaya-Arenas, P. & Kaplan, I. Quantifying pesticide exposure risk for monarch caterpillars on milkweeds bordering agricultural land. Front. Ecol. Evol. https://doi.org/10.3389/fevo.2019.00223 (2019).
Article 

Google Scholar 
Olaya-Arenas, P., Hauri, K., Scharf, M. E. & Kaplan, I. Larval pesticide exposure impacts monarch butterfly performance. Sci. Rep. 10(1), 1–12 (2020).Article 

Google Scholar 
Cameron, S. A. et al. Patterns of widespread decline in North American bumble bees. PNAS 108(2), 662–667 (2011).ADS 
MathSciNet 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Epstein, L. Fifty years since silent spring. Annu. Rev. Phytopathol. 52, 377–402 (2014).CAS 
PubMed 
Article 

Google Scholar 
Rayor, L. S. Effects of monarch larval host plant chemistry and body size on Polistes wasp predation. In The Monarch Butterfly Biology and Conservation (eds Oberhauser, K. S. & Solensky, M. J.) 39–46 (Cornell University Press, 2004).
Google Scholar 
Baker, A. M. & Potter, D. A. Invasive paper wasp turns urban pollinator gardens into ecological traps for monarch butterfly larvae. Sci. Rep. 10(1), 1–7 (2020).Article 

Google Scholar 
Castellanos, I. & Barbosa, P. Dropping from host plants in response to predators by a polyphagous caterpillar. J. Lepid. Soc. 65(4), 270–272 (2011).
Google Scholar 
Kessler, S. C. et al. Bees prefer foods containing neonicotinoid pesticides. Nature 521(7550), 74–76 (2015).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Liao, L.-H., Wu, W.-Y. & Berenbaum, M. R. Behavioral responses of honey bees (Apis mellifera) to natural and synthetic xenobiotics in food. Sci. Rep. 7(1), 1–8 (2017).Article 

Google Scholar 
Musser, R. O. et al. Caterpillar saliva beats plant defences. Nature 416(6881), 599–600 (2002).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Schmidt, J. & Smith, J. Host examination walk and oviposition site selection of Trichogramma minutum: Studies on spherical hosts. J. Insect Behav. 2(2), 143–171 (1989).Article 

Google Scholar 
Ramos, R. S. et al. Investigation of the lethal and behavioral effects of commercial insecticides on the parasitoid wasp Copidosoma truncatellum. Chemosphere 191, 770–778 (2018).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Chareonviriyaphap, T. et al. Pesticide avoidance behavior in Anopheles albimanus, a malaria vector in the Americas. J. Am. Mosq. Control Assoc. 13(2), 171–183 (1997).CAS 
PubMed 

Google Scholar 
Nansen, C., Baissac, O., Nansen, M., Powis, K. & Baker, G. Behavioral avoidance-will physiological insecticide resistance level of insect strains affect their oviposition and movement responses?. PLoS ONE 11(3), e0149994 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Martini, X., Kincy, N. & Nansen, C. Quantitative impact assessment of spray coverage and pest behavior on contact pesticide performance. Pest Manag. Sci. 68(11), 1471–1477 (2012).CAS 
PubMed 
Article 

Google Scholar 
Bull, D. & Coleman, R. Effects of pesticides on Trichogramma spp. Southwest. Entomol. Suppl. 8, 156–168 (1985).CAS 

Google Scholar 
Thubru, D., Firake, D. & Behere, G. Assessing risks of pesticides targeting lepidopteran pests in cruciferous ecosystems to eggs parasitoid, Trichogramma brassicae (Bezdenko). Saudi J. Biol. Sci. 25(4), 680–688 (2018).CAS 
PubMed 
Article 

Google Scholar 
Selwood, K. & Zimmer, H. Refuges for biodiversity conservation: A review of the evidence. Biol. Conserv. 245, 108502 (2020).Article 

Google Scholar 
Chmiel, J. A., Daisley, B. A., Pitek, A. P., Thompson, G. J. & Reid, G. Understanding the effects of sublethal pesticide exposure on honey bees: A role for probiotics as mediators of environmental stress. Front. Ecol. Evol. 8, 22 (2020).Article 

Google Scholar 
Chittka, L., Williams, N., Rasmussen, H. & Thomson, J. Navigation without vision: Bumblebee orientation in complete darkness. Proc. R. Soc. B 266(1414), 45–50 (1999).PubMed Central 
Article 

Google Scholar 
Young, M. W. & Kay, S. A. Time zones: A comparative genetics of circadian clocks. Nat. Rev. Genet. 2(9), 702–715 (2001).CAS 
PubMed 
Article 

Google Scholar 
Mallet, J. Gregarious roosting and home range in Heliconius butterflies. Natl. Geogr. Res. 2(2), 198–215 (1986).
Google Scholar 
Chang, Y.-M. et al. Roosting site usage, gregarious roosting and behavioral interactions during roost-assembly of two Lycaenidae butterflies. Zool. Stud. 59, e10 (2020).PubMed 
PubMed Central 

Google Scholar 
Vulinec, K. Collective security aggregation by insects as a defence. In Insect Defences. Adaptive Mechanisms of Prey and Predators (eds Evans, D. L. & Schmidt, J. O.) 251–288 (State University of New York, 1990).
Google Scholar 
Salcedo, C. Environmental elements involved in communal roosting in Heliconius butterflies (Lepidoptera: Nymphalidae). Environ. Entomol. 39(3), 907–911 (2010).PubMed 
Article 

Google Scholar 
Giordano, B. V., McGregor, B. L., Runkel, A. E. IV. & Burkett-Cadena, N. D. Distance diminishes the effect of deltamethrin exposure on the monarch butterfly, Danaus plexippus. J. Am. Mosq. Control Assoc. 36(3), 181–188 (2020).PubMed 
Article 

Google Scholar 
Matzrafi, M. Climate change exacerbates pest damage through reduced pesticide efficacy. Pest Manag. Sci. 75(1), 9–13 (2019).CAS 
PubMed 
Article 

Google Scholar 
Hewitt, A. Spray drift: Impact of requirements to protect the environment. Crop Prot. 19(8–10), 623–627 (2000).Article 

Google Scholar 
Nail, K. R., Stenoien, C. & Oberhauser, K. S. Immature monarch survival: Effects of site characteristics, density and time. Ann. Entomol. Soc. 108(5), 680–690 (2015).Article 

Google Scholar 
Payne, C. C. & Mertens, P. P. Cytoplasmic polyhedrosis viruses. In The Reoviridae (ed. Joklik, K.) 425–504 (Springer, 1983).Chapter 

Google Scholar 
Zalucki, M. P. et al. It’s the first bites that count: Survival of first-instar monarchs on milkweeds. Austral. Ecol. 26(5), 547–555 (2001).Article 

Google Scholar 
Salvato, M. Influence of mosquito control chemicals on butterflies (Nymphalidae, Lycaenidae, Hesperiidae) of the lower Florida keys. J. Lepid. Soc. 55(1), 8–14 (2001).
Google Scholar 
Frey, D. F. & Leong, K. L. Can microhabitat selection or differences in ‘catchability’ explain male-biased sex ratios in overwintering populations of monarch butterflies?. Anim. Behav. 45(5), 1025 (1993).Article 

Google Scholar 
Macgregor, C. J. & Scott-Brown, A. S. Nocturnal pollination: An overlooked ecosystem service vulnerable to environmental change. Emerg. Top. Life Sci. 4(1), 19–32 (2020).PubMed 
PubMed Central 
Article 

Google Scholar  More