Aim The aim of this study is to model the past, current and future distribution of J. phoenicea s.s., J. turbinata and J. canariensis, based on bioclimatic variables using a maximum entropy model (MaxEnt) in the Mediterranean and Macaronesian regions. Location Mediterranean and Macaronesian Taxon Cupressaceae, Juniperus Methods Data on the occurrence of the J. phoenicea complex was obtained from the GBIF, the literature, herbaria, and the authors’ field notes. The bioclimatic variables were obtained from the WorldClim database (http://worldclim.org/) and Paleoclim (http://www.paleoclim.org/). The climate data related to species localities were used for predictions of niches by implementation of MaxEnt and we evaluated the model with ENMeval. Results The potential niches of Juniperus phoenicea during the LIG, LGM and MH covered 30%, 10% and almost 100%, respectively, of the current potential niche. Climate warming could reduce potential niches by 30% and 90% in scenarios RCP2.6 and RCP8.5, respectively. The potential niches of Juniperus turbinata had a broad circum-Mediterranean and Canarian distribution during the LIG and the MH, extending its distribution during the LGM when it was found in more areas than at present; the predicted warming in scenario RCP2.6 and RCP8.5 could reduce the current potential niche by 30% and 50%, respectively. The model did not find suitable niches for J. canariensis during the LIG and the LGM, but during the MH its potential niche was 30% larger than at present. The climate warming scenario RCP2.6 indicates a reduction of the potential niche by 30%, while RCP8.5 does so by almost 60%. Main conclusions This research can provide information to increase the protection of the juniper forest and to try to counteract the phenomenon of local extinctions caused by anthropic pressure and climate changes.

The impact of rising global temperatures on survival and reproduction is putting many species at risk of extinction. In particular, it has recently been shown that thermal effects on reproduction, especially limits to male fertility, can underpin species distributions in insects. However, the physiological factors influencing fertility at high temperatures are poorly understood. Key factors that affect somatic thermal tolerance such as hardening, the ability to phenotypically increase thermal tolerance after a mild heat shock, and the differential impact of temperature on different life stages, are largely unexplored for thermal fertility tolerance. Here, we examine the impact of high temperatures on male fertility in the cosmopolitan fruit fly Drosophila virilis. We first determined whether temperature stress at either the pupal or adult life-history stage impacts fertility. We then tested the capacity for heat-hardening to mitigate heat-induced sterility. We found that thermal stress reduces fertility in different ways in pupae and adults. Pupal heat stress delays sexual maturity, whereas males heated as adults can reproduce initially following heat stress, but lose the ability to produce offspring. We also found evidence that while heat-hardening in D. virilis can improve high temperature survival, there is no significant protective impact of this same hardening treatment on fertility. These results suggest that males may be unable to prevent the costs of high temperature stress on fertility through heat-hardening which limits a species' ability to quickly and effectively reduce fertility loss in the face of short-term high temperature events.

Time series data are often observed in ecological monitoring. Frequently such data exhibit nonlinear trends over time potentially due to complex relationships between observed and auxiliary variables, and there may also be sudden declines over time due to major disturbances. This poses substantial challenges for modelling such data and also for model-based adaptive monitoring. We propose novel methods for finding adaptive designs for monitoring when historical data show such nonlinear patterns and sudden declines over time. This work is motivated by a coral reef monitoring program that has been established at Scott Reef; a coral reef off the Western coast of Australia. Data collected for monitoring the health of Scott Reef are considered, and semiparametric and interrupted time series modelling approaches are adopted to describe how these data vary over time. New methods are then proposed that enable adaptive monitoring designs to be found based on such modelling approaches. These methods are then applied to find future monitoring designs at Scott Reef and form a set of recommendations for future monitoring. Through applying the proposed methods, it was found that future information gain is expected to be similar across a variety of different sites, suggesting that no particular location needed to be prioritised at Scott Reef for the next monitoring phase. In addition, it was found that omitting some sampling sites/reef locations was possible without substantial loss in expected information gain, depending upon the disturbances that were observed. The resulting adaptive designs are used to provide recommendations for future monitoring in this region, and for reefs where changes to the current monitoring practices are being sought. Furthermore, as the methods used and developed throughout this study are generic in nature, this research has the potential to improve ecological monitoring more broadly where complex data are being collected over time.

Rising ocean temperatures are threatening marine species and populations worldwide, and ectothermic taxa are particularly vulnerable to warming. Echinoderms are an ecologically important phylum of marine ectotherms and shifts in their population dynamics can have profound impacts on the marine environment. The effects of warming on echinoderms are highly variable across controlled laboratory-based studies. Accordingly, synthesis of these studies will facilitate the better understanding of broad patterns in responses of echinoderms to ocean warming. Herein, a meta-analysis incorporating the results of 85 studies (710 individual responses) is presented, exploring the effects of warming on various performance predictors. The mean responses of echinoderms to all magnitudes of warming were compared across multiple biological responses, ontogenetic life stages, taxonomic classes, and regions, facilitated by multivariate linear mixed effects models. Further models were conducted which only incorporated responses to warming greater than the projected end-of-century mean annual temperatures at the collection sites. This meta-analysis provides evidence that ocean warming will generally accelerate metabolic rate (+ 32%) and reduce survival (– 35%) in echinoderms, and echinoderms from sub-tropical (– 9%) and tropical (– 8%) regions will be the most vulnerable. The relatively high vulnerability of echinoderm larvae to warming (– 20%) indicates that this life stage may be a significant developmental bottleneck in the near-future, likely reducing successful recruitment into populations. Furthermore, asteroids appear to be the class of echinoderms that are most negatively affected by elevated temperature (– 30%). When considering only responses to magnitudes of warming representative of end-of-century climate change projections, the negative impacts on asteroids, tropical species and juveniles were exacerbated (– 51%, – 34%, and – 40% respectively). The results of these analyses enable better predictions of how keystone and invasive echinoderm species may perform in a warmer ocean, and the possible consequences for populations, communities, and ecosystems.

Earth systems are nearing a global tipping point, beyond which the dynamics of biological systems will become unstable. One major driver of instability is species invasion, especially by organisms that act as “ecosystem engineers” through their modification of abiotic and biotic factors. In a mosaic landscape of non-invaded and invaded habitat, ecosystems modified through invasion may serve as “sink” habitat. To understand how native organisms respond to habitat that is becoming increasingly modified, it is essential to examine biological communities within invaded and non-invaded habitat, identifying compositional shifts between native and non-native taxa as well as measuring how modification has affected interactions among community members. Using dietary metabarcoding, our study examines the response of a native Hawaiian generalist predator to habitat modification by comparing biotic interactions across metapopulations of spiders collected in native forest and sites invaded by kahili ginger. Our study shows that, although there are shared components of the dietary community, spiders in invaded habitat are eating a less consistent and more diverse diet consisting of more non-native arthropods which are rarely or entirely undetected in spiders collected from native forest. Additionally, the frequency of novel interactions with parasites was significantly higher in invaded sites, reflected by the frequency and diversity of non-native Hymenoptera parasites and entomopathogenic fungi. The study highlights the role of habitat modification driven by an invasive plant in altering community structure and biotic interactions, appearing to serve as a “sink” for native arthropods and thereby threatening the stability of the ecosystem.

To predict suitable growing regions for Leonurus japonicus and to provide scientific sopport for the habitat conservation and the exploitation and utilization of germplasm resources under climate change conditions, this study combined niche and priority conservation models to assess the future potential distribution of L. japonicus in China. To this end, distribution points and samples of L. japonicus were gathered through online and field surveys. The Maxent model with optimized parameters was used for predicting the suitable habitats of L. japonicus at different stages, and the Marxan model was used to determine the priority of protected areas. The results showed that the highest temperature in the hottest month, the lowest temperature in the coldest month, the precipitation in the wettest month, the precipitation in the driest month, and altitude were the main environmental factors influencing the distribution of L. japonicus. Under the three climate change scenarios, the centroid of the suitable area of L. japonicus migrated northward, and the migration position tended to expand further northwest. In the future, there would be no significant niche differentiation of L. japonicus; the Marxan results showed that priority protected areas for L. japonicus were in southwestern central China, Lingnan, southern east China, and Guizhou. Overall, the results of this research can provide a strategy for the determination of priority protection areas for Leonurus japonicus in China.

Seasonal dietary shifts of animals are important ecological adaptation strategies. An increasing number of studies have shown that seasonal dietary shifts can influence or even determine the composition of gut microbiota. The turpan wonder gecko Teratoscincus roborowskii lives in extreme desert environments, which have flexible dietary shift to fruit-eating in warm seasons. But the impact of such shifts on gut microbiota is poorly understood. Here, 16SrRNA sequencing and LC-MS metabolomics we used to examine the changes of gut microbiota composition and metabolic pattern of T. roborowskii. The results demonstrated that the gut microbes of T. roborowskii had significant seasonal changes, the diversity and abundance of gut microbiota in autumn were higher than those of in spring. Firmicutes, Bacteroidetes and Proteobacteria were the core gut microbes of T. roborowskii. Verrucomicrobia and Proteobacteria exhibit dynamic pattern of ebb and flow between spring and autumn.The composition and structure of gut microbes in different seasons perform specific metabolic functions, and this change may be an important adaptation for T. roborowskii to cope with dietary shifts and improve energy acquisition. Our study will provide a theoretical basis for exploring the adaptive evolution to special frugivorous behavior of the T. roborowskii, which is an important supplement to the study of the gut microbiology of desert lizards.

1. The fixation index, FIS has been a staple measure to detect selection or departures from random mating in populations. However, current Next Generation Sequencing (NGS) cannot easily estimate Fis, in multi-locus gene families, which contain multiple loci having similar or identical arrays of variant sequences of ≥1 kilobase, which differ at multiple positions. In these families, high-quality short-read NGS data typically identify variants, but not the genomic location, which is required to calculate Fis (based on locus-specific observed and expected heterozygosity). Thus, to assess assortative mating, or selection on heterozygotes, from NGS of multi-locus gene families, we need a method that does not require knowledge of which variants are allelic at which locus in the genome. 2. We developed such a method. Like Fis, our novel measure, 1His, is based on the principle that positive assortative mating, or selection against heterozygotes, reduces within-individual variability relative to the population. 3. We demonstrate high accuracy of 1His on a wide-range of simulated scenarios, and two datasets from natural populations of penguins and dolphins. 4. 1His is important because multi-locus gene families are often involved in assortative mating, or selection on heterozygotes. 1His is particularly useful for multi-locus gene families such as toll-like receptors, the major-histocompatibility-complex in animals, homeobox genes in fungi and self-incompatibility genes in plants.

Plant and soil biodiversity can have significant effects on herbivore resistance mediated by plant metabolites. Here, we disentangled the independent effects of plant diversity and soil legacy on constitutive and herbivore-induced plant metabolomes of three plant species in two complementary microcosm experiments. First, we grew plants in sterile soil with three different plant diversity levels. Second, single plant species were grown on soil with different plant diversity-induced soil legacies. We infested a subset of all plants with Spodoptera exigua larvae, a generalist leaf-chewing herbivore, and assessed foliar and root metabolomes. Neither plant diversity nor soil legacy had significant effects on overall foliar, root, or herbivore-induced metabolome composition. Herbivore-induced metabolomes, however, differed from those of control plants. We also detected 139 significantly regulated metabolites by comparing plants grown in monocultures with conspecifics growing in plant or soil legacy mixtures. Moreover, plant-plant and plant-soil interactions regulated 141 metabolites in herbivore-induced plants. Taken together, plant diversity and soil legacy independently alter the concentration and induction of plant metabolites, thus affecting the plant's defensive capability. This is a first step towards disentangling plant and soil biodiversity effects on herbivore resistance, thereby improving our understanding of the mechanisms that govern ecosystem functioning.

1. High latitude ecosystems are experiencing the most rapid warming on earth, expected to trigger a diverse array of ecological responses. Climate warming affects the ecophysiology of fish, and fish close to the cold end of their thermal distribution are expected to increase somatic growth from increased temperatures and a prolonged growth season, which in turn affects maturation schedules, reproduction and survival, boosting population growth. Accordingly, fish species living in ecosystems close to their northern range edge should increase in numerical importance and possibly displace cold-water adapted species. 2. We aim to document if and how population level effects of warming mediated by individual level responses to increased temperatures, shift community structure and composition in high latitude ecosystems. 3. We studied 11 cool-water adapted freshwater fish populations in communities dominated by cold-water adapted species to investigate changes in the relative importance of cool-water fish during the last 30 years of rapid warming in high latitude lakes. In addition, we studied the individual level responses to warming to clarify the potential mechanisms underlying the population effects. 4. Our long-term series‘ (1991-2020) reveal a marked increase in numerical importance of the cool-water fish species, perch, in ten out of eleven populations, and in most fish communities the cool-water species is now dominant. Moreover, we show that climate warming affects population level processes via direct and indirect temperature effects on the individuals. Specifically, the increase in abundance arises from increased survival of 0+ individuals, faster juvenile growth and ensuing earlier maturation, all boosted by climate warming. 5. The speed and magnitude of the response to warming in these high latitude fish communities strongly suggest that cold-water fish will be displaced by fish adapted to warmer water. Consequently, management should focus on climate adaptation limiting future introductions and invasions of cool-water fish and mitigating harvesting pressure on cold-water fish.

With many plant-pollinator interactions undergoing change, we require a better understanding of how the addition of new interacting partners, such as antagonists, can affect plant reproduction. One such group of antagonistic floral visitors, nectar robbers, can deplete plants of nectar rewards without contributing to pollination. The addition of nectar robbing to the floral visitor assemblage could therefore have costs to the plant´s reproductive output. We focus on a recent plant colonist, Digitalis purpurea, a plant that in its native range is rarely robbed, but experiences intense nectar robbing in areas it has been introduced to. Here, we test the costs to reproduction following experimental nectar robbing. To identify any changes in the behaviour of the principal pollinators in response to nectar robbing, we measured visitation rates, visit duration, proportion of flowers visited and rate of rejection of inflorescences. To find the effects of robbing on fitness, we used proxies for female and male components of reproductive output, by measuring the seeds produced per fruit and the pollen export respectively. Nectar robbing significantly reduced the rate of visitation and lengths of visits by bumblebees. Additionally, bumblebees visited a lower proportion of flowers on an inflorescence that had robbed flowers. We found that flowers in the robbed treatment produced significantly fewer seeds per fruit on average but did not export fewer pollen grains. Our finding that robbing leads to reduced seed production could be due to fewer and shorter visits to flowers leading to less effective pollination. We discuss the potential consequences of new pollinator environments, such as exposure to nectar robbing, for plant reproduction.

Lacking of systematic evaluations in soil quality and microbial community recovery after different amendments addition limits optimization of amendments combination in coal mine-soils. We performed a short-term incubation experiment over 12 weeks to assess the effects of three amendments (biochar: C; nitrogen fertilizer at three levels: N-N1~N3; microbial agent at two levels: M-M1~M2) based on C/N ratio (regulated by biochar and N level: 35:1, 25:1, 12.5:1) on soil quality and microbial community in the Qilian Mountains, China. Over the incubation period, soil pH and MBC/MBN were significantly lower than unamended treatment in N addition and C+M+N treatments, respectively. Soil organic carbon (SOC), total nitrogen (TN), available nitrogen (AN), available phosphorus (AP), available potassium (AK), microbial biomass carbon (MBC) and microbial biomass nitrogen (MBN) contents had a significant increase in all amended treatments (P<0.001). Higher AP, AK, MBC, MBN and lower MBC/MBN were observed in N2-treated soil(corresponding to C/N ratio of 25:1). Meanwhile, N2-treated soil significantly increased species richness and diversity of soil bacterial community (P<0.05). Principal coordinate analysis further showed that soil bacterial community compositions were significantly separated by N level. C-M-N treatments (especially at N2 and N1 levels) significantly increased the relative abundance (>1%) of the bacterial phyla Bacteroidetes and Firmicutes, and decreased the relative abundance of fungal phyla Chytridiomycota (P<0.05). Redundancy analysis illustrated the importance of soil nutrients in explaining variability in bacteria community composition (74.73%) than fungal (35.0%). Our results indicated that N and M addition based on biochar can improve soil quality by neutralizing soil pH and increasing soil nutrient contents, and the appropriate C/N ratio (25:1: biochar+N2-treated soil) can better promote mass, richness and diversity of soil bacterial community. Our study provided a new insight for achieving restoration of damaged habitats by changing microbial structure, diversity and mass by regulating C/N ratio of amendments

Target-site insensitive mutations and overexpression of detoxification genes are two major mechanisms conferring insecticide resistance. Many molecular assays were applied to detect these two kinds of resistance genetic markers in insect populations. Unfortunately, these assays are time-consuming and have high false-positive rates. RNA-Seq data, which contains information on the variation within expressed regions of the genome and expression information of detoxification genes, provides us a valuable resource to detect resistance-associated markers. At present, there is no corresponding method at present. Here, we collected 66 reported resistance mutations of four main insecticide targets (AChE, VGSC, RyR, and nAChR) of 82 insect species. Next, we obtained 403 sequences of the four target genes and 12,665 sequences of three kinds of detoxification genes including P450, GST, and CCE. Here, we developed a Perl program, FastD, to detect insecticide target-site insensitive mutations and overexpressed detoxification genes from RNA-Seq data, and constructed a web server for FastD (http://www.insect-genome.com/fastd). FastD program was then applied to detect two kinds of resistant markers in five populations of two insects, Plutella xylostella and Aphis gossypii. Results showed that RyR mutation G4946E was detected in all P. xylostella populations, with higher frequencies in two resistant populations, ZZ (66.1%) and CHR (94.55%), than a susceptible population CHS (2.32%). CYP6a2 was overexpressed 10.82-fold in ZZ population. As to A. gossypii, nAChR mutation R81T was detected in resistant population KR with 49.85% frequency, but not in susceptible population NS. CYP6CY22 and CYP6CY13 were overexpressed 39.61- and 22.04-fold respectively in KR population. FastD is a program using RNA-Seq data to detect two types of resistance markers to estimate resistance level of insect populations. Generally, resistance level estimated by FastD were consistent with previous reports, confirming the reliability of this program in predicting population resistance at omics-level.

Resource fluctuation is a major driver of animal movement, influencing strategic choices such as residency vs nomadism, or social dynamics. The Arctic tundra is characterized by strong seasonality: resources are abundant during the short summers but scarce in winters. Therefore, expansion of boreal-forest species onto the tundra raises questions on how they cope with winter-resource scarcity. We examined a recent incursion by red foxes (Vulpes vulpes) onto the coastal tundra of western Hudson Bay, an area historically occupied by Arctic foxes (Vulpes lagopus) that lacks access to anthropogenic foods, and compared seasonal shifts in space use of the two species. We used 4 years of telemetry data following 8 red foxes and 11 Arctic foxes to test the hypothesis that the movement strategies of both species are primarily driven by temporal variability of resources. We also predicted that the harsh tundra conditions in winter affect red foxes more than Arctic foxes, which are adapted to this environment. Dispersal was the most frequent winter movement strategy in both fox species, despite its association with high mortality (winter mortality was 9.4 times higher in dispersers than residents). Red foxes consistently dispersed towards the boreal forest, whereas Arctic foxes primarily used sea ice to disperse. Home range size of red and Arctic foxes did not differ in summer, but resident red foxes substantially increased their home range size in winter, whereas home range size of resident Arctic foxes did not change seasonally. As climate changes, abiotic constraints on some species may relax, but associated declines in prey communities may lead to local extirpation of many predators, notably by favoring dispersal during resource scarcity.

Comparing genome scans among species is a powerful approach for investigating the patterns left by evolutionary processes. In particular, this offers a way to detect candidate genes that drive convergent evolution. We compared genome scan results to investigate if patterns of genetic diversity and divergence are shared among divergent species within the stickleback order (Gasterosteiformes): the threespine stickleback (Gasterosteus aculeatus), ninespine stickleback (Pungitius pungitus) and tubesnout (Aulorhynchus flavidus). Populations were sampled from the southern and northern edges of each species’ range, to identify patterns associated with latitudinal changes in genetic diversity. Weak correlations in genetic diversity (FST and expected heterozygosity) and three different patterns in the genomic landscape were found among these species. Additionally, no candidate genes for convergent evolution were detected. This is a counterexample to the growing number of studies that have shown overlapping genetic patterns, demonstrating that genome scan comparisons can be noisy due to the effects of several interacting evolutionary forces.

A new coccidian species, Isospora elliotae n. sp., from the Australian magpie Gymnorhina tibicen (Latham, 1801) in Western Australia is described and characterised morphologically and molecularly. Microscopic analysis of a faecal sample identified subspheroidal oöcysts (n = 20), 20–22 × 18–20 (20.7 × 18.7); length/width (L/W) ratio 1.05–1.14 (1.10). Wall bi-layered, 1.0–1.3 (1.2) thick, outer layer smooth, c.2/3 of total thickness. Micropyle and oocyst residuum absent, but usually two polar granules are present. Sporocysts (n = 28) ovoidal, 12–13 × 9–11 (12.6 × 9.7); L/W ratio 1.22–1.35 (1.30). Stieda body present, flattened to half-moon-shaped, c. 0.5 deep × 2.0 wide; sub-Stieda indistinct or barely discernible, c. 1.0 deep × 2.5 wide; para-Stieda body absent; sporocyst residuum present, composed of granules dispersed among the sporozoites. Sporozoites vermiform, with anterior and posterior refractile bodies and nucleus. Segments of three gene loci (18S rRNA, 28S rRNA and COI) were sequenced and I. elliotae n. sp. exhibited 99.8% genetic similarity to Isospora sp. MAH-2013a (KF648870) followed by 99.7% genetic similarity to Isospora neochmiae Yang, Brice & Ryan, 2016 (KT224380) at the 18S rRNA gene locus. It shared 97.0% genetic similarity with an unnamed Isospora sp. (AY283852) at the 28S rRNA gene locus and it also shared the highest genetic similarity of 99.8% with the unnamed Isospora sp. from an American crow (OL999120) at the COI gene locus. Based on morphological and molecular data, this isolate is a new species named as I. elliotae n. sp.

1. Wetlands belong to the globally most threatened habitats, and organisms depending on them are of conservation concern. Wetland destruction and quality loss may affect negatively also boreal breeding ducks in which habitat selection often needs balancing between important determinants of habitat suitability. In Finland duck population trajectories are habitat-specific, while the reasons behind are not known. 2. In this research, the balance of nest predation risk and invertebrate food abundance in boreal breeding ducks was studied in Finland at 45 lakes and ponds in 2017 and 2018. Nest predation experiments were conducted with artificial nests followed by wildlife cameras during seven days. Invertebrates were sampled from the study water bodies using emergence and activity traps. Duck pairs and broods were also surveyed from these and 18 additional water bodies. 3. The wildlife camera results indicate that predation risk was higher in the water bodies surrounded by agricultural land than forestland. Ponds (seasonal, beaver and man-made) had lower nest predation risk and they were also more invertebrate-rich habitats than permanent lakes. In addition, artificial nests further away from water bodies had higher survival than shoreline nests. Habitat use of duck pairs (prior to nesting) was not associated with invertebrate food, but duck broods preferred habitats rich in food. 4. High nest predation pressure in shorelines of especially agricultural landscapes may contribute the declining population trends of ducks in Finland. Controlling predators would be an important conservation action to improve duck breeding success. This research underlines the benefits of the availability of different water body types for the breeding ducks. There is an urgent need to pay attention to protecting seasonal ponds, while the lack of flooded waters may be mitigated by favouring beavers or man-made ponds.

Understanding population change across long time scales and at fine spatiotemporal resolutions is important for confronting a broad suite of conservation challenges. However, this task is hampered by a lack of quality long-term census data for multiple species collected across large geographic regions. Here, we used century-long (1919-2018) data from the Audubon Christmas Bird Count (CBC) survey to assess population changes in over 300 avian species in North America and evaluate their temporal non-stationarity. To estimate population sizes across the entire century, we employed a Bayesian hierarchical model that accounts for species detection probabilities, variable sampling effort, and missing data. We evaluated population trends using generalized additive models (GAMs) and assessed temporal non-stationarity in the rate of population change by extracting the first derivatives from the fitted GAM functions. We then summarized the population dynamics across species, space, and time using a non-parametric clustering algorithm that categorized individual population trends into four distinct trend clusters. We found that species varied widely in their population trajectories, with over 90% of species showing a considerable degree of spatial and/or temporal non-stationarity, and many showing strong shifts in the direction and magnitude of population trends throughout the past century. Species were roughly equally distributed across the four clusters of population trajectories, though grassland, forest, and desert specialists more commonly showed declining trends. Interestingly, for many species, region-wide population trends often differed from those observed at individual sites, suggesting that conservation decisions need to be tailored to fine spatial scales. Together, our results highlight the importance of considering spatial and temporal non-stationarity when assessing long-term population changes. More generally, we demonstrate the promise of novel statistical techniques for improving the utility and extending the temporal scope of existing citizen science datasets.