Scale gaps in landscape phenology: challenges and opportunities

Pursuit and escape drive fine-scale movement variation during migration in a temperate alpine ungulate

Climate change reduces snowpack, advances snowmelt phenology, drives summer warming, alters growing season precipitation regimes, and consequently modifies vegetation phenology in mountain systems. Elevational migrants track spatial variation in seasonal plant growth by moving between ranges at different elevations during spring, so climate-driven vegetation change may disrupt historic benefits of migration. Elevational migrants can furthermore cope with short-term environmental variability by undertaking brief vertical movements to refugia when sudden adverse conditions arise. We uncover drivers of fine-scale vertical movement variation during upland migration in an endangered alpine specialist, Sierra Nevada bighorn sheep (Ovis canadensis sierrae) using a 20-year study of GPS collar data collected from 311 unique individuals. We used integrated step-selection analysis to determine factors that promote vertical movements and drive selection of destinations following vertical movements. Our results reveal that relatively high temperatures consistently drive uphill movements, while precipitation likely drives downhill movements. Furthermore, bighorn select destinations at their peak annual biomass and maximal time since snowmelt. These results indicate that although Sierra Nevada bighorn sheep seek out foraging opportunities related to landscape phenology, they compensate for short-term environmental stressors by undertaking brief up- and downslope vertical movements. Migrants may therefore be impacted by future warming and increased storm frequency or intensity, with shifts in annual migration timing, and fine-scale vertical movement responses to environmental variability.

Defining the roles of local precipitation and anthropogenic water sources in driving the abundance of Aedes aegypti, an emerging disease vector in urban, arid landscapes

Understanding drivers of disease vectors' population dynamics is a pressing challenge. For short-lived organisms like mosquitoes, landscape-scale models must account for their highly local and rapid life cycles. Aedes aegypti, a vector of multiple emerging diseases, has become abundant in desert population centers where water from precipitation could be a limiting factor. To explain this apparent paradox, we examined Ae. aegypti abundances at > 660 trapping locations per year for 3 years in the urbanized Maricopa County (metropolitan Phoenix), Arizona, USA. We created daily precipitation layers from weather station data using a kriging algorithm, and connected localized daily precipitation to numbers of mosquitoes trapped at each location on subsequent days. Precipitation events occurring in either of two critical developmental periods for mosquitoes were correlated to suppressed subsequent adult female presence and abundance. LASSO models supported these analyses for female presence but not abundance. Precipitation may explain 72% of Ae. aegypti presence and 90% of abundance, with anthropogenic water sources supporting mosquitoes during long, precipitation-free periods. The method of using kriging and weather station data may be generally applicable to the study of various ecological processes and patterns, and lead to insights into microclimates associated with a variety of organisms' life cycles.

Varying effects of tree cover on relationships between satellite-observed vegetation greenup date and spring temperature across Eurasian boreal forests

Earth observation satellites have facilitated the quantification of how vegetation phenology responds to climate warming on large scales. However, satellite image pixels may contain a mixture of multiple vegetation types or species with diverse phenological responses to climate variability. It is unclear how these mixed pixels affect the statistical relationships between satellite-derived vegetation phenology and climate factors. Here, we aim to investigate the impacts of percent tree cover (PTC), a measure of mixed pixel, on the statistical relationships between satellite-derived vegetation greenup date (GUD) and spring air temperature across Eurasian boreal forests at a 0.05° spatial resolution. We estimated GUD using Moderate Resolution Imaging Spectroradiometer (MODIS) time series data. The responses of GUD to interannual variation in spring temperature (April to May) during 2001-2020 were characterized by correlation coefficient (RTAM) and sensitivity (STAM). We then evaluated the local impacts of PTC on spatial variations in RTAM and STAM using partial correlation analysis through spatial moving windows. Our results indicate that, for most areas, forests with higher PTC were associated with stronger RTAM and STAM. Moreover, PTC had stronger local impacts on RTAM and STAM than mean annual temperature and temperature seasonality for 37.3% and 27.4% of the moving windows, respectively. These impacts were spatially varying and different among forest types. Specifically, deciduous broadleaf forests and deciduous needleleaf forests tend to have a higher proportion of these impacts compared to other forest types. Our findings demonstrate the nonnegligible effects of PTC on the statistical responses of GUD to temperature variability at coarse spatial resolution (0.05°) across Eurasian boreal forests.

Urban environments provide new perspectives for forecasting vegetation phenology responses under climate warming

Given that already-observed temperature increase within cities far exceeds the projected global temperature rise by the end of the century, urban environments often offer a unique opportunity for studying ecosystem response to future warming. However, the validity of thermal gradients in space serving as a substitute for those in time is rarely tested. Here, we investigated vegetation phenology dynamics in China's 343 cities and empirically test whether phenological responses to spatial temperature rise in urban settings can substitute for those to temporal temperature rise in their natural counterparts based on satellite-derived vegetation phenology and land surface temperature from 2003 to 2018. We found prevalent advancing spring phenology with "high confidence" and delaying autumn phenology with "medium confidence" under the context of widespread urban warming. Furthermore, we showed that space cannot substitute for time in predicting phenological shifts under climate warming at the national scale and for most cities. The thresholds of ~11°C mean annual temperature and ~600 mm annual precipitation differentiated the magnitude of phenological sensitivity to temperature across space and through time. Below those thresholds, there existed stronger advanced spring phenology and delayed autumn phenology across the spatial urbanization gradients than through time, and vice versa. Despite the complex and diverse relationships between phenological sensitivities across space and through time, we found that the directions of the temperature changes across spatial gradients were converged (i.e., mostly increased), but divergent through temporal gradients (i.e., increased or decreased without a predominant direction). Similarly, vegetation phenology changes more uniformly over space than through time. These results suggested that the urban environments provide a real-world condition to understand vegetation phenology response under future warming.

Citizen science data reveal regional heterogeneity in phenological response to climate in the large milkweed bug, Oncopeltus fasciatus

Regional populations of geographically widespread species may respond to different environmental factors across the species' range, generating divergent effects of climate change on life-history phenology. Using thousands of citizen science observations extracted from iNaturalist and associated with corresponding temperature, precipitation, elevation, and daylength information, we examined the drivers of adult mating and of nymphal phenology, development, and group size for populations of the large milkweed bug, Oncopeltus fasciatus, in different ecoregions. Research-grade iNaturalist images were correctly identified 98.3% of the time and yielded more than 3000 observations of nymphal groups and 1000 observations of mating adults spanning 18 years. Mating phenology showed distinct regional patterns, ranging from year-round mating in California to temporally restricted mating in the Great Lakes Northeastern Coast ecoregion. Relative temperature increases of 1°C for a given daylength expanded the mating season by more than a week in western ecoregions. While increases in relative temperature delayed mating phenology in all ecoregions, greater winter precipitation advanced mating in the California ecoregion. In the eastern ecoregions, nymphal phenology was delayed by increases in summer rainfall but was advanced by relative temperature increases, whereas in western regions, relative temperature increases delayed nymphal phenology. Furthermore, accumulated growing degree days (AGDD) was a poor predictor of developmental progression, as we found a positive but weak correlation between AGDD and age structure only for the Appalachian Southeast North America and the Great Lakes Northern Coast ecoregions. These complex phenological responses of O. fasciatus are just one example of how populations may be differentially susceptible to a diversity of climatic effects; using data across a species' whole distribution is critical for exposing regional variations, especially for species with large, continental-scale ranges. This study demonstrates the potential of photodocumented biodiversity data to aid in the monitoring of life history, host plant-insect interactions, and climate responsiveness.

Climate lags and genetics determine phenology in quaking aspen (Populus tremuloides)

Spatiotemporal patterns of phenology may be affected by mosaics of environmental and genetic variation. Environmental drivers may have temporally lagged impacts, but patterns and mechanisms remain poorly known. We combine multiple genomic, remotely sensed, and physically modeled datasets to determine the spatiotemporal patterns and drivers of canopy phenology in quaking aspen, a widespread clonal dioecious tree species with diploid and triploid cytotypes. We show that over 391 km² of southwestern Colorado: greenup date, greendown date, and growing season length vary by weeks and differ across sexes, cytotypes, and genotypes; phenology has high phenotypic plasticity and heritabilities of 31-61% (interquartile range); and snowmelt date, soil moisture, and air temperature predict phenology, at temporal lags of up to 3 yr. Our study shows that lagged environmental effects are needed to explain phenological variation and that the effect of cytotype on phenology is obscured by its correlation with topography. Phenological patterns are consistent with responses to multiyear accumulation of carbon deficit or hydraulic damage.

Weather anomalies more important than climate means in driving insect phenology

Studies of long-term trends in phenology often rely on climatic averages or accumulated heat, overlooking climate variability. Here we test the hypothesis that unusual weather conditions are critical in driving adult insect phenology. First, we generate phenological estimates for Lepidoptera (moths and butterflies) across the Eastern USA, and over a 70 year period, using natural history collections data. Next, we assemble a set of predictors, including the number of unusually warm and cold days prior to, and during, the adult flight period. We then use phylogenetically informed linear mixed effects models to evaluate effects of unusual weather events, climate context, species traits, and their interactions on flight onset, offset and duration. We find increasing numbers of both warm and cold days were strong effects, dramatically increasing flight duration. This strong effect on duration is likely driven by differential onset and termination dynamics. For flight onset, impact of unusual climate conditions is dependent on climatic context, but for flight cessation, more unusually cold days always lead to later termination particularly for multivoltine species. These results show that understanding phenological responses under global change must account for unusual weather events, especially given they are predicted to increase in frequency and severity.

Antecedent climatic conditions spanning several years influence multiple land-surface phenology events in semi-arid environments

Ecological processes are complex, often exhibiting non-linear, interactive, or hierarchical relationships. Furthermore, models identifying drivers of phenology are constrained by uncertainty regarding predictors, interactions across scales, and legacy impacts of prior climate conditions. Nonetheless, measuring and modeling ecosystem processes such as phenology remains critical for management of ecological systems and the social systems they support. We used random forest models to assess which combination of climate, location, edaphic, vegetation composition, and disturbance variables best predict several phenological responses in three dominant land cover types in the U.S. Northwestern Great Plains (NWP). We derived phenological measures from the 25-year series of AVHRR satellite data and characterized climatic predictors (i.e., multiple moisture and/or temperature based variables) over seasonal and annual timeframes within the current year and up to 4 years prior. We found that antecedent conditions, from seasons to years before the current, were strongly associated with phenological measures, apparently mediating the responses of communities to current-year conditions. For example, at least one measure of antecedent-moisture availability [precipitation or vapor pressure deficit (VPD)] over multiple years was a key predictor of all productivity measures. Variables including longer-term lags or prior year sums, such as multi-year-cumulative moisture conditions of maximum VPD, were top predictors for start of season. Productivity measures were also associated with contextual variables such as soil characteristics and vegetation composition. Phenology is a key process that profoundly affects organism-environment relationships, spatio-temporal patterns in ecosystem structure and function, and other ecosystem dynamics. Phenology, however, is complex, and is mediated by lagged effects, interactions, and a diversity of potential drivers; nonetheless, the incorporation of antecedent conditions and contextual variables can improve models of phenology.

The ecological implications of intra- and inter-species variation in phenological sensitivity

Plant-pollinator mutualisms rely upon the synchrony of interacting taxa. Climate change can disrupt this synchrony as phenological responses to climate vary within and across species. However, intra- and interspecific variation in phenological responses is seldom considered simultaneously, limiting our understanding of climate change impacts on interactions among taxa across their ranges. We investigated how variation in phenological sensitivity to climate can alter ecological interactions simultaneously within and among species using natural history collections and citizen science data. We focus on a unique system, comprising a wide-ranged spring ephemeral with varying color morphs (Claytonia virginica) and its specialist bee pollinator (Andrena erigeniae). We found strongly opposing trends in the phenological sensitivities of plants vs their pollinators. Flowering phenology was more sensitive to temperature in warmer regions, whereas bee phenology was more responsive in colder regions. Phenological sensitivity varied across flower color morphs. Temporal synchrony between flowering and pollinator activity was predicted to change heterogeneously across the species' ranges in the future. Our work demonstrates the complexity and fragility of ecological interactions in time and the necessity of incorporating variation in phenological responses across multiple axes to understand how such interactions will change in the future.

Consistent trait-temperature interactions drive butterfly phenology in both incidental and survey data

Data availability limits phenological research at broad temporal and spatial extents. Butterflies are among the few taxa with broad-scale occurrence data, from both incidental reports and formal surveys. Incidental reports have biases that are challenging to address, but structured surveys are often limited seasonally and may not span full flight phenologies. Thus, how these data source compare in phenological analyses is unclear. We modeled butterfly phenology in relation to traits and climate using parallel analyses of incidental and survey data, to explore their shared utility and potential for analytical integration. One workflow aggregated “Pollard” surveys, where sites are visited multiple times per year; the other aggregated incidental data from online portals: iNaturalist and eButterfly. For 40 species, we estimated early (10%) and mid (50%) flight period metrics, and compared the spatiotemporal patterns and drivers of phenology across species and between datasets. For both datasets, inter-annual variability was best explained by temperature, and seasonal emergence was earlier for resident species overwintering at more advanced stages. Other traits related to habitat, feeding, dispersal, and voltinism had mixed or no impacts. Our results suggest that data integration can improve phenological research, and leveraging traits may predict phenology in poorly studied species.

New directions in tropical phenology

Earth's most speciose biomes are in the tropics, yet tropical plant phenology remains poorly understood. Tropical phenological data are comparatively scarce and viewed through the lens of a 'temperate phenological paradigm' expecting phenological traits to respond to strong, predictably annual shifts in climate (e.g., between subfreezing and frost-free periods). Digitized herbarium data greatly expand existing phenological data for tropical plants; and circular data, statistics, and models are more appropriate for analyzing tropical (and temperate) phenological datasets. Phylogenetic information, which remains seldom applied in phenological investigations, provides new insights into phenological responses of large groups of related species to climate. Consistent combined use of herbarium data, circular statistical distributions, and robust phylogenies will rapidly advance our understanding of tropical - and temperate - phenology.

Multiple UAV Flights across the Growing Season Can Characterize Fine Scale Phenological Heterogeneity within and among Vegetation Functional Groups

Grasslands and shrublands exhibit pronounced spatial and temporal variability in structure and function with differences in phenology that can be difficult to observe. Unpiloted aerial vehicles (UAVs) can measure vegetation spectral patterns relatively cheaply and repeatably at fine spatial resolution. We tested the ability of UAVs to measure phenological variability within vegetation functional groups and to improve classification accuracy at two sites in Montana, U.S.A. We tested four flight frequencies during the growing season. Classification accuracy based on reference data increased by 5–10% between a single flight and scenarios including all conducted flights. Accuracy increased from 50.6 to 61.4% at the drier site, while at the more mesic/densely vegetated site, we found an increase of 59.0 to 64.4% between a single and multiple flights over the growing season. Peak green-up varied by 2–4 weeks within the scenes, and sparse vegetation classes had only a short detectable window of active phtosynthesis; therefore, a single flight could not capture all vegetation that was active across the growing season. The multi-temporal analyses identified differences in the seasonal timing of green-up and senescence within herbaceous and sagebrush classes. Multiple UAV measurements can identify the fine-scale phenological variability in complex mixed grass/shrub vegetation.

Phenological displacement is uncommon among sympatric angiosperms

Interactions between species can influence successful reproduction, resulting in reproductive character displacement, where the similarity of reproductive traits - such as flowering time - among close relatives growing together differ from when growing apart. Evidence for the overall prevalence and direction of this phenomenon, and its stability under environmental change, remains untested across large scales. Using the power of crowdsourcing, we gathered phenological information from over 40 000 herbarium specimens, and investigated displacement in flowering time across 110 animal-pollinated species in the eastern USA. Overall, flowering time displacement is not common across large scales. However, displacement is generally greater among species pairs that flower close in time, regardless of direction. Furthermore, with climate change, the flowering times of closely related species are predicted, on average, to shift further apart by the mid-21^st century. We demonstrate that the degree and direction of phenological displacement among co-occurring closely related species pairs varies tremendously. However, future climate change may alter the differences in reproductive timing among many of these species pairs, which may have significant consequences for species interactions and gene flow. Our study provides one promising path towards understanding how the phenological landscape is structured and may respond to future environmental change.

Observing the Observers: How Participants Contribute Data to iNaturalist and Implications for Biodiversity Science

The availability of citizen science data has resulted in growing applications in biodiversity science. One widely used platform, iNaturalist, provides millions of digitally vouchered observations submitted by a global user base. These observation records include a date and a location but otherwise do not contain any information about the sampling process. As a result, sampling biases must be inferred from the data themselves. In the present article, we examine spatial and temporal biases in iNaturalist observations from the platform's launch in 2008 through the end of 2019. We also characterize user behavior on the platform in terms of individual activity level and taxonomic specialization. We found that, at the level of taxonomic class, the users typically specialized on a particular group, especially plants or insects, and rarely made observations of the same species twice. Biodiversity scientists should consider whether user behavior results in systematic biases in their analyses before using iNaturalist data.