More stories

  • in

    A global inventory of animal diversity measured in different grazing treatments

    Synthesis and data extractionData were collected using a literature search of Web of Science for peer-reviewed journal articles published between 1970 and November 2019. We conducted two sets of searches to capture grazing with discrete comparisons (e.g., grazed/ungrazed, moderate vs. heavy intensity grazing) and a range of grazing intensities. The search terms used for each were as follows 1) (graz* OR livestock) AND (exclosure* OR exclusion OR exclude* OR ungrazed OR retire* OR fallow* OR fence* OR paddock*), 2) (“grazing intensity” OR “grazing gradient” OR “stocking rate” OR “rotation*grazing”). Our synthesis includes domesticated and wild grazer species, with the latter defined as an undomesticated species naturally occurring in the study area during the study. Wild grazers are typically native species to the region (e.g., the American bison in Western North America) but can include non-native species that are naturalized in the area (e.g., feral horses on Sable Island).We excluded any study that did not test the effect of grazing animals. A grazer was defined using the definition provided by the authors of the respective study to account for the proportion of forage types in a herbivore’s diet that varies between seasons and habitats. For example, we included animals where their diet is assumed to come from all (e.g., cattle, sheep), most (e.g., wapiti, kangaroos), or some (e.g., deer species) grass species. However, within the included studies, these animals were classified as grazers as most of their diet was grass for the duration of the study. For added clarity about the herbivore composition in each study, we extracted a list of any herbivores listed in the paper regardless of foraging type or if any data was provided.We only included studies that measured animal diversity or abundance as a response variable and included data we could extract or contact the author to obtain9. We included any study with a grazing treatment and included observations within these studies of any grazed and ungrazed sites. All studies with grazing included a comparison to either ungrazed sites, different grazing practices (e.g., cattle vs. sheep), and/or differences in intensities (e.g., heavy/light, extensive/intensive). Studies that only measured plants or soil biota were excluded because syntheses of grazing effects on these groups have already been conducted7,11,12, and our goal was to provide a robust inventory of animal diversity. However, if a study included plants, lichens, or fungi in addition to animals, we included this data. Studies discussing marine grazing or aquatic systems were also excluded. From these preliminary filters, we identified 3,489 published manuscripts. We reviewed these 3,489 published articles and found 245 studies that surveyed animals in grazed sites. In total, we extracted 16,105 observations for over 1,200 species.We extracted 28 variables that focus on management systems, assemblages of grazer species, ecosystem characteristics, and survey type (Table 1). The latitudes, longitudes, and elevations of each study were included when provided for use with geospatial data. In addition, we included variables about the study site’s disturbance history, including last time grazed, if a flood event or fire had occurred, if fertilization was used, if the area was open or fenced off, and if the area was publicly or privately owned. Furthermore, the timeline for the study (i.e., the years the authors initiated and completed the study) was also provided. Study initiation was described by the authors and could include when the grazing treatment started, another treatment was applied, and/or animal surveys began. These timeline columns can be useful in identifying long-term studies and differentiating single grazing events or multi-year experiments. Finally, we generalized the characteristics of the ecosystem of the sites used in each study based on the climate and dominant vegetation.Table 1 The attributes and description of the metadata.csv file that lists the general characteristics of each study.Full size tableWithin the grazing data, we included information about the grazer when provided, including any measurement of the intensity of grazing (e.g., animals per hectare, the height of residual vegetation). We also provided two columns that detailed whether the study tested grazing effects using a discrete comparison or gradient of intensities (Table 2). The value for the target specimens extracted may represent either a single observation or a summarized statistic (e.g., mean animals per site). We identify unique observations as “count” and summarized statistics by the metric used, such as mean, median, standard deviation (column stat in grazingData.csv). When possible, we also included any record of other grazers that co-occurred with the observed grazer species. The data for these variables were extracted from the papers by a single researcher who read through each paper and filled in available data on the mentioned variables.Table 2 The attributes and description of the grazingData.csv file that has the extracted data from each study.Full size tableWe extracted information about the target specimen, site, year, experimental replicate, and response estimate (Table 2). We included multiple categorizations of the target species to assist future users in synthesizing similar taxa (Table 2). When a species name or genus was provided, we conducted a search query (see detailedTaxa.r) through the global biodiversity information facility (GBIF.org) to determine the taxonomic classification of the species, including kingdom, phylum, order, class, and family. When a species name was not included, we provided the lowest taxonomic resolution available. We also included a broader classification of ‘higherTaxon’ to distinguish plants, fungi, vertebrates, and invertebrates. These columns may help group similar species together for community-level analyses. Lastly, we included the characteristic of the plant community (i.e., planted or self-assembled, tilled, and its vegetation class) when plant data was reported.Patterns among studiesMost of the studies took place in the United States (26%), Australia (9%), and the United Kingdom (7%) (Fig. 1). As expected, most studies were conducted in grasslands (n = 206), followed by forests (n = 92) and shrublands (n = 82) (Fig. 2). We included publications from the entire range of years (i.e., 1970–2019), but most were published after 2000 (76%). The number of sites in a study and the study duration showed a bimodal distribution with a long tail (Fig. 3). Most studies included one to eight different sites, and few were conducted longer than five years (Fig. 3). A few studies were highly replicated, while many were limited in their replication (Fig. 3).Fig. 1The locations of studies that measured the response of animals to domestic or wild grazing.Full size imageFig. 2The number of grazing studies conducted in ecosystems around the world. We generalized the characteristics of the ecosystem of the sites used in each study based on the climate and dominant vegetation community. We separated grassland communities into those that were (a) semi-natural without recent cultivation or seeding (self-assembled), (b) recently cultivated or had supplemental seeding (planted/cultivated), and (c) a combination of both. In most grasslands, the cultivation history was unclear.Full size imageFig. 3The number of independent sites surveyed and the duration of each study. Most studies were conducted at either a single site or with some replication (e.g., 6–8 sites). Similarly, most studies were either conducted in one year ( >30%) or over a few years (e.g., 3–6 years). Very few studies (32) or lasted longer than 15 years.Full size imageSite and management data were not reported in all studies, as found in other reviews of grazing impacts on ecological processes10. Of the studies that mentioned the ownership status of the land used, 46% were on private land, 42% were on public land, and 12% had a history of both public and private ownership. Most studies included binary comparisons (56%) of grazed vs. ungrazed plots or sites, though some also included a discrete (22%) or a continuous estimate of grazing intensities (18%).Of the studies that reported plant community origin, 76% were self-assembled, 17% were planted communities, and the remaining included sites were a combination of the two. Domesticated grazers as the focal herbivore made up 67% of the studies, with 12% of the studies having wild grazers as the focal herbivore, and 21% having both present. Domesticated livestock were the most frequently surveyed grazers including cattle (n = 164), sheep (n = 83), and horses (n = 21), but studies are included that examined wild grazers, such as kangaroos (n = 6), elephants (n = 5), and pronghorn (n = 5) (Fig. 4).Fig. 4The frequency in which a study reported herbivores. We included any mention of herbivores regardless of being a grazer, browser, granivore, or other class. This list was obtained by the text within the manuscript and is different than the representation of species in the database (i.e., the measured species).Full size image More

  • in

    Exceptional longevity in northern peripheral populations of Wels catfish (Siluris glanis)

    Roff, D. A. The Evolution of Life Histories (Chapman & Hall, 1992).
    Google Scholar 
    Stearns, S. C. The Evolution of Life Histories (Oxford University Press, 1992).
    Google Scholar 
    Tibblin, P. et al. Evolutionary divergence of adult body size and juvenile growth in sympatric subpopulations of a top predator in aquatic ecosystems. Am. Nat. 186, 98–110 (2015).PubMed 

    Google Scholar 
    Voituron, Y., de Fraipont, M., Issartel, J., Guillaume, O. & Clobert, J. Extreme lifespan of the human fish (Proteus anguinus): A challenge for ageing mechanisms. Biol. Lett. 7, 105–107 (2011).PubMed 

    Google Scholar 
    Longhurst, A. Murphy’s law revisited: Longevity as a factor in recruitment to fish populations. Fish. Res. 56, 125–131 (2002).
    Google Scholar 
    Schaffer, W. M. Optimal reproductive effort in fluctuating environments. Am. Nat. 108, 783–790 (1974).
    Google Scholar 
    Beamish, R. J., McFarlane, G. A. & Benson, A. Longevity overfishing. Prog. Oceanogr. 68, 289–302 (2006).ADS 

    Google Scholar 
    Conti, B. Considerations on temperature, longevity and aging. Cell. Mol. Life Sci. 65, 1626–1630 (2008).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Inness, C. L. W. & Metcalfe, N. B. The impact of dietary restriction, intermittent feeding and compensatory growth on reproductive investment and lifespan in a short-lived fish. Proc. R. Soc. Lond. B Biol. Sci. 275, 1703–1708 (2008).
    Google Scholar 
    Liu, R. K. & Walford, R. L. Increased growth and life-span with lowered ambient temperature in the annual fish, Cynolebias adloffi. Nature 212, 1277–1278 (1966).ADS 

    Google Scholar 
    Trip, E. D., Clements, K. D., Raubenheimer, D. & Choat, J. H. Temperature-related variation in growth rate, size, maturation and life span in a marine herbivorous fish over a latitudinal gradient. J. Anim. Ecol. 83, 866–875 (2014).PubMed 

    Google Scholar 
    Munch, S. B. & Salinas, S. Latitudinal variation in lifespan within species is explained by the metabolic theory of ecology. Proc. Natl. Acad. Sci. U.S.A. 106, 13860–13864 (2009).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Britton, J. R., Pegg, J., Sedgwick, R. & Page, R. Investigating the catch returns and growth rate of wels catfish, Silurus glanis, using mark-recapture. Fish. Man. Ecol. 14, 263–268 (2007).
    Google Scholar 
    Hamel, M. J. et al. Range-wide age and growth characteristics of shovelnose sturgeon from mark–recapture data: Implications for conservation and management. Can. J. Fish. Aquat. Sci. 72, 71–82 (2015).
    Google Scholar 
    Hamel, M. J. et al. Using mark–recapture information to validate and assess age and growth of long-lived fish species. Can. J. Fish. Aquat. Sci. 71, 559–566 (2014).
    Google Scholar 
    Casale, P., Mazaris, A. D., Freggi, D., Vallini, C. & Argano, R. Growth rates and age at adult size of loggerhead sea turtles (Caretta caretta) in the Mediterranean Sea, estimated through capture-mark-recapture records. Sci. Mar. 73, 589–595 (2009).
    Google Scholar 
    IUCN (International Union for Conservation of Nature) 2008. Siluris glanis. The IUCN Red List of Threatened Species. Version 2021-3 (2010). https://www.iucnredlist.org. (Accessed 25 February 2021).Copp, G. H. et al. Voracious invader or benign feline? A review of the environmental biology of European catfish Silurus glanis in its native and introduced ranges. Fish. Fish. 10, 252–282 (2009).
    Google Scholar 
    Palm, S., Vinterstare, J., Nathanson, J. E., Triantafyllidis, A. & Petersson, E. Reduced genetic diversity and low effective size in peripheral northern European catfish Silurus glanis populations. J. Fish. Biol. 95, 1407–1421 (2019).PubMed 

    Google Scholar 
    Jensen, A., Lillie, M., Bergstrom, K., Larsson, P. & Hoglund, J. Whole genome sequencing reveals high differentiation, low levels of genetic diversity and short runs of homozygosity among Swedish wels catfish. Heredity 127, 79–91 (2021).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Cucherousset, J. et al. Ecology, behaviour and management of the European catfish. Rev. Fish. Biol. Fish. 28, 177–190 (2017).
    Google Scholar 
    Kuzishchin, K. V., Gruzdeva, M. A. & Pavlov, D. S. Traits of biology of European Wels Catfish Silurus glanis from the Volga-Ahtuba water system, the Lower Volga. J. Ichthyol. 58, 833–844 (2019).
    Google Scholar 
    Alp, A., Kara, C., Üçkardeş, F., Carol, J. & García-Berthou, E. Age and growth of the European catfish (Silurus glanis) in a Turkish Reservoir and comparison with introduced populations. Rev. Fish. Biol. Fish. 21, 283–294 (2010).
    Google Scholar 
    Carol, J., Benejam, L. B. & García-Berthou, E. Growth and diet of European catfish (Silurus glanis) in early and late invasion stages. Fund. Appl. Limnol. 174, 317–328 (2009).
    Google Scholar 
    Severov, Y. A. Size–age structure, growth rate, and fishery of European Catfish Silurus glanis in the lower Kama Reservoir. J. Ichthyol. 60, 118–121 (2020).
    Google Scholar 
    Lessmark, O. Malprovfiske i Möckeln 2006. Länsstyrelsens rapportserie (2006).Lessmark, O. Malprovfiske i Möckeln 2007. Länsstyrelsens rapportserie (2007).Harka, A. Studies on the growth of the sheatfish (Silurus glanis L.) in River Tisza. Aquac. Hung. (Szarvas) 4, 135–144 (1984).
    Google Scholar 
    Edwards, J. E. et al. Advancing research for the management of long-lived species: A case study on the Greenland shark. Front. Mar. Sci. https://doi.org/10.3389/fmars.2019.00087 (2019).Article 

    Google Scholar 
    Pikitch, E. K., Doukakis, P., Lauck, L., Chakrabarty, P. & Erickson, D. L. Status, trends and management of sturgeon and paddlefish fisheries. Fish. Fish. 6, 233–265 (2005).
    Google Scholar 
    Pironon, S. et al. Geographic variation in genetic and demographic performance: New insights from an old biogeographical paradigm. Biol. Rev. 92, 1877–1909 (2017).PubMed 

    Google Scholar 
    Antonovics, J., McKane, A. J. & Newman, T. J. Spatiotemporal dynamics in marginal populations. Am. Nat. 167, 16–27 (2006).CAS 
    PubMed 

    Google Scholar 
    Alp, A., Kara, C. & Büyükcapar, H. M. Reproductive biology in a Native European Catfish, Siluris glanis L., 1758, population in Menzelet Resevoir. Turk. J. Vet. Ani. Sci. 28, 613 (2004).
    Google Scholar 
    Boulêtreau, S. & Santoul, F. The end of the mythical giant catfish. Ecosphere 7(11), e01606. https://doi.org/10.1002/ecs2.1606 (2016).Article 

    Google Scholar 
    Bergmann, C. Ober die verhaltnisse der warmeokonomie der thiere zu ihrer grosse. Gottinger Studien 3, 595–708 (1847).
    Google Scholar 
    Blanck, A. & Lamouroux, N. Large-scale intraspecific variation in life-history traits of European freshwater fish. J. Biogeogr. 34, 862–875 (2007).
    Google Scholar 
    Charnov, E. L., Turner, T. F. & Winemiller, K. O. Reproductive constraints and the evolution of life histories with indeterminate growth. Proc. Natl. Acad. Sci. U.S.A. 98, 9460–9464 (2001).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Ricklefs, R. E. Embryo development and ageing in birds and mammals. Proc. R. Soc. B 273, 2077–2082 (2006).PubMed 
    PubMed Central 

    Google Scholar 
    Lee, W. S., Monaghan, P. & Metcalfe, N. B. Experimental demonstration of the growth rate-lifespan trade-off. Proc. R. Soc. B 280, 20122370 (2013).PubMed 
    PubMed Central 

    Google Scholar 
    Rennie, M. D., Kraft, C., Sprules, W. G. & Johnson, T. B. Factors affecting the growth and condition of lake whitefish (Coregonus clupeaformis). Can. J. Fish. Aquat. Sci. 66, 2096–2108 (2009).
    Google Scholar 
    Prats, J., Val, R., Armengol, J. & Dolz, J. Temporal variability in the thermal regime of the lower Ebro River (Spain) and alteration due to anthropogenic factors. J. Hydrol. 387, 105–118 (2010).ADS 

    Google Scholar 
    Kale, S. & Sönmez, A. Y. Climate change effects on annual streamflow of Filyos River (Turkey). J. Water Clim. Change 11, 420–433 (2020).
    Google Scholar 
    Britton, J. R., Cucherousset, J., Davies, G. D., Godard, M. J. & Copp, G. H. Non-native fishes and climate change: Predicting species responses to warming temperatures in a temperate region. Freshw. Biol. 55, 1130–1141 (2010).
    Google Scholar 
    Garcia, V. B., Lucifora, L. O. & Myers, R. A. The importance of habitat and life history to extinction risk in sharks, skates, rays and chimaeras. Proc. R. Soc. B 275, 83–89 (2008).PubMed 

    Google Scholar 
    Jackson, J. B. C. et al. Historical overfishing and the recent collapse of coastal ecosystems. Science 293, 629–637 (2001).CAS 
    PubMed 

    Google Scholar 
    Kuparinen, A. & Merilä, J. Detecting and managing fisheries-induced evolution. TREE 22, 652–659 (2007).PubMed 

    Google Scholar 
    Swedish University of Agricultural Sciences (SLU). National Data Host Lakes and Watercourses, and National Data Host Agricultural Land (Swedish University of Agricultural Sciences, 2021).
    Google Scholar 
    Emåförbundet. Vattenflöden och Nivåer (n.d.). http://www.eman.se/sv/vattenhushallning/vattenfloden-och-nivaer/historik/. (Accessed 12 May 2021)Fabens, A. J. Properties and fitting of the Von Bertalanffy growth curve. Growth 29, 265–289 (1965).CAS 
    PubMed 

    Google Scholar 
    R Core Team. R: A Language and Environment for Statistical Computing. (R Foundation for Statistical Computing, 2021). https://www.R-project.org/. (Accessed 13 April 2021)Bokor, Z. et al. Survival and growth rates of wels catfish (Siluris glanis Linnaeus, 1758) larvae originating from fertilization with cryopreserved or fresh sperm. J. Appl. Ichthyol. 31, 164–168 (2015).
    Google Scholar 
    du Sert, N. P. et al. The ARRIVE guidelines 2.0: Updated guidelines for reporting animal research. PLoS Biol. 18, e3000410 (2020).
    Google Scholar 
    Horoszewicz, L. & Backiel, T. Growth of Wels (Silurus glanis L.) in the Vistula river and the Zegrzyñski reservoir. Arch. Polish Fish. 11, 115–121 (2003).
    Google Scholar  More

  • in

    Island biogeography and human practices drive ecological connectivity in mosquito species richness in the Lakshadweep Archipelago

    MacArthur, R. H. & Wilson, E. O. The theory of island biogeography (Princeton University Press, 1967).
    Google Scholar 
    MacArthur, R. H. & Wilson, E. O. An equilibrium theory of insular zoogeography. Evolution 17, 373–387 (1968).
    Google Scholar 
    Caraballo, H. Emergency department management of mosquito-borne illness: malaria, dengue, and west nile virus. Emerg. Med. Pract. 16(5), 1–2 (2014).MathSciNet 
    PubMed 

    Google Scholar 
    Rejmánková, E., Grieco, J., Achee, N., Roberts, DR. Ecology of larval habitats. In: Manguin S, editor. Anopheles mosquitoes: new insights into malaria vectors 9th. InTech; Rijeka: pp. 397–446. (2013).Sharma, M., Quader, S., Guttal, V. & Isvaran, K. The enemy of my enemy: multiple interacting selection pressures lead to unexpected anti-predator responses. Oecologia 192(1), 1–12 (2020).ADS 
    PubMed 

    Google Scholar 
    Yee, D. A., Kesavaraju, B. & Juliano, S. A. Interspecific differences in feeding behavior and survival under food-limited conditions for larval Aedes albopictus and Aedes aegypti (Diptera: Culicidae). Ann. Entomol. Soc. Am. 97, 720–728 (2006).
    Google Scholar 
    Messina, J. P. et al. The current and future global distribution and population at risk of dengue. Nat. Microbiol. 4, 1508–1515 (2019).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Rose, N. H. et al. Climate and urbanization drive mosquito preference for humans. Curr. Biol. 30, 3570-3579.e6 (2020).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Day, J. F. Mosquito oviposition behavior and vector control. Insects 7(4), 65 (2016).PubMed Central 

    Google Scholar 
    McBride, C. S. Genes and odors underlying the recent evolution of mosquito preference for humans. Curr. Biol. 26, R41–R46 (2016).MathSciNet 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Southerst, R. W. Global change and human vulnerability to vector-borne diseases. Clin. Microbiol. Rev. 17, 136–173 (2004).
    Google Scholar 
    Vitousek, P. M. Nutrient cycling and limitation: Hawai‘i as a model system (Princeton University Press, 2004).
    Google Scholar 
    Grant, P. R. & Grant, B. R. How and why species multiply: the radiation of darwin’s finches (Princeton University Press, 2011).
    Google Scholar 
    Cliff, A. D. & Haggett, P. The epidemiological significance of islands. Health Place. 1, 199–209 (1995).
    Google Scholar 
    Arrhenius, O. Species and area. J. Ecol. 9(1), 95–99 (1921).
    Google Scholar 
    Preston, F. W. Time and space and the variation of species. Ecology 41(4), 611–627 (1960).
    Google Scholar 
    Rosenzweig, M. L. Species diversity in space and time (Cambridge University Press, 1995).
    Google Scholar 
    Drakare, S. et al. The imprint of the geographical, evolutionary and ecological context on species-area relationships. Ecol. Lett. 9: 215 227. (2006).Kotiaho, J., Kaitala, V., Komonen, A. & Päivinen, J. Predicting the risk of extinction from shared ecological characteristics. Proc. Natl. Acad. Sci. USA 102, 1963–1967 (2005).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Bataille, A. et al. Natural colonization and adaptation of a mosquito species in Galápagos and its implications for disease threats to endemic wildlife. Proc. Nat. Acad. Sci. 106(25), 10230–10235 (2009).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Sinka, M. E. et al. A new malaria vector in Africa: predicting the expansion range of Anopheles stephensi and identifying the urban populations at risk. Proc. Nat. Acad. Sci. 117(40), 24900–24908 (2020).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Powell, J.R. Genetic variation in insect vectors: death of typology? Insects. 11;9(4):139. (2018).Whittaker, R. H. Communities and ecosystems (Macmillan, 1975).
    Google Scholar 
    Nekola, J. C. & White, P. S. The distance decay of similarity in biogeography and ecology. J. Biogeogr. 26, 867–878 (1999).
    Google Scholar 
    Green, J. L. et al. Spatial scaling of microbial eukaryote diversity. Nature 432, 747–750 (2004).ADS 
    CAS 
    PubMed 

    Google Scholar 
    Horner-Devine, M. C., Lage, M. & Hughes, J. B. Bohannan BJ A taxa-area relationship for bacteria. Nature 432, 750–753 (2004).ADS 
    CAS 
    PubMed 

    Google Scholar 
    Martiny, J, B. H., Eisen, J.A., Penn, K., Allison, S.D., Horner-Devine, M.C. Drivers of bacterial beta-diversity depend on spatial scale. Proc. Natl. Acad. Sci. USA 108(19):7850−4. (2011).Segre, H., Ron, R., de Malach, N., Henkin, Z., Mandel, M., Kadmon, R. Competitive exclusion, beta diversity, and deterministic vs. stochastic drivers of community assembly. Ecol. Lett., 17(11):1400−8. (2014).Ishtiaq, F. et al. Biogeographical patterns of blood parasite lineage diversity in avian hosts from southern Melanesian islands. J. Biogeogr. 37, 120–132 (2010).
    Google Scholar 
    Barrera, R., Amador, M. & MacKay, A. J. Population dynamics of Aedes aegypti and dengue as influenced by weather and human behavior in San Juan. Puerto Rico. PLoS Negl. Trop. Dis. 5(12), e1378. https://doi.org/10.1371/journal.pntd.0001378 (2011).Article 
    PubMed 

    Google Scholar 
    Campbell, K. M., Lin, C. D., Iamsirithaworn, S. & Scott, T. W. The complex relationship between weather and dengue virus transmission in Thailand. Am. J. Trop. Med. Hyg. 89, 1066–1080. https://doi.org/10.4269/ajtmh.13-0321 (2013).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Evans, M. V. et al. Microclimate and larval habitat density predict adult Aedes albopictus abundance in Urban Areas. Am. J. Trop. Med. Hyg. 101(2), 362–370 (2019).PubMed 
    PubMed Central 

    Google Scholar 
    Mustak, M. S. et al. The peopling of Lakshadweep Archipelago. Sci. Rep. 9, 6968 (2019).ADS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Sharma, S. K. & Hamzakoya, K. K. Geographical spread of Anopheles stephensi, vector of urban malaria, Aedes aegypti vector of Dengue/DHF, in the Arabian sea islands of Lakshadweep. India. Dengue Bull. 25, 88–91 (2001).
    Google Scholar 
    Sharma RS, Ali, MKS, Dhillon GPS. Epidemiological and entomological aspects of an outbreak of chikungunya in Lakshadweep islands, India, during 2007. Dengue Bull., 178–185 (2008).Subramaniam, H., Ramoo, H. & Sumanam, S. D. Filariasis survey in the Laccadive, minicoy and amindivi Islands. Madras state. Indian J. Malariol. 12, 115–127 (1958).CAS 
    PubMed 

    Google Scholar 
    Roy, R. G., Joy, C. T., Hussain, C. M. & Mohamed, I. K. Malaria in Lakshadweep Islands. Indian J. Med. Res. 67, 924–925 (1978).CAS 
    PubMed 

    Google Scholar 
    Ali, S. M. K. et al. Study on the ecoepidemiology of chikungunya in UT of Lakshadweep. J. Commun. Dis. 41(2), 81–92 (2009).
    Google Scholar 
    Samuel, P. P., Krishnamoorthi, R., Hamzakoya, K. K. & Aggarwal, C. S. Entomo-epidemiological investigations on chikungunya outbreak in the Lakshadweep Islands. Indian Ocean. Indian J. Med. Res. 129(4), 442–445 (2009).PubMed 

    Google Scholar 
    Jayalakshmi, K. & Mathiarasan, L. Prevalence of disease vectors in Lakshadweep Islands during post-monsoon season. J. Vector Borne Dis. 55, 189–196 (2018).
    Google Scholar 
    Su, C. L. et al. Molecular epidemiology of Japanese encephalitis virus in mosquitoes in Taiwan during 2005–2012. PLoS Negl. Trop. Dis. 8, e3122 (2014).PubMed 
    PubMed Central 

    Google Scholar 
    Muslim, A. et al. Armigeres subalbatus incriminated as a vector of zoonotic Brugia pahangi filariasis in suburban Kuala Lumpur. Peninsular Malaysia. Parasites Vectors 6, 219 (2013).PubMed 

    Google Scholar 
    Wilke, A. B. B. et al. Community composition and year-round abundance of vector species of mosquitoes make Miami-Dade County, Florida a receptive gateway for arbovirus entry to the United States. Sci. Rep. 9, 8732 (2019).ADS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Medeiros-Sousa, A. R., Fernandes, A., Ceretti-Junior, W., Wilke, A. B. B. & Marrelli, M. T. Mosquitoes in urban green spaces: using an island biogeographic approach to identify drivers of species richness and composition. Sci. Rep. 7, 17826 (2017).ADS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Lum, J. K., Kaneko, A., Taleo, G., Amos, M. & Reiff, D. M. Genetic diversity and gene flow of humans, Plasmodium falciparum, and Anopheles farauti s.s. of Vanuatu. inferred malaria dispersal and implications for malaria control. Acta Trop. 103, 102–107 (2007).CAS 
    PubMed 

    Google Scholar 
    Marques, T. C. et al. Mosquito (Diptera: Culicidae) assemblages associated with Nidularium and Vriesea bromeliads in Serra do Mar, Atlantic Forest, Brazil. Parasites Vectors 5, 41 (2012).PubMed 
    PubMed Central 

    Google Scholar 
    Laporta, G. Z. & Sallum, M. A. M. Coexistence mechanisms at multiple scales in mosquito assemblages. BMC Ecol. 14, 30 (2014).PubMed 
    PubMed Central 

    Google Scholar 
    Koenraadt, C. J. & Takken, W. Cannibalism and predation among larvae of the Anopheles gambiae complex. Med. Vet. Entomol. 17(1), 61–66 (2003).CAS 
    PubMed 

    Google Scholar 
    Chathuranga, W. G. D., Karunaratne, S. H. P. P., Priyanka, W. A. & De Silva, P. Predator–prey interactions and the cannibalism of larvae of Armigeres subalbatus (Diptera: Culicidae). J. Asia-Pac. Entomol. 23, 124–131 (2020).
    Google Scholar 
    Focks, D. A. & Chadee, D. D. Pupal survey: an epidemiologically significant surveillance method for Aedes aegypti: an example using data from Trinidad. Am. J. Trop. Med. Hyg. 56(2), 159–167 (1997).CAS 
    PubMed 

    Google Scholar 
    Lounibos, L. P., Bargielowski, I., Carrasquilla, M. C. & Nishimura, N. Coexistence of Aedes aegypti and Aedes albopictus (Diptera: Culicidae) in Peninsular Florida two decades after competitive displacements. J. Med. Entomol. 53, 1385–1390 (2016).PubMed 

    Google Scholar 
    Juliano, S. A. Species interactions among larval mosquitoes: context dependence across habitat gradients. Annu. Rev. Entomol. 54, 37–56 (2009).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Bargielowski, I.E., Lounibos, L.P., Carrasquilla, M.C. Evolution of resistance to satyrization through reproductive character displacement in populations of invasive dengue vectors. Proc. Natl. Acad. Sci. 19:110(8):2888–92. (2013).Chadee, D. D. Dengue cases and Aedes aegypti indices in Trinidad. West Indies. Acta Trop. 112(2), 174–180 (2009).CAS 
    PubMed 

    Google Scholar 
    XX. https://www.census2011.co.in/census/state/lakshadweep.htmlChristophers, S. R. The fauna of British India, including Ceylon and Burma; Diptera: Family Culicidae; Tribe Anophelini Vol. 4 (Taylor & Francis, 1933).
    Google Scholar 
    Barraud, P.J. The fauna of British India, including Ceylon and Burma. Diptera V. Family Culicidae. Tribes Megarhinini and Culicini. London: Taylor and Francis p. 463. (1934).Walther, B. A., Cotgreave, P., Price, R. D., Gregory, R. D. & Clayton, D. H. Sampling effort and parasite species richness. Parasitol. Today 11, 306–310 (1995).CAS 
    PubMed 

    Google Scholar 
    Chao, A. Non-parametric estimation of the number of classes in a population. Scand. J. Stat. 11, 265–270 (1984).
    Google Scholar 
    Oksanen, J. et al. Vegan: community ecology package. R Package Version 2(10), 2013 (2015).
    Google Scholar 
    R Core Team. R Development Core Team. R A Lang. Environ. Stat. Comput. 55, 275–286 (2016).McFadden, D. Conditional logit analysis of qualitative choice behavior. Front. Econ. 1, 105–142 (1974).
    Google Scholar 
    Burnham, K. P., Anderson, D. R. & Huyvaert, K. P. AIC model selection and multimodel inference in behavioral ecology: some background, observations, and comparisons. Behav. Ecol. Sociobiol. 65, 23–35 (2011).
    Google Scholar 
    Bray, J. R. & Curtis, J. T. An ordination of the upland forest communities of southern Wisconsin. Ecol. Monograph. 27, 325–349 (1957).
    Google Scholar 
    Sokal, R. R. & Rohlf, F. J. Biometry: the principles and practice of statistics in biological research 3rd edn. (Freeman, 1995).MATH 

    Google Scholar 
    Fortin, M. J. & Dale, M. R. T. Spatial analysis: a guide for ecologists 1–30 (Cambridge University Press, 2005).
    Google Scholar 
    Hartig, F. DHARMa: residual diagnostics for hierarchical (multi-level/mixed) regression models. http://florianhartig.github.io/DHARMa/. (2019).Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).
    Google Scholar 
    World Health Organization, Guidelines for dengue surveillance and mosquito control. Western Pacific Education in Action Series No.8 (WHO, Geneva, 1995) More

  • in

    A species diversity dataset of beetles by three passive acquisition methods in Tei Tong Tsai (Hong Kong)

    Study sitesThe sample site Tei Tong Tsai is located within the Island District (112°5’ E, 22°5’ N Hong Kong, China) and connected to Lantau Country Park. The rich woods in Tei Tong Tsai provide a suitable environment for insects to survive, with rich biodiversity. Weather records (Supplement 1) for May 2019 show that the highefst temperature was 27.2 °C, the lowest was 15.7 °C, the average was 21.7 °C; and the annual average rainfall was 297.8 mm. The suitable temperature and rainfall have created a suitable ecological environment and high biodiversity, establishing Tei Tong Tsai as a prime location for studying beetle diversity. In May 2019, a 13 sample sites were selected for beetle collection (Fig. 1). All latitude and longitude formats were converted to degrees, minutes, and seconds.Fig. 1Sampling points for the three passive acquisition methods used in the Tei Tong Tsai sampling site (indicated by red dots).Full size imageExperimental protocolIn this study, three passive collection methods were used for beetle collection. FIT is an efficient collecting method for insects with strong flying abilities and was first developed and used abroad14. MT and PT collect insects that are not strong flyers and live on the surface. A flight interception trap, a malaise trap, and 10 pitfall traps were set up to collect beetles in each sample site. Samples were selected to cover ecological environments at different longitudes, latitudes, altitudes, and distances from water sources. Reasonable sampling distances (depending on the terrain, with an interval between 100 and 200 m) were set up between sample sites to fully cover Tei Tong Tsai’s habitats. Due to the topography, the distance between the 10th and 11th sample points was about 350 m. The distance between two other close sample points were in the range of 100–200 m. All three traps were based on the original device to maximize the advantages and achieve better collection results.Collection devices. The flight interception trap (Fig. 2a) mainly comprises an interceptor screen (plastic net, PVC plastic glass, or plexiglas) and an insect specimen receiver (PVC), which is an efficient collection device for intercepting and collecting insects with strong flight ability. The detailed installation steps include the following: Firstly, punch two holes on the long side of the PVC screen with a hole puncher spaced about 30 cm apart; then, fix the screen to a bamboo pole with silk, install the specimen receiver, fix all three, bolt the rope, and fix it in the air with a thick rope (the sink is about 0.5–1 m from the ground). After installation, relevant drugs were placed inside the specimen receiver to poison the insects. The drugs used depend on the purpose of the study. For morphological studies, saline (5 mmol/L NaCl solution) or water with detergent is used. By contrast, DNA molecular studies use a mixture of 2% SDS (sodium dodecyl sulfate) and EDTA (ethylene diamine tetraacetic acid, 0.1 mol/L, PH = 8) or highly concentrated alcohol, which effectively controls the degradation of DNA. Currently, high-concentration alcohol, SDS and EDTA mixtures are commonly used. The device is widely applicable and can be installed in almost any habitat; however, it is best installed along the insects’ flight paths, including roads, rivers, or creeks between valleys. In this experiment, we improved this device by increasing the size of the water trough considering the actual situation of the sample site. Also, to properly conduct the molecular experiments, the reagents we used were a mixture of SDS and EDTA. Therefore, the improved device was more suitable for diverse habitats, and the insect species collected were abundant, reflecting good collection practices14.Fig. 2Three passive acquisition methods: (a) flight interception trap; (b) malaise trap; (c) pitfall trap.Full size imageMalaise traps (Fig. 2b) are large tent-like structures constructed from thin mesh. They are among the most commonly used static non-attractant insect traps and insect collection devices. Invented by Malaise (1937) and later improved upon by Townes and Sharkey, these traps are important tools for insect collection and monitoring worldwide15. The malaise trap used at the Tei Tong Tsai Country Park was the Townes type, which is generally set up in forest areas with rich habitats and relatively stable ground. The material is usually meshed mosquito netting fabricated into a tent-shaped insect interception field. The insects hit the net vertically, continue to fly upward, and are gradually led into the trap by the tilted top. The drug in the trap is usually anhydrous ethanol, which intercepts beetles with weak flying abilities16,17.The pitfall trap (Fig. 2c) is an effective method for capturing surface beetles; it is simple to use, easy to carry, and a common device for collection in the wild. The PT is created by digging a pit into the ground with the same depth as a wide-mouth plastic cup (20 cm high, 10 cm in diameter); The upper edge of the cup must be flushed with the soil surface, and a mixture of absolute ethanol is poured inside to collect flightless beetles14. About one-quarter of the way from the top, small holes are punched above the wide-mouth cup to prevent the loss of specimens from rainwater filling the cups. The 10 sets of traps in this experiment were not evenly distributed, but they were all in suitable habitats.Specimen samplingThe sampling site for this study was Tei Tong Tsai, and the sampling period was from 1st May to 28th May (2019). FIT, and PTs were collected once every two days. Due to the small number of beetles collected by MT, mt was collected only once. All beetles were picked out and arranged separately after collection, added to anhydrous ethanol, preserved, and labeled. The beetles collected by the three passive acquisition methods were picked according to morphological species.Specimen identificationThe taxonomic status for the family level of all samples was determined based on the relevant literature18,19,20,21. Relevant experts completed further identification (Supplement 2).All the specimens collected in this study are currently in the zoological museum of the Institute of Zoology, Chinese Academy of Sciences (Beijing, China).Specimen photographyBeetles were poured from the bottle and arranged separately according to the general species. Firstly, we used tweezers or a brush to place the beetles on unbreakable and unwrinkled paper (as far as possible with the backside upwards to keep them tight and neat, reducing the space left, and considering the label in the photograph). Simultaneously, we captured multiple photos according to the size and species of insect for the large specimens in the tube, adjusted the light near them to brighten the background, placed graph paper next to the beetles as a reference scale, then adjusted our Olympus camera settings to the appropriate photographing parameters. Finally, we inserted the photographed beetles and matching labels back into the tube and added anhydrous ethanol for preservation (Fig. 3). The labels were set in the photos as 2019 DTZ-FIT/MT/PTX-5XX-5XX (-N), in which 2019 represents the collection time, DTZ represents Tei Tong Tsai, FIT/MT/PT signifies the collection method, X represents the number of sampling points, 5XX-5XX represents sampling time, and N represents the photo number. If a sample site had many insects on the same date and required more than one photo, n was used to represent the number of photos. See the Supplement 3 for the complete document.Fig. 3Examples of beetles collected from three passive acquisition methods: overall photos of beetles collected by (a) FIT, (b) PT, and (c) MT. On the bottom right corner shows scale in each photo.Full size imageAfter the morphological data of the samples were collected, their Latin name and collection information were recorded in a table. Each passive acquisition method corresponded to a table, and each table was divided into 13 sheets according to 13 sampling points. The collection time was listed horizontally on each sheet, and the beetles’ species names were listed vertically (were named in the morphological species order as 1, 2, 3, …, N). The number of beetles was recorded in the corresponding position and the Supplement 4 file.Finally, data show the beetles’ biodiversity collected from each sampling site. Firstly, we summarized the data from each sampling point after completing the data statistics. Afterward, we counted the number of beetle individuals collected under the different passive acquisition methods at different points (Fig. 4). In Fig. 4, red, blue, and green represent the number of beetle individuals collected by MT, PT, and FIT, respectively. Fig. 4 shows that MT collected fewer beetles than FIT and PT. Secondly, the data of 13 sampling points in each collecting method were summarized to obtain the total number of families and species collected by each method (Fig. 5). A graph created in Excel 2016 displays the collection method as the horizontal coordinate and the number as the vertical coordinate. In the graph, red represents the number of families, and blue represents the number of species. Fig. 5 shows that FIT collected more beetle species and individuals than PT and MT, and MT collected the least. Thirdly, all data from the 13 sampling points and the three collection methods were summarized. The number of species collected in all families was counted. Families with more than ten species were selected (a total of 11 families) for data presentation (Fig. 6). Finally, a graphic was drawn in Excel 2016. Fig. 6 shows that the number of species in Staphylinidae, Curculionidae, and Chrysomelidae accounted for a large number, and the diversity was relatively high.Fig. 4Data table of numbers of individual beetles collected by different methods at 13 sampling points. The red, blue, and green columns represent the number of beetles collected by MT, PT, and FIT, respectively.Full size imageFig. 5The number of beetles collected by different passive acquisition methods. Horizontal coordinates represent collection methods. The red column and blue column represent the number of beetles collected on the family level and species level, respectively.Full size imageFig. 6Families with more than ten species (a total of 11 families) were selected for presentation. The sample sizes of each groups were also shown.Full size image More

  • in

    Ranking threats to biodiversity and why it doesn’t matter

    The difficulties inherent in ranking global threats are due to them being context-dependent, which result from conditions and the nature of the threats themselves differing among locations, habitats, and taxa (Fig. 1). Current high-risk hotspots from habitat loss and overexploitation are primarily located in the tropics, whereas Europe is documented as a threat hotspot for pollution6. On islands, biological invasions mainly threaten biodiversity in the Pacific and Atlantic Oceans, while islands in the Indian Ocean and near the coasts of Asia are mostly threatened by overexploitation and agriculture3. Climate change affects species more at higher latitudes and altitudes because species are constrained by the physical environment (geographic barriers and mountain tops) to follow their optimal isotherms.Fig. 1: Divergence of global threat rankings across different references and international agencies.IPBES, WWF, and IUCN established global rankings of the five threats responsible for the current biodiversity crisis (B: central, yellow panel). However, the relative importance of each threat depends on the taxon, system, species’ characteristics, time, and/or the metric considered, resulting in divergences. Global biodiversity threats are represented by colors and symbols, given in the top panel. This figure encapsulates results combined from different studies detailed in Supplementary Table 1 with their associated references.Full size imageThe relative importance of threats also depends on the taxon considered. At the global scale, vertebrates are primarily threatened by habitat loss, overexploitation, and then biological invasions. But even within the vertebrates rankings differ — birds and mammals are mainly affected by overexploitation, while amphibians have a higher probability of succumbing to habitat loss6. Because of species-specific traits and adaptations, some species are likely to respond differently to global threats even within a clade. Large-bodied vertebrates are more likely to be threatened by overexploitation, whereas small-bodied vertebrates are more prone to habitat loss or pollution (Fig. 1). Threat ranking also depends on the habitat under consideration. Marine mammals are more threatened by overexploitation and pollution than terrestrial mammals for which habitat loss is the primary threat (Fig. 1). On islands, habitat loss is secondary to the pressures of biological invasions in freshwater systems, but the former is more important for terrestrial vertebrates and plants3. Another source of uncertainty is that most studies examining threats are based on well-studied taxa such as terrestrial vertebrates, which only represent a small subset of the tree of life. For instance, only 0.2% of fungi, 1.7% of invertebrates, and 10% of described plants are assessed in the IUCN update of 20197, potentially underestimating the intensity of some threats and biasing conservation priorities for these groups. Similarly, there is a bias of research effort towards regions with high-income countries, while research from low or middle-income countries is generally underrepresented8. This may give the false impression of absence of threats in some regions of the world.Likewise, period-specific global threat ranks are subject to the vagaries of temporal dynamics (Fig. 1). However, distinguishing past, current, and future threats is essential for current or future conservation interventions. Historically, overexploitation caused most of the Pleistocene megafauna extinctions, likely exacerbated by climate change. As agricultural practices intensified, habitat loss played a major role in extinctions. As humans later colonized islands, biological invasions caused the extinction of hundreds of species worldwide3. In contrast, climate change is only predicted to become major in the near future9. In fact, the effects of recent threats might be masked by delayed species’ responses, especially in under-studied regions, resulting in a large extinction debt. For instance, the severity of biological invasions often causes native species to decline rapidly to local extinction, while other threats such as habitat loss might affect species more slowly. In both cases, the eventual extinctions are ultimately if similar magnitude. More

  • in

    A global 0.05° dataset for gross primary production of sunlit and shaded vegetation canopies from 1992 to 2020

    Cox, P. & Jones, C. Climate change – Illuminating the modern dance of climate and CO2. Science 321, 1642–1644 (2008).CAS 
    PubMed 
    Article 

    Google Scholar 
    Gilmanov, T. G. et al. Gross primary production and light response parameters of four Southern Plains ecosystems estimated using long-term CO2-flux tower measurements. Glob. Biogeochem. Cycle 17, 1071 (2003).ADS 
    Article 
    CAS 

    Google Scholar 
    Running, S. W. Climate change – Ecosystem disturbance, carbon, and climate. Science 321, 652–653 (2008).CAS 
    PubMed 
    Article 

    Google Scholar 
    Sun, Z. et al. Spatial pattern of GPP variations in terrestrial ecosystems and its drivers: Climatic factors, CO2 concentration and land-cover change, 1982–2015. Ecol. Inform. 46, 156–165 (2018).CAS 
    Article 

    Google Scholar 
    Running, S. W. et al. A global terrestrial monitoring network integrating tower fluxes, flask sampling, ecosystem modeling and EOS satellite data. Remote Sens. Environ. 70, 108–127 (1999).ADS 
    Article 

    Google Scholar 
    Madani, N. et al. The Impacts of Climate and Wildfire on Ecosystem Gross Primary Productivity in Alaska. J. Geophys. Res.-Biogeosci. 126, e2020JG006078 (2021).ADS 
    Article 

    Google Scholar 
    Morales, P. et al. Comparing and evaluating process-based ecosystem model predictions of carbon and water fluxes in major European forest biomes. Glob. Change Biol. 11, 2211–2233 (2005).ADS 
    Article 

    Google Scholar 
    Tramontana, G., Ichii, K., Camps-Valls, G., Tomelleri, E. & Papale, D. Uncertainty analysis of gross primary production upscaling using Random Forests, remote sensing and eddy covariance data. Remote Sens. Environ. 168, 360–373 (2015).ADS 
    Article 

    Google Scholar 
    Canadell, J. G. et al. Carbon metabolism of the terrestrial biosphere: A multitechnique approach for improved understanding. Ecosystems 3, 115–130 (2000).CAS 
    Article 

    Google Scholar 
    Fletcher, B. J. et al. Photosynthesis and productivity in heterogeneous arctic tundra: consequences for ecosystem function of mixing vegetation types at stand edges. J. Ecol. 100, 441–451 (2012).CAS 
    Article 

    Google Scholar 
    Liu, L., Guan, L. & Liu, X. Directly estimating diurnal changes in GPP for C3 and C4 crops using far-red sun-induced chlorophyll fluorescence. Agr. Forest Meteorol. 232, 1–9 (2017).ADS 
    Article 

    Google Scholar 
    Xu, X. et al. Long-term trend in vegetation gross primary production, phenology and their relationships inferred from the FLUXNET data. J. Environ. Manage. 246, 605–616 (2019).PubMed 
    Article 

    Google Scholar 
    Baldocchi, D. D. How eddy covariance flux measurements have contributed to our understanding of Global Change Biology. Glob. Change Biol. 26, 242–260 (2020).ADS 
    Article 

    Google Scholar 
    He, L., Chen, J. M., Liu, J., Belair, S. & Luo, X. Assessment of SMAP soil moisture for global simulation of gross primary production. J. Geophys. Res.-Biogeosci. 122, 1549–1563 (2017).Article 

    Google Scholar 
    Wang, S., Ibrom, A., Bauer-Gottwein, P. & Garcia, M. Incorporating diffuse radiation into a light use efficiency and evapotranspiration model: An 11-year study in a high latitude deciduous forest. Agr. Forest Meteorol. 248, 479–493 (2018).ADS 
    Article 

    Google Scholar 
    Wang, S. et al. Recent global decline of CO2 fertilization effects on vegetation photosynthesis. Science 370, 1295–1300 (2020).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    Yu, G., Fu, Y., Sun, X., Wen, X. & Zhang, L. Recent progress and future directions of ChinaFLUX. Sci. China Ser. D-Earth Sci. 49, 1–23 (2006).ADS 
    Article 

    Google Scholar 
    McCallum, I. et al. Improved light and temperature responses for light-use-efficiency-based GPP models. Biogeosciences 10, 6577–6590 (2013).ADS 
    Article 

    Google Scholar 
    Stocker, B. D. et al. Drought impacts on terrestrial primary production underestimated by satellite monitoring. Nature Geoscience 12, 264‐+ (2019).ADS 
    Article 
    CAS 

    Google Scholar 
    Cheng, S. J. et al. Variations in the influence of diffuse light on gross primary productivity in temperate ecosystems. Agr. Forest Meteorol. 201, 98–110 (2015).ADS 
    Article 

    Google Scholar 
    Zhang, M. et al. Effects of cloudiness change on net ecosystem exchange, light use efficiency, and water use efficiency in typical ecosystems of China. Agr. Forest Meteorol. 151, 803–816 (2011).ADS 
    Article 

    Google Scholar 
    Oliphant, A. J. et al. The role of sky conditions on gross primary production in a mixed deciduous forest. Agr. Forest Meteorol. 151, 781–791 (2011).ADS 
    Article 

    Google Scholar 
    Urban, O. et al. Ecophysiological controls over the net ecosystem exchange of mountain spruce stand. Comparison of the response in direct vs. diffuse solar radiation. Glob. Change Biol. 13, 157–168 (2007).ADS 
    Article 

    Google Scholar 
    Zhou, H. et al. Large contributions of diffuse radiation to global gross primary productivity during 1981–2015. Glob. Biogeochem. Cycle 35, e2021GB006957 (2021).ADS 
    CAS 
    Article 

    Google Scholar 
    Guanter, L. et al. Retrieval and global assessment of terrestrial chlorophyll fluorescence from GOSAT space measurements. Remote Sens. Environ. 121, 236–251 (2012).ADS 
    Article 

    Google Scholar 
    Guanter, L. et al. Global and time-resolved monitoring of crop photosynthesis with chlorophyll fluorescence. Proc. Natl. Acad. Sci. USA 111, E1327–E1333 (2014).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Liu, L. & Cheng, Z. Detection of vegetation light-use efficiency based on solar-induced chlorophyll fluorescence separated from canopy radiance spectrum. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 3, 306–312 (2010).ADS 
    Article 

    Google Scholar 
    MacBean, N. et al. Strong constraint on modelled global carbon uptake using solar-induced chlorophyll fluorescence data (vol 8, 1973, 2018). Sci. Rep. 8, 10420 (2018).ADS 
    PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Meroni, M. et al. Remote sensing of solar-induced chlorophyll fluorescence: Review of methods and applications. Remote Sens. Environ. 113, 2037–2051 (2009).ADS 
    Article 

    Google Scholar 
    Zheng, T. & Chen, J. M. Photochemical reflectance ratio for tracking light use efficiency for sunlit leaves in two forest types. ISPRS-J. Photogramm. Remote Sens. 123, 47–61 (2017).ADS 
    Article 

    Google Scholar 
    Damm, A. et al. Remote sensing of sun-induced fluorescence to improve modeling of diurnal courses of gross primary production (GPP). Glob. Change Biol. 16, 171–186 (2010).ADS 
    Article 

    Google Scholar 
    Lee, J. E. et al. Simulations of chlorophyll fluorescence incorporated into the Community Land Model version 4. Glob. Change Biol. 21, 3469–3477 (2015).ADS 
    Article 

    Google Scholar 
    Pinto, F. et al. Sun-induced chlorophyll fluorescence from high-resolution imaging spectroscopy data to quantify spatio-temporal patterns of photosynthetic function in crop canopies. Plant Cell Environ. 39, 1500–1512 (2016).CAS 
    PubMed 
    Article 

    Google Scholar 
    Porcar-Castell, A. et al. Linking chlorophyll a fluorescence to photosynthesis for remote sensing applications: mechanisms and challenges. J. Exp. Bot. 65, 4065–4095 (2014).CAS 
    PubMed 
    Article 

    Google Scholar 
    Xie, X., Li, A., Jin, H., Yin, G. & Nan, X. Derivation of temporally continuous leaf maximum carboxylation rate (V-cmax) from the sunlit leaf gross photosynthesis productivity through combining BEPS model with light response curve at tower flux sites. Agr. Forest Meteorol. 259, 82–94 (2018).ADS 
    Article 

    Google Scholar 
    Chen, J. M., Liu, J., Leblanc, S. G., Lacaze, R. & Roujean, J. L. Multi-angular optical remote sensing for assessing vegetation structure and carbon absorption. Remote Sens. Environ. 84, 516–525 (2003).ADS 
    Article 

    Google Scholar 
    Chen, J. M. et al. Effects of foliage clumping on the estimation of global terrestrial gross primary productivity. Glob. Biogeochem. Cycle 26, GB1019 (2012).ADS 
    Article 
    CAS 

    Google Scholar 
    Running, S. W., Thornton, P. E., Nemani, R. & Glassy, J. M. in Methods in Ecosystem Science. Ch.3 (Springer, New York, NY. Press, 2000).Wu, C., Munger, J. W., Niu, Z. & Kuang, D. Comparison of multiple models for estimating gross primary production using MODIS and eddy covariance data in Harvard Forest. Remote Sens. Environ. 114, 2925–2939 (2010).ADS 
    Article 

    Google Scholar 
    Makela, A. et al. Developing an empirical model of stand GPP with the LUE approach: analysis of eddy covariance data at five contrasting conifer sites in Europe. Glob. Change Biol. 14, 92–108 (2008).ADS 
    Article 

    Google Scholar 
    McCallum, I. et al. Satellite-based terrestrial production efficiency modeling. Carbon Balanc. Manag. 4, 8–8 (2009).Article 

    Google Scholar 
    Wang, H. et al. Deriving maximal light use efficiency from coordinated flux measurements and satellite data for regional gross primary production modeling. Remote Sens. Environ 114, 2248–2258 (2010).ADS 
    Article 

    Google Scholar 
    Yu, R. An improved estimation of net primary productivity of grassland in the Qinghai-Tibet region using light use efficiency with vegetation photosynthesis model. Ecol. Model. 431, 109121 (2020).Article 

    Google Scholar 
    Yuan, W. et al. Deriving a light use efficiency model from eddy covariance flux data for predicting daily gross primary production across biomes. Agr. Forest Meteorol. 143, 189–207 (2007).ADS 
    Article 

    Google Scholar 
    Beer, C. et al. Terrestrial gross carbon dioxide uptake: global distribution and covariation with climate. Science 329, 834–838 (2010).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    Running, S. W. et al. A continuous satellite-derived measure of global terrestrial primary production. Bioscience 54, 547–560 (2004).Article 

    Google Scholar 
    Zhang, Y. et al. Development of a coupled carbon and water model for estimating global gross primary productivity and evapotranspiration based on eddy flux and remote sensing data. Agr. Forest Meteorol. 223, 116–131 (2016).ADS 
    Article 

    Google Scholar 
    He, M. et al. Development of a two-leaf light use efficiency model for improving the calculation of terrestrial gross primary productivity. Agr. Forest Meteorol. 173, 28–39 (2013).ADS 
    Article 

    Google Scholar 
    Zhou, Y. et al. Global parameterization and validation of a two-leaf light use efficiency model for predicting gross primary production across FLUXNET sites. J. Geophys. Res.-Biogeosci. 121, 1045–1072 (2016).Article 

    Google Scholar 
    Friedlingstein, P. et al. Uncertainties in CMIP5 Climate Projections due to Carbon Cycle Feedbacks. J. Clim. 27, 511–526 (2014).ADS 
    Article 

    Google Scholar 
    Raich, J. W. et al. Potential net primary productivity in South-America – application of a global-model. Ecol. Appl. 1, 399–429 (1991).CAS 
    PubMed 
    Article 

    Google Scholar 
    Li, J. et al. An algorithm differentiating sunlit and shaded leaves for improving canopy conductance and vapotranspiration estimates. J. Geophys. Res.-Biogeosci. 124, 807–824 (2019).ADS 
    Article 

    Google Scholar 
    Chen, J. M., Liu, J., Cihlar, J. & Goulden, M. L. Daily canopy photosynthesis model through temporal and spatial scaling for remote sensing applications. Ecol. Model. 124, 99–119 (1999).CAS 
    Article 

    Google Scholar 
    Keenan, T. F. et al. Recent pause in the growth rate of atmospheric CO2 due to enhanced terrestrial carbon uptake. Nat. Commun. 7, 13428 (2016).ADS 
    CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Huang, M. et al. Air temperature optima of vegetation productivity across global biomes. Nat. Ecol. Evol. 3, 772–779 (2019).PubMed 
    PubMed Central 
    Article 

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

    Google Scholar 
    Korson, L., Drosthan, W. & Millero, F. J. Viscosity of water at various temperatures. J. Phys. Chem. 73, 34–39 (1969).CAS 
    Article 

    Google Scholar 
    Olofsson, P., Van Laake, P. E. & Eklundh, L. Estimation of absorbed PAR across Scandinavia from satellite measurements Part I: Incident PAR. Remote Sens. Environ. 110, 252–261 (2007).ADS 
    Article 

    Google Scholar 
    González, J. A. & Calbó, J. Modelled and measured ratio of PAR to global radiation under cloudless skies. Agr. Forest Meteorol. 110, 319–325 (2002).ADS 
    Article 

    Google Scholar 
    Zhang, X., Zhang, Y. & Zhoub, Y. Measuring and modelling photosynthetically active radiation in Tibet Plateau during April–October. Agr. Forest Meteorol. 102, 207–212 (2000).ADS 
    Article 

    Google Scholar 
    Yang, Y., Xiao, P., Feng, X. & Li, H. Accuracy assessment of seven global land cover datasets over China. ISPRS-J. Photogramm. Remote Sens. 125, 156–173 (2017).ADS 
    Article 

    Google Scholar 
    Liu, Y., Liu, R. & Chen, J. M. GLOBMAP global Leaf Area Index since 1981. Zenodo https://doi.org/10.5281/zenodo.4700264 (2019).Vermote, E. MOD09A1 MODIS/Terra Surface Reflectance 8-Day L3 Global 500m SIN Grid V006. NASA EOSDIS Land Processes DAAC https://doi.org/10.5067/MODIS/MOD09A1.006 (2015).Deng, F., Chen, J. M., Plummer, S., Chen, M. & Pisek, J. Algorithm for global leaf area index retrieval using satellite imagery. IEEE Trans. Geosci. Remote Sens. 44, 2219–2229 (2006).ADS 
    Article 

    Google Scholar 
    Vermote, E. NOAA CDR Program. NOAA Climate Data Record (CDR) of AVHRR Leaf Area Index (LAI) and Fraction of Absorbed Photosynthetically Active Radiation (FAPAR), Version 5. LAI. NOAA National Centers for Environmental Information https://doi.org/10.7289/V5TT4P69 (2019).He, L., Chen, J. M., Pisek, J., Schaaf, C. & Strahler, A. Global clumping index map derived from the MODIS BRDF product. Remote Sens. Environ. 119, 118–130 (2012).ADS 
    Article 

    Google Scholar 
    Liu, R. G. & Liu, Y. Generation of new cloud masks from MODIS land surface reflectance products. Remote Sens. Environ. 133, 21–37 (2013).ADS 
    CAS 
    Article 

    Google Scholar 
    Chen, J. M., Deng, F. & Chen, M. Locally adjusted cubic-spline capping for reconstructing seasonal trajectories of a satellite-derived surface parameter. IEEE Trans. Geosci. Remote Sens. 44, 2230–2238 (2006).ADS 
    Article 

    Google Scholar 
    Harris, I.C. CRU JRA: Collection of CRU JRA forcing datasets of gridded land surface blend of Climatic Research Unit (CRU) and Japanese reanalysis (JRA) data. Centre for Environmental Data Analysis http://catalogue.ceda.ac.uk/uuid/863a47a6d8414b6982e1396c69a9efe8 (2019).Li, X., Liang, H. & Cheng, W. Evaluation and comparison of light use efficiency models for their sensitivity to the diffuse PAR fraction and aerosol loading in China. Int. J. Appl. Earth Obs. Geoinf. 95, 102269 (2021).
    Google Scholar 
    Duan, Q. Y., Sorooshian, S. & Gupta, V. Effective and efficient global optimization for conceptual rain full-runoff models. Water Resour. Res. 28, 1015–1031 (1992).ADS 
    Article 

    Google Scholar 
    Gu, L. H. et al. Advantages of diffuse radiation for terrestrial ecosystem productivity. J. Geophys. Res.-Atmos. 107, 4050 (2002).ADS 

    Google Scholar 
    Bi, W. & Zhou, Y. A global 0.05° dataset for gross primary production of sunlit and shaded vegetation canopies (1992–2020). Dryad https://doi.org/10.5061/dryad.dfn2z352k (2022).Ogutu, B. O. & Dash, J. Assessing the capacity of three production efficiency models in simulating gross carbon uptake across multiple biomes in conterminous USA. Agr. Forest Meteorol. 174, 158–169 (2013).ADS 
    Article 

    Google Scholar 
    Cai, W. et al. Large differences in terrestrial vegetation production derived from satellite-based light use efficiency models. Remote Sens. 6, 8945–8965 (2014).ADS 
    Article 

    Google Scholar 
    Anav, A. et al. Spatiotemporal patterns of terrestrial gross primary production: a review. Rev. Geophys. 53, 785–818 (2015).ADS 
    Article 

    Google Scholar 
    Li, X. & Xiao, J. Mapping photosynthesis solely from solar-induced chlorophyll fluorescence: A global, fine-resolution dataset of gross primary production derived from OCO-2. Remote Sens. 11, 2563 (2019).ADS 
    Article 

    Google Scholar 
    Alemohammad, S. H. et al. Water, Energy, and Carbon with Artificial Neural Networks (WECANN): a statistically based estimate of global surface turbulent fluxes and gross primary productivity using solar-induced fluorescence. Biogeosciences 14, 4101–4124 (2017).ADS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Joiner, J. et al. Estimation of terrestrial global gross primary production (GPP) with satellite data-driven models and eddy covariance flux data. Remote Sens. 10, 1346 (2018).ADS 
    Article 

    Google Scholar 
    Wang, S., Zhang, Y., Ju, W., Qiu, B. & Zhang, Z. Tracking the seasonal and inter-annual variations of global gross primary production during last four decades using satellite near-infrared reflectance data. Sci. Total Environ. 755, 142569 (2021).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    Zheng, Y. et al. Improved estimate of global gross primary production for reproducing its long-term variation, 1982–2017. Earth Syst. Sci. Data 12, 2725–2746 (2020).ADS 
    Article 

    Google Scholar 
    Running, S., Mu, Q. & Zhao, M. MOD17A2H MODIS/Terra Gross Primary Productivity 8-Day L4 Global 500m SIN Grid V006. NASA EOSDIS Land Processes DAAC https://doi.org/10.5067/MODIS/MOD17A2H.006 (2015).Ciais, P. et al. A three-dimensional synthesis study of delta O-18 in atmospheric CO2 .1. Surface fluxes. J. Geophys. Res.-Atmos. 102, 5857–5872 (1997).ADS 
    CAS 
    Article 

    Google Scholar 
    Zhang, Y., Joiner, J., Gentine, P. & Zhou, S. Reduced solar-induced chlorophyll fluorescence from GOME-2 during Amazon drought caused by dataset artifacts. Glob. Change Biol. 24, 2229–2230 (2018).ADS 
    Article 

    Google Scholar 
    Xie, X. et al. Assessment of five satellite-derived LAI datasets for GPP estimations through ecosystem models. Sci. Total Environ. 690, 1120–1130 (2019).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    Fang, H., Wei, S., Jiang, C. & Scipal, K. Theoretical uncertainty analysis of global MODIS, CYCLOPES, and GLOBCARBON LAI products using a triple collocation method. Remote Sens. Environ. 124, 610–621 (2012).ADS 
    Article 

    Google Scholar 
    Camacho, F., Cemicharo, J., Lacaze, R., Baret, F. & Weiss, M. GEOV1: LAI, FAPAR essential climate variables and FCOVER global time series capitalizing over existing products. Part 2: Validation and intercomparison with reference products. Remote Sens. Environ. 137, 310–329 (2013).ADS 
    Article 

    Google Scholar 
    Prince, S. D. & Goward, S. N. Global primary production: A remote sensing approach. J. Biogeogr. 22, 815–835 (1995).Article 

    Google Scholar 
    Verma, S. B. et al. Annual carbon dioxide exchange in irrigated and rainfed maize-based agroecosystems. Agr. Forest Meteorol. 131, 77–96 (2005).ADS 
    Article 

    Google Scholar 
    Yan, H. et al. Improved global simulations of gross primary product based on a new definition of water stress factor and a separate treatment of C3 and C4 plants. Ecol. Model. 297, 42–59 (2015).CAS 
    Article 

    Google Scholar 
    Jiang, S. et al. Comparison of satellite-based models for estimating gross primary productivity in agroecosystems. Agr. Forest Meteorol. 297, 108253 (2021).ADS 
    Article 

    Google Scholar 
    Yang, X. et al. Solar-induced chlorophyll fluorescence that correlates with canopy photosynthesis on diurnal and seasonal scales in a temperate deciduous forest. Geophys. Res. Lett. 42, 2977–2987 (2015).ADS 
    CAS 
    Article 

    Google Scholar 
    Zhou, H. et al. Responses of gross primary productivity to diffuse radiation at global FLUXNET sites. Atmos. Environ. 244, 117905 (2021).CAS 
    Article 

    Google Scholar 
    Han, J. et al. Effects of diffuse photosynthetically active radiation on gross primary productivity in a subtropical coniferous plantation vary in different timescales. Ecol. Indic. 115, 106403 (2020).Article 

    Google Scholar 
    Grant, I. F., Prata, A. J. & Cechet, R. P. The impact of the diurnal variation of albedo on the remote sensing of the daily mean albedo of grassland. J. Appl. Meteorol. 39, 231–244 (2000).ADS 
    Article 

    Google Scholar 
    Singarayer, J. S., Ridgwell, A. & Irvine, P. Assessing the benefits of crop albedo bio-geoengineering. Environ. Res. Lett. 4, 045110 (2009).ADS 
    Article 

    Google Scholar 
    Tang, S. et al. LAI inversion algorithm based on directional reflectance kernels. J. Environ. Manage. 85, 638–648 (2007).CAS 
    PubMed 
    Article 

    Google Scholar  More

  • in

    Impacts of larval host plant species on dispersal traits and free-flight energetics of adult butterflies

    Ehrlich, P. R. & Raven, P. H. Butterflies and plants: A study in coevolution. Evolution 18, 586 (1964).Article 

    Google Scholar 
    Raguso, R. A. et al. The raison d’être of chemical ecology. Ecology 96, 617–630 (2015).PubMed 
    Article 

    Google Scholar 
    Kariyat, R. R. & Portman, S. L. Plant–herbivore interactions: Thinking beyond larval growth and mortality. Am. J. Bot. 103, 789–791 (2016).PubMed 
    Article 

    Google Scholar 
    Raubenheimer, D. & Simpson, S. J. Nutritional ecology and foraging theory. Curr. Opin. Insect Sci. 27, 38–45 (2018).PubMed 
    Article 

    Google Scholar 
    Goehring, L. & Oberhauser, K. S. Effects of photoperiod, temperature, and host plant age on induction of reproductive diapause and development time in Danaus plexippus. Ecol. Entomol. 27, 674–685 (2002).Article 

    Google Scholar 
    Hahn, D. A. Larval nutrition affects lipid storage and growth, but not protein or carbohydrate storage in newly eclosed adults of the grasshopper Schistocerca americana. J. Insect Physiol. 51, 1210–1219 (2005).CAS 
    PubMed 
    Article 

    Google Scholar 
    Portman, S. L., Kariyat, R. R., Johnston, M. A., Stephenson, A. G. & Marden, J. H. Cascading effects of host plant inbreeding on the larval growth, muscle molecular composition, and flight capacity of an adult herbivorous insect. Funct. Ecol. 29, 328–337 (2015).Article 

    Google Scholar 
    Johnson, C. G. Physiological factors in insect migration by flight. Nature 198, 423–427 (1963).Article 

    Google Scholar 
    Harrison, R. G. Dispersal polymorphisms in insects. Annu. Rev. Ecol. Syst. 11, 95–118 (1980).Article 

    Google Scholar 
    Zera, A. J. & Denno, R. F. Physiology and ecology of dispersal polymorphism in insects. Annu. Rev. Entomol. 42, 207–231 (1997).CAS 
    PubMed 
    Article 

    Google Scholar 
    Marden, J. H. et al. Weight and nutrition affect pre-mRNA splicing of a muscle gene associated with performance, energetics and life history. J. Exp. Biol. 211, 3653–3660 (2008).CAS 
    PubMed 
    Article 

    Google Scholar 
    Raguso, R. A., Ojeda-Avila, T., Desai, S., Jurkiewicz, M. A. & Arthur Woods, H. The influence of larval diet on adult feeding behaviour in the tobacco hornworm moth, Manduca sexta. J. Insect Physiol. 53, 923–932 (2007).CAS 
    PubMed 
    Article 

    Google Scholar 
    Cease, A. J. et al. Nutritional imbalance suppresses migratory phenotypes of the Mongolian locust (Oedaleus asiaticus). R. Soc. Open Sci. 4, https://doi.org/10.1098/rsos.161039 (2017).Reichstein, T., Von Euw, J., Parsons, J. A. & Rothschild, M. Heart poisons in the monarch butterfly. Science 161, 861–866 (1968).CAS 
    PubMed 
    Article 

    Google Scholar 
    Brower, L. P., Ryerson, W. N., Coppinger, L. L. & Glazier, S. C. Ecological chemistry and the palatability spectrum. Science 161, 1349–1351 (1968).CAS 
    PubMed 
    Article 

    Google Scholar 
    Young, A. M. An evolutionary-ecological model of the evolution of migratory behavior in the Monarch Butterfly, and its absence in the Queen Butterfly. Acta Biotheor. 31, 219–237 (1982).Article 

    Google Scholar 
    Agrawal, A. A. Monarchs and Milkweed: A Migrating Butterfly, a Poisonous Plant, and Their Remarkable Story of Coevolution. (Princeton University Press, 2017).Batalden, R. V. & Oberhauser, K. S. Potential changes in eastern north American monarch migration in response to an introduced Milkweed, Asclepias curassavica. in Monarchs in a Changing World: Biology and Conservation of an Iconic Butterfly 215–224 (2015).Tyler Flockhart, D. T. et al. Tracking multi-generational colonization of the breeding grounds by monarch butterflies in eastern North America. Proc. R. Soc. B Biol. Sci. 280, 20131087 (2013).Saunders, S. P., Ries, L., Oberhauser, K. S., Thogmartin, W. E. & Zipkin, E. F. Local and cross-seasonal associations of climate and land use with abundance of monarch butterflies Danaus plexippus. Ecography. 41, 278–290 (2018).Article 

    Google Scholar 
    Pleasants, J. M. & Oberhauser, K. S. Milkweed loss in agricultural fields because of herbicide use: Effect on the monarch butterfly population. Insect Conserv. Divers. 6, 135–144 (2013).Article 

    Google Scholar 
    Borders, B. & Lee-Mäder, B. B. Project milkweed. in Monarchs in a Changing World: Biology and Conservation of an Iconic Butterfly. pp.190-196 (Cornell University press, 2015).Agrawal, A. A., Petschenka, G., Bingham, R. A., Weber, M. G. & Rasmann, S. Toxic cardenolides: Chemical ecology and coevolution of specialized plant-herbivore interactions. N. Phytologist 194, 28–45 (2012).CAS 
    Article 

    Google Scholar 
    Malcolm, S. B. Milkweeds, monarch butterflies and the ecological significance of cardenolides. Chemoecology 5–6, 101–117 (1994).Article 

    Google Scholar 
    Pocius, V. M., Debinski, D. M., Bidne, K. G., Hellmich, R. L. & Hunter, F. K. Performance of early Instar Monarch Butterflies (Danaus plexippus L.) on nine Milkweed species native to Iowa. J. Lepid. Soc. 71, 153–161 (2017).
    Google Scholar 
    Ali, J. G. & Agrawal, A. A. Specialist versus generalist insect herbivores and plant defense. Trends Plant Sci. 17, 293–302 (2012).CAS 
    PubMed 
    Article 

    Google Scholar 
    Zalucki, M. P., Brower, L. P. & Alonso-M, A. Detrimental effects of latex and cardiac glycosides on survival and growth of first-instar monarch butterfly larvae Danaus plexippus feeding on the sandhill milkweed Asclepias humistrata. Ecol. Entomol. 26, 212–224 (2001).Article 

    Google Scholar 
    Agrawal, A. A., Hastings, A. P., Patrick, E. T. & Knight, A. C. Specificity of herbivore-induced hormonal signaling and defensive traits in five closely related milkweeds (Asclepias spp.). J. Chem. Ecol. 40, 717–729 (2014).CAS 
    PubMed 
    Article 

    Google Scholar 
    Agrawal, A. A., Ali, J. G., Rasmann, S. & Fishbein, M. Macroevolutionary trends in the defense of milkweeds against monarchs. Monarch. a Chang. World Biol. Conserv. Iconic Insect. Cornell University Press, Ithaca, NY. pp. 47–59 (2011).Pocius, V. M. et al. Milkweed matters: Monarch butterfly (Lepidoptera: Nymphalidae) survival and development on nine midwestern milkweed species. Environ. Entomol. 46, 1098–1105 (2017).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Petschenka, G. et al. Stepwise evolution of resistance to toxic cardenolides via genetic substitutions in the na+/k+-atpase of milkweed butterflies (lepidoptera: Danaini). Evolution (N. Y). 67, 2753–2761 (2013).CAS 

    Google Scholar 
    Agrawal, A. A. et al. Cardenolides, toxicity, and the costs of sequestration in the coevolutionary interaction between monarchs and milkweeds. Proc. Natl Acad. Sci. USA 118, e2024463118 (2021).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Marden, J. H. Variability in the size, composition, and function of insect flight muscles. Annu. Rev. Physiol. 62, 157–178 (2000).CAS 
    PubMed 
    Article 

    Google Scholar 
    Bicudo, J. E. P. W., Buttemer, W. A., Chappell, M. A., Pearson, J. T. & Bech, C. Ecological and Environmental Physiology of Birds. Ecological and Environmental Physiology of Birds 3 (Oxford University Press, 2010).Bailey, E. Biochemistry of Insect Flight. in Insect Biochemistry and Function. pp. 89–176 (Springer, 1975).Dudley, R. The biomechanics of insect flight: form, function, evolution. Annals of the Entomological Society of America 93 (Princeton University Press, 2000).Solensky, M. J. Overview of monarch migration. in The Monarch Butterfly: Biology and Conservation 79–83 (2004).Urquhart, F. A. & Urquhart, N. R. Monarch butterfly (Danaus plexippus L.) overwintering population in Mexico (Lep. Danaidae). Atalanta 7, 56–61 (1976).
    Google Scholar 
    Brower, L. P. Understanding and misunderstanding the migration of the monarch butterfly (Nymphalidae) in North America: 1857–1995. J. – Lepid. Soc. 49, 304–385 (1995).
    Google Scholar 
    Fisher, K. E., Adelman, J. S. & Bradbury, S. P. Employing Very High Frequency (VHF) radio telemetry to recreate monarch butterfly flight paths. Environ. Entomol. 49, 312–323 (2020).PubMed 
    Article 

    Google Scholar 
    Reppert, S. M. & de Roode, J. C. Demystifying monarch butterfly migration. Curr. Biol. 28, R1009–R1022 (2018).CAS 
    PubMed 
    Article 

    Google Scholar 
    Zhu, H., Gegear, R. J., Casselman, A., Kanginakudru, S. & Reppert, S. M. Defining behavioral and molecular differences between summer and migratory monarch butterflies. BMC Biol. 7, 1–14 (2009).Heinze, S. & Reppert, S. M. Anatomical basis of sun compass navigation I: The general layout of the monarch butterfly brain. J. Comp. Neurol. 520, 1599–1628 (2012).PubMed 
    Article 

    Google Scholar 
    Zhan, S. et al. The genetics of monarch butterfly migration and warning colouration. Nature 514, 317–321 (2014).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Soule, A. J., Decker, L. E. & Hunter, M. D. Effects of diet and temperature on monarch butterfly wing morphology and flight ability. J. Insect Conserv. 24, 961–975 (2020).Article 

    Google Scholar 
    Decker, L. E., Soule, A. J., de Roode, J. C. & Hunter, M. D. Phytochemical changes in milkweed induced by elevated CO2 alter wing morphology but not toxin sequestration in monarch butterflies. Funct. Ecol. 33, 411–421 (2019).Article 

    Google Scholar 
    Heinrich, B. Temperature regulation of the sphinx moth, Manduca sexta. I. Flight energetics and body temperature during free and tethered flight. J. Exp. Biol. 54, 141–152 (1971).CAS 
    PubMed 
    Article 

    Google Scholar 
    Nicolson, S. W. & Louw, G. N. Simultaneous measurement of evaporative water loss, oxygen consumption, and thoracic temperature during flight in a carpenter bee. J. Exp. Zool. 222, 287–296 (1982).Article 

    Google Scholar 
    Rothe, U. & Nachtigall, W. Flight of the honey bee IV. J. Comp. Physiol. B 158, 711–718 (1989).Article 

    Google Scholar 
    Nachtigall, W., Hanauer-Thieser, U. & Mörz, M. Flight of the honey bee VII: Metabolic power versus flight speed relation. J. Comp. Physiol. B 165, 484–489 (1995).Article 

    Google Scholar 
    Niven, J. E. & Scharlemann, J. P. W. Do insect metabolic rates at rest and during flight scale with body mass? Biol. Lett. 1, 346–349 (2005).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Zalucki, M. P., Parry, H. R. & Zalucki, J. M. Movement and egg laying in Monarchs: To move or not to move, that is the equation. Austral. Ecol. 41, 154–167 (2016).Article 

    Google Scholar 
    Marden, J. H. & Chai, Peng Aerial predation and butterfly design: How palatability, mimicry, and the need for evasive flight constrain mass allocation. Am. Nat. 138, 15–36 (1991).Article 

    Google Scholar 
    Levin, E., Lopez-Martinez, G., Fane, B. & Davidowitz, G. Hawkmoths use nectar sugar to reduce oxidative damage from flight. Science 355, 733–735 (2017).CAS 
    PubMed 
    Article 

    Google Scholar 
    Petschenka, G. & Agrawal, A. A. Milkweed butterfly resistance to plant toxins is linked to sequestration, not coping with a toxic diet. Proc. R. Soc. B Biol. Sci. 282, 20151865 (2015).Petschenka, G. & Agrawal, A. A. How herbivores coopt plant defenses: Natural selection, specialization, and sequestration. Curr. Opin. Insect Sci. 14, 17–24 (2016).PubMed 
    Article 

    Google Scholar 
    Tan, W. H., Tao, L., Hoang, K. M., Hunter, M. D. & de Roode, J. C. The effects of milkweed induced defense on parasite resistance in monarch butterflies, Danaus plexippus. J. Chem. Ecol. 44, 1040–1044 (2018).CAS 
    PubMed 
    Article 

    Google Scholar 
    Brower, L. P. & Glazier, S. C. Localization of heart poisons in the monarch butterfly. Science 188, 19–25 (1975).CAS 
    PubMed 
    Article 

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

    Google Scholar 
    Zalucki, M. P., Malcolm, S. B., Hanlon, C. C. & Paine, T. D. First-instar monarch larval growth and survival on milkweeds in Southern California: Effects of latex, leaf hairs and cardenolides. Chemoecology 22, 75–88 (2012).Article 

    Google Scholar 
    Ziegler, R. & Van Antwerpen, R. Lipid uptake by insect oocytes. Insect Biochem. Mol. Biol. 36, 264–272 (2006).Beenakkers, A. M. T., Van der Horst, D. J. & Van Marrewijk, W. J. A. Insect flight muscle metabolism. Insect Biochem. 14, 243–260 (1984).CAS 
    Article 

    Google Scholar 
    Beall, G. The fat content of a butterfly, Danaus Plexippus Linn., as affected by migration. Ecology 29, 80–94 (1948).Article 

    Google Scholar 
    James, D. G. Phenology of weight, moisture and energy reserves of Australian monarch butterflies, Danaus plexippus. Ecol. Entomol. 9, 421–428 (1984).Article 

    Google Scholar 
    Briegel, H. Metabolic relationship between female body size, reserves, and fecundity of Aedes aegypti. J. Insect Physiol. 36, 165–172 (1990).Article 

    Google Scholar 
    Hines, W. J. W. & Smith, M. J. H. Some aspects of intermediary metabolism in the desert locust (Schistocerca gregaria Forskål). J. Insect Physiol. 9, 463–468 (1963).CAS 
    Article 

    Google Scholar 
    Inagaki, S. & Yamashita, O. Metabolic shift from lipogenesis to glycogenesis in the last instar larval fat body of the silkworm, Bombyx mori. Insect Biochem. 16, 327–331 (1986).CAS 
    Article 

    Google Scholar 
    Venkatesh, K. & Morrison, P. E. Studies of weight changes and amount of food ingested by the stable fly, stomoxys calcitrans (Diptera: Muscidae). Can. Entomol. 112, 141–149 (1980).Article 

    Google Scholar 
    Arrese, E. L. & Soulages, J. L. Insect fat body: Energy, metabolism, and regulation. Annu. Rev. Entomol. 55, 207–225 (2010).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Mevi-Schütz, J. & Erhardt, A. Larval nutrition affects female nectar amino acid preference in the map butterfly (Araschnia levana). Ecology 84, 2788–2794 (2003).Article 

    Google Scholar 
    Wassenaar, L. I. & Hobson, K. A. Natal origins of migratory monarch butterflies at wintering colonies in Mexico: New isotopic evidence. Proc. Natl Acad. Sci. USA 95, 15436–15439 (1998).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Majewska, A. A. & Altizer, S. Exposure to Non-Native Tropical Milkweed Promotes Reproductive Development in Migratory Monarch Butterflies. Insects 10, 253 (2019).Howard, E., Aschen, H. & Davis, A. K. Citizen science observations of monarch butterfly overwintering in the Southern United States. Psyche: A Journal of Entomology 2010, https://doi.org/10.1155/2010/689301 (2010).Satterfield, D. A., Maerz, J. C. & Altizer, S. Loss of migratory behaviour increases infection risk for a butterfly host. Proc. R. Soc. B Biol. Sci. 282, 20141734 (2015).Petschenka, G. et al. Relative selectivity of plant cardenolides for Na+/K+-ATPases from the monarch butterfly and non-resistant insects. Front. Plant Sci. 9, 1424 (2018).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Jones, P. L., Petschenka, G., Flacht, L. & Agrawal, A. A. Cardenolide intake, sequestration, and excretion by the monarch butterfly along gradients of plant toxicity and larval ontogeny. J. Chem. Ecol. 45, 264–277 (2019).CAS 
    PubMed 
    Article 

    Google Scholar 
    Tao, L., Hoang, K. M., Hunter, M. D. & de Roode, J. C. Fitness costs of animal medication: antiparasitic plant chemicals reduce fitness of monarch butterfly hosts. J. Anim. Ecol. 85, 1246–1254 (2016).PubMed 
    Article 

    Google Scholar 
    Lederhouse, R. C. The effect of female mating frequency on egg fertility in the black swallowtail, Papilio polyxenes asterius (Papilionidae). J. Lepid. Soc. 35, 266–277 (1981).
    Google Scholar 
    Jones, R. E., Hart, J. R. & Bull, G. D. Temperature, size and egg production in the Cabbage Butterfly, Pieris rapae L. Aust. J. Zool. 30, 159–168 (1982).Article 

    Google Scholar 
    Haukioja, E. & Neuvonen, S. The relationship between size and reproductive potential in male and female Epirrita autumnata (Lep., Geometridae). Ecol. Entomol. 10, 267–270 (1985).Article 

    Google Scholar 
    Altizer, S. M., Oberhauser, K. S. & Brower, L. P. Associations between host migration and the prevalence of a protozoan parasite in natural populations of adult monarch butterflies. Ecol. Entomol. 25, 125–139 (2000).Article 

    Google Scholar 
    Masters, A. R., Malcolm, S. B. & Brower, L. P. Monarch butterfly (Danaus plexippus) thermoregulatory behavior and adaptations for overwintering in Mexico. Ecology 69, 458–467 (1988).Article 

    Google Scholar 
    Kammer, A. E. Thoracic temperature, shivering, and flight in the monarch butterfly, Danaus plexippus (L.). Z. Vgl. Physiol. 68, 334–344 (1970).Article 

    Google Scholar 
    Pendar, H. & Socha, J. J. Estimation of instantaneous gas exchange in flow-through respirometry systems: A modern revision of bartholomew’s ztransform method. PLoS One 10, e0139508 (2015).PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Lighton, J. R. B. Measuring Metabolic Rates: A Manual for Scientists. (Oxford University Press, 2008).Alonso-Mejía, A., Rendon-Salinas, E., Montesinos-Patiño, E. & Brower, L. P. Use of lipid reserves by monarch butterflies overwintering in Mexico: Implications for conservation. Ecol. Appl. 7, 934–947 (1997).Article 

    Google Scholar 
    Diaz, R., Overholt, W. A., Hahn, D. & Samayoa, A. C. Diapause induction in Gratiana boliviana (Coleoptera: Chrysomelidae), a biological control agent of tropical soda apple in Florida. Ann. Entomol. Soc. Am. 104, 1319–1326 (2011).Article 

    Google Scholar 
    Tschinkel, W. R. Sociometry and sociogenesis of colonies of the fire ant Solenopsis invicta during one annual cycle. Ecol. Monogr. 63, 425–457 (1993).Article 

    Google Scholar 
    Fink, L. S. & Brower, L. P. Birds can overcome the cardenolide defence of monarch butterflies in Mexico. Nature 291, 67–70 (1981).CAS 
    Article 

    Google Scholar 
    Ali, J. G. & Agrawal, A. A. Trade-offs and tritrophic consequences of host shifts in specialized root herbivores. Funct. Ecol. 31, 153–160 (2017).Article 

    Google Scholar 
    Woodson, R. E. The North American Species of Asclepias L. Ann. Mo. Bot. Gard. 41, 1 (1954).Article 

    Google Scholar 
    NRCS USDA. The PLANTS Database. National Plant Data Center. http://plants.usda.gov (2006).Agrawal, A. A., Salminen, J. P. & Fishbein, M. Phylogenetic trends in phenolic metabolism of milkweeds (Asclepias): Evidence for escalation. Evolution (N. Y). 63, 663–673 (2009).CAS 

    Google Scholar 
    Pocius, V. M. et al. Monarch butterflies show differential utilization of nine midwestern milkweed species. Front. Ecol. Evol. 6, 169 (2018).Pocius, V. M., Debinski, D. M., Pleasants, J. M., Bidne, K. G. & Hellmich, R. L. Monarch butterflies do not place all of their eggs in one basket: Oviposition on nine Midwestern milkweed species. Ecosphere 9, e02064 (2018).Article 

    Google Scholar 
    Ladner, D. T. & Altizer, S. Oviposition preference and larval performance of North American monarch butterflies on four Asclepias species. Entomol. Exp. Appl. 116, 9–20 (2005).Article 

    Google Scholar 
    Borders, B. A guide to the native milkweeds of Oregon. Xerces Soc. Invertebr. Conserv. www.xerces.org, 5, 12-23 (2012). More

  • in

    Intolerant baboons avoid observer proximity, creating biased inter-individual association patterns

    All research methods included in this study were performed in accordance with the relevant guidelines and regulations, under ZA/LP/81996 research permit, with ethical approval from the Animal Welfare Ethical Review Board (AWERB) at Durham University. The authors confirm the study was carried out in compliance with ARRIVE guidelines.All inter-individual association data was collected between June 2018 and June 2019 on a wild habituated group of Afro-montane chacma baboons in the western Soutpansberg Mountains, South Africa (central coordinates S29.44031°, E23.02217°) (for study site description see2). The study group was habituated circa 2005 and was the focus of intermittent research attention until 2014. The study area experienced long-term anthropogenic activities (local farming, forestry, and residences) prior to 2005, as such, consistent interactions with humans have been ongoing with this population for some time. From 2007 onwards numerous researchers were able to collect expansive datasets on the study group (e.g. Refs.17,18), indicating that habituation was at a typical level found elsewhere (also validated by AA and RH, who had researched chacma baboons elsewhere). From 2014 the group received full day (dawn until dusk) follows 3–4 days a week, with occasional gaps of up to 5 weeks in duration. These gaps did not appear to effect habituation levels, likely due to the presence of other researchers at the field site who always tried to act benignly when encountering the habituated group. The follow schedule was designed so that the study group retained as much of their natural interactions with predators as possible by ensuring the baboons spent significant time without observers who may influence the frequency and nature of predator–prey interactions19.The study site was located in a private nature reserve and the study group was not hunted during observation gaps or engaged in any conflict with humans, other than occasionally being scared (chasing, yelling, throwing stones etc.) from a small plantation by local workers, usually resulting in alarm barks and fleeing responses. However, the study group appeared adept at recognising the differences between researchers and these threats20. The majority of the study group’s home-range typically overlapped with the core area of the Lajuma Research Centre, and as a result, interactions with staff living in the area, unfamiliar researchers, and tourists were frequent. However, the baboons had not engaged in ‘raiding’ residences, threatening humans, or any other potentially negative symptom of habituation before the end of this study.Sampling methodology for proximity associations30-s focal sampling was used to collect proximity associations between all group members (excluding infants). All data was collected between June and December 2018 and January and June 2019; the majority of 2018s data was collected during the wet season, whilst most of 2019s data was collected during the dry season. To account for time of day, each day was split into four time-periods that were seasonally adjusted ensuring each period accounted for 25% of the current day length. A randomly ordered list of individuals was produced for each day, the first individual identified from the top 15 (approx. 20% of group size) individuals on the list was sampled immediately. Individuals could only be sampled once per time period per day, and a maximum of twice total per day. All individuals received at least 14 focal observations per time period (56 total) across the study period (see below for how we handled uneven sampling for some individuals). A video camera was used by AA (the only observer to collect this data) to record all focal observations (Panasonic HC-W580 Camcorder). At the end of the 30-s focal observation the identities of all neighbouring conspecifics within 5 m, 2.5 m, 1 m, and touching the focal animal were recorded (audibly by AA). We chose the end of the focal observation to record this data as this was most likely to reflect the conditions during the focal, i.e., the observer had been in proximity for at least 30 s.Neighbour information was extracted from video footage and entered manually by AA and AW. Data was split into separate years to reflect an observation gap of several weeks and to understand whether there was consistency in the hypothesized effects through time and to reflect underlying differences in environmental conditions during the two study periods; during the dry season fruits and seeds are scarce and day lengths are several hours shorter than in the wet season such that day journey lengths are often shorter in the dry season and animals are much more sedentary which could impact inter-individual spacings. In 2018 each individual was sampled between 28 and 30 times; 28 focals were randomly selected from each individual to make sampling even. For 2019 there were between 25 and 27 focals per individual; 25 samples of each individual were randomly selected. Observations were undertaken at a range of distances. For both years the median end observer distance was 4.5 m; data was thus split into close focal observations of less than or equal to 4.5 m (2018: n = 918, 2019: n = 809), and observations greater than 4.5 m (2018: n = 902 2019: n = 816). See supporting information Table S1 for summary statistics of the observation distances of each individual.We did not make any attempt to record our focal data evenly across the various habitats at our field site (see Supporting information text S1 for complete habitat descriptions) as our previous research indicated there was little difference in general spatial cohesion/inter-individual proximity patterns across these habitats (see Supporting information text S2 and Table S2). As a result, we considered it unlikely that there were fundamental differences in inter-individual association patterns across habitats, or that observers struggled to reliably detect or identify neighbours in dense habitats. We do acknowledge, however, that there will always be an element of bias with such methods, as observations were avoided, aborted, or excluded if visual obstructions (e.g., cliffs, rocks, walls, buildings, very dense vegetation etc.) prohibited accurate assessments; the observations used in the current study are from occasions when these factors were not an issue.During this study the group contained between 85 and 92 individuals. Age-sex class was defined according to secondary sexual characteristics (e.g., testes descending/enlarging, sexual swelling, canine eruption) and changes in pelage throughout juvenile development (see Supporting information text S3 for full descriptions). All 65 non-infant individuals that were present during 2017 (when displacement tolerances were calculated) and still remaining in the group by the end of 2019 were used in this study (4 individuals from the prior FID study were no longer present). There were a high number of births between 2018 and 2019, but none were independent by the time either of our sampling periods begun in 2018 or 2019. There was no immigration of foreign individuals, but two individuals disappeared, both during the 2018 focal sampling period. As a result, we had a very consistent pool of individuals to sample from during this study. We removed all data associated with the two individuals who disappeared as their occurrences as neighbours would have been poorly sampled (due to missing more than half the study) relative to the rest of the group which would have led to statistical biases21.Flight initiation distance procedureIndividual displacement tolerance estimates were previously quantified in our previous research2 using a flight initiation distance (FID) procedure22 that was completed between October 2017 and April 2018, prior and independent to the commencement of proximity association focal sampling in June 2018. Individual baboons were approached by an observer, and the distance at which the animal displaced away from the observer measured (see Supporting information Table S2 for summary statistics). This procedure was repeated 24 times for each individual baboon, with approaches spread evenly across two observers differing in familiarity. At the beginning of each approach we also recorded several behavioural, social, and environmental factors that could have hypothetically influenced an individual’s FID2 including whether the animal was engaged (e.g., digging or grooming) or not engaged (e.g., resting, chewing food, being groomed), habitat type (open/closed: see Supporting text S1), whether the animal was on the ground or sat on a low branch or rock within 50 cm of the ground, the number of conspecifics within 5 m of the focal animal, and whether there had been any external events within the preceding 5 min (e.g., alarm calls, aggressions, encountering another group or predator). During the approach, we also recorded the visual orientation distance (the distance at which the focal animal directed its line of vision towards the head of the approaching observer) and whether one of the focal animal’s neighbours had displaced/fled before the focal animal. Although all but neighbour flee first and external events showed some importance for predicting looking (see Table S4), FID was found to be distinct amongst individuals and repeatable within each individual, evidence that displacement tolerance may be an individual level trait2. Full details of methods, statistical analysis, and results (including comparison to the original model) for this updated model are in Supporting information text S4, with model summary results for the previous and updated models in Tables S3 and S4.The notion of an observer approaching a habituated primate may be considered atypical or likely to result in habituation/sensitization effects or agonistic behaviours being directed towards the approaching observers. However, our previous study2 showed that almost all approaches resulted in the animal passively relocating (98.85%), a very benign response identical to the behaviours of subordinate baboons displacing away from dominant conspecifics. This suggests that in this group, observers may be considered equivalent to a high-level social threat2. Throughout observation periods on habituated animals, observers are likely to approach or displace animals either incidentally or accidentally multiple times throughout the day, especially during lengthy focal observations. As such, the approach methodology is unlikely to represent a stimulus outside of the norm for our study animals. This may explain why displacement responses were so passive and why there was no evidence of habituation or sensitization effects across the group or individually through a range of temporal periods2 or after life-threatening events20. As a result, our situation was possible without risk of causing stress or anxiety in the study subjects, eliciting agonistic behaviours towards observers, or interfering with their prior habituation levels.Statistical analysisInfluence of tolerance and observer distance on inter-individual association patternsQuantifying displacement toleranceTo quantify displacement tolerance towards observers we extracted the individual conditional modes from the updated FID model using the ranef function in brms. Conditional modes are often referred to as Best Linear Unbiased Predictors (BLUPs) and are the difference between the predicted mean population-level response for a given set of treatments (i.e., population-level effects) and the predicted responses for each individual, and therefore infer the extent to which each individual differs from the population mean. The conditional modes and their associated standard deviations can be found in supporting information Table S5.To validate that the conditional modes from the updated model were both representative of the individual’s flight responses and in line with the estimates produced from our previous study2 we performed additional tests. Firstly, we performed a Pearson’s correlation between the conditional modes from the updated model and the conditional modes from the previous article. Individual tolerance estimates were consistent (r(63) = 0.915, p  More