Co‐occurrence models fail to infer underlying patterns of avoidance and aggregation when closure is violated

Advances in multi‐species monitoring have prompted an increase in the use of multi‐species occupancy analyses to assess patterns of co‐occurrence among species, even when data were collected at scales likely violating the assumption that sites were closed to changes in the occupancy state for the target species. Violating the closure assumption may lead to erroneous conclusions related to patterns of co‐occurrence among species. Occurrence for two hypothetical species was simulated under patterns of avoidance, aggregation, or independence, when the closure assumption was either met or not. Simulated populations were sampled at two levels (N = 250 or 100 sites) and two scales of temporal resolution for surveys. Sample data were analyzed with conditional two‐species occupancy models, and performance was assessed based on the proportion of simulations recovering the true pattern of co‐occurrence. Estimates of occupancy were unbiased when closure was met, but biased when closure violations occurred; bias increased when sample size was small and encounter histories were collapsed to a large‐scale temporal resolution. When closure was met and patterns of avoidance and aggregation were simulated, conditional two‐species models tended to correctly find support for non‐independence, and estimated species interaction factors (SIF) aligned with predicted values. By contrast, when closure was violated, models tended to incorrectly infer a pattern of independence and power to detect simulated patterns of avoidance or aggregation that decreased with smaller sample size. Results suggest that when the closure assumption is violated, co‐occurrence models often fail to detect underlying patterns of avoidance or aggregation, and incorrectly identify a pattern of independence among species, which could have negative consequences for our understanding of species interactions and conservation efforts. Thus, when closure is violated, inferred patterns of independence from multi‐species occupancy should be interpreted cautiously, and evidence of avoidance or aggregation is likely a conservative estimate of true pattern or interaction.     

Context‐dependency in carnivore co‐occurrence across a multi‐use conservation landscape

Carnivore intraguild dynamics depend on a complex interplay of environmental affinities and interspecific interactions. Context‐dependency is commonly expected with varying suites of interacting species and environmental conditions but seldom empirically described. In South Africa, decentralized approaches to conservation and the resulting multi‐tenure conservation landscapes have markedly altered the environmental stage that shapes the structure of local carnivore assemblages. We explored assemblage‐wide patterns of carnivore spatial (residual occupancy probability) and temporal (diel activity overlap) co‐occurrence across three adjacent wildlife‐oriented management contexts—a provincial protected area, a private ecotourism reserve, and commercial game ranches. We found that carnivores were generally distributed independently across space, but existing spatial dependencies were context‐specific. Spatial overlap was most common in the protected area, where species occur at higher relative abundances, and in game ranches, where predator persecution presumably narrows the scope for spatial asymmetries. In the private reserve, spatial co‐occurrence patterns were more heterogeneous but did not follow a dominance hierarchy associated with higher apex predator densities. Pair‐specific variability suggests that subordinate carnivores may alternate between pre‐emptive behavioral strategies and fine‐scale co‐occurrence with dominant competitors. Consistency in species‐pairs diel activity asynchrony suggested that temporal overlap patterns in our study areas mostly depend on species'' endogenous clock rather than the local context. Collectively, our research highlights the complexity and context‐dependency of guild‐level implications of current management and conservation paradigms; specifically, the unheeded potential for interventions to influence the local network of carnivore interactions with unknown population‐level and cascading effects.     

Space‐for‐time is not necessarily a substitution when monitoring the distribution of pelagic fishes in the San Francisco Bay‐Delta

Occupancy models are often used to analyze long‐term monitoring data to better understand how and why species redistribute across dynamic landscapes while accounting for incomplete capture. However, this approach requires replicate detection/non‐detection data at a sample unit and many long‐term monitoring programs lack temporal replicate surveys. In such cases, it has been suggested that surveying subunits within a larger sample unit may be an efficient substitution (i.e., space‐for‐time substitution). Still, the efficacy of fitting occupancy models using a space‐for‐time substitution has not been fully explored and is likely context dependent. Herein, we fit occupancy models to Delta Smelt (Hypomesus transpacificus) and Longfin Smelt (Spirinchus thaleichthys) catch data collected by two different monitoring programs that use the same sampling gear in the San Francisco Bay‐Delta, USA. We demonstrate how our inferences concerning the distribution of these species changes when using a space‐for‐time substitution. Specifically, we found the probability that a sample unit was occupied was much greater when using a space‐for‐time substitution, presumably due to the change in the spatial scale of our inferences. Furthermore, we observed that as the spatial scale of our inferences increased, our ability to detect environmental effects on system dynamics was obscured, which we suspect is related to the tradeoffs associated with spatial grain and extent. Overall, our findings highlight the importance of considering how the unique characteristics of monitoring programs influences inferences, which has broad implications for how to appropriately leverage existing long‐term monitoring data to understand the distribution of species.     

Multi‐surveyor capture‐mark‐recapture as a powerful tool for butterfly population monitoring in the pre‐imaginal stage

For many elusive insect species, which are difficult to cover by standard monitoring schemes, innovative survey methods are needed to gain robust data on abundance and population trends. We suggest a monitoring of overwintering larvae for the endangered nymphalid butterfly Limenitis reducta. We tested different removal and capture‐mark‐recapture (CMR) approaches in a field study in the “Alb‐Donau” region, Germany. Classical removal and CMR studies require movement of the organisms under study, but in our approach, we replaced movement of the study organisms by random movement of multiple different surveyors. We tested the validity of the approach by comparing detection frequencies from our field data with simulated detections. Our results indicate that multi‐surveyor removal/CMR techniques are suitable for estimating abundance of overwintering L. reducta larvae. Depending on surveyor experience, the average detection probability ranged between 16% for novices and 35% for experts. The uncertainty of population estimates increased with a decrease in personnel expenditure. Estimated larval densities on a spruce clear‐cut varied between one and three individuals per 100 m², probably related to habitat conditions. We suggest a CMR approach with three to four trained surveyors for the monitoring of L. reducta populations in the overwintering stage. Compared with previous sampling methods, our approach is a powerful tool with clear advantages: long survey period, estimates of the absolute population size accompanied by uncertainty measures, and estimates of overwinter mortality. The proposed method can be adapted and used for several different butterfly species, other insect taxa with specific immobile life stages, and some sessile organisms, for example, elusive plants, fungi, or corals.     

Concordant geographic and genetic structure revealed by genotyping‐by‐sequencing in a New Zealand marine isopod

Population genetic structure in the marine environment can be influenced by life‐history traits such as developmental mode (biphasic, with distinct adult and larval morphology, and direct development, in which larvae resemble adults) or habitat specificity, as well as geography and selection. Developmental mode is thought to significantly influence dispersal, with direct developers expected to have much lower dispersal potential. However, this prediction can be complicated by the presence of geophysical barriers to dispersal. In this study, we use a panel of 8,020 SNPs to investigate population structure and biogeography over multiple spatial scales for a direct‐developing species, the New Zealand endemic marine isopod Isocladus armatus. Because our sampling range is intersected by two well‐known biogeographic barriers (the East Cape and the Cook Strait), our study provides an opportunity to understand how such barriers influence dispersal in direct developers. On a small spatial scale (20 km), gene flow between locations is extremely high, suggestive of an island model of migration. However, over larger spatial scales (600 km), populations exhibit a clear pattern of isolation‐by‐distance. Our results indicate that I. armatus exhibits significant migration across the hypothesized barriers and suggest that large‐scale ocean currents associated with these locations do not present a barrier to dispersal. Interestingly, we find evidence of a north‐south population genetic break occurring between Māhia and Wellington. While no known geophysical barrier is apparent in this area, it coincides with the location of a proposed border between bioregions. Analysis of loci under selection revealed that both isolation‐by‐distance and adaption may be contributing to the degree of population structure we have observed here. We conclude that developmental life history largely predicts dispersal in the intertidal isopod I. armatus. However, localized biogeographic processes can disrupt this expectation, and this may explain the potential meta‐population detected in the Auckland region.     

An open spatial capture–recapture model for estimating density,movement, and population dynamics from line‐transect surveys

The purpose of many wildlife population studies is to estimate density, movement, or demographic parameters. Linking these parameters to covariates, such as habitat features, provides additional ecological insight and can be used to make predictions for management purposes. Line‐transect surveys, combined with distance sampling methods, are often used to estimate density at discrete points in time, whereas capture–recapture methods are used to estimate movement and other demographic parameters. Recently, open population spatial capture–recapture models have been developed, which simultaneously estimate density and demographic parameters, but have been made available only for data collected from a fixed array of detectors and have not incorporated the effects of habitat covariates. We developed a spatial capture–recapture model that can be applied to line‐transect survey data by modeling detection probability in a manner analogous to distance sampling. We extend this model to a) estimate demographic parameters using an open population framework and b) model variation in density and space use as a function of habitat covariates. The model is illustrated using simulated data and aerial line‐transect survey data for North Atlantic right whales in the southeastern United States, which also demonstrates the ability to integrate data from multiple survey platforms and accommodate differences between strata or demographic groups. When individuals detected from line‐transect surveys can be uniquely identified, our model can be used to simultaneously make inference on factors that influence spatial and temporal variation in density, movement, and population dynamics.     

Benefits and limitations of a new genome‐based PCR‐RFLP genotyping assay (GB‐RFLP): A SNP‐based detection method for identification of species in extremely young adaptive radiations

High‐throughput DNA sequencing technologies make it possible now to sequence entire genomes relatively easily. Complete genomic information obtained by whole‐genome resequencing (WGS) can aid in identifying and delineating species even if they are extremely young, cryptic, or morphologically difficult to discern and closely related. Yet, for taxonomic or conservation biology purposes, WGS can remain cost‐prohibitive, too time‐consuming, and often constitute a “data overkill.” Rapid and reliable identification of species (and populations) that is also cost‐effective is made possible by species‐specific markers that can be discovered by WGS. Based on WGS data, we designed a PCR restriction fragment length polymorphism (PCR‐RFLP) assay for 19 Neotropical Midas cichlid populations (Amphilophus cf. citrinellus), that includes all 13 described species of this species complex. Our work illustrates that identification of species and populations (i.e., fish from different lakes) can be greatly improved by designing genetic markers using available “high resolution” genomic information. Yet, our work also shows that even in the best‐case scenario, when whole‐genome resequencing information is available, unequivocal assignments remain challenging when species or populations diverged very recently, or gene flow persists. In summary, we provide a comprehensive workflow on how to design RFPL markers based on genome resequencing data, how to test and evaluate their reliability, and discuss the benefits and pitfalls of our approach.     

Temperature‐induced multi‐species cohort effects in sympatric snakes

In reptiles, reproductive maturity is often determined by size rather than age. Consequently, growth early in life may influence population dynamics through effects on generation time and survival to reproduction. Because reproductive phenology and pre‐ and post‐natal growth are temperature dependent, environmental conditions may induce multi‐species cohort effects on body size in sympatric reptiles. I present evidence of this using 10 years of neonatal size data for three sympatric viviparous snakes, Dekay''s Brown snakes (Storeria dekayi), Red‐bellied Snakes (S. occipitomaculata), and Common Garter snakes (Thamnophis sirtalis). End‐of‐season neonatal size varied in parallel across species such that snout–vent length was 36%–61% greater and mass was 65%–223% greater in years when gestating females could achieve higher April–May (vs. June–July or August–September) operative temperatures. Thus, temperature had a larger impact during follicular enlargement and ovulation than during gestation or post‐natal growth. Multi‐species cohort effects like these may affect population dynamics and the magnitude of these effects may increase with climate change.     

Niche theory‐based modeling of assembly processes of viral communities in bats

Understanding the assembly processes of symbiont communities, including viromes and microbiomes, is important for improving predictions on symbionts’ biogeography and disease ecology. Here, we use phylogenetic, functional, and geographic filters to predict the similarity between symbiont communities, using as a test case the assembly process in viral communities of Mexican bats. We construct generalized linear models to predict viral community similarity, as measured by the Jaccard index, as a function of differences in host phylogeny, host functionality, and spatial co‐occurrence, evaluating the models using the Akaike information criterion. Two model classes are constructed: a “known” model, where virus–host relationships are based only on data reported in Mexico, and a “potential” model, where viral reports of all the Americas are used, but then applied only to bat species that are distributed in Mexico. Although the “known” model shows only weak dependence on any of the filters, the “potential” model highlights the importance of all three filter types—phylogeny, functional traits, and co‐occurrence—in the assemblage of viral communities. The differences between the “known” and “potential” models highlight the utility of modeling at different “scales” so as to compare and contrast known information at one scale to another one, where, for example, virus information associated with bats is much scarcer.     

Camera trap reveals the co‐occurrence patterns of two sympatric muntjac species in southern Anhui Province,China: No spatial segregation

The competitive relationship and coexistence pattern among close related species have long been one of the hot issues in ecological research. Interspecies interactions can exert important influences on the local distribution of rare species. Black muntjac Muntiacus crinifrons is an endemic species to eastern China, currently restricted to limited regions. In contrast, Chinese muntjac Muntiacus reevesi is the most common and widespread deer in southern China. Both species co‐occur in southern Anhui and western Zhejiang Province. Little is known about the interaction of these two sympatric‐related species. In this study, to investigate the site use determinants and co‐occurrence pattern of the two sympatric muntjac species, we conducted a camera trap survey across about 250 km² in mountainous area of southern Anhui Province, China. We adopted a multistep approach to incorporate habitat preferences while modeling occupancy and detection. We found that the two species did not separate along elevation gradient (range from 400 m to 1,400 m) as described in previous studies. Results of single‐species occupancy models indicated that elevation had positive effects on the site use of both species, while slope had an opposite influence on their site use. Positive effects of elevation on the site use implied that both species try to avoid human interference at low elevations. Significant negative effect of slope on the site use of black muntjac suggested that the species prefer habitat with gentle slope and avoided steep. Co‐occurrence models and species interaction factors provided evidence that the two muntjac species had an independent occupancy (ψ ^{BM CM} = ψ ^{BM cm}, SIF = 1) and exhibited a positive species interaction in detection probability (p ^BM < r ^{BM CM}). Combined with the results of previous studies, we suggested that it was fine differentiation in microhabitats and food resources utilization rather spatial or temporal segregation that allowed the two species co‐occurrence. The site use determinants revealed in our study would be useful for the habitat conservation and restoration for the rare black muntjac, and the co‐occurrence pattern of the two sympatric muntjac species could provide useful information for deep understanding of the coexistence mechanism among forest‐dwelling ungulates.     

Multi‐modal adaptor‐clathrin contacts drive coated vesicle assembly

Clathrin‐coated pits are formed by the recognition of membrane and cargo by the AP2 complex and the subsequent recruitment of clathrin triskelia. A role for AP2 in coated‐pit assembly beyond initial clathrin recruitment has not been explored. Clathrin binds the β2 subunit of AP2, and several binding sites have been identified, but our structural knowledge of these interactions is incomplete and their functional importance during endocytosis is unclear. Here, we analysed the cryo‐EM structure of clathrin cages assembled in the presence of β2 hinge‐appendage (β2HA). We find that the β2‐appendage binds in at least two positions in the cage, demonstrating that multi‐modal binding is a fundamental property of clathrin‐AP2 interactions. In one position, β2‐appendage cross‐links two adjacent terminal domains from different triskelia. Functional analysis of β2HA‐clathrin interactions reveals that endocytosis requires two clathrin interaction sites: a clathrin‐box motif on the hinge and the “sandwich site” on the appendage. We propose that β2‐appendage binding to more than one triskelion is a key feature of the system and likely explains why assembly is driven by AP2.     

High‐resolution imaging sheds new light on a multi‐tier symbiotic partnership between a “walking” solitary coral,a sipunculan,and a bivalve from East Africa

Marine symbioses are integral to the persistence of ecosystem functioning in coral reefs. Solitary corals of the species Heteropsammia cochlea and Heterocyathus aequicostatus have been observed to live in symbiosis with the sipunculan worm Aspidosiphon muelleri muelleri, which inhabits a cavity within the coral, in Zanzibar (Tanzania). The symbiosis of these photosymbiotic corals enables the coral holobiont to move, in fine to coarse unconsolidated substrata, a process termed as “walking.” This allows the coral to escape sediment cover in turbid conditions which is crucial for these light‐dependent species. An additional commensalistic symbiosis of this coral‐worm holobiont is found between the Aspidosiphon worm and the cryptoendolithic bivalve Jousseaumiella sp., which resides within the cavity of the coral skeleton. To understand the morphological alterations caused by these symbioses, interspecific relationships, with respect to the carbonate structures between these three organisms, are documented using high‐resolution imaging techniques (scanning electron microscopy and µCT scanning). Documenting multi‐layered symbioses can shed light on how morphological plasticity interacts with environmental conditions to contribute to species persistence.     

Changes in vegetation structure and composition of a lowland mire over a sixty‐five‐year interval

1. Mires are characterized by plant communities of high conservation and societal value, which have experienced a major decline in area in many parts of the world, particularly Europe. Evidence suggests that they may be particularly vulnerable to changes in climate and nutrient addition. Although they have been the focus of extensive paleoecological research, few attempts have been made to examine the dynamics of mire vegetation during the current era of anthropogenic environmental change.
2. To assess long‐term change in the spatial structure and composition of a lowland mire community, in 2016 we resurveyed plots first surveyed in 1951. Measures of species richness and composition were compared between the two surveys, and changes in community composition were related to plant traits.
3. Overall, mean species richness declined by 26%. The area of occupancy declined in 37% of species, which were primarily oligotrophic species typical of nutrient‐poor bog communities. Conversely, occupancy increased in 21% of species, especially those that were more tolerant of higher nutrient availability. These changes were associated with variation in plant functional traits, as indicated by an increase mean Ellenberg trait values for nitrogen and mean temperature, and a decline in values for precipitation. These results suggest that eutrophication and climate change have been key drivers of floristic change on this site.
4. Synthesis. This investigation provides a rare assessment of the dynamics of a mire community over a multi‐decadal interval. Results indicate that substantial change has occurred in the composition of the community, and the distribution of species within it. The investigation provides evidence of the impact of environmental change on the composition and structure of a lowland mire community, and highlights challenges for its future conservation.

     

White‐tailed deer exploit temporal refuge from multi‐predator and human risks on roads

Although most prey have multiple predator species, few studies have quantified how prey respond to the temporal niches of multiple predators which pose different levels of danger. For example, intraspecific variation in diel activity allows white‐tailed deer (Odocoileus virginianus) to reduce fawn activity overlap with coyotes (Canis latrans) but finding safe times of day may be more difficult for fawns in a multi‐predator context. We hypothesized that within a multi‐predator system, deer would allocate antipredation behavior optimally based on combined mortality risk from multiple sources, which would vary depending on fawn presence. We measured cause‐specific mortality of 777 adult (>1‐year‐old) and juvenile (1–4‐month‐old) deer and used 300 remote cameras to estimate the activity of deer, humans, and predators including American black bears (Ursus americanus), bobcats (Lynx rufus), coyotes, and wolves (Canis lupus). Predation and vehicle collisions accounted for 5.3 times greater mortality in juveniles (16% mortality from bears, coyotes, bobcats, wolves, and vehicles) compared with adults (3% mortality from coyotes, wolves, and vehicles). Deer nursery groups (i.e., ≥1 fawn present) were more diurnal than adult deer without fawns, causing fawns to have 24–38% less overlap with carnivores and 39% greater overlap with humans. Supporting our hypothesis, deer nursery groups appeared to optimize diel activity to minimize combined mortality risk. Temporal refuge for fawns was likely the result of carnivores avoiding humans, simplifying diel risk of five species into a trade‐off between diurnal humans and nocturnal carnivores. Functional redundancy among multiple predators with shared behaviors may partially explain why white‐tailed deer fawn predation rates are often similar among single‐ and multi‐predator systems.     

Habitat selection in a fluctuating ground squirrel population: Density‐dependence and fitness consequences

Investigating individual‐based habitat settlement decisions is a central theme in ecology, yet studies that quantify density‐dependent habitat selection or tie fitness to resource selection decisions remain rare. We quantified habitat selection in golden‐mantled ground squirrels (Callospermophilus lateralis) across two spatial scales (home‐range placement, and occurrence within the home range) by using 11 consecutive years of data on individual space use, and we used resource selection functions and multilevel modeling to address how habitat preferences may be influenced by density or linked to fitness outcomes. Squirrels preferred dry meadow over other habitat types (wet meadow, aspen, spruce, and willow) at both spatial scales. Squirrels were more likely to use dry meadow that contained shorter vegetation and vision‐enhancing prominences such as rocks (“perches”). The use of dry meadow at each scale was not influenced by changes in density. The use of dry meadow did not lead to increased litter size, pre‐hibernation mass, or survival. However, squirrels that experienced a greater number of perches or lower local densities had higher survival rates. Our results suggest that a lack of visual obstruction, probably facilitating detection of predators, drives habitat selection in this system. Surprisingly, squirrels maintained their preference for dry meadow as density increased, and they experienced reduced survival as a result. This work furthers our understanding about the causes and consequences of changes in habitat use, informing wildlife management and conservation.     

Diauxic lags explain unexpected coexistence in multi‐resource environments

How the coexistence of species is affected by the presence of multiple resources is a major question in microbial ecology. We experimentally demonstrate that differences in diauxic lags, which occur as species deplete their own environments and adapt their metabolisms, allow slow‐growing microbes to stably coexist with faster‐growing species in multi‐resource environments despite being excluded in single‐resource environments. In our focal example, an Acinetobacter species (Aci2) competitively excludes Pseudomonas aurantiaca (Pa) on alanine and on glutamate. However, they coexist on the combination of both resources. Experiments reveal that Aci2 grows faster but Pa has shorter diauxic lags. We establish a tradeoff between Aci2’s fast growth and Pa’s short lags as their mechanism for coexistence. We model this tradeoff to accurately predict how environmental changes affect community composition. We extend our work by surveying a large set of competitions and observe coexistence nearly four times as frequently when the slow‐grower is the fast‐switcher. Our work illustrates a simple mechanism, based entirely on supplied‐resource growth dynamics, for the emergence of multi‐resource coexistence.     

CellRegMap: a statistical framework for mapping context‐specific regulatory variants using scRNA‐seq

Single‐cell RNA sequencing (scRNA‐seq) enables characterizing the cellular heterogeneity in human tissues. Recent technological advances have enabled the first population‐scale scRNA‐seq studies in hundreds of individuals, allowing to assay genetic effects with single‐cell resolution. However, existing strategies to analyze these data remain based on principles established for the genetic analysis of bulk RNA‐seq. In particular, current methods depend on a priori definitions of discrete cell types, and hence cannot assess allelic effects across subtle cell types and cell states. To address this, we propose the Cell Regulatory Map (CellRegMap), a statistical framework to test for and quantify genetic effects on gene expression in individual cells. CellRegMap provides a principled approach to identify and characterize genotype–context interactions of known eQTL variants using scRNA‐seq data. This model‐based approach resolves allelic effects across cellular contexts of different granularity, including genetic effects specific to cell subtypes and continuous cell transitions. We validate CellRegMap using simulated data and apply it to previously identified eQTL from two recent studies of differentiating iPSCs, where we uncover hundreds of eQTL displaying heterogeneity of genetic effects across cellular contexts. Finally, we identify fine‐grained genetic regulation in neuronal subtypes for eQTL that are colocalized with human disease variants.     

Using a multi‐model ensemble approach to determine biodiversity hotspots with limited occurrence data in understudied areas: An example using freshwater mussels in México

Species distribution models (SDMs) are an increasingly important tool for conservation particularly for difficult‐to‐study locations and with understudied fauna. Our aims were to (1) use SDMs and ensemble SDMs to predict the distribution of freshwater mussels in the Pánuco River Basin in Central México; (2) determine habitat factors shaping freshwater mussel occurrence; and (3) use predicted occupancy across a range of taxa to identify freshwater mussel biodiversity hotspots to guide conservation and management. In the Pánuco River Basin, we modeled the distributions of 11 freshwater mussel species using an ensemble approach, wherein multiple SDM methodologies were combined to create a single ensemble map of predicted occupancy. A total of 621 species‐specific observations at 87 sites were used to create species‐specific ensembles. These predictive species ensembles were then combined to create local diversity hotspot maps. Precipitation during the warmest quarter, elevation, and mean temperature were consistently the most important discriminatory environmental variables among species, whereas land use had limited influence across all taxa. To the best of our knowledge, our study is the first freshwater mussel‐focused research to use an ensemble approach to determine species distribution and predict biodiversity hotspots. Our study can be used to guide not only current conservation efforts but also prioritize areas for future conservation and study.     

Remotely sensed between‐individual functional trait variation in a temperate forest

1. Trait‐based ecology holds the promise to explain how plant communities work, for example, how functional diversity may support community productivity. However, so far it has been difficult to combine field‐based approaches assessing traits at the level of plant individuals with limited spatial coverage and approaches using remote sensing (RS) with complete spatial coverage but assessing traits at the level of vegetation pixels rather than individuals. By delineating all individual‐tree crowns within a temperate forest site and then assigning RS‐derived trait measures to these trees, we combine the two approaches, allowing us to use general linear models to estimate the influence of taxonomic or environmental variation on between‐ and within‐species variation across contiguous space.
2. We used airborne imaging spectroscopy and laser scanning to collect individual‐tree RS data from a mixed conifer‐angiosperm forest on a mountain slope extending over 5.5 ha and covering large environmental gradients in elevation as well as light and soil conditions. We derived three biochemical (leaf chlorophyll, carotenoids, and water content) and three architectural traits (plant area index, foliage‐height diversity, and canopy height), which had previously been used to characterize plant function, from the RS data. We then quantified the contributions of taxonomic and environmental variation and their interaction to trait variation and partitioned the remaining within‐species trait variation into smaller‐scale spatial and residual variation. We also investigated the correlation between functional trait and phylogenetic distances at the between‐species level. The forest consisted of 13 tree species of which eight occurred in sufficient abundance for quantitative analysis.
3. On average, taxonomic variation between species accounted for more than 15% of trait variation in biochemical traits but only around 5% (still highly significant) in architectural traits. Biochemical trait distances among species also showed a stronger correlation with phylogenetic distances than did architectural trait distances. Light and soil conditions together with elevation explained slightly more variation than taxonomy across all traits, but in particular increased plant area index (light) and reduced canopy height (elevation). Except for foliage‐height diversity, all traits were affected by significant interactions between taxonomic and environmental variation, the different responses of the eight species to the within‐site environmental gradients potentially contributing to the coexistence of the eight abundant species.
4. We conclude that with high‐resolution RS data it is possible to delineate individual‐tree crowns within a forest and thus assess functional traits derived from RS data at individual level. With this precondition fulfilled, it is then possible to apply tools commonly used in field‐based trait ecology to partition trait variation among individuals into taxonomic and potentially even genetic variation, environmental variation, and interactions between the two. The method proposed here presents a promising way of assessing individual‐based trait information with complete spatial coverage and thus allowing analysis of functional diversity at different scales. This information can help to better understand processes shaping community structure, productivity, and stability of forests.