Drivers Of Species Distribution Research Articles

Thermal Tolerance and Species Distributions: Interactions Between Latitude, Elevation and Arboreality in Ants

ABSTRACTAimGlobal warming has highlighted the importance of understanding the role of thermal tolerance as a driver of species distributions, especially for ectotherms. Here we analyse interactions between latitude, elevation and arboreality as predictors of geographic patterns of thermal tolerance in ants.LocationThe collected data are distributed globally.MethodsWe first tested the effect of latitude, elevation and its interactive effect on ant CTmax and CTmin. Second, we tested whether CTmax and CTmin are phylogenetic clustering. Finally, we tested whether CTmax and CTmin can be explained by nesting microhabitat (ground vs. tree‐nesting species) and whether the probability of occurrence of tree‐nesting species along thermal gradients helps explain the global pattern in ant CTmax.ResultsCTmax and CTmin displayed high and low phylogenetic signals respectively and therefore showed different responses to geographic gradients. Notably, we found that CTmax was higher in higher latitudes. This was explained by a lack of elevational turnover at high latitudes among tree‐nesting species, which are exposed to warmer microclimates and have higher CTmax compared with ground‐nesting species. CTmin decreased with elevation at low latitudes, but did not vary with elevation at higher latitudes.Main ConclusionsOur findings highlight the important influence of arboreality on the macroecology of thermal tolerance, substantially modifying traditional notions of variation along latitudinal and elevation gradients.

A matter of scale: Identifying the best spatial and temporal scale of environmental variables to model the distribution of a small cetacean.

The importance of scale when investigating ecological patterns and processes is recognised across many species. In marine ecosystems, the processes that drive species distribution have a hierarchical structure over multiple nested spatial and temporal scales. Hence, multi-scale approaches should be considered when developing accurate distribution models to identify key habitats, particularly for populations of conservation concern. Here, we propose a modelling procedure to identify the best spatial and temporal scale for each modelled and remotely sensed oceanographic variable to model harbour porpoise (Phocoena phocoena) distribution within the Irish Exclusive Economic Zone. Harbour porpoise sightings were recorded during dedicated line-transect aerial surveys conducted in the summers of 2016, 2021 and 2022. Binary generalised additive models were used to assess the relationships between porpoise presence and oceanographic variables at different spatial (5-40 km) and temporal (daily, monthly and across survey period) scales. Selected variables included sea surface temperature, thermal fronts, chlorophyll-a, sea surface height, mixed layer depth and salinity. A total of 30,514 km was covered on-effort with 216 harbour porpoise sightings recorded. Overall, the best spatial scale corresponded to the coarsest resolution considered in this study (40 km), while porpoise presence showed stronger association with oceanographic variables summarised at a longer temporal scale. Habitat models including covariates at coarse spatial and temporal scales may better reflect the processes driving availability and abundance of resources at these large scales. These findings support the hypothesis that a multi-scale approach should be applied when investigating species distribution. Identifying suitable spatial and temporal scale would improve the functional interpretation of the underlying relationships, particularly when studying how a small marine predator interacts with its environment and responds to climate and ecosystem changes.

Spatial distribution of the range-expanding species Seriola fasciata (Bloch, 1793) in Mediterranean Sea: From past to future

The Mediterranean Sea is a highly susceptible area to climate change, that facilitates the introduction of warm-affinity exotic species, contributing to the expansion of their biogeographical range. One such thermophilic species is the Atlantic fish Seriola fasciata, which has colonised this area over the past three decades. The present study analyzed its spatial distribution in the Mediterranean Sea to identify aggregation areas and dynamics over time, and the environmental predictors influencing its presence.The utilized statistical tools and the Species Distribution Model proved effective in identifying specific spatial and temporal distribution patterns, as well as discerning some environmental variables influencing the species presence, with distinctions recorded between juveniles and adults.S. fasciata was observed to be established in the central Mediterranean, with Fishing Aggregating Devices potentially influencing its presence, particularly of juveniles. Sea floor temperature and habitats emerged as the primary factors driving species distribution. An aggregation area in the Levant Sea, conducive mainly for the adults, was identified and is expected to intensify over time.These findings contribute valuable insights into a relatively understudied species and its presence in the Mediterranean Sea, where climate change is affecting marine biodiversity.

Spatially‐nested hierarchical species distribution models to overcome niche truncation in national‐scale studies

Spatial truncation in species distribution models (SDMs) might cause niche truncation and model transferability issues, particularly when extrapolating models to non‐analog environmental conditions. While broad calibration extents reduce truncation issues, they usually overlook local ecological factors driving species distributions at finer resolution. Spatially‐nested hierarchical SDMs (HSDMs) address truncation by merging (a) a global model calibrated with broadly extended, yet typically low‐resolution, basic, and imprecise data; and (b) a regional model calibrated with spatially restricted but more precise and reliable data. This study aimed to examine HSDMs' efficacy to overcome spatial truncation in national‐scale studies. We compared two hierarchical strategies (‘covariate', which uses the global model output as a covariate for the regional model, and ‘multiply', which calculates the geometric mean of the global and regional models) and a non‐hierarchical strategy. The three strategies were compared in terms of niche truncation, environmental extrapolation, model performance, species' predicted distributions and shifts, and trends in species richness. We examined the consistency of the results over two study areas (Spain and Switzerland), 108 tree species, and four future climate scenarios. Only the non‐hierarchical strategy was susceptible to niche truncation, and environmental extrapolation issues. Hierarchical strategies, particularly the ‘covariate' one, presented greater model accuracy than non‐hierarchical strategies. The non‐hierarchical strategy predicted the highest overall values and the lowest decreases over time in species distribution ranges and richness. Differences between strategies were more evident in Switzerland, which was more affected by niche truncation issues. Spain was more negatively affected by climate change and environmental extrapolation. The ‘covariate' strategy exhibited higher model performance than the ‘multiply' one. However, uncertainties regarding model temporal transferability advocate for adopting and further examining multiple hierarchical approaches. This research underscores the importance of adopting spatially‐nested hierarchical SDMs given the compromised reliability of non‐hierarchical approaches due to niche truncation and extrapolation issues.

Overcoming the lack of distribution data for range‐restricted habitat specialist frogs

AbstractWith over 40% of species threatened with extinction, the distributional range of most amphibians is still unknown, generating uncertainty whether species are naturally rare, under‐sampled, or difficult to detect. We implemented a modelling approach that uses bromeliads as surrogates to predict the distribution of habitat specialist frogs that lack distributional data. We aimed to predict and survey potential new sites for the occurrence of a rare and microendemic bromeliad‐dwelling frog from the Atlantic Rainforest of Brazil. We used Maxent to predict suitable areas based on climate and topographic profiles, combined with 21 occurrence records of bromeliads within which our target frog species (Crossodactylodes itambe) is strictly restricted. The bromeliad‐based models identified four areas potentially suitable, and subsequent surveys revealed an entirely new species of a rare bromeliad‐dwelling frog in one of the areas. We demonstrate that using easy‐to‐survey surrogate species that have a strong relationship with species that are hard to detect has enormous potential to reveal crucial information on the potential ranges and distribution of cryptic taxa. Such a surrogate modelling approach could be extended to other habitat specialist species that lack distribution data, such as amphibians restricted to specific refugia and/or reproduction sites. This could improve the targeting of surveys and increase understanding of the patterns of rarity and the drivers of species distribution, especially for areas with high endemicity and range‐restricted frogs.

Data-centric species distribution modeling: Impacts of modeler decisions in a case study of invasive European frog-bit.

Species distribution models (SDMs) are widely utilized to guide conservation decisions. The complexity of available data and SDM methodologies necessitates considerations of how data are chosen and processed for modeling to enhance model accuracy and support biological interpretations and ecological applications. We built SDMs for the invasive aquatic plant European frog-bit using aggregated and field data that span multiple scales, data sources, and data types. We tested how model results were affected by five modeler decision points: the exclusion of (1) missing and (2) correlated data and the (3) scale (large-scale aggregated data or systematic field data), (4) source (specimens or observations), and (5) type (presence-background or presence-absence) of occurrence data. Decisions about the exclusion of missing and correlated data, as well as the scale and type of occurrence data, significantly affected metrics of model performance. The source and type of occurrence data led to differences in the importance of specific explanatory variables as drivers of species distribution and predicted probability of suitable habitat. Our findings relative to European frog-bit illustrate how specific data selection and processing decisions can influence the outcomes and interpretation of SDMs. Data-centric protocols that incorporate data exploration into model building can help ensure models are reproducible and can be accurately interpreted in light of biological questions.

Guidelines for the use of spatially varying coefficients in species distribution models

AbstractAimSpecies distribution models (SDMs) are increasingly applied across macroscales using detection‐nondetection data. These models typically assume that a single set of regression coefficients can adequately describe species–environment relationships and/or population trends. However, such relationships often show nonlinear and/or spatially varying patterns that arise from complex interactions with abiotic and biotic processes that operate at different scales. Spatially varying coefficient (SVC) models can readily account for variability in the effects of environmental covariates. Yet, their use in ecology is relatively scarce due to gaps in understanding the inferential benefits that SVC models can provide compared to simpler frameworks.InnovationHere we demonstrate the inferential benefits of SVC SDMs, with a particular focus on how this approach can be used to generate and test ecological hypotheses regarding the drivers of spatial variability in population trends and species–environment relationships. We illustrate the inferential benefits of SVC SDMs with simulations and two case studies: one that assesses spatially varying trends of 51 forest bird species in the eastern United States over two decades and a second that evaluates spatial variability in the effects of five decades of land cover change on grasshopper sparrow (Ammodramus savannarum) occurrence across the continental United States.Main conclusionsWe found strong support for SVC SDMs compared to simpler alternatives in both empirical case studies. Factors operating at fine spatial scales, accounted for by the SVCs, were the primary divers of spatial variability in forest bird occurrence trends. Additionally, SVCs revealed complex species–habitat relationships with grassland and cropland area for grasshopper sparrow, providing nuanced insights into how future land use change may shape its distribution. These applications display the utility of SVC SDMs to help reveal the environmental factors that drive species distributions across both local and broad scales. We conclude by discussing the potential applications of SVC SDMs in ecology and conservation.

Accounting for dispersal and intraspecific variation in forecasts of species distribution under climate change

Abstract Dispersal and local adaptation play an important role in driving species distributions under climate change. Many studies aim to estimate relationships between species occurrences and environmental variables to predict range shift and biodiversity loss, ignoring dispersal and intraspecific variation contributing to complex spatial and temporal dynamics. We accounted for dispersal and intraspecific variation in forecasts of species distribution under climate change with species distribution models (SDMs) using two cold‐adapted, low‐dispersal Platycerus species (Coleoptera: Lucanidae), each with distinct subspecies distributions, as focus species. The results showed that the subspecies‐level model performed significantly better than the species‐level model when considering dispersal constraints in SDMs. Whether or not dispersal or intraspecific variation is accounted for, the predicted species range of Platycerus albisomni is expected to decrease in the future. For Platycerus takakuwai, accounting for dispersal constraints in SDMs indicated that its potential distribution area would increase at the subspecies level under climate change, but decrease at the species level. These divergent results show that SDMs at the subspecies level can detect impacts of climate change that may be overlooked in species‐level models. Therefore, models that consider intraspecific variation and dispersal constraints may provide a more realistic perspective on the impacts of climate change. Because accurate mapping of potential habitats is needed for conservation purposes, demographic studies should include dispersal explicitly and explore how and when intraspecific variation in dispersal affects local population dynamics. This approach could help evaluate species' habitat shifts, thus enabling suitable conservation strategies.

How substrate type and elevation drive woody communities and influence species ecological strategies in the Brazilian savanna at local scale

AbstractEnvironmental heterogeneity can influence the formation of different vegetation types, even at small spatial scales. However, little is known about the factors that drive communities at local scales and species' mechanisms to adapt to different environments. Here, we investigated the influence of local environmental conditions and spatiality on the floristic variation of woody vegetation and ecological strategies of shrub‐tree species in the Brazilian savanna. We sampled adjacent sites in three substrate types—deep soil at high altitude (in Red Oxisol, with vegetation locally called Typical Cerrado—TCRO), deep soil, sandy and at low elevation (Typical Cerrado in Entisol Quartzipsamment—TCEQ), and shallow soil with rocky outcrops at high elevation (in Entisol Lithic, with vegetation known locally as Rupestrian Cerrado—RCEL). At each environment, we inventoried the vegetation, analyzed the physical–chemical soil properties, and recorded the elevation and geographic coordinates. For the most representative species of each site and the most common one among the environments, we evaluated the functional attributes (leaf nutrient concentration, specific leaf area, and maximum population heights). We showed that substrate type and elevation explain most of the floristic variation at the local scale and that the low nutritional demand constitutes the main strategy used by the species in RCEL, while in TCEQ and TCRO, species are more efficient in resource acquisition. Our results revealed that local ecological factors and processes drive species distribution. Therefore, variations in vegetation, even at the local scale must be considered to conserve biodiversity effectively.

Drivers of Species Distribution and Niche Dynamics for Ornamental Plants Originating at Different Latitudes

Human activities provide migration opportunities for many ornamental plants and make them become a new potential invasion risk, threatening the local ecosystem. However, ornamental plants come from a wide range of sources, and there is still a lack of understanding on the distribution driving factors, ecological niche dynamics and invasion ability of ornamental plants based on the origin of different latitudes to evaluate their potential invasion risks. In this study, an ensemble of ecological niche model and a niche dynamic model were used to analyze the invasion potential of herbaceous and woody ornamental plants originating from different latitudes. The results showed that there were significant differences in environmental factors driving the distribution of plants originating from different latitudes, and climate-related factors were the primary driving force for each plant in the native and introduced regions. Urban land was the most influential factor in the introduced areas of most plants, potentially reflecting the importance of human activities in the distribution of ornamental plants. Additionally, only woody plants originating from mid-latitudes showed greater diffusivity than those originating in high latitudes and low latitudes, and the niche widths of all the herbaceous plants in the introduced regions nearly exceeded those in the native regions. This phenomenon was observed only in woody plants with mid-latitude origins. The niche similarity of all plant species between the introduced and native regions was high, indicating that all species in the introduced regions inherited niche characteristics from plants in the native regions.

Glycerol metabolism supports oral commensal interactions.

During oral biofilm development, interspecies interactions drive species distribution and biofilm architecture. To understand what molecular mechanisms determine these interactions, we used information gained from recent biogeographical investigations demonstrating an association of corynebacteria with streptococci. We previously reported that Streptococcus sanguinis and Corynebacterium durum have a close relationship through the production of membrane vesicle and fatty acids leading to S. sanguinis chain elongation and overall increased fitness supporting their commensal state. Here we present the molecular mechanisms of this interspecies interaction. Coculture experiments for transcriptomic analysis identified several differentially expressed genes in S. sanguinis. Due to its connection to fatty acid synthesis, we focused on the glycerol-operon. We further explored the differentially expressed type IV pili genes due to their connection to motility and biofilm adhesion. Gene inactivation of the glycerol kinase glpK had a profound impact on the ability of S. sanguinis to metabolize C. durum secreted glycerol and impaired chain elongation important for their interaction. Investigations on the effect of type IV pili revealed a reduction of S. sanguinis twitching motility in the presence of C. durum, which was caused by a decrease in type IV pili abundance on the surface of S. sanguinis as determined by SEM. In conclusion, we identified that the ability to metabolize C. durum produced glycerol is crucial for the interaction of C. durum and S. sanguinis. Reduced twitching motility could lead to a closer interaction of both species, supporting niche development in the oral cavity and potentially shaping symbiotic health-associated biofilm communities.

Decline in cold-water coral growth promotes molluscan diversity: A paleontological perspective from a cold-water coral mound in the western Mediterranean Sea

Framework-forming scleractinian cold-water corals (CWCs) act as ecosystem engineers, building and supporting biodiversity hotspots in the deep sea worldwide. While spatial patterns and drivers of species distributions have been evaluated on modern CWC reefs, little is known about how reef diversity is affected by habitat variability over geologic time – the scale at which CWC reefs initiate, thrive, and decline. Using three CWC reef sediment cores as species diversity archives, we investigated temporal trends of molluscan diversity over the last ~13 kyr from a CWC mound in the Alboran Sea (western Mediterranean Sea) to evaluate (a) how spatial patterns of CWC-associated diversity are recorded in reef sediments, (b) the potential of CWC reefs as biodiversity hotspots when coral growth is flourishing and when it is not, and (c) which palaeoceanographic conditions or habitat characteristics may be driving biodiversity. Our results reveal that at the ecosystem scale ecological differences between CWC habitats are more pronounced than ecological signatures of molluscan assemblages associated with intervals of CWC framework (flourishing growth) or non-framework (negligible CWC growth). However, within habitats, significant differences emerge between these assemblages with lower molluscan diversity associated with flourishing CWC growth. Significant negative correlations between molluscan diversity and palaeoceanographic conditions conducive for CWC growth (high food availability, strong hydrodynamics, optimal bottom-water temperatures and salinities, and high aggradation rates indicative of flourishing CWC growth also imply that CWC growth and relevant environmental conditions contribute to reduced molluscan diversity. Additionally, high coral volume, used here as a proxy for habitat structural complexity, is positively correlated with molluscan diversity just as high habitat complexity is in living CWC reefs. Altogether, these patterns detected over geologic time resemble those observed spatially across living CWC reefs today – where competition with resources, particularly food, prevents high reef biodiversity in the immediate vicinity of dense living CWC colonies. Overall, our study demonstrates that (1) ecological paradigms of living CWCs are preserved in their sedimentary record, (2) flourishing CWC growth and conditions promoting CWC growth drive habitat-scale diversity patterns, and (3) a geological approach can be applied to study long-term diversity dynamics in CWC ecosystems.

Responses to saltwater exposure vary across species, populations and life stages in anuran amphibians.

To predict the impacts of environmental change on species, we must first understand the factors that limit the present-day ranges of species. Most anuran amphibians cannot survive at elevated salinities, which may drive their distribution in coastal locations. Previous research showed that coastal Hyla cinerea are locally adapted to brackish habitats in North Carolina, USA. Although Hyla squirella and Hyla chrysoscelis both inhabit coastal wetlands nearby, they have not been observed in saline habitats. We take advantage of naturally occurring microgeographic variation in coastal wetland occupancy exhibited by these congeneric tree frog species to explore how salt exposure affects oviposition site choice, hatching success, early tadpole survival, plasma osmolality and tadpole body condition across coastal and inland locations. We observed higher survival among coastal H. cinerea tadpoles than inland H. cinerea, which corroborates previous findings. But contrary to expectations, coastal H. cinerea had lower survival than H. squirella and H. chrysoscelis, indicating that all three species may be able to persist in saline wetlands. We also observed differences in tadpole plasma osmolality across species, locations and salinities, but these differences were not associated with survival rates in salt water. Instead, coastal occupancy may be affected by stage-specific processes like higher probability of total clutch loss as shown by inland H. chrysoscelis or maladaptive egg deposition patterns as shown by inland H. squirella. Although we expected salt water to be the primary filter driving species distributions along a coastal salinity gradient, it is likely that the factors dictating anuran ranges along the coast involve stage-, species- and location-specific processes that are mediated by ecological processes and life history traits.

Multiple genetic sources facilitate the northward range expansion of an intertidal oyster along China's coast.

Coastal artificial structures on the former mudflats provide available habitats for the rocky intertidal species which can establish new populations in these emerging habitats over their former distribution range limits. As a former southern species, the oyster Crassostrea sikamea has become a pioneer and rapidly invaded the artificial shorelines in northern China. We used a seascape genomics approach to investigate the population structure and genetic sources of C. sikamea on the coastal artificial structures, which is crucial for understanding the genetic mechanisms driving species distribution range expansion and invasion pathway of intertidal species. Five C. sikamea populations, including two artificial substrate populations (WGZ and ZAP), one oyster reef population (LS), and two natural rocky shore populations (ZS and XM), were measured using single nucleotide polymorphism (SNPs) obtained from double digest restriction-site associated DNA sequencing (ddRAD-Seq). Redundancy analyses (RDA) were implemented for investigating the relationship between local temperature variables and the temperature adaptability of C. sikamea. Genetic diversity, direction and strength of gene flow, and population structure all revealed that the LS and ZS populations were the genetic sources for the oyster populations on the emerging northern coastal artificial structures. Results of RDA showed that there were different adaptive potentials for northern and southern populations to local temperature variables and the oyster reef population which frequently suffers from heat stress owned high heat adaptability. The ZS population as a genetic source nearby the Yangtze River estuary provided mass larvae for the northern populations, and the other genetic source, the heat-tolerant LS population, in the oyster reef played an important role in the post-settlement success by providing preadapted genotypes. These results highlight the importance of multiple sources with divergent adaptative capabilities for biological invasion, and also emphasize the importance of the oyster reef in coastal biodiversity and conservation.

Species assemblage networks identify regional connectivity pathways among hydrothermal vents in the Northwest Pacific.

The distribution of species among spatially isolated habitat patches supports regional biodiversity and stability, so understanding the underlying processes and structure is a key target of conservation. Although multivariate statistics can infer the connectivity processes driving species distribution, such as dispersal and habitat suitability, they rarely explore the structure. Methods from graph theory, applied to distribution data, give insights into both connectivity pathways and processes by intuitively formatting the data as a network of habitat patches. We apply these methods to empirical data from the hydrothermal vent habitats of the Northwest Pacific. Hydrothermal vents are "oases" of biological productivity and endemicity on the seafloor that are imminently threatened by anthropogenic disturbances with unknown consequences to biodiversity. Here, we describe the structure of species assemblage networks at hydrothermal vents, how local and regional parameters affect their structure, and the implications for conservation. Two complementary networks were formed from an extensive species assemblage dataset: a similarity network of vent site nodes linked by weighted edges based on their pairwise assemblage similarity and a bipartite network of species nodes linked to vent site nodes at which they are present. Using these networks, we assessed the role of individual vent sites in maintaining network connectivity and identified biogeographic sub-regions. The three sub-regions and two outlying sites are separated by their spatial arrangement and local environmental filters. Both networks detected vent sites that play a disproportionately important role in regional pathways, while the bipartite network also identified key vent sites maintaining the distinct species assemblages of their sub-regions. These regional connectivity pathways provide insights into historical colonization routes, while sub-regional connectivity pathways are of value when selecting sites for conservation and/or estimating the multivent impacts from proposed deep-sea mining.

Recommendations for quantifying and reducing uncertainty in climate projections of species distributions.

Projecting the future distributions of commercially and ecologically important species has become a critical approach for ecosystem managers to strategically anticipate change, but large uncertainties in projections limit climate adaptation planning. Although distribution projections are primarily used to understand the scope of potential change-rather than accurately predict specific outcomes-it is nonetheless essential to understand where and why projections can give implausible results and to identify which processes contribute to uncertainty. Here, we use a series of simulated species distributions, an ensemble of 252 species distribution models, and an ensemble of three regional ocean climate projections, to isolate the influences of uncertainty from earth system model spread and from ecological modeling. The simulations encompass marine species with different functional traits and ecological preferences to more broadly address resource manager and fishery stakeholder needs, and provide a simulated true state with which to evaluate projections. We present our results relative to the degree of environmental extrapolation from historical conditions, which helps facilitate interpretation by ecological modelers working in diverse systems. We found uncertainty associated with species distribution models can exceed uncertainty generated from diverging earth system models (up to 70% of total uncertainty by 2100), and that this result was consistent across species traits. Species distribution model uncertainty increased through time and was primarily related to the degree to which models extrapolated into novel environmental conditions but moderated by how well models captured the underlying dynamics driving species distributions. The predictive power of simulated species distribution models remained relatively high in the first 30 years of projections, in alignment with the time period in which stakeholders make strategic decisions based on climate information. By understanding sources of uncertainty, and how they change at different forecast horizons, we provide recommendations for projecting species distribution models under global climate change.

Standardized NEON organismal data for biodiversity research

AbstractUnderstanding patterns and drivers of species distribution and abundance, and thus biodiversity, is a core goal of ecology. Despite advances in recent decades, research into these patterns and processes is currently limited by a lack of standardized, high‐quality, empirical data that span large spatial scales and long time periods. The NEON fills this gap by providing freely available observational data that are generated during robust and consistent organismal sampling of several sentinel taxonomic groups within 81 sites distributed across the United States and will be collected for at least 30 years. The breadth and scope of these data provide a unique resource for advancing biodiversity research. To maximize the potential of this opportunity, however, it is critical that NEON data be maximally accessible and easily integrated into investigators' workflows and analyses. To facilitate its use for biodiversity research and synthesis, we created a workflow to process and format NEON organismal data into the ecocomDP (ecological community data design pattern) format that were available through the ecocomDP R package; we then provided the standardized data as an R data package (neonDivData). We briefly summarize sampling designs and data wrangling decisions for the major taxonomic groups included in this effort. Our workflows are open‐source so the biodiversity community may: add additional taxonomic groups; modify the workflow to produce datasets appropriate for their own analytical needs; and regularly update the data packages as more observations become available. Finally, we provide two simple examples of how the standardized data may be used for biodiversity research. By providing a standardized data package, we hope to enhance the utility of NEON organismal data in advancing biodiversity research and encourage the use of the harmonized ecocomDP data design pattern for community ecology data from other ecological observatory networks.

Differential nutrient limitation and tree height control leaf physiology, supporting niche partitioning in tropical dipterocarp forests

Abstract Revealing the mechanisms of environmental niche partitioning within lowland tropical forests is important for understanding the drivers of current species distributions and potential vulnerability to environmental change. Tropical forest structure and species composition change across edaphic gradients in Borneo over short distances. However, our understanding of how edaphic conditions affect tree physiology and whether these relationships drive niche partitioning within Bornean forests remains incomplete. This study evaluated how leaf physiological function changes with nutrient availability across a fine‐scale edaphic gradient and whether these relationships vary according to tree height. Furthermore, we tested whether intraspecific leaf trait variation allows generalist species to populate a wider range of environments. We measured leaf traits of 218 trees ranging in height from 4 to 66 m from 13 dipterocarp species within four tropical forest types (alluvial, mudstone, sandstone and kerangas) occurring along an <5 km edaphic gradient in North Borneo. The traits measured included saturating photosynthesis (Asat), maximum photosynthetic capacity (Vcmax), leaf dark respiration (Rleaf), leaf mass per area (LMA), leaf thickness, minimum stomatal conductance (gdark) and leaf nutrient concentrations (N, P, Ca, K and Mg). Across all species, leaf traits varied consistently in response to soil nutrient availability across forest types except Rleaf_mass, [Mg]leaf and [Ca]leaf. Changes in photosynthesis and respiration rates were related to different leaf nutrients across forest types, with greater nutrient‐use efficiency in more nutrient‐poor environments. Generalist species partially or fully compensated reductions in mass‐based photosynthesis through increasing LMA in more nutrient‐poor environments. Leaf traits also varied with tree height, except Vcmax_mass, but only in response to height‐related modifications of leaf morphology (LMA and leaf thickness). These height–trait relationships did not vary across the edaphic gradient, except for Asat, [N]leaf, [P]leaf and [K]leaf. Our results highlight that modification of leaf physiological function and morphology act as important adaptations for Bornean dipterocarps in response to edaphic and vertical environmental gradients. Meanwhile, multiple nutrients appear to contribute to niche partitioning and could drive species distributions and high biodiversity within Bornean forest landscapes. Read the free Plain Language Summary for this article on the Journal blog.

Larix species range dynamics in Siberia since the Last Glacial captured from sedimentary ancient DNA

Climate change is expected to cause major shifts in boreal forests which are in vast areas of Siberia dominated by two species of the deciduous needle tree larch (Larix). The species differ markedly in their ecosystem functions, thus shifts in their respective ranges are of global relevance. However, drivers of species distribution are not well understood, in part because paleoecological data at species level are lacking. This study tracks Larix species distribution in time and space using target enrichment on sedimentary ancient DNA extracts from eight lakes across Siberia. We discovered that Larix sibirica, presently dominating in western Siberia, likely migrated to its northern distribution area only in the Holocene at around 10,000 years before present (ka BP), and had a much wider eastern distribution around 33 ka BP. Samples dated to the Last Glacial Maximum (around 21 ka BP), consistently show genotypes of L. gmelinii. Our results suggest climate as a strong determinant of species distribution in Larix and provide temporal and spatial data for species projection in a changing climate.

Topographic depressions provide potential microrefugia for ground-dwelling arthropods

Species can survive periods of unfavorable conditions in small areas that are protected from climate-related disturbances, such as increasing temperature and severe drought. These areas are known as “microrefugia” and are increasingly recognized by conservationists. Although some studies suggest that the in situ survival of invertebrate species may be mediated by topographically complex environments, there is little information about the main environmental factors that drive species distributions within such areas. Here, we investigate the spatial patterns and species trait composition (moisture preference, body size, dispersal capacity, and feeding habit) of five groups of ground-dwelling arthropods—spiders, woodlice, ants, ground beetles, and rove beetles—in topographic depressions (i.e., “solution dolines” or “solution sinkholes”) and on the surrounding plateau within a forested karst landscape and analyze the microhabitat conditions that affect these arthropod assemblages. We found that dolines have the capacity to maintain characteristic arthropod assemblages—including species that may be particularly vulnerable to climate change (e.g., species associated with moist habitats)—and thus, they may contribute to the landscape-scale biodiversity of karst landscapes. We also found that doline bottoms have the potential to maintain permanently moist conditions not only in spring and autumn but also during drier periods of the year. This ability of dolines may indicate the presence of potential hydrologic microrefugia. Furthermore, dolines displayed specific sets of species traits (e.g., more small-bodied spiders, more carnivorous ground beetles, and more rove beetles with high dispersal capacity occurred in dolines than on the plateau), highlighting that they may facilitate the persistence of some species and traits that are purged from the surrounding landscape. Future studies may reveal the long-term ecological consequences of different climatic and anthropogenic factors on the distribution and functional traits of arthropod taxa within microrefugia and on the refugial capacity of these safe havens under a warming climate.