1. Introduced large herbivores have partly filled ecological gaps formed in the late Pleistocene, when many of the Earth's megafauna were driven extinct. However, extant predators are generally considered incapable of exerting top-down influences on introduced megafauna, leading to unusually strong disturbance and herbivory relative to native herbivores.2. We report on the first documented predation of juvenile feral donkeys Equus africanus asinus by cougars Puma concolor in the Mojave and Sonoran Deserts of North America. We then investigated how cougar predation corresponds with differences in feral donkey behaviour and associated effects on desert wetlands.3. Focusing on a feral donkey population in the Death Valley National Park, we used camera traps and vegetation surveys to compare donkey activity patterns and impacts between wetlands with and without cougar predation. 4. Donkeys were primarily diurnal at wetlands with cougar predation, thereby avoiding cougars. However, donkeys were active throughout the day and night at sites without predation. Donkeys were ~87% less active (measured as hours of activity a day) at wetlands with predation (p < 0.0001). Sites with predation had reduced donkey disturbance and herbivory, including ~46% fewer access trails, 43% less trampled bare ground and 192% more canopy cover (PERMANOVA, R 2 = 0.22, p = 0.0003). 5. Our study is the first to reveal a trophic cascade involving cougars, feral equids and vegetation. Cougar predation appears to rewire an ancient food web, with diverse implications for modern ecosystems. Our results suggest that protecting apex predators could have important implications for the ecological effects of introduced megafauna.

Advances in climate science have rendered obsolete the gridded observation data widely used in downstream applications. Novel climate reanalysis products outperform legacy data products in accuracy, temporal resolution, and provision of uncertainty metrics. Consequently, there is an urgent need to develop a workflow through which to integrate these improved data into biological analyses. The ERA5 product family (ERA5 and ERA5-Land) are the latest and most advanced global reanalysis products created by the European Center for Medium-range Weather Forecasting (ECMWF). These data products offer up to 83 essential climate variables (ECVs) at hourly intervals for the time-period of 1981 to today with preliminary back-extensions being available for 1950-1981. Spatial resolutions range from 30x30km (ERA5) to 11x11km (ERA5-Land) and can be statistically downscaled to study-requirements at finer spatial resolutions. Kriging is one such method to interpolate data to finer resolutions and has the advantages that one can leverage additional covariate information and obtain the uncertainty associated with the downscaling. The KrigR R-package enables users to (1) download ERA5(-Land) climate reanalysis data for a user-specified region, and time-period, (2) aggregate these climate products to desired temporal resolutions and metrics, (3) acquire topographical co-variates, and (4) statistically downscale spatial data to a user-specified resolution using co-variate data via kriging. KrigR can execute all these tasks in a single function call, thus enabling the user to obtain any of 83 (ERA5) / 50 (ERA5-Land) climate variables at high spatial and temporal resolution with a single R-command. Additionally, KrigR contains functionality for computation of bioclimatic variables and aggregate metrics from the variables offered by ERA5(-Land). This R-package provides an easy-to-implement workflow for implementation of state-of-the-art climate data while avoiding issues of storage limitations at high temporal and spatial resolutions by providing data according to user-needs rather than in global data sets. Consequently, KrigR provides a toolbox to obtain a wide range of tailored climate data at unprecedented combinations of high temporal and spatial resolutions thus enabling the use of world-leading climate data through the R-interface and beyond.

There is an increasing need for high spatial and temporal resolution climate data for the wide community of researchers interested in climate change and its consequences. Currently, there is a large mismatch between the spatial resolutions of global climate model and reanalysis datasets (at best around 0.25° and 0.1° respectively) and the resolutions needed by many end-users of these datasets, which are typically on the scale of 30 arcseconds (∼900 m). This need for improved spatial resolution in climate datasets has motivated several groups to statistically downscale various combinations of observational or reanalysis datasets. However, the variety of downscaling methods and inputs used makes it difficult to reconcile the resultant differences between these high-resolution datasets. Here we make use of the KrigR R-package to statistically downscale the world-leading ERA5(-Land) reanalysis data using kriging. We show that kriging can accurately recover spatial heterogeneity of climate data given strong relationships with co-variates; that by preserving the uncertainty associated with the statistical downscaling, one can investigate and account for confidence in high-resolution climate data; and that the statistical uncertainty provided by KrigR can explain much of the difference between widely used high resolution climate datasets (CHELSA, TerraClimate, and WorldClim2) depending on variable, timescale, and region. This demonstrates the advantages of using KrigR to generate customized high spatial and/or temporal resolution climate data.

The ecological sciences have joined the big data revolution. However, despite exponential growth in data availability, broader interoperability amongst datasets is still needed to unlock the potential of open access. The interface of demography and functional traits is well-positioned to benefit from said interoperability. Trait-based ecological approaches have been criticised because of their inability to predict fitness components, the core of demography; likewise, demographic approaches are data-hungry, and so using traits as ecological shortcuts to understanding and forecasting population viability could offer great value. Here, we introduce MOSAIC, an open-access trait database that unlocks the demographic potential stored in the COMADRE, COMPADRE, and PADRINO open-access databases. MOSAIC data have been digitised and curated through a combination of existing datasets and additional taxonomic and/or trait records sourced from primary literature. In its first release, MOSAIC (v. 1.0.0) includes 14 trait fields for 300 animal and plant species: biomass, height, growth determination, regeneration, sexual dimorphism, mating system, hermaphrodism, sequential hermaphrodism, dispersal capacity, type of dispersal, mode of dispersal, dispersal classes, volancy, and aquatic habitat dependency. MOSAIC also includes species-level phylogenies for 1,359 species and population-specific climate data where locations are recorded. Using MOSAIC, we highlight a taxonomic mismatch of widely used trait databases with existing structured population models. Despite millions of trait records available in open-access databases, taxonomic overlap between open-access demographic and trait databases is <5%. We identify where traits of interest to ecologists can benefit from database integration and start to quantify traits that are poorly quantified (e.g., growth determination, modularity). The MOSAIC database evidences the importance of improving interoperability in open-access efforts in ecology as well as the need for complementary digitisation to fill targeted taxonomic gaps. In addition, MOSAIC highlights emerging challenges associated with the disparity between locations where different trait records are sourced.