Mapping human pressures on biodiversity across the planet uncovers anthropogenic threat complexes

Diana E. Bowler1,2,3,4 | Anne D. Bjorkman5,6 | Maria Dornelas7 | Isla H. Myers-Smith8 | Laetitia M. Navarro2,9 | Aidin Niamir1 | Sarah R. Supp10 | Conor Waldock11 | Marten Winter2,12 | Mark Vellend13 | Shane A. Blowes2,14 | Katrin Böhning-Gaese1 | Helge Bruelheide2,9 | Robin Elahi15 | Laura H. Antão7,16,17 | Jes Hines2,18 | Forest Isbell19 | Holly P. Jones20 | Anne E. Magurran7 | Juliano Sarmento Cabral2,21 | Amanda E. Bates11,22


| INTRODUC TI ON
Human activities are reshaping biological communities and impacting ecosystem functioning across the Earth (Barnosky et al., 2011;Dornelas et al., 2014;Isbell et al., 2017). Meeting the global challenge of the conservation and sustainable use of nature requires not only quantifying biodiversity change, but also identifying the underlying causes of change (Isbell et al., 2017;Tittensor et al., 2014). Climate change, habitat change, exploitation, pollution and invasive alien species have been recognized as the most important and widespread direct anthropogenic causes of biodiversity change (IPBES, 2019;IPCC, 2013;Pereira, Navarro, & Martins, 2012). These five main drivers have been linked with changes in multiple dimensions of biodiversity, including genetic diversity, species' population sizes, community richness and ecosystem functioning (Pereira et al., 2012). The impacts of these drivers on a biological community in a given region critically depend on the amount of exposure to each driver, which is described by its local magnitude or change (such as the strength of climate change or intensity of land-use). An important, but so far underexplored, step towards understanding the global patterns of biodiversity change is characterizing the exposure patterns of biological communities to different types of environmental change.
Global maps of pressures such as the terrestrial human footprint (Sanderson et al., 2002;Venter et al., 2016), marine pressures  and river threats  characterize the geographic hotspots of anthropogenic threats to biodiversity. These maps have estimated that at least 75% of terrestrial land has been exposed to some sort of land-use change  with only 23.2% classified as wilderness . Similarly, the marine realm has been exposed to multiple pressures  with 13.2% of the oceans regarded as wilderness (Jones et al., 2018;IPBES, 2019). The intensity of the terrestrial human footprint has been linked with spatial variation in ecological processes; for instance reduced animal movement was found in areas with a higher human footprint across different species (Tucker et al., 2018). However, these global maps show the summed pressure of different drivers related to human activities and ignore any relationships among them. Areas identified as having high human pressure could be underlaid by different combinations of drivers with varying intensity, each of which may have contrasting impacts on biodiversity. Hence, while useful, these cumulative threat maps do not help understand the relative contribution of each driver to global biodiversity change, nor the potential for interactive impacts among drivers.
The relative importance of different drivers for biodiversity change and ecosystem services is a key component of both policy-oriented assessments such as IPBES framework (Diaz et al., 2015) and conservation targets such as CBD Aichi Biodiversity Targets (Tittensor et al., 2014). While many studies have examined the relationships between different anthropogenic drivers and changes in biodiversity and ecosystems, they usually focus on only one or two drivers, such as climate change and land use change (Sirami et al., 2017). Multiple drivers may be associated with each other, either coincidentally (due to shared causes) or causally (when one driver affects the intensity of another) (Geary, Nimmo, Doherty, Ritchie, & Tulloch, 2019). Moreover, when drivers co-occur, their impacts on communities may be additive, or interact synergistically (total impact stronger when together) or antagonistically (total impact weaker when together) (Cote, Darling, & Brown, 2016;Sirami et al., 2017). Nonetheless, few studies have examined the effects of multiple drivers on biodiversity (Mazor et al., 2018;Sirami et al., 2017). Hence, unpacking the spatial patterns of exposure of different drivers, and assessing the extent to which they jointly act on communities, is a key area of research.
Here, we analysed the spatial relationships among variables related to the main anthropogenic drivers of biodiversity change, and show how they overlap in different biogeographical regions across the entire surface of the world. The motivation for our study was primarily to advance the study of the attribution of biodiversity terrestrial realm (e.g. in the tundra and boreal forest with high climate change but low human use and pollution), whereas the opposite was true in the marine realm (e.g. in the Indo-Pacific with high climate change and high fishing).

We show that different regions of the world can be defined by Anthropogenic
Threat Complexes (ATCs), distinguished by different sets of drivers with varying intensities. We identify 11 ATCs that can be used to test hypotheses about patterns of biodiversity and ecosystem change, especially about the joint effects of multiple drivers. 5. Our global analysis highlights the broad conservation priorities needed to mitigate the impacts of anthropogenic change, with different priorities emerging on land and in the ocean, and in different parts of the world. change, especially regarding potential drivers of long-term changes in species' populations and communities (Bowler et al., 2017;Dornelas et al., 2014Dornelas et al., , 2019. Hence, our questions were focused towards the relevance of the spatial patterns of drivers for biodiversity change rather than on explaining the spatial patterns in drivers themselves. We determined the extent of driver overlap to assess the potential for different drivers to either act alone on communities, rendering their specific fingerprints easier to detect, or for drivers to act in combination, with the potential for interactive effects. For many drivers, it can be hypothesized that exposure patterns may be inter-linked due to related local or regional human activities, driven by local human population density (Ellis, Goldewijk, Siebert, Lightman, & Ramankutty, 2010;Geary et al., 2019). In contrast, climate change is expected to be distributed differently to other variables because it is an outcome of processes at regional and global scales (IPCC, 2013).
Towards these aims, we quantified the strengths of the spatial relationships among the intensities of variables related to different anthropogenic drivers: climate change, habitat conversion and exploitation (grouped together as 'human use'), pollution and potential for alien species immigration (Table 1). On the basis of these relationships, we defined 'Anthropogenic Threat Complexes' that typify the combinations of drivers impacting different regions of the world. Studies mapping drivers of biodiversity change have so far considered the terrestrial and marine realms separately. By employing a standardized analysis for both realms, our study also highlights similarities and differences in anthropogenic environmental changes across the world, including across realms.

| Approach to data selection
By drivers, we refer to anthropogenic drivers only and not natural drivers, such as inter-annual climatic variation, nor biological processes, such as dispersal and demographical rates, that more proximately lead to changes in species' abundances. Our analysis required datasets that were (a) global; (b) spatially gridded at high-resolution and (c) publicly available. We selected datasets on variables related to the main five drivers that were included in previous studies of global environmental change in the context of biodiversity change (Table 1; Sanderson et al., 2002). We further searched for data on other relevant variables following the IUCN threats categories (Table S1; . As we focused on a land versus ocean comparison, we did not specifically consider freshwater threats . Variables representing drivers of biodiversity change were in most cases based on data between 1990 and 2010 (except climate change and forest loss; see below). In total, we had datasets for 16 driver-related variables ( Table 1). The terrestrial datasets came from various sources (Table S2). Most of the marine datasets came from the landmark study of . All variables were continuous quantitative variables. Although the specific variables differ among realms, we aligned each variable to the dominant drivers that are common across both realms (Table 1). Because habitat conversion and exploitation were difficult to classify separately across terrestrial and marine ecosystems, we combined both into a single 'human use' category.

| Climate change
Climate change has multiple components (IPCC, 2013), hence we characterized climate change by several variables using global spatiotemporal gridded temperature data for the terrestrial  and marine realm . To calculate temperature trends, we used data between 1950, proposed as the start of the Anthropocene (Waters et al., 2016), and 2010. Temperature trends were estimated by fitting simple linear regression models to annual temperature means of each grid cell and extracting the coefficient for the effect of year.
Temperature divergence was inferred from the t-statistic of this linear regression, that is the trend divided by its standard error, representing the significance of the trend. Velocity of climate change was calculated as the ratio between the temporal temperature trend and the local spatial gradient in temperature and reflects the speed at which an organism would have to move to track its current climatic conditions (Loarie et al., 2009). Trends of extreme temperatures were calculated by whichever was largest of the temporal trends in mean temperature of the warmest or coolest month. We further included aridity trend for the terrestrial realm and ocean acidification for the marine realm. The aridity trend was estimated by taking monthly and annual datasets on potential evapotranspiration and precipitation, and calculating the monthly:annual ratio (Zorner, Trabucco, Bossio, & Verchot, 2008), followed by the temporal trend of the annual monthly average of this ratio. We focusing on aridity trend since reduced water availability can exacerbate the impacts of rising temperatures (Mitchell et al., 2018). Ocean acidification, from the Halpern layers, was

| Human use
In the terrestrial realm, we collated human use variables related to the proportional coverage of different types of land conversion or use: cropland, pasture land, cattle density, urban area and forest loss. Data on crop land, pasture land and urban/built land-cover were taken from different databases, primarily based on satellite imagery-crop land , pasture land  and urban land . We additionally included information on the amount of forest loss since deforestation itself is a recognized threat (Barlow et al., 2016). Forest loss, based on FAO (Food and Agriculture Association of the United Nations) wood harvest statistics, was calculated as the loss of primary forest for the same time frame as our climate change statistics, that is between 1950 and 2010 (Hurtt et al., 2011). We also included data on cattle density, which was based on subnational livestock data that were statistically downscaled using multiple predictors . In the marine realm, human use variables were based on different commercial fishing activities separated by gear types (e.g. dredging or cast nets), which determine their selectivity and impact on the surrounding seascape . These fishing types were pelagic low-bycatch, pelagic high-bycatch, demersal habitat-modifying, demersal non-habitat-modifying low-bycatch and demersal non-habitat-modifying high bycatch. These data were based on FAO and other commercial catch data sources, available at a 0.5° resolution.  subsequently standardized the catch data by average ocean productivity (based on a Vertically Generalized Production Model) and disaggregated the data to 1 km 2 resolution to align with their other datasets .

| Pollution
Nitrogen from both fossil fuel combustion and agriculture is one is the biggest pollutants impacting biodiversity (De Schrijver et al., 2011). We included data on nitrogen pollution for the terrestrial realm in the form of atmospheric nitrogen  and fertilizer use , and for the marine realm as fertilizer use, based on watershed boundaries, amount used on nearby agricultural land and a river plume model . We also included data on pesticide use in both realms , another important component of agricultural intensification negatively affecting biodiversity (Geiger et al., 2010). National estimates of fertilizer use and pesticide were downscaled by the data providers according to land-cover maps; thus, these datasets were not fully independent of the cropland data and were expected to be correlated, but they still represent the best current knowledge of the spatial distribution of these variables. We also included a layer estimating the extent of ocean pollution, based on the distribution of shipping lanes . Finally, we included night-time light pollution detected by satellite imagery , also included in previous terrestrial and marine threat maps Venter et al., 2016).

| Alien species potential
The spread of alien species is among the greatest threats to biodiversity and ecosystem services (Blackburn, Bellard, & Ricciardi, 2019).
Alien species are defined as species that are introduced into areas beyond their historical range, usually through human transport, accidentally or incidentally (Hulme, 2009;Seebens et al., 2015). We used information on transport infrastructure related to human movement and trade that depict possible species transportation pathways and vectors (Davidson, Scianni, Minton, & Ruiz, 2018;Hulme, 2009).
Specifically, we used spatial datasets of transport connectivity (including data on road and rail networks and navigable rivers) in the terrestrial realm and cargo volume at ports in the marine realm (Table   S1). While these do not represent the only invasion pathways, they are commonly accepted proxies for human-mediated propagule pressure, which is known to be among the most important determinant of alien species establishment (Hulme, 2009;Seebens et al., 2015).
We note two limitations of our alien species data. First, our proxy was for alien species in general and not specifically for invasive and harmful alien species. Only a proportion of alien species become invasive, but this proportion differs among taxa (Jeschke & Pyšek, 2018). Although the most harmful alien species have been identified for some regions (Katsanevakis et al., 2014), this is an ongoing process at the global level (Pagad, Genovesi, Carnevali, Schigel, & McGeoch, 2018). Hence, we refer to this variable as alien species potential rather than specifically as invasive alien species. Second, we only used a proxy for alien species propagule pressure and not direct measures of their distribution. Our analysis required high-resolution global gridded maps, but high-resolution alien species distribution data are not available for neither realm. However, information on alien species distribution was available at a regional, subnational and national levels for some taxonomic groups, including birds  and plants (van Kleunen et al., 2019). We used these datasets, representing taxa with low and high mobility, to assess the validity of our proxy by testing the correlation between alien species richness and mean connectivity at the spatial scale of the distribution data. We found a significant rank correlation for both datasets (birds, ρ = 0.42; plants, ρ = 0.46, see Figure S1 for more details), suggesting our proxy represents a reasonable estimate of alien species potential.

| Human population
We also included 'human population density' as a separate driver (CIESIN, 2017) accounting for the effects of human activities not falling into the other categories . Although we recognize that the relationship between local human population density and the local human activities negatively impacting on biodiversity can strongly vary among countries, for example with economy size. By including human population density, we could also test the relationship between human population density and the other drivers.

| Justification for layer exclusion
We did not use data for some variables that were previously included in the terrestrial human footprint or the Halpern layers.
The human footprint includes data on roads, railways and navigable waterways . Although we did not separately include these data, these data were taken into account in the terrestrial connectivity variable (for alien species potential). In the marine realm, we excluded a shipping lane variable since the ocean pollution variable was already based on the distribution of shipping lanes . Additional available marine layers that we did not use were: UV radiation, oils rigs (based on night lights, already included), inorganic pollution (highly correlated with other land-based coastal pollutants that were already included) and artisanal fishing (data poor and mostly modelled; Halpern et al., 2008).

| Geographical region data
To understand how the driver distributions compared with classic biome distributions that are already used in biodiversity research, we also extracted geographical region data on biomes. Data on the spatial distribution of terrestrial biomes (including deserts, temperate broadleaf forest and boreal forest among others) were taken from WWF (Olson et al., 2001). Marine regions were defined by combining coastal/shelf region polygon data-MEOW (Spalding et al., 2007) and ocean polygon data (natur alear thdata.com). We did not use marine ecosystem data as used by others  because the ecosystems spatially overlapped in our 2-D global raster grid, when, in reality, these different ecosystems occur at different depths in the water column.

| Data processing
We harmonized each dataset to a standard global grid. The resolutions of the original datasets were approximately at a 100 km 2 grid (or 1°) or finer resolution; hence, we aggregated all datasets to a standard grid of 100 km 2 grid cells by taking the mean value of the grid cells. Atmospheric nitrogen deposition was only available at a coarser resolution (see Table S1); however, we disaggregated this also to 100 km. Datasets were bound between latitudes of −58 and 78 (due to poorer data at extreme latitudes) and reprojected onto a common equal-area map projection (Eckert IV; ESPG = 54,012). Missing values in some of the human activity datasets (primarily in pesticide use) occurred in remote regions (e.g. very high latitudes and deserts) with likely absent or low variable values (since they were in areas with no crop cover) and were imputed as zero. Greenland was excluded due to missing data in several of the datasets.
Because each dataset comprised data in different units (e.g. temperature data in °C and fertilizer data in kg/ha), it was not possible to directly compare their absolute values. Instead, we ranked and scaled (between 0 and 1) values of each dataset for ease of interpretation ( Figure S2 shows the distributions of the original values of each variable and Figure S3 shows global maps of the ranked and scaled data). This processing also reduced the large

| Data analysis
To examine the relationships among the intensities of the 16 different variables, we calculated Spearman's rank correlation coefficients (ρ) for each pairwise combination of variables across all grid cells in each realm. We chose this statistic because it is equivalent to the commonly used Pearson's correlation on ranked data, which was consistent with our aforementioned data processing. We used Dutilleul's modified t test to account for spatial autocorrelation before testing the significance of the correlations (Dutilleul, Pelletier, & Alpargu, 2008). We also used Moran's I and correlograms to assess the extent of spatial autocorrelation within each variable and understand the spatial scale at which each variable tends to act. For the marine realm, correlations were also examined separately for grid cells whose centroid overlapped with oceanic or coastal regions. To assess the importance of the drivers in different parts of the world, we compared the average intensity of each of the five main drivers (climate change, human use, human population, pollution and alien species potential) for different geographical regions. For this, we first calculated the mean driver value for each grid cell, that is, averaging across associated driver-related variables (Table 1). We then plotted the distribution of these mean values across all grid cells within each terrestrial biome and marine region.
Since we found many variables were correlated, we used a clustering technique to group together the variables with similar spatial patterns and examine where these groups can be found across the world. Although, driver combinations vary in a continuous manner, we chose a clustering method that produces discrete grouping to provide the simplest description of the main patterns in the data. We used k-medoid clustering, with the partitioning around the medoids algorithm and Euclidean distances (Maechler, Rousseeuw, Struyf, Hubert, & Hornik, 2018), for clustering grid cells according to their extent of exposure of all the variables. To focus the clustering on the main axes of variation, we first used PCA to reduce the number of variables in driver groups with multiple variables, that is climate change, human use and pollution.
We performed the clustering with two or three PCA components (whichever captured at least 75% of the variation) for each of these drivers along with the variables for human population density and potential alien species. We selected the number of clusters by comparing the changes in cluster silhouette width with increasing cluster number; however, we limited the cluster number to <10.
Two terrestrial clusters were subsequently grouped since they only slightly differed in one variable and had similar mean driver intensities. To smooth the maps, we used a moving window to assign each cell the mode of its 3 × 3 cell neighbourhood. Finally, to repack the datasets into cumulative driver maps, we summed the number of variables for which each grid cell was in the upper 10% of values (based on values greater than zero). Analyses were run in R v. 3.4.1 (R Core Team, 2018), using the packages raster (Hijmans, 2017), SpatialPack (Vallejos, Osorio, & Bevilacqua, 2018) and cluster (Maechler et al., 2018).

| Sensitivity analyses
To examine the effect of the grain size, we repeated the data processing steps except harmonized the datasets to a global grid of 500 km resolution and repeated the correlation analysis (similar results were obtained-see Figure S4).

| RE SULTS
We found that anthropogenic drivers of biodiversity change were more spatially coupled in the terrestrial than in the marine realm ( Figure 1; Figure S7). On land, 40% of the possible pairwise relationships between variables (excluding climate change-related variables) showed positive correlation strengths of at least 0.7. Thus, terrestrial areas with high intensities of one variable also tended to have high intensities of other variables. Moreover, correlations were found between different types of terrestrial drivers; for instance high crop land-cover and cattle density were associated with high pollution, high alien potential and high human population density. Conversely, in the marine realm, we found fewer correlations-only 15% of the possible pairwise relationships (excluding climate change-related variables) showed a strong positive correlation (>0.7)-and these relationships were mostly within, rather than between, different driver types; for instance among different types of human use (e.g. different types of demersal fishing; Figure 1). Across all variables, oceanic regions showed fewer correlations compared to coastal regions ( Figure S8). Spatial autocorrelation was present in all variables and tended to reach greater distances in the marine human-uses and climate-change variables ( Figures S9, S10), and shorter distances in the coastal-based marine variables. However, the correlations among drivers in both realms were statistically significant (all p < .05 for links shown in Figure 1) after accounting for autocorrelation. In neither realm were there strong negative correlations among variables ( Figure S11 shows the full correlation matrix). Terrestrial biomes exposed to some of the strongest climate change, such as the tundra, boreal forest and desert, have experienced the lowest intensities of human use (Figure 3; Figure S12).
In contrast, other terrestrial biomes, such as tropical dry broadleaf forest and tropical coniferous forest, have had high intensities of human use, pollution and invasions but lower intensities of climate change. Temperate broadleaf and mixed forest is distinct in the terrestrial realm by experiencing higher than average intensities of all drivers (Figure 3). In the marine realm, central and western Indo-Pacific were most exposed to multiple drivers, including rapid climate change and multiple human uses, especially fishing ( Figure 3).
The temperate Northern Atlantic, which includes the North Sea, has been also strongly exposed to multiple drivers, whereas low mean intensities of multiple drivers emerge in the Pacific Ocean.
The cluster analysis defined multiple regions for each realm (five terrestrial regions and six marine regions) that had different exposure patterns to the different driver variables (Figures 4, 5). These exposure patterns can be regarded as 'anthropogenic threat complexes' (ATC) that characterize the typical combinations of environmental change emerging at large scales. ATCs I and VI (dark blue-grey regions in Figure 4) represent terrestrial and marine areas with high exposure to multiple drivers, especially non-climatic drivers such as human use, pollution and alien potential (Figure 4). At the other extreme, ATC XI represents areas (light orange) that were exposed to low intensities of most variables and was only found in the marine realm ( Figure 4). On land, climate change variables were strongest within ATCs IV and V (orange regions), despite variation in aridity trend ( Figure 5). In the marine realm, climate change variables were generally strongest within ATCs IX and X (orange regions), except ocean acidification ( Figure 5).
The largest terrestrial ATC was V (34% of terrestrial grid cells), covering regions exposed to medium-high climate change and lower intensities of other drivers. The largest marine ATC was IX (23% of marine grid cells), including regions exposed to both medium human use, especially fishing, and medium-high climate change (Figures 4, 5).
Here we produced a truly global map ( Figure 6) highlighting the areas exposed to high intensities of multiple drivers and connect the ATCs to previous cumulative human impact maps produced separately for the terrestrial (Sanderson et al., 2002;Venter et al., 2016) F I G U R E 1 Strong and positive relationships among anthropogenic drivers of biodiversity change. Each grey link between two variables represents a significant positive correlation with strength >0.7 of the variable intensities across 100 km 2 grids covering each realm. We find a higher number of correlations between drivers in the terrestrial versus the marine realm F I G U R E 2 Relationships between climate change and non-climatic drivers (human use, pollution and alien species). The length of each bar represents the correlation coefficient between each variable and temperature change (upper bar, air or sea surface temperature-SST) or VOCC (lower bar, velocity of climate change). We find negative correlations in the terrestrial realm and positive correlations in the marine realm. *denotes statistical significance after accounting for spatial autocorrelation F I G U R E 3 Regions of the terrestrial and marine realms are exposed to distinct combinations of drivers. The violin plots show the distribution of values for each driver in each region. Violins with a median greater than the global median of each driver (the dashed zero line; median of biome values) are coloured in a darker shade. Regions are presented in declining order of total driver exposure: regions at the top are exposed to the greater intensities of drivers and those at the bottom are exposed to the least. Names of the terrestrial regions were shortened for presentation purposes. Figure S12 shows the full distributions for each individual driver variable in each region as well as the full names of the terrestrial regions and marine realms .
Terrestrial regions with the highest cumulative intensities across multiple drivers tended to be within ATC I, especially parts of Europe, India as well as Northern and Eastern China ( Figure 6). In the marine realm, the highest cumulative intensities were found within ATCs VI and IX, especially in the North Sea of the Northern Atlantic, East China Sea and central Indo-Pacific. Terrestrial regions with the lowest cumulative intensities of multiple drivers were commonly within ATCs III and V, covering parts of South America, Africa, Russia and Australia. In the marine realm, the lowest cumulative intensities were within ATCs X and XI, in the southern Atlantic Ocean and Pacific Ocean.

| D ISCUSS I ON
Our findings suggest that many drivers of biodiversity change are unlikely to act alone, but rather to jointly impact biological communities. In the terrestrial realm, drivers that tend to be simultaneously present in a region were associated with agricultural activities (including land conversion to cropland, use of land for cattle rearing, use of pesticides/fertilizer), human population density (urban cover and light pollution) and potential for alien species.
These drivers have especially high spatial overlap in the temperate broadleaf and mixed forests of eastern North America, Europe and F I G U R E 4 Geography of the Anthropocene. Different geographical regions of the world are exposed to different Anthropogenic Threat Complexes (labelled I-XI). In the world map, each colour in each realm denotes a region exposed to the same ATC. The bar charts below are a legend for the colours of the map that describe the magnitudes of climate change and non-climatic drivers in each ATC. Similar ATCs across realms are shaded with the same colour. The bars represent the lower and upper quartiles of the intensities of variables within each driver group (climate change or non-climatic drivers that included human use, pollution and alien species potential) China. In the marine realm, drivers that tend to spatially overlap include different types of pollutants and the potential for alien species, with the coastal regions of the Indo-Pacific and Northern Atlantic (especially the North Sea) showing the highest overlap.
Taken together with recent studies that highlight the lack of research on interactive effects of drivers (Mazor et al., 2018;Sirami et al., 2017), our findings emphasize the widespread relevance of further research on understanding how multiple drivers combine to impact biological communities.
Correlations among drivers were more frequent in the terrestrial realm than in the marine realm. Spatial relationships among drivers in the terrestrial realm were likely based on different activities supporting proximal human populations, leading to simultaneous pressures from land-use change and pollution (Ellis et al., 2010  . Different histories of human activities in the terrestrial and marine realm likely also plays a role. Human exploitation in the terrestrial realm is more advanced, spanning millennia, whereas technological advances have only recently allowed greater exploitation of the world's oceans (Knapp et al., 2017). ships between a driver and a metric of biodiversity change may in fact be caused by a correlated driver. Although spatial heterogeneity at smaller spatial scales (e.g. neighbouring sites that differ in only one driver, such as land-use change) can estimate the local effect of drivers such as habitat conversion (Newbold et al., 2015), correlated large-scale drivers affect regional species pools and hence still influence local community dynamics (Harrison & Cornell, 2008). Analysis of the relationship between ATCs and species' population trends or diversity hotspots for conservation may help identify the most globally harmful combinations of drivers.
While we limited our analysis to coarse patterns, the ATCs could be potentially resolved to finer spatial scales to allow for more regional analysis. This could be used, for instance, to inform the design of biodiversity observatories that treat driver exposure as a natural experiment. Observatories could be selected along a given driver gradient, keeping all but this driver constant. Alternatively, regions where specific drivers emerge could be selected for targeted research to isolate the particular impact of one driver, or a combination of drivers, over others.
Large-scale information on drivers is most immediately relevant for global conservation policy (IPBES, 2019; Tittensor et al., 2014). For biodiversity conservation at specific locations, local data on the magnitudes of different types of anthropogenic pressures are essential to guide management plans. However, even in this latter context, there are a number of advantages of having F I G U R E 6 Regions of the world exposed to high intensities of multiple drivers. The main map shows the number of the 16 driver variables for which each grid cell was in the highest 10% of values within each realm. Regions in the darkest orange are exposed to high intensities of multiple variables, whereas those in off-white are exposed to lower intensities (i.e. within the 90% quantile of values) of all variables. The smaller plots below show the same for each of the separate drivers. Larger versions of the driver plots are presented in Figure  S13 knowledge on the large-scale patterns of drivers. First, large-scale patterns allow local management to be modified according to the broader anthropogenic land-or seascape context (Harrison & Cornell, 2008). While landscape ecology and habitat connectivity already inform conservation planning, knowledge of all drivers present in the region could influence local decision making by enabling prediction of possible future threats and biodiversity change at local scales, for example large-scale land-use change affecting local community changes through dispersal (Hansen & DeFries, 2007). Second, by characterizing regions of the world in terms of the nature of anthropogenic environmental change, the ATCs suggest how information and data might be pooled and synthesized across regions, even across the marine and terrestrial realms. Regions exposed to the same ATC, regardless of location, would benefit from exchanging knowledge about prioritization strategies and management of the multiple drivers, as well as implementing cross-border and inter-regional strategies to minimize their impact (Bonebrake et al., 2019).
Quantifying exposure to environmental change is the first step towards determining which species, in which places, are most impacted by human activities. However, the realized outcome of different drivers on biodiversity will ultimately depend on a combination of both the magnitude of exposure to drivers and species' sensitivities to environmental change (Foden et al., 2013). While exposure characterizes the amount of environmental change (e.g. temperature change in °C), sensitivity refers to the tolerance of the biological community to that driver (e.g. abundance (or individual fitness) change per °C temperature change). We intentionally focused on exposure patterns so that our results were not species-specific and therefore potentially relevant for any taxa or ecosystem. Unlike exposure, sensitivities vary among taxa according to characteristics such as their life history, traits and niche breadth (Sunday et al., 2015) and therefore need to be examined separately for different taxa. Hence, despite similar exposure patterns, we can expect a diversity of biodiversity responses within each ATC due to variation in species' sensitivities.
Our findings are clearly dependent on the included layers-we selected layers based on their past inclusion in global threat maps and their linkage to IUCN threat categories. However, many of the recognized IUCN threats to biodiversity (Joppa et al., 2016; were not available as high-resolution spatial datasets at a global scale, which was necessary for our analysis. These threats include energy production and mining, hunting and other forms of human disturbance . Data availability was especially limited for the marine realm and the threat from invasive alien species globally. Rather than use proxy variables of species transport and thus of propagule pressure, spatially explicit maps of the number of invasive alien species would have improved our analysis. The Global Register of Introduced and Invasive Species (www.griis.org) is working towards better knowledge of the distribution of invasive species (Pagad et al., 2018). Our study also only focused on spatial patterns and did not explicitly consider the consequences of spatiotemporal changes in human activities over time. Ongoing projects such as the Copernicus project (http://www.coper nicus.eu/) will greatly increase the availability of high-resolution spatiotemporal datasets on different variables in the coming years, enabling better attribution of biodiversity change to the underlying drivers.
Our macroecological approach to mapping the anthropogenic drivers of biodiversity change contributes to the development of broad conservation policy targeted towards the mitigation of specific driver complexes. Development of conservation strategies that simultaneously attempt to tackle multiple drivers are likely to be more efficient in the long-term (Bates et al., 2019;Bonebrake et al., 2019). The ATC framework emphasizes the fact that multiple drivers usually affect local and regional biological communities and discourages the prevalent simplistic focus on one or two drivers. Much more research should focus on understanding the joint impacts of multiple drivers and how different drivers have interactive effects on biological communities. Our findings help direct some of these future studies by identifying which drivers most commonly overlap, and in which regions of the world. Finally, by taking a cross-realm approach, we hope to encourage information exchange across regions of the world that are exposed to similar suites of drivers, regardless of environmental realm, and the development of joined-up conservation policies across national borders and the terrestrial-marine interface.

ACK N OWLED G EM ENTS
This paper arose from discussions at the sChange Workshop (www. idiv.de/schange), which was supported by sDiv, the Synthesis Centre of iDiv, the German Centre for Integrative Biodiversity Research Halle-Jena-Leipzig. We thank all other participants of this workshop for the stimulating discussion and compilation of data that directly fed into this project. We thank the Associate Editor, Julian

CO N FLI C T O F I NTE R E S T
The authors have no conflicts of interest to declare.

DATA AVA I L A B I L I T Y S TAT E M E N T
The Data Sources section shows the sources of each dataset and