Oceanographic connectivity explains the intra-specific diversity of mangrove forests at global scales.

The distribution of mangrove intra-specific biodiversity can be structured by historical demographic processes that enhance or limit effective population sizes. Oceanographic connectivity (OC) may further structure intra-specific biodiversity by preserving or diluting the genetic signatures of historical changes. Despite its relevance for biogeography and evolution, the role of oceanographic connectivity in structuring the distribution of mangrove's genetic diversity has not been addressed at global scale. Here we ask whether connectivity mediated by ocean currents explains the intra-specific diversity of mangroves. A comprehensive dataset of population genetic differentiation was compiled from the literature. Multigenerational connectivity and population centrality indices were estimated with biophysical modeling coupled with network analyses. The variability explained in genetic differentiation was tested with competitive regression models built upon classical isolation-by-distance (IBD) models considering geographic distance. We show that oceanographic connectivity can explain the genetic differentiation of mangrove populations regardless of the species, region, and genetic marker (significant regression models in 95% of cases, with an average R-square of 0.44 ± 0.23 and Person's correlation of 0.65 ± 0.17), systematically improving IBD models. Centrality indices, providing information on important stepping-stone sites between biogeographic regions, were also important in explaining differentiation (R-square improvement of 0.06 ± 0.07, up to 0.42). We further show that ocean currents produce skewed dispersal kernels for mangroves, highlighting the role of rare long-distance dispersal events responsible for historical settlements. Overall, we demonstrate the role of oceanographic connectivity in structuring mangrove intra-specific diversity. Our findings are critical for mangroves' biogeography and evolution, but also for management strategies considering climate change and genetic biodiversity conservation.

Exposure of protected areas in Central America to extreme weather events.

Central America and the Caribbean are regularly battered by megadroughts, heavy rainfall, heat waves, and tropical cyclones. Although 21st-century climate change is expected to increase the frequency, intensity, and duration of these extreme weather events (EWEs), their incidence in regional protected areas (PAs) remains poorly explored. We examined historical and projected EWEs across the region based on 32 metrics that describe distinct dimensions (i.e., intensity, duration, and frequency) of heat waves, cyclones, droughts, and rainfall and compared trends in PAs with trends in unprotected lands. From the early 21st century onward, exposure to EWEs increased across the region, and PAs were predicted to be more exposed to climate extremes than unprotected areas (as shown by autoregressive model coefficients at p < 0.05 significance level). This was particularly true for heat waves, which were projected to have a significantly higher average (tested by Wilcoxon tests at p < 0.01) intensity and duration, and tropical cyclones, which affected PAs more severely in carbon-intensive scenarios. PAs were also predicted to be significantly less exposed to droughts and heavy rainfall than unprotected areas (tested by Wilcoxon tests at p < 0.01). However, droughts that could threaten connectivity between PAs are increasingly common in this region. We estimated that approximately 65% of the study area will experience at least one drought episode that is more intense and longer lasting than previous droughts. Collectively, our results highlight that new conservation strategies adapted to threats associated with EWEs need to be tailored and implemented promptly. Unless urgent action is taken, significant damage may be inflicted on the unique biodiversity of the region.

Ciclones, olas de calor, sequías y lluvias intensas son eventos comunes en Centroamérica y el Caribe, cuya frecuencia, intensidad y duración se espera aumente durante el siglo XXI a causa del cambio climático. Sin embargo, en la actualidad, se desconoce cuál será la incidencia de estos eventos meteorológicos extremos (EME) dentro de las áreas protegidas. En este estudio examinamos la exposición histórica y futura a los extremos climáticos y comparamos el grado de exposición dentro y fuera de las áreas protegidas de toda la región por medio de 32 métricas que describen distintas dimensiones (intensidad, duración y frecuencia) de las olas de calor, los ciclones, las sequías y las precipitaciones. Los resultados indican que a medida que aumente el número de EME, las áreas protegidas estarán más expuestas a los extremos climáticos que las áreas no protegidas. Esto es especialmente cierto en el caso de las olas de calor, que, según las proyecciones, tendrán una intensidad y una duración medias significativamente mayores, y de los ciclones tropicales, que afectarán más gravemente a las zonas protegidas en los escenarios intensivos en carbono. Nuestros resultados también indican que las zonas protegidas estarán significativamente menos expuestas a sequías o lluvias torrenciales que las zonas no protegidas. Sin embargo, las sequías que podrían amenazar la conectividad entre áreas protegidas son cada vez más frecuentes en esta región. Se estima que aproximadamente el 65% del área de estudio experimentará al menos un episodio de sequía más intenso y duradero que las sequías anteriores. En conjunto, nuestros resultados ponen de relieve la necesidad de diseñar y aplicar con prontitud nuevas estrategias de conservación adaptadas a las amenazas asociadas a los EWE. A menos que se tomen medidas urgentes, la biodiversidad única de la región podría sufrir daños considerables.

Climate correlates of bluetongue incidence in southern Portugal.

Model forecasts of the spatiotemporal occurrence dynamics of diseases are necessary and can help understand and thus manage future disease outbreaks. In our study, we used ecological niche modelling to assess the impact of climate on the vector suitability for bluetongue disease, a disease affecting livestock production with important economic consequences. Specifically, we investigated the relationship between the occurrence of bluetongue outbreaks and the environmental suitability of each of the four vector species studied. We found that the main vector for bluetongue disease, Culicoides imicola, a typically tropical and subtropical species, was a strong predictor for disease outbreak occurrence in a region of southern Portugal from 2004 to 2021. The results highlight the importance of understanding the climatic factors that might influence vector presence to help manage infectious disease impacts. When diseases impact economically relevant species, the impacts go beyond mortality and have important economic consequences.

Response of an Afro-Palearctic bird migrant to glaciation cycles.

Migration allows animals to exploit spatially separated and seasonally available resources at a continental to global scale. However, responding to global climatic changes might prove challenging, especially for long-distance intercontinental migrants. During glacial periods, when conditions became too harsh for breeding in the north, avian migrants have been hypothesized to retract their distribution to reside within small refugial areas. Here, we present data showing that an Afro-Palearctic migrant continued seasonal migration, largely within Africa, during previous glacial-interglacial cycles with no obvious impact on population size. Using individual migratory track data to hindcast monthly bioclimatic habitat availability maps through the last 120,000 y, we show altered seasonal use of suitable areas through time. Independently derived effective population sizes indicate a growing population through the last 40,000 y. We conclude that the migratory lifestyle enabled adaptation to shifting climate conditions. This indicates that populations of resource-tracking, long-distance migratory species could expand successfully during warming periods in the past, which could also be the case under future climate scenarios.

Human disturbances affect the topology of food webs.

Networks describe nodes connected by links, with numbers of links per node, the degree, forming a range of distributions including random and scale-free. How network topologies emerge in natural systems still puzzles scientists. Based on previous theoretical simulations, we predict that scale-free food webs are favourably selected by random disturbances while random food webs are selected by targeted disturbances. We assume that lower human pressures are more likely associated with random disturbances, whereas higher pressures are associated with targeted ones. We examine these predictions using 351 empirical food webs, generally confirming our predictions. Should the topology of food webs respond to changes in the magnitude of disturbances in a predictable fashion, consistently across ecosystems and scales of organisation, it would provide a baseline expectation to understand and predict the consequences of human pressures on ecosystem dynamics.

Potential for invasion of traded birds under climate and land-cover change.

Humans have moved species away from their native ranges since the Neolithic, but globalization accelerated the rate at which species are being moved. We fitted more than half million distribution models for 610 traded bird species on the CITES list to examine the separate and joint effects of global climate and land-cover change on their potential end-of-century distributions. We found that climate-induced suitability for modelled invasive species increases with latitude, because traded birds are mainly of tropical origin and much of the temperate region is 'tropicalizing.' Conversely, the tropics are becoming more arid, thus limiting the potential from cross-continental invasion by tropical species. This trend is compounded by forest loss around the tropics since most traded birds are forest dwellers. In contrast, net gains in forest area across the temperate region could compound climate change effects and increase the potential for colonization of low-latitude birds. Climate change has always led to regional redistributions of species, but the combination of human transportation, climate, and land-cover changes will likely accelerate the redistribution of species globally, increasing chances of alien species successfully invading non-native lands. Such process of biodiversity homogenization can lead to emergence of non-analogue communities with unknown environmental and socioeconomic consequences.

Exploring the Effects of Geopolitical Shifts on Global Wildlife Trade.

International wildlife trade is a major driver of species extinction and biological invasions. Anticipating environmental risks requires inferences about trade patterns, which are shaped by geopolitics. Although the future cannot be predicted, scenarios can help deal with the uncertainty of future geopolitical dynamics. We propose a framework for generating and analyzing scenarios based on four geopolitical storylines, distinguished by combinations of international trade barrier strength and domestic law enforcement degree across countries supplying and demanding wildlife. We then use historical data on bird trade to classify countries into geopolitical profiles and confirm that trade barriers and law enforcement allow predicting bird trade patterns, supporting our scenarios' plausibility and enabling projections for future global bird trade. Our framework can be used to examine the consequences of geopolitical changes for wildlife trade and to advise policy and legislation. Reducing demand for wildlife and ameliorating global inequality are key for curbing trade related risks.

Discriminating climate, land-cover and random effects on species range dynamics.

Species are reportedly shifting their distributions poleward and upward in several parts of the world in response to climate change. The extent to which other factors might play a role driving these changes is still unclear. Land-cover change is a major cause of distributional changes, but it cannot be discarded that distributional dynamics might be at times caused by other mechanisms (e.g. dispersal, ecological drift). Using observed changes in the distribution of 82 breeding birds in Great Britain between three time periods 1968-72 (t1 ), 1988-91 (t2 ) and 2007-2011 (t3 ), we examine whether observed bird range shifts between t1 -t2 and t1 -t3 are best explained by climate change or land-cover change, or whether they are not distinguishable from what would be expected by chance. We found that range shifts across the rear edge of northerly distributed species in Great Britain are best explained by climate change, while shifts across the leading edge of southerly distributed species are best explained by changes in land-cover. In contrast, at the northern and southern edges of Great Britain, range dynamics could not be distinguished from that expected by chance. The latter observation could be a consequence of boundary effects limiting the direction and magnitude of range changes, stochastic demographic mechanisms neither associated with climate nor land-cover change or with complex interactions among factors. Our results reinforce the view that comprehensive assessments of climate change effects on species range shifts need to examine alternative drivers of change on equal footing and that null models can help assess whether observed patterns could have arisen by chance alone.

Thermal tolerance and the importance of microhabitats for Andean frogs in the context of land use and climate change.

Global warming is having impacts across the Tree of Life. Understanding species' physiological sensitivity to temperature change and how they relate to local temperature variation in their habitats is crucial to determining vulnerability to global warming. We ask how species' vulnerability varies across habitats and elevations, and how climatically buffered microhabitats can contribute to reduce their vulnerability. We measured thermal sensitivity (critical thermal maximum-CTmax ) of 14 species of Pristimantis frogs inhabiting young and old secondary, and primary forests in the Colombian Andes. Exposure to temperature stress was measured by recording temperature in the understorey and across five microhabitats. We determined frogs' current vulnerability across habitats, elevations and microhabitats accounting for phylogeny and then ask how vulnerability varies under four warming scenarios: +1.5, +2, +3 and +5°C. We found that CTmax was constant across species regardless of habitat and elevation. However, species in young secondary forests are expected to become more vulnerable because of increased exposure to higher temperatures. Microhabitat variation could enable species to persist within their thermal temperature range as long as regional temperatures do not surpass +2°C. The effectiveness of microhabitat buffering decreases with a 2-3°C increase, and is almost null under a 5°C temperature increase. Microhabitats will provide thermal protection to Andean frog communities from climate change by enabling tracking of suitable climates through short distance movement. Conservation strategies, such as managing landscapes by preserving primary forests and allowing regrowth and reconnection of secondary forest would offer thermally buffered microhabitats and aid in the survival of this group.

Para determinar la vulnerabilidad de las especies al calentamiento global es indispensable considerar la tolerancia fisiológica de las especies al cambio de temperatura y las condiciones ambientales a las que están expuestas. En este estudio exploramos la vulnerabilidad de especies a través de diferentes hábitats y altitudes y examinamos si ciertos microhábitats contribuyen a reducir la vulnerabilidad al calentamiento global. Medimos la tolerancia térmica (CTmax ) de catorce especies de ranas Pristimantis en bosques secundarios jóvenes y viejos, y bosques primarios en los Andes tropicales. Registramos la temperatura a la que estas especies están expuestas en el sotobosque así como dentro de cinco microhábitats. Usando CTmax y las temperaturas a las que están expuestas, determinamos la vulnerabilidad de las especies en diferentes hábitats, elevaciones y microhábitats. También preguntamos cómo cambiará esta vulnerabilidad si la temperatura incrementa: 1.5°C, 2°C, 3°C y 5°C. CTmax fue constante en todos los hábitats y elevaciones. Las especies de bosques secundarios jóvenes son más vulnerables pues están expuestas a temperaturas más altas. Al utilizar microhábitats, las especies estarán protegidas si el aumento de temperatura no supera los + 2°C. Todos los microhábitats seguirán proporcionando refugio térmico si la temperatura aumenta 1.5°C, pero esta protección térmica disminuirá si la temperatura aumenta 2-3°C y será casi nula con un aumento de temperatura de 5°C. Los microhábitats proporcionarán protección térmica a la comunidad de ranas de los Andes contra el cambio climático. Estrategias de conservación, como la regeneración natural y la reconexión de bosques secundarios y la preservación de bosques primarios, ayudaría a la supervivencia de las ranas al tener microhábitats que ofrecen refugio térmico.

Intraspecific variation in lizard heat tolerance alters estimates of climate impact.

Research addressing the effects of global warming on the distribution and persistence of species generally assumes that population variation in thermal tolerance is spatially constant or overridden by interspecific variation. Typically, this rationale is implicit in sourcing one critical thermal maximum (CTmax ) population estimate per species to model spatiotemporal cross-taxa variation in heat tolerance. Theory suggests that such an approach could result in biased or imprecise estimates and forecasts of impact from climate warming, but limited empirical evidence in support of those expectations exists. We experimentally quantify the magnitude of intraspecific variation in CTmax among lizard populations, and the extent to which incorporating such variability can alter estimates of climate impact through a biophysical model. To do so, we measured CTmax from 59 populations of 15 Iberian lizard species (304 individuals). The overall median CTmax across all individuals from all species was 42.8°C and ranged from 40.5 to 48.3°C, with species medians decreasing through xeric, climate-generalist and mesic taxa. We found strong statistical support for intraspecific differentiation in CTmax by up to a median of 3°C among populations. We show that annual restricted activity (operative temperature > CTmax ) over the Iberian distribution of our study species differs by a median of >80 hr per 25-km2 grid cell based on different population-level CTmax estimates. This discrepancy leads to predictions of spatial variation in annual restricted activity to change by more than 20 days for six of the study species. Considering that during restriction periods, reptiles should be unable to feed and reproduce, current projections of climate-change impacts on the fitness of ectotherm fauna could be under- or over-estimated depending on which population is chosen to represent the physiological spectra of the species in question. Mapping heat tolerance over the full geographical ranges of single species is thus critical to address cross-taxa patterns and drivers of heat tolerance in a biologically comprehensive way.

The effect of multiple biotic interaction types on species persistence.

No species can persist in isolation from other species, but how biotic interactions affect species persistence is still a matter of inquiry. Is persistence more likely in communities with higher proportion of competing species, or in communities with more positive interactions? How do different components of community structure mediate this relationship? We address these questions using a novel simulation framework that generates realistic communities with varying numbers of species and different proportions of biotic interaction types within and across trophic levels. We show that when communities have fewer species, persistence is more likely if positive interactions-such as mutualism and commensalism-are prevalent. In species-rich communities, the disproportionate effect of positive interactions on persistence is diluted and different combinations of biotic interaction types can coexist without affecting persistence significantly. We present the first theoretical examination of how multiple-interaction networks with varying architectures relate to local species persistence, and provide insight about the underlying causes of stability in communities.

Projected climate changes threaten ancient refugia of kelp forests in the North Atlantic.

Intraspecific genetic variability is critical for species adaptation and evolution and yet it is generally overlooked in projections of the biological consequences of climate change. We ask whether ongoing climate changes can cause the loss of important gene pools from North Atlantic relict kelp forests that persisted over glacial-interglacial cycles. We use ecological niche modelling to predict genetic diversity hotspots for eight species of large brown algae with different thermal tolerances (Arctic to warm temperate), estimated as regions of persistence throughout the Last Glacial Maximum (20,000 YBP), the warmer Mid-Holocene (6,000 YBP), and the present. Changes in the genetic diversity within ancient refugia were projected for the future (year 2100) under two contrasting climate change scenarios (RCP2.6 and RCP8.5). Models predicted distributions that matched empirical distributions in cross-validation, and identified distinct refugia at the low latitude ranges, which largely coincide among species with similar ecological niches. Transferred models into the future projected polewards expansions and substantial range losses in lower latitudes, where richer gene pools are expected (in Nova Scotia and Iberia for cold affinity species and Gibraltar, Alboran, and Morocco for warm-temperate species). These effects were projected for both scenarios but were intensified under the extreme RCP8.5 scenario, with the complete borealization (circum-Arctic colonization) of kelp forests, the redistribution of the biogeographical transitional zones of the North Atlantic, and the erosion of global gene pools across all species. As the geographic distribution of genetic variability is unknown for most marine species, our results represent a baseline for identification of locations potentially rich in unique phylogeographic lineages that are also climatic relics in threat of disappearing.

How complex should models be? Comparing correlative and mechanistic range dynamics models.

Criticism has been levelled at climate-change-induced forecasts of species range shifts that do not account explicitly for complex population dynamics. The relative importance of such dynamics under climate change is, however, undetermined because direct tests comparing the performance of demographic models vs. simpler ecological niche models are still lacking owing to difficulties in evaluating forecasts using real-world data. We provide the first comparison of the skill of coupled ecological-niche-population models and ecological niche models in predicting documented shifts in the ranges of 20 British breeding bird species across a 40-year period. Forecasts from models calibrated with data centred on 1970 were evaluated using data centred on 2010. We found that more complex coupled ecological-niche-population models (that account for dispersal and metapopulation dynamics) tend to have higher predictive accuracy in forecasting species range shifts than structurally simpler models that only account for variation in climate. However, these better forecasts are achieved only if ecological responses to climate change are simulated without static snapshots of historic land use, taken at a single point in time. In contrast, including both static land use and dynamic climate variables in simpler ecological niche models improve forecasts of observed range shifts. Despite being less skilful at predicting range changes at the grid-cell level, ecological niche models do as well, or better, than more complex models at predicting the magnitude of relative change in range size. Therefore, ecological niche models can provide a reasonable first approximation of the magnitude of species' potential range shifts, especially when more detailed data are lacking on dispersal dynamics, demographic processes underpinning population performance, and change in land cover.

Additive threats from pathogens, climate and land-use change for global amphibian diversity.

Amphibian population declines far exceed those of other vertebrate groups, with 30% of all species listed as threatened by the International Union for Conservation of Nature. The causes of these declines are a matter of continued research, but probably include climate change, land-use change and spread of the pathogenic fungal disease chytridiomycosis. Here we assess the spatial distribution and interactions of these primary threats in relation to the global distribution of amphibian species. We show that the greatest proportions of species negatively affected by climate change are projected to be found in Africa, parts of northern South America and the Andes. Regions with the highest projected impact of land-use and climate change coincide, but there is little spatial overlap with regions highly threatened by the fungal disease. Overall, the areas harbouring the richest amphibian faunas are disproportionately more affected by one or multiple threat factors than areas with low richness. Amphibian declines are likely to accelerate in the twenty-first century, because multiple drivers of extinction could jeopardize their populations more than previous, mono-causal, assessments have suggested.

Consequences of climate change on the tree of life in Europe.

Many species are projected to become vulnerable to twenty-first-century climate changes, with consequent effects on the tree of life. If losses were not randomly distributed across the tree of life, climate change could lead to a disproportionate loss of evolutionary history. Here we estimate the consequences of climate change on the phylogenetic diversities of plant, bird and mammal assemblages across Europe. Using a consensus across ensembles of forecasts for 2020, 2050 and 2080 and high-resolution phylogenetic trees, we show that species vulnerability to climate change clusters weakly across phylogenies. Such phylogenetic signal in species vulnerabilities does not lead to higher loss of evolutionary history than expected with a model of random extinctions. This is because vulnerable species have neither fewer nor closer relatives than the remaining clades. Reductions in phylogenetic diversity will be greater in southern Europe, and gains are expected in regions of high latitude or altitude. However, losses will not be offset by gains and the tree of life faces a trend towards homogenization across the continent.

The effects of phenotypic plasticity and local adaptation on forecasts of species range shifts under climate change.

Species are the unit of analysis in many global change and conservation biology studies; however, species are not uniform entities but are composed of different, sometimes locally adapted, populations differing in plasticity. We examined how intraspecific variation in thermal niches and phenotypic plasticity will affect species distributions in a warming climate. We first developed a conceptual model linking plasticity and niche breadth, providing five alternative intraspecific scenarios that are consistent with existing literature. Secondly, we used ecological niche-modeling techniques to quantify the impact of each intraspecific scenario on the distribution of a virtual species across a geographically realistic setting. Finally, we performed an analogous modeling exercise using real data on the climatic niches of different tree provenances. We show that when population differentiation is accounted for and dispersal is restricted, forecasts of species range shifts under climate change are even more pessimistic than those using the conventional assumption of homogeneously high plasticity across a species' range. Suitable population-level data are not available for most species so identifying general patterns of population differentiation could fill this gap. However, the literature review revealed contrasting patterns among species, urging greater levels of integration among empirical, modeling and theoretical research on intraspecific phenotypic variation.

Predictors of contraction and expansion of area of occupancy for British birds.

Geographical range dynamics are driven by the joint effects of abiotic factors, human ecosystem modifications, biotic interactions and the intrinsic organismal responses to these. However, the relative contribution of each component remains largely unknown. Here, we compare the contribution of life-history attributes, broad-scale gradients in climate and geographical context of species' historical ranges, as predictors of recent changes in area of occupancy for 116 terrestrial British breeding birds (74 contractors, 42 expanders) between the early 1970s and late 1990 s. Regional threat classifications demonstrated that the species of highest conservation concern showed both the largest contractions and the smallest expansions. Species responded differently to climate depending on geographical distribution-northern species changed their area of occupancy (expansion or contraction) more in warmer and drier regions, whereas southern species changed more in colder and wetter environments. Species with slow life history (larger body size) tended to have a lower probability of changing their area of occupancy than species with faster life history, whereas species with greater natal dispersal capacity resisted contraction and, counterintuitively, expansion. Higher geographical fragmentation of species' range also increased expansion probability, possibly indicating a release from a previously limiting condition, for example through agricultural abandonment since the 1970s. After accounting statistically for the complexity and nonlinearity of the data, our results demonstrate two key aspects of changing area of occupancy for British birds: (i) climate is the dominant driver of change, but direction of effect depends on geographical context, and (ii) all of our predictors generally had a similar effect regardless of the direction of the change (contraction versus expansion). Although we caution applying results from Britain's highly modified and well-studied bird community to other biogeographic regions, our results do indicate that a species' propensity to change area of occupancy over decadal scales can be explained partially by a combination of simple allometric predictors of life-history pace, average climate conditions and geographical context.

Integrating multiple lines of evidence into historical biogeography hypothesis testing: a Bison bison case study.

One of the grand goals of historical biogeography is to understand how and why species' population sizes and distributions change over time. Multiple types of data drawn from disparate fields, combined into a single modelling framework, are necessary to document changes in a species's demography and distribution, and to determine the drivers responsible for change. Yet truly integrated approaches are challenging and rarely performed. Here, we discuss a modelling framework that integrates spatio-temporal fossil data, ancient DNA, palaeoclimatological reconstructions, bioclimatic envelope modelling and coalescence models in order to statistically test alternative hypotheses of demographic and potential distributional changes for the iconic American bison (Bison bison). Using different assumptions about the evolution of the bioclimatic niche, we generate hypothetical distributional and demographic histories of the species. We then test these demographic models by comparing the genetic signature predicted by serial coalescence against sequence data derived from subfossils and modern populations. Our results supported demographic models that include both climate and human-associated drivers of population declines. This synthetic approach, integrating palaeoclimatology, bioclimatic envelopes, serial coalescence, spatio-temporal fossil data and heterochronous DNA sequences, improves understanding of species' historical biogeography by allowing consideration of both abiotic and biotic interactions at the population level.

Data error propagation in stacked bioclimatic envelope models.

Stacking is the process of overlaying inferred species potential distributions for multiple species based on outputs of bioclimatic envelope models (BEMs). The approach can be used to investigate patterns and processes of species richness. If data limitations on individual species distributions are inevitable, but how do they affect inferences of patterns and processes of species richness? We investigate the influence of different data sources on estimated species richness gradients in China. We fitted BEMs using species distributions data for 334 bird species obtained from (1) global range maps, (2) regional checklists, (3) museum records and surveys, and (4) citizen science data using presence-only (Mahalanobis distance), presence-background (MAXENT), and presence-absence (GAM and BRT) BEMs. Individual species predictions were stacked to generate species richness gradients. Here, we show that different data sources and BEMs can generate spatially varying gradients of species richness. The environmental predictors that best explained species distributions also differed between data sources. Models using citizen-based data had the highest accuracy, whereas those using range data had the lowest accuracy. Potential richness patterns estimated by GAM and BRT models were robust to data uncertainty. When multiple data sets exist for the same region and taxa, we advise that explicit treatments of uncertainty, such as sensitivity analyses of the input data, should be conducted during the process of modeling.

Power laws in species' biotic interaction networks can be inferred from co-occurrence data.

Inferring biotic interactions from species co-occurrence patterns has long intrigued ecologists. Yet recent research revealed that co-occurrences may not reliably represent pairwise biotic interactions. We propose that examining network-level co-occurrence patterns can provide valuable insights into community structure and assembly. Analysing ten bipartite networks of empirically sampled biotic interactions and associated species spatial distribution, we find that approximately 20% of co-occurrences correspond to actual interactions. Moreover, the degree distribution shifts from exponential in co-occurrence networks to power laws in networks of biotic interactions. This shift results from a strong interplay between species' biotic (their interacting partners) and abiotic (their environmental requirements) niches, and is accurately predicted by considering co-occurrence frequencies. Our work offers a mechanistic understanding of the assembly of ecological communities and suggests simple ways to infer fundamental biotic interaction network characteristics from co-occurrence data.