Species Richness-Environment Relationships of European Arthropods at Two Spatial Grains: Habitats and Countries

We study how species richness of arthropods relates to theories concerning net primary productivity, ambient energy, water-energy dynamics and spatial environmental heterogeneity. We use two datasets of arthropod richness with similar spatial extents (Scandinavia to Mediterranean), but contrasting spatial grain (local habitat and country). Samples of ground-dwelling spiders, beetles, bugs and ants were collected from 32 paired habitats at 16 locations across Europe. Species richness of these taxonomic groups was also determined for 25 European countries based on the Fauna Europaea database. We tested effects of net primary productivity (NPP), annual mean temperature (T), annual rainfall (R) and potential evapotranspiration of the coldest month (PETmin) on species richness and turnover. Spatial environmental heterogeneity within countries was considered by including the ranges of NPP, T, R and PETmin. At the local habitat grain, relationships between species richness and environmental variables differed strongly between taxa and trophic groups. However, species turnover across locations was strongly correlated with differences in T. At the country grain, species richness was significantly correlated with environmental variables from all four theories. In particular, species richness within countries increased strongly with spatial heterogeneity in T. The importance of spatial heterogeneity in T for both species turnover across locations and for species richness within countries suggests that the temperature niche is an important determinant of arthropod diversity. We suggest that, unless climatic heterogeneity is constant across sampling units, coarse-grained studies should always account for environmental heterogeneity as a predictor of arthropod species richness, just as studies with variable area of sampling units routinely consider area.


Introduction
Since the early observations of Darwin and von Humboldt, ecologists have attempted to explain why species diversity increases towards the Equator. Relationships of biodiversity with net primary productivity (NPP) are frequently suggested as potential explanations [1][2][3][4], and NPP is usually the best correlate of biodiversity [5]. Productivity-diversity relationships are assumed to be driven by the NPP of an ecosystem as a result of increased provision of vital resources [1]. For example, climates which are highly deficient in water (such as desert) or energy (such as arctic) have both low NPP and low species richness. There is, however, little consensus on the mechanisms underpinning increases in diversity from intermediate to high levels of NPP [2,6]. Furthermore, the shape of productivity-diversity relationships is typically dependent on the spatial grain of the analysis. For example, linear increases in species richness in large sampling units (such as degree grids) contrast with the multiple ways in which species richness within local habitats responds to NPP [7][8][9]. Thus, although productivity-diversity relationships are widespread [5], the underlying mechanisms still need to be resolved [10].
In other cases however, biodiversity is more strongly correlated to ambient energy than to NPP [11][12][13][14]. In contrast to NPP, ambient energy does not include water availability and can be expressed as annual mean temperature (T). Possible mechanisms for increasing species richness with ambient energy include tropical niche conservatism [13], dispersal limitation after glaciation [14] and metabolic theory [11,15]. Based on the evolutionary origin of many taxa in tropical climates, their occurrence in cooler climates depends on the evolution of coldtolerance. Within larger taxonomic groups, communities in warm climates often include many basal taxa, while communities in temperate to cold climates are increasingly restricted to few derived taxa. This results in a positive relationship of species richness and T. Analyses of the phylogenetic structure of communities provide evidence for this mechanism, notably in butterflies [13]. Another effect of historical climate on contemporary richness patterns has been described for European dung beetles [14]. Their limit of thermal tolerance during the last glacial maximum marks a strong change in current richness with low species numbers north of this limit. Thus, limited colonization of areas with historically unsuitable climate can contribute to current correlations between richness and ambient energy. Metabolic theory provides an additional mechanism for higher species richness in warmer climates based on the influence of temperature on metabolic rates and on rates of speciation [11]. Metabolic theory predicts a direct, monotonic relationship of species richness with T, whereby the slope of log e -transformed species richness with inverse T is predicted to be 20.65 [15].
A third framework for broad-scale patterns of species richness is biological relativity to water-energy dynamics. It is based on the dependence of all life on the availability of water in a liquid form [16][17]. This framework suggests that species richness is proportional to (a) the availability of liquid water (annual rainfall, R), and (b) the lowest monthly value of potential evapotranspiration (PET min ). This ''interim general model'' explains almost 80% of the variation in species richness of trees and shrubs in eastern and southern Africa [16]. At higher latitudes (Europe and North America) the effect of rainfall on species richness of trees dominates (r = 0.64), and effects of PET min or other energy variables are non-significant [18].
The final determinant of species richness patterns considered here is spatial environmental heterogeneity. Its importance can be expected to increase with increasing size of sampling units (spatial grain [2,19]). For example, only limited spatial climatic heterogeneity can be expected along a 21 m transect within one habitat type (the small sampling grain in our study, see Methods). In contrast, large differences in internal spatial climatic heterogeneity exist among countries (the large sampling grain in our study). For example, T differs by 13.3uC between the warmest and the coldest 10610 arc minute square within Switzerland -more than ten times the difference within Denmark, which has a similar surface area ( Table 1). Effects of spatial environmental heterogeneity on species richness can be explained with niche theory, which assumes different environmental preferences and tolerances among species [20,21]. Among environmental variables, climate is central to the distribution and persistence of species worldwide [22]. In particular, the distribution of numerous taxa is influenced by T, for example plants [23], beetles [14,24], spiders [25] and birds [26]. Thus, heterogeneity of environmental conditions (such as T) within large sampling units can enhance overall species richness by providing suitable conditions for larger numbers of species with different ecological niches. Increased species richness due to an increased importance of spatial environmental heterogeneity will consequently shift the focus from alpha diversity (local species richness) to beta diversity (turnover of species in space [27]).
Here, we explore species richness-environment relationships of European arthropods. We combine the results of a continent-wide standardised sampling programme of local ground-dwelling arthropod communities (local grain) with existing coarse-grained country inventories (country grain) of comparable spatial extent. At the local grain ground-dwelling ants (Formicidae), beetles (Coleoptera), bugs (Hemiptera) and spiders (Araneae) were sampled in 32 habitats at 16 locations across Europe, ranging from boreal to Mediterranean in climate. At the country grain, inventories of 25 European countries were obtained for the same groups [28]. We used these data to test each of the above hypotheses by first comparing the explanatory power of productivity, ambient energy, the interim general model, and the best possible statistical model (drawn from all variables) for biodiversity within local habitats (alpha diversity). Secondly, we compared the potential for environmental heterogeneity to explain species turnover across locations (beta diversity). Thirdly, we compared the explanatory power of productivity, ambient energy, the interim general model, spatial environmental heterogeneity and the best possible statistical model (drawn from all variables) for biodiversity within countries (gamma diversity).

Ethics
-Field sites were selected and established within the EU FP6 ALARM project to form a long-lasting research network -Each field site had a site manager, responsible for contacts to local authorities and/or land owners -In most cases the land belonged to regional research stations -More detailed descriptions can be found in reference [29] -Protected areas or rare habitats were not included into this field site network -We did not include protected species.
The employed pitfall traps capture invertebrates, with no protected species affected in the habitats we sampled. No permissions are needed to use pitfall traps outside of protected areas.

Data
Ground-dwelling arthropods were sampled in 32 habitats at 16 locations across Europe (Table 2, Fig. 1E). Thirteen locations were part of the ALARM field site network [30], and three sites (Bern, Silkeborg, Wien) were added to fill geographic gaps. As far as possible, one forest as an example of a near-natural habitat and one cereal field as an example of an intensive agricultural habitat were sampled in each location. When unavailable, other nearnatural habitats (scrubland or extensive grassland) and other intensive agricultural habitats (intensive grassland or olive grove) were sampled instead ( Table 2). Trapping took place in 2006 and started five days after the beginning of the vegetation period (the onset of growth in the majority of plant species) in each location [31]. In each habitat, eight pitfall traps of 7 cm diameter were placed along a transect and separated by 3 m from each other. The traps were filled with 0.1 L of a 4% formaldehyde solution, to which sodium dodecyl sulphate was added as detergent. Three sampling periods of two weeks were separated by pauses of two weeks. Adult arthropods were identified to species level by specialists. The following groups were considered: ants (Hymenoptera: Formicidae), beetles (Coleoptera: Carabidae, Curculionoidea, and Staphylinidae), bugs (Hemiptera: Auchenorrhyncha and Heteroptera), and spiders (Araneae). In addition to the analysis of taxonomic groups, we divided the studied arthropods into trophic groups according to the dominant feeding type in the respective family. Herbivore families were all Curculionoidea, all Auchenorrhyncha and the heteropteran families Berytidae, Cydnidae, Lygaeidae, Miridae, Pentatomidae, Piesmatidae, Plataspidae, Pyrrhocoridae, Rhopalidae, Scutelleridae and Tingidae. Spiders, ants and the remaining beetle and bug families were carnivores. Detritivores could not be analysed because this feeding type did not dominate in any of the sampled families. Species numbers within 25 European countries were taken from the Fauna Europaea database [28] (Table 1, Fig. 1F). Countries smaller than 30,000 km 2 were excluded, as were all countries for which the known number of arthropod species lay below the 95% confidence interval of the species -log area relationship, indicating effects of insularity (Republic of Ireland) or incomplete knowledge of the arthropod fauna (Ukraine, Belarus). We further excluded islands such as the Balearic Islands, Corsica, Greek Islands, Sardinia and Sicily from the respective mainland areas.
Environmental variables were extracted from a European gridded data set with a monthly time step and a spatial resolution of 10610 arc minutes, which corresponds approximately to 16 km [32]. The primary climatic variables temperature (T, Fig. 1B) and precipitation were constructed through interpolation from station observations [33]. Annual rainfall (R; Fig. 1C) was calculated as the sum of precipitation in all months with an average temperature .0uC [18]. NPP (Fig. 1A) was estimated by running the LPJ-GUESS ecosystem model [34,35] with the same climate data, and parameterized for the potential natural vegetation of Europe [36,37]. LPJ-GUESS and the closely related LPJ-DGVM [34] have formerly been shown to reproduce observed variations in NPP across various types of vegetation and climates [38][39][40]. PET min (Fig. 1D) represents potential evapotranspiration of the coldest month of each year, and was calculated using the Thornthwaite equation, which only requires knowledge of air temperature [41]. We used long-term annual means of NPP, T, R and PET min from 1971 to 2000. For analysing species richness within habitats and species turnover across habitats, we used NPP, T, R and PET min values of the grid in which the habitats were located. For analysing species richness within countries, explanatory variables were averaged across all grids of the country. In addition, spatial heterogeneity in NPP, T, R and PET min were calculated as ranges for each country by subtracting the minimal from the maximal value, respectively (i.e. difference between the grid cells with the highest and lowest value).

Analysis
Local species richness and species richness in European countries was analysed using linear models with standardised explanatory terms (mean = zero, standard deviation = 1) in the statistical environment R version 2.12.0 [42]. To account for possible differences in sampling efficiency between locations, we used the number of captured individuals N as a covariate in the analyses of local species richness. We accounted for possible effects of spatial autocorrelation of the habitats within the locations and among the locations and countries with generalised least squares [43] with spatial simultaneous autoregressive error models [44]. Models were based on neighbourhood matrices that linked the two habitats within a location and each location with at least one other location for the local grain analyses and allowed each country to be in the neighbourhood of at least one other country, i.e. at a Environmental variables: area (in 10 3 km 2 , islands excluded), averages of net primary productivity (NPP; in kg C per m 2 per year) and annual mean temperature (T; in uC), annual rainfall (R; in mm) and minimal potential evapotranspiration (PET min ; in mm), and the spatial heterogeneity in primary productivity (NPP range ), annual mean temperature (T range ), annual rainfall (R range ), minimal potential evapotranspiration (PET minrange ). Species richness for spiders (S a ), beetles (S c ), bugs (S h ), ants (S f ), herbivores (S her ) and carnivores (S car ). doi:10.1371/journal.pone.0045875.t001   Table 1 maximum distance of about 850 km from centre to centre. For this we used the package spdep [45]. In addition, we calculated Moran's I correlograms for the residuals of models with and without correction for spatial autocorrelation to assess if the tested theories miss important spatially structured environmental variables. Missing crucial spatially structured environmental variables will lead to significant residual spatial autocorrelation of the uncorrected models. Separate models of local species richness were calculated according to productivity-diversity relationships, ambient energy and the interim general model, plus one ''Best'' model in which all explanatory variables relevant for the different theories (NPP, ) entered the initial models. We included interactions of all linear terms with habitat to test if there are different responses in the different habitat types. We identified the minimal adequate models by a backwards variable selection procedure according to the second order Akaike information criterion (AICc). Linear terms were always kept in the model when the respective quadratic term increased the model fit. In cases of high collinearity (Pearsson r .0.5) of linear terms (see Tables 3, 4) we calculated separate models always containing only one of these terms, and the best model was chosen based on the AICc model selection criteria. This restricted the models to only one energy variable (either T or PET min ). It also reduced the risk of overfitting, which is considerable given the low numbers of replicates (N = 16 locations and N = 25 countries).
As relationships of species richness with NPP can be either linear or hump-shaped [3,7,10], we allowed the quadratic term of NPP to remain in the productivity-diversity relationship models if that resulted in lower AICc values. Ambient energy models were calculated using untransformed species richness and T. To test predictions made by metabolic theory, we calculated the slope of log e (species richness) with 1/[0.0000862(273+T)] for comparison with the predicted slope of 20.65 [15,46]. We used the first version of the interim general model (IGM1), where species richness is explained by a linear term of R plus a linear and quadratic term of PET min [16]. At the country grain, R and PET min had a Pearson correlation coefficient of 0.6. Nevertheless, we also tested the full model including both variables for means of completeness. We also performed an influence analysis using Cook's distance. If data points had a Cook's distance .0.5, indicating disproportional weight in the regression analysis, then the effect of excluding those data points from the model was examined.
Relationships of species turnover with environmental variables were analysed using Mantel tests [47]. We used presence-absence data of the trapped species per location (natural and disturbed habitat combined). Community dissimilarities were calculated as Morisita-Horn distances and related to Euclidean environmental distances between all possible pairs of sites. Each environmental variable was tested separately, and separate tests were calculated for spiders, beetles, bugs, ants, all groups combined, herbivores and carnivores using the function mantel in the package vegan (default settings [48]). The significance was based on Monte Carlo tests with 999 permutations.
Species richness of arthropods in European countries was analysed in a similar way as local species richness, with the following additions. We corrected for the area of the countries by including log 10 (area) as an additional explanatory variable in all models [49]. In addition to productivity-diversity relationships, ambient energy and interim general models, we calculated a model containing spatial environmental heterogeneity. Variables considered were the ranges of NPP, T, R and PET min . However, all these range variables were highly intercorrelated (Table 4) and thus it was not possible to include them simultaneously in one model. Therefore, we calculated separate models always containing one of these terms and selected the best model according to AICc. Again, we calculated ''Best'' models in which all explanatory variables relevant for the different theories entered the set of initial models containing only one of the range variables characterising spatial environmental heterogeneity.

Species Richness in Local Habitats
The samples contained 33223 individuals of our focal taxa that comprised 83 ant species (Formicidae), 444 beetle species (Coleoptera), 185 bug species (Hemiptera) and 354 spider species (Araneae). Relationships between species richness and environmental variables at the local grain were highly variable ( Table 5). All arthropods combined and carnivores considered separately showed a hump-shaped relationship with NPP, while lacking a significant effect of ambient energy or variables from the interim general model. The ambient energy model was best for spiders, while beetles conformed most to the energy term of the interim general model. Bugs and herbivores were not significantly affected by any environmental variable. Ants showed a negative relationship to R plus an interactive effect of T and habitat type. Ant species richness increased with T in near-natural habitats, but did not change significantly with T in intensive agriculture ( Fig. 2A). A more detailed examination of the response of ants to climate and habitat type has been given elsewhere [50]. Habitat type had a significant effect only on spiders, with higher species richness in near-natural than in intensive agricultural habitats (Table 5). In contrast to ants, there was a significant interactive effect of T and habitat type on spider richness, whereby species richness increased with T in agricultural habitats but did not significantly change with T in near-natural habitats (Fig. 2B)

Species Turnover
Species turnover across the 16 locations was most strongly correlated with differences in T (Table 6). Correlations were highest for spiders, carnivores, beetles and all groups combined, followed by ants and herbivores. Only bugs showed no significant relationship of species turnover with environmental variables. In the remaining groups, correlations of species turnover with differences in T were at least 49% stronger than with any other environmental variable. Correlations of species turnover with differences in NPP and/or PET min were significant, but substantially less strong than those with differences in T.

Species Richness in Countries
All tested environmental variables showed some significant effects on species richness within countries (Table 7). With respect to productivity-diversity relationships, species richness of spiders, beetles, bugs and all groups combined increased with NPP. The ambient energy-models revealed increased species richness of beetles, bugs, ants, herbivores, carnivores and all groups combined with T. However, for all groups combined this relationship was no longer significant when an overly influential data point (Portugal) was omitted from the analysis. As for metabolic theory, the slope of log e (species richness) versus T 21 came close to the predicted value of 20.65 for ants (20. With respect to biological relativity to water-energy dynamics, the full interim general model including both variables (R and linear and quadratic terms of PET min ) had constantly higher AICc values than simplified models ( Table 7). The reduced interim general model for spiders revealed an unexpected negative response of species richness to rainfall ( Table 7). The remaining groups showed hump-shaped relationships with PET min in accordance with biological relativity to water-energy dynamics. Portugal was overly influential in the interim general model for ants, and no significant model remained after its removal. There were consistent positive relationships of arthropod richness with spatial environmental heterogeneity (Table 7). T range gave a better model fit than NPP range , R range and PET minrange in all cases.
Models with free variable selection always combined variables from several theories ( Table 7). They were statistically superior to any single theory according to their higher explanatory power and lower AICc values (DAICc .6.1). Residuals showed significant spatial autocorrelation in the majority of single-theory models, but in none of the models with free variable selection ( Table 7, ''Best''). This suggests that the models with free variable selection included the majority of relevant variables while single theories tended to miss crucial information. In accordance with a high importance of spatial environmental heterogeneity, species richness increased with T range in all models with free variable selection (Table 7).

Discussion
Although our results from the local habitat samples were variable with respect to environmental effects on species richness, there were numerous significant effects of environmental variables on species richness of the same groups at the country grain. This suggests that productivity-diversity relationships, ambient energy, the interim general model and spatial environmental heterogeneity all contribute to the explanation of arthropod species richness of European countries. However, model selection according to AICc identified T range as the strongest predictor of arthropod richness across all studied groups. Independent support for a strong role of spatial heterogeneity in T comes from the significant relationship of species turnover across locations with differences in T. If species turnover across locations is driven by T, then countries with a high T range will contain higher beta diversity and consequently more species in total than countries with more uniform temperatures. In the following, we will discuss the different theories for broad-scale gradients in species richness and what can be concluded from our data.

Productivity-diversity Relationships
The observed increase of species richness in countries with NPP is in accordance with the majority of studies on broad-scale relationships of species richness with climate [5]. In contrast, relationships of species richness with NPP at the local grain were hump-shaped and restricted to carnivores and to the sum of all arthropod species. This accords with a generally reduced effect size [51], and with a transition from monotonous to hump-shaped productivity-diversity relationships towards small spatial grain [7]. The differences between grains cannot be explained by differences in gradient length, because NPP varied only slightly more among locations than among countries (locations: 0.39 g C m 22 a 21 in Lesvos to 0.66 g C m 22 a 21 near Kraków; countries: 0.41 g C m 22 a 21 in Norway to 0.65 g C m 22 a 21 in France; Tables 1,2). Nevertheless, the NPP gradient was relatively short at both grains. When gradients include areas with very low NPP, stronger effects at the local grain would be expected. Thus, any conclusions with respect to small-scale productivity-diversity relationships from our data should be made with caution.

Ambient Energy
We found significant effects of T on species richness at both local and country grains. The increase of ant species richness with T in both near-natural local habitats and in countries accords with ambient energy theories. In contrast, spider richness increased with T in intensive agricultural habitats but not at the country grain. This suggests that the richness pattern of spiders in agricultural habitats is not indicative of other habitat types and thus of limited relevance for their overall species richness in countries. Beetles, bugs, carnivores and herbivores showed significant positive relationships with T at the country, but not at the local grain. T can influence species via its effect on NPP. Low temperatures limit terrestrial NPP in temperate to arctic climates [5], and large parts of our study area lie in the temperate to boreal region. However, the latitudinal gradient of NPP was unimodal in our study, with decreasing NPP from temperate to Mediterranean climate (Fig. 1A). At the country scale, all groups except spiders and beetles were significantly affected by T in addition to NPP, suggesting direct effects of ambient energy on species richness. Across countries, the slope of log e (species richness) versus T 21 accorded with metabolic theory only for ants. Species richness of ants has been found earlier to conform with metabolic theory [12,52]. Our results suggest that this may be an exception rather than the rule among terrestrial arthropods. Based on widespread nonlinearity, geographic and taxonomic dependence of temperature-richness relationships, metabolic theory has been more generally questioned [46]. Deviations from metabolic theory can be due to violations of its assumptions [53]. For example, the assumption of body size invariance with temperature is violated by the significant increase of spider body size across Europe with temperature [54]. Tests of tropical niche conservatism and dispersal limitation after glaciation require phylogenetic analyses that exceed the scope of the current investigation [13,14]. These historical climatic explanations predict more basal taxa in warm climates and high richness of few derived taxa in cooler climate. Thus, the richness of higher taxonomic categories such as families should increase more strongly towards warm climate than the number of species. Such a pattern is present in our spider data: while species richness shows no significant relationship with ambient energy at the country grain, the number of spider families increases with T across countries (t 1,23 = 4.9, p,0.001). This indicates niche conservatism in warmer climates and encourages more detailed phylogenetic exploration of the distribution of European arthropods.

Temperature, Species Turnover and Spatial Environmental Heterogeneity
In contrast to its variable effect on species richness within local habitats (alpha diversity), differences in T had strong effects on species turnover across locations (beta diversity; Table 6). This role of T confirms that it represents an important niche dimension of European arthropods [24,25]. The number of available niches in a given area thus correlates to the range of temperatures present in that area. Our results are consistent with the ideas that species richness is enhanced by (i) elevational range and (ii) habitat heterogeneity in an area [19,55,56]. Elevation is a main driver of temperature variation in mountains, leading to correlations between elevational range and spatial heterogeneity in temperature (r = 0.92 for the countries studied here). The occurrence of the same ecosystems at similar temperatures across the world that have contrasting elevations demonstrates that temperature is more crucial for biodiversity than elevation per se [57]. Being difficult to measure, habitat heterogeneity is often determined by the number of distinguishable vegetation types present in an area [56]. In near-  Table 7. Models of species richness in N = 25 European countries. natural situations, vegetation types are in turn determined by the environmental preferences of their constituent plant species, including their temperature preference [58]. Thus, spatial environmental heterogeneity, elevational range and habitat heterogeneity are interrelated, and we presume that climate often has the most direct influence on biodiversity. In our study, variability in T had a dominant effect on species turnover and gamma diversity. Apart from the existence of more niches along large temperature gradients, climatic heterogeneity can also buffer species extinctions by allowing species confronted with climatic fluctuations to relocate to suitable climatic refuges [59]. The dominant role of T range for arthropod species richness in European countries is in accordance with both enhanced niche availability and reduced extinction during climatic fluctuations.

Biological Relativity to Water-energy Dynamics
Biological relativity to water-energy dynamics is expressed in the interim general model [16,17]. It predicts increasing species richness with R (water term) and a unimodal relationship of species richness with PET min (energy term). We found only partial support for IGM, since reduced models always resulted in lower AICc values compared to the full model (Tables 5,7). Numbers of spider species per country decreased with R in the interim general model, but increased with R in the model with free variable selection. Dominant effects of the energy term in the interim general models at the country grain were replaced by group-specific positive, negative (ants) or absent (herbivores) effects of R in the models with free variable selection (Table 7). Given these inconsistencies, the interim general model provided no robust explanation of arthropod richness in our study. Nevertheless, the numerous significant relationships with R and PET min suggest that the interim general model may apply to European arthropods, but that larger datasets are necessary to disentangle its components. The dominant effect of the water term in the models with free variable selection is in accordance with Hawkins et al. [18], who found that the energy term of the interim general model becomes dispensable in temperate to arctic climate.

Herbivores versus Carnivores
The strength of latitudinal diversity gradients has been found to increase across trophic levels [51]. In our study, the major difference between herbivores and carnivores was at the local grain, where carnivores showed a hump-shaped relationship to NPP and herbivores no significant relationship at all. At the country grain, the results for herbivores and carnivores were similar to each other and to those of bugs and all arthropods combined. This suggests that the observed differences between arthropod taxa are due to their different phylogeny or other lifehistory traits rather than caused by their trophic position.

Countries versus Local Habitats
The stronger and more consistent effects in countries versus locations are in accordance with the general decrease of speciesrichness environment relationships towards small spatial grain [51,56,60]. Ecological processes are scale-dependent [61], and effects of the studied broad-scale environmental conditions may affect regional species pools rather than local assemblages. However, species pools can affect local species richness [62], especially in mobile organisms such as the studied arthropod groups. In addition, some of the mechanisms to explain broadscale patterns in species richness include local processes such as resource partitioning (productivity-diversity relationships), metabolism (metabolic theory), or cold tolerance (tropical niche conservatism). Thus, contrasts such as the variable role of ambient Table 7. Cont. Differences in data quality could have contributed to the stronger effects in countries versus locations [10]. First, our pitfall traps sampled only ground-dwelling arthropods over a limited time period. The 33223 sampled individuals of 1066 species already represent a major identification effort. Nevertheless, sampling intensity per habitat corresponded only to a minimal effort that is expected to encompass around 75% of all species attainable at the respective site with pitfall traps [63]. Even more problematic can be differences in sampling efficiency between sites due, for example, to weather or habitat structure [64]. We reduced these differences by applying strictly standardized sampling methods and by including the number of individuals captured per habitat as a factor in the models of species richness. By using (log) individual numbers as a factor, we assume that true abundances are similar across habitats and that observed differences in individual numbers are due to variation in sampling efficiency. However, true abundances may differ. Results at the local grain would change strongly if individual numbers were excluded from the models -the only consistency being increased ant richness with ambient energy and variable effects between groups (results not shown). This highlights that the difficulty to obtain large, representative arthropod samples from defined areas remains a main obstacle in community ecology. Accordingly, conclusions from the results at the local grain should be drawn with care. In contrast to these sampling issues at the local grain, species inventories of countries are the results of many decades of research and have reached asymptotes in almost all cases [65]. A second difference in data quality between locations and countries relates to the environmental data. Interpolated climatic variables will result in relatively accurate values across large areas such as countries, but have only limited accuracy at the grain of local habitats. Sources of error include spatially and temporally unpredictable factors such as rainfall, as well as anthropogenic effects on productivity, especially in disturbed habitats. Thus, sampling error at the locations is probably much higher than at the country grain, underlining the need for additional high-quality inventories of invertebrates along environmental gradients of large spatial extent.

Main Conclusions
Our study supports the scale-dependence of species richnessenvironment relationships. While relationships of local species richness with environmental variables were contingent on the arthropod group, species richness patterns at the country grain were more consistent and partly supported all tested theories. Niche theory provides a plausible link between the two grains: On the one hand, differences in temperature were the best correlate of species turnover across locations. On the other hand, spatial heterogeneity in annual mean temperature had the strongest effects on arthropod diversity within European countries. These two independent findings suggest that temperature is an important niche dimension and that countries with wider ranges in annual mean temperature provide a greater breadth of niche space and so can support larger numbers of arthropod species. Unless environmental heterogeneity is constant across sampling units (thinkable e.g. in marine environments), we strongly suggest that studies with large sampling units take into account environmental heterogeneity, just as studies with variable area of sampling units nowadays routinely consider area [49].