Occurrence-based diversity estimation reveals macroecological and conservation knowledge gaps for global woody plants.

Incomplete sampling of species' geographic distributions has challenged biogeographers for many years to precisely quantify global-scale biodiversity patterns. After correcting for the spatial inequality of sample completeness, we generated a global species diversity map for woody angiosperms (82,974 species, 13,959,780 occurrence records). The standardized diversity estimated more pronounced latitudinal and longitudinal diversity gradients than the raw data and improved the spatial prediction of diversity based on environmental factors. We identified areas with potentially high species richness and rarity that are poorly explored, unprotected, and threatened by increasing human pressure: They are distributed mostly at low latitudes across central South America, Central Africa, subtropical China, and Indomalayan islands. These priority areas for botanical exploration can help to efficiently fill spatial knowledge gaps for better describing the status of biodiversity and improve the effectiveness of the protected area network for global woody plant conservation.

Global raster dataset on historical coastline positions and shelf sea extents since the Last Glacial Maximum.

Motivation: Historical changes in sea level caused shifting coastlines that affected the distribution and evolution of marine and terrestrial biota. At the onset of the Last Glacial Maximum (LGM) 26 ka, sea levels were >130 m lower than at present, resulting in seaward-shifted coastlines and shallow shelf seas, with emerging land bridges leading to the isolation of marine biota and the connection of land-bridge islands to the continents. At the end of the last ice age, sea levels started to rise at unprecedented rates, leading to coastal retreat, drowning of land bridges and contraction of island areas. Although a growing number of studies take historical coastline dynamics into consideration, they are mostly based on past global sea-level stands and present-day water depths and neglect the influence of global geophysical changes on historical coastline positions. Here, we present a novel geophysically corrected global historical coastline position raster for the period from 26 ka to the present. This coastline raster allows, for the first time, calculation of global and regional coastline retreat rates and land loss rates. Additionally, we produced, per time step, 53 shelf sea rasters to present shelf sea positions and to calculate the shelf sea expansion rates. These metrics are essential to assess the role of isolation and connectivity in shaping marine and insular biodiversity patterns and evolutionary signatures within species and species assemblages. Main types of variables contained: The coastline age raster contains cells with ages in thousands of years before present (bp), representing the time since the coastline was positioned in the raster cells, for the period between 26 ka and the present. A total of 53 shelf sea rasters (sea levels <140 m) are presented, showing the extent of land (1), shelf sea (0) and deep sea (NULL) per time step of 0.5 kyr from 26 ka to the present. Spatial location and grain: The coastline age raster and shelf sea rasters have a global representation. The spatial resolution is scaled to 120 arcsec (0.333° × 0.333°), implying cells of c. 3,704 m around the equator, 3,207 m around the tropics (±30°) and 1,853 m in the temperate zone (±60°). Time period and temporal resolution: The coastline age raster shows the age of coastline positions since the onset of the LGM 26 ka, with time steps of 0.5 kyr. The 53 shelf sea rasters show, for each time step of 0.5 kyr, the position of the shelf seas (seas shallower than 140 m) and the extent of land. Level of measurement: Both the coastline age raster and the 53 shelf sea rasters are provided as TIFF files with spatial reference system WGS84 (SRID 4326). The values of the coastline age raster per grid cell correspond to the most recent coastline position (in steps of 0.5 kyr). Values range from 0 (0 ka, i.e., present day) to 260 (26 ka) in bins of 5 (0.5 kyr). A value of "no data" is ascribed to pixels that have remained below sea level since 26 ka. Software format: All data processing was done using the R programming language.

Multiple drivers of the COVID-19 spread: The roles of climate, international mobility, and region-specific conditions.

Following its initial appearance in December 2019, coronavirus disease 2019 (COVID-19) quickly spread around the globe. Here, we evaluated the role of climate (temperature and precipitation), region-specific COVID-19 susceptibility (BCG vaccination factors, malaria incidence, and percentage of the population aged over 65 years), and human mobility (relative amounts of international visitors) in shaping the geographical patterns of COVID-19 case numbers across 1,020 countries/regions, and examined the sequential shift that occurred from December 2019 to June 30, 2020 in multiple drivers of the cumulative number of COVID-19 cases. Our regression model adequately explains the cumulative COVID-19 case numbers (per 1 million population). As the COVID-19 spread progressed, the explanatory power (R2) of the model increased, reaching > 70% in April 2020. Climate, host mobility, and host susceptibility to COVID-19 largely explained the variance among COVID-19 case numbers across locations; the relative importance of host mobility and that of host susceptibility to COVID-19 were both greater than that of climate. Notably, the relative importance of these factors changed over time; the number of days from outbreak onset drove COVID-19 spread in the early stage, then human mobility accelerated the pandemic, and lastly climate (temperature) propelled the phase following disease expansion. Our findings demonstrate that the COVID-19 pandemic is deterministically driven by climate suitability, cross-border human mobility, and region-specific COVID-19 susceptibility. The identification of these multiple drivers of the COVID-19 outbreak trajectory, based on mapping the spread of COVID-19, will contribute to a better understanding of the COVID-19 disease transmission risk and inform long-term preventative measures against this disease.

Past and future decline of tropical pelagic biodiversity.

A major research question concerning global pelagic biodiversity remains unanswered: when did the apparent tropical biodiversity depression (i.e., bimodality of latitudinal diversity gradient [LDG]) begin? The bimodal LDG may be a consequence of recent ocean warming or of deep-time evolutionary speciation and extinction processes. Using rich fossil datasets of planktonic foraminifers, we show here that a unimodal (or only weakly bimodal) diversity gradient, with a plateau in the tropics, occurred during the last ice age and has since then developed into a bimodal gradient through species distribution shifts driven by postglacial ocean warming. The bimodal LDG likely emerged before the Anthropocene and industrialization, and perhaps ∼15,000 y ago, indicating a strong environmental control of tropical diversity even before the start of anthropogenic warming. However, our model projections suggest that future anthropogenic warming further diminishes tropical pelagic diversity to a level not seen in millions of years.

Integrating multiple sources of ecological data to unveil macroscale species abundance.

The pattern of species abundance, represented by the number of individuals per species within an ecological community, is one of the fundamental characteristics of biodiversity. However, despite their obvious significance in ecology and biogeography, there is still no clear understanding of these patterns at large spatial scales. Here, we develop a hierarchical modelling approach to estimate macroscale patterns of species abundance. Using this approach, estimates of absolute abundance of 1248 woody plant species at a 10-km-grid-square resolution over East Asian islands across subtropical to temperate biomes are obtained. We provide two examples of the basic and applied use of the estimated species abundance for (1) inference of macroevolutionary processes underpinning regional biodiversity patterns and (2) quantitative community-wide assessment of a national red list. These results highlight the potential of the elucidation of macroscale species abundance that has thus far been an inaccessible but critical property of biodiversity.

A geometric approach to scaling individual distributions to macroecological patterns.

Understanding macroecological patterns across scales is a central goal of ecology and a key need for conservation biology. Much research has focused on quantifying and understanding macroecological patterns such as the species-area relationship (SAR), the endemic-area relationship (EAR) and relative species abundance curve (RSA). Understanding how these aggregate patterns emerge from underlying spatial pattern at individual level, and how they relate to each other, has both basic and applied relevance. To address this challenge, we develop a novel spatially explicit geometric framework to understand multiple macroecological patterns, including the SAR, EAR, RSA, and their relationships. First, we provide a general theory that can be used to derive the asymptotic slopes of the SAR and EAR, and demonstrates the dependency of RSAs on the shape of the sampling region. Second, assuming specific shapes of the sampling region, species geographic ranges, and individual distribution patterns therein based on theory of stochastic point processes, we demonstrate various well-documented macroecological patterns can be recovered, including the tri-phasic SAR and various RSAs (e.g., Fisher's logseries and the Poisson lognormal distribution). We also demonstrate that a single equation unifies RSAs across scales, and provide a new prediction of the EAR. Finally, to demonstrate the applicability of the proposed model to ecological questions, we provide how beta diversity changes with spatial extent and its grain over multiple scales. Emergent macroecological patterns are often attributed to ecological and evolutionary mechanisms, but our geometric approach still can recover many previously observed patterns based on simple assumptions about species geographic ranges and the spatial distribution of individuals, emphasizing the importance of geometric considerations in macroecological studies.

Spatially explicit approach to estimation of total population abundance in field surveys.

Population abundance is fundamental in ecology and conservation biology, and provides essential information for predicting population dynamics and implementing conservation actions. While a range of approaches have been proposed to estimate population abundance based on existing data, data deficiency is ubiquitous. When information is deficient, a population estimation will rely on labor intensive field surveys. Typically, time is one of the critical constraints in conservation, and management decisions must often be made quickly under a data deficient situation. Hence, it is important to acquire a theoretical justification for survey methods to meet a required estimation precision. There is no such theory available in a spatially explicit context, while spatial considerations are critical to any field survey. Here, we develop a spatially explicit theory for population estimation that allows us to examine the estimation precision under different survey designs and individual distribution patterns (e.g. random/clustered sampling and individual distribution). We demonstrate that clustered sampling decreases the estimation precision when individuals form clusters, while sampling designs do not affect the estimation accuracy when individuals are distributed randomly. Regardless of individual distribution, the estimation precision becomes higher with increasing total population abundance and the sampled fraction. These insights provide theoretical bases for efficient field survey designs in information deficiency situations.