Scale-dependent environmental filtering of ground-dwelling predators in winter wheat and adjacent set-aside areas in Hungary

Agricultural intensification may act as an environmental filter shaping invertebrate assemblages at multiple spatial scales. However, it is not fully understood which scale is the most influential. Therefore, we utilized a hierarchical approach to examine the effect of local management (inorganic fertilization and soil properties; within-field scale), habitat type (winter wheat field and set-aside field; between-field scale) and landscape complexity (landscape scale) on assemblage structure and functional diversity of two important groups of natural enemies, carabids and spiders, in a cultivated lowland landscape in Hungary. Environmental filtering affected natural enemies at different spatial scales; likely as a result of enemies’ different dispersal ability and sensitivity to fertilizer use. Carabids were strongly affected at the within-field scale: positively by soil pH, negatively by soil organic matter and fertilization. At the between-field scale, carabids had higher activity density in the set-aside fields than in the winter wheat fields and simple landscapes enhanced carabids diversity, species richness and activity density at the landscape scale. Spiders were more abundant and species-rich in the set-aside fields than in the winter wheat fields. Although highly mobile (macropterous) carabids might disperse to arable crops from greater distances, while spiders possibly depended more on the proximity of set-aside fields, the winter wheat fields (where pest control should be delivered) were utilized mostly by common agrobiont species. Increasing crop heterogeneity within arable fields could be a potential option to increase the diversity of carabids and spiders in the studied region.


Introduction
Reduced landscape heterogeneity, as a result of progressive agricultural intensification, is the major driver of biodiversity loss in many countries across Europe (Kleijn et al. 2009;Fahrig et al. 2011). It also acts as an environmental filter (Duflot et al. 2014;Gámez-Virués et al. 2015) that reduces the capacity of landscapes to support large species-pools of arthropods (Ewers and Didham 2006;Hendrickx et al. 2007;Tscharntke et al. 2012). Environmental filtering increases a particular kind of species turnover when the most specialized species are replaced with generalist ones (biotic homogenization: Olden et al. 2004;Clavel et al. 2011;Duflot et al. 2014;Gámez-Virués et al. 2015). Such homogenized assemblages in the ecosystem can provide a limited range of ecosystem services due to their reduced functional diversity (Olden et al. 2004;Cadotte and Tucker 2017). Since these species share the same generalized traits, some resources may become inaccessible, while competition for another set of resources will increase. Additionally, more specialized Electronic supplementary material The online version of this article (doi:https ://doi.org/10.1007/s1084 1-020-00249 -9) contains supplementary material, which is available to authorized users. arthropods with low dispersal ability cannot reach complementary resources and they are out-competed by more mobile species (Hendrickx et al. 2007;Tscharntke et al. 2012).
Environmental filtering caused by agricultural intensification has multiple components and affects arthropod assemblages at different spatial scales (e.g. Clough et al. 2005;Gallé et al. 2019). Croplands are subject to periodical withinfield disturbances, such as tillage, fertilization and pesticide applications, with long-term detrimental consequences on arthropods (e.g. Batáry et al. 2008;Sądej et al. 2012). In some cases, fields are left bare during winter (Kleijn et al. 2009). The availability of resources in croplands is therefore highly dynamic, resulting in strong spatiotemporal variation in habitat utilization (Kromp 1999;Rand et al. 2006). When crops are disturbed, arthropods, as well as other farmland animals, may depend on perennial non-crop habitats adjacent to arable fields (between-field scale). Various types of non-crop habitats (e.g. field boundaries, fallows, grasslands, woodlands and hedgerows) provide refuge, alternative prey sources, nesting and overwintering sites (e.g. Landis et al. 2000;Clough et al. 2005;Kovács-Hostyánszki et al. 2011;Kovács-Hostyánszki and Báldi 2012). During the growing season, these habitats may act as sources for arthropods for crop recolonization (Clough et al. 2005;Tscharntke et al. 2005;Bianchi et al. 2006;Schellhorn et al. 2014). Hence, the potential species pool filtered by landscape simplification might be further influenced at local scales by habitat type and crop management resulting in locally distinct arthropod assemblages (Duflot et al. 2014).
Based on the aforementioned studies, it is evident that agricultural intensification and environmental filtering significantly shape the structure of invertebrate assemblages in croplands. Overall, this resulted in a decrease in arthropod diversity (Hendrickx et al. 2007;Gámez-Virués et al. 2015). In order to counteract or, at least, mitigate the negative effect of environmental filtering, the European Community introduced specific agri-environmental schemes (AES, Batáry et al. 2015). In Hungary, set-aside fields were introduced as part of the national AES program in 2002 following accession to the European Community (Ángyán et al. 2003). Set-aside fields were designed specifically for the protection of the great bustard (Otis tarda) and other protected farmland bird species, but their establishment has also been beneficial for non-target groups of animals and plants (Kovács-Hostyánszki et al. 2011;Kovács-Hostyánszki and Báldi 2012;Tóth et al. 2018). In consequence, AES can also act to enhance ecosystem services, such as pest control or pollination. This ecological intensification (Bommarco et al. 2013) can sustain agricultural production with the least possible adverse effects on the environment by complementing or replacing artificial inputs with ecosystem services to increase crop productivity (Kleijn et al. 2019).
Nevertheless, the effectiveness of specific AES habitats at alleviating environmental filtering and delivering benefiting ecosystem services is still poorly understood because the environmental impact assessment of AES focuses on broad measures (Batáry et al. 2015). Moreover, it is also relatively unknown how arthropod assemblages and especially their functional traits are filtered at multiple spatial scales because studies addressing functional-trait based approaches are still rare (Gámez-Virués et al. 2015, but see Aviron et al. 2005 andGallé et al. 2019).
These knowledge gaps prompted us to examine the assemblage structure of natural enemies using a hierarchical approach, from local to landscape scales. As model groups, we selected ground beetles (Coleoptera: Carabidae, hereafter carabids) and spiders (Araneae) because of their significant contribution to the biodiversity (including functional diversity) in agroecosystems (Greenstone 1999;Kromp 1999). Although carabids and spiders are important ecosystem service providers (Bommarco et al. 2013;Gagic et al. 2017;Martin et al. 2019), the effect of environmental filtering on them can vary depending on their functional traits. Utilizing a functional-trait based approach may provide better insight into the various processes in dynamic of agricultural landscapes than only considering taxonomic differences, such as species richness and/or abundance (Knapp et al. 2019;Magura and Lövei 2019;Martin et al. 2019). Whereas some functional traits, such as body size, wing morphology (as a proxy for dispersal) and habitat affinity (a proxy for fidelity) reflect environmental filtering through dispersal and their connection with feeding is indirect, hunting strategy is directly linked with the potential impact on pest populations (Schmidt-Entling and Döbeli 2009;Duflot et al. 2014;Anjum-Zubair et al. 2015;Gámez-Virués et al. 2015;Gallé et al. 2018Gallé et al. , 2019Martin et al. 2019).
In this study, we investigated how the local management in the form of inorganic fertilization and soil properties (within-field scale), habitat type including the newlyestablished AES habitat (between-field scale) and landscape complexity (landscape scale) affected carabid and spider assemblages in Hungary. In particular, we focused on their species richness, abundance and functional diversity. In the latter case, species were categorized according to their functional traits (body size, habitat affinity, wing morphology or hunting strategy), and these groups were analyzed separately to distinguish between the effects attributable to a specific group. We aimed to test the following hypotheses: (i) The abundance, species richness and functional diversity of both carabids and spiders would be higher in croplands with improved soil fertility (i.e. neutral pH, high amount of soil organic matter, hereafter SOM) and without inorganic fertilization. Contrarily, these assemblage patterns would be lower in croplands, where intensive use of fertilizer attempted to compensate the SOM deficit. (ii) Newly-established set-aside fields, as the most common AES habitat adopted in the study region, adjacent to croplands provide a more suitable environment for carabids and spiders and would act as source habitats for recolonization in arable crops. Thus we expect that presence of set-asides improves either the taxonomic or the functional diversity of these pest control providers. (iii) The high amount of semi-natural habitats (e.g. forests, hedgerows and grasslands) around arable lands leads to higher complexity at the landscape scale. This improved spatial heterogeneity in landscape composition is expected to support high taxonomic and functional diversity of carabids and spiders through improved dispersal and colonization conditions compared to simple landscapes.

Experimental area
Our study was conducted in the lowland of Heves region in north-eastern Hungary. The area is characterized by arable fields (dominated by winter wheat, spring and winter barley, and maize) and areas under the AES program. In general, the number of available sites within the study region was restricted. We selected seven winter wheat fields with a size of 5-10 ha (for GPS coordinates, see ESM- Table 1). Each winter wheat field had an adjacent, newly-established (1-3 years old) set-aside field pair (size 1.98-5.43 ha). Set-aside fields are the most common AES habitats adopted in Heves region (Kovács-Hostyánszki and Báldi 2012). Originally, these set-aside fields used to be arable fields and have been taken out of production for a 5-year contract period. At the beginning of this contract period, they were sown by a legume-grass seed mixture composed of Medicago sativa, grasses of Lolium sp. and Festuca sp. to ensure an organic matter enrichment (Ángyán et al. 2003). Apart from legumes and grass, several wild plants germinated from the seed bank. The vegetation cover was similar for all set-aside fields; they were mown once in a year in the second half of June, leaving the hay on the ground. No insecticides or herbicides were used during the contract period.
In each of seven winter wheat fields, we assigned a plot of 45 × 20 m in the central section adjacent to a set-aside field to avoid any edge effects from other areas. Farmers were not allowed to use insecticides within this experimental plot, while normal conventional management was applied to the rest of the field. Each experimental plot was divided into two subplots of 20 × 20 m with a 5 m separation strip to avoid any interference between treatments (Fig. 1). One subplot was treated with NPK fertilizer at the usual rate of 95 kg N/ ha in mid-April 2014, while the other half, without fertilizer input, served as control. Since fertilization treatment was not applied to the set-aside field, our sampling design was not fully factorial between habitats and fertilizer use. The distance between study site pairs (i.e. pairs of winter wheat field and adjacent set-aside field) was between 1.21 and 10.08 km (mean: 4.09 km).
In each winter wheat subplot, 15 soil samples were taken randomly when the wheat plants reached the stem elongation stage 10 + at the BBCH scale (Zadoks et al. 1974). Soil samples were stored at 4 °C until SOM and pH were measured. Prior to analyses, samples from each experimental field were sieved (4 mm mesh size), blended and then combined for analyses. Soil analyses were conducted according to Hungarian standards (no. 08-0452 and 08-0206-2) and the literature (Mason 1983). SOM content varied from 2.45 to 4.09%, and the pH range was 4.74-6.59.
In 1 km radius around the experimental winter wheat fields, we estimated the complexity of landscape using a GIS database implemented by local authorities. Landscapes with > 20% of semi-natural habitats (semi-natural grasslands, forests, hedgerows and other non-crop elements) within this radius were classified as complex, while those with < 20% of such habitats were classified as simple (Batáry et al. 2008Tscharntke et al. 2011). In addition, the average distance between fields was greater than 2 km, thus exceptional partial overlap of landscapes did not bias the spatial independence. Four of the seven landscape windows around the seven winter wheat fields were complex, and three were simple (see ESM- Table 1).

Arthropod sampling and functional traits
In each winter wheat subplot, we installed four pitfall traps to sample carabids and spiders at 10 and 20 m (two traps 10 m apart in each distance, Fig. 1) from the inner edge of the wheat fields adjacent to the set-aside field. Four pitfall traps were also placed in set-aside fields using the same sampling arrangement. In total, we installed 84 traps: 28 in fertilized winter wheat (four traps per subplot, seven subplots in total), 28 in control winter wheat and 28 in set-aside fields. Pitfall traps were plastic cups of 114 mm diameter filled with approximately 250 ml of 50% propylene glycol solution saturated with NaCl and a drop of odorless detergent to reduce surface tension. Green plastic roofs were installed to protect traps from litter and rain.
All traps were open for 10 days at the beginning of June 2014 (between 3rd and 13th June). The collected material was sorted in the laboratory and spiders and carabids were identified to species level using keys in Loksa (1969Loksa ( , 1972, Heimer and Nentwig (1991), Nentwig et al. (2019), and Hůrka (1996), respectively. Even though our sampling season was rather short, it coincided with the flowering and milky ripening stages of winter wheat one month before harvest, and with the late spring activity peaks of ground-dwelling arthropods (Anjum-Zubair et al. 2015;Gagic et al. 2017), making it suitable for our comparative purposes.
Activity density (counts from pitfall traps are more properly described as activity density than abundance) and species richness of both predatory groups per trap were calculated. For each carabid species, we recorded the following functional traits: body size (small: < 10 mm, medium: 10.1-15 mm or large: > 15.1 mm), flight ability (macropterous or brachypterous) and habitat affinity (generalist, open or forest habitat specialist) according to Hůrka (1996). For spiders, body size (small: < 5 mm, medium: 5.1-10 mm or large: > 10.1 mm), hunting strategy (running spiders, stalkers, and sheet-or aerial-web builders) and habitat affinity (generalist, open or forest habitat specialist) was obtained from Uetz et al. (1999), Buchar and Růžička (2002) and Nentwig et al. (2019).

Statistical analyses
The expected species richness of carabids and spiders was estimated in set-aside fields and winter wheat fields with the Chao estimator (Chao 1987) using the specaccum function in the 'vegan' package (Oksanen et al. 2017). The standard deviations were generated from 10,000 reshufflings of the sample order. The diversity of carabid and spider assemblages at the within-field (fertilization treatment), the between-field (habitat type), and landscape scale (landscape complexity) was compared by the Rényi one-parametric diversity index family. This approach does not consider only one numerical value (Tóthmérész 1995), but calculates an assemblage diversity profile by plotting several diversities against a scale parameter (Lövei 2005;Ricotta 2005;Elek and Tóthmérész 2010). When the scale parameter, α = 0, the value of the Rényi diversity equals the logarithm of species richness. At this stage, the index is very sensitive to the presence of rare species. When α ≈ 1, the Rényi diversity corresponds to the Shannon diversity. As the scale parameter increases, the Rényi Fig. 1 Spatial arrangement of pitfall traps in winter wheat subplots and set-aside fields diversity becomes more sensitive to the relative abundance of the more common species. At α = 2, the plotted value equals the quadratic (Simpson) diversity. Towards α ≈ + ∞, values approach the inverse of the Berger-Parker dominance index, which is determined by the abundance of the most common species only. If the Rényi diversity profiles of two assemblages intersect, the assemblages can still be compared, but the comparison should be given for different section of profiles in relation to rarity and dominance of the species set (Lövei 2005). For the calculation of the Rényi diversity, we used the renyi function from the 'vegan' package (Oksanen et al. 2017).
We applied separate analyses for activity density and species richness of carabids and spiders and selected functional traits of both groups. For functional trait analyses, a subset was created for each category in a given functional trait. For instance, the trait of body size of carabids was included as a factor with three trait levels: small, medium-sized and large carabids. However, only categories with at least 200 individuals were used for further analyses to achieve robustness of the data; thus, small, medium-sized, macropterous and open-habitat carabids and small, medium-sized, running, stalking, open-habitat and generalist spiders were included in the analyses. First, we tested the predictive power of all studied explanatory variables and possible correlations using the random forest algorithm from the 'party' package (Strobl et al. 2008). Then, we used linear mixed-effect models, where species richness of carabids and spiders and activity density of running spiders and the whole assemblage were squareroot transformed; while the rest of response variables were log transformed to normalize data distribution. We used the lme function from the 'nlme' package (Pinheiro et al. 2017) and included treatment (two levels: control or fertilizer), habitat type (two levels: winter wheat field or set-aside field), pH (continuous), SOM (continuous) and landscape complexity (two levels: simple or complex) as explanatory variables and study site ID (seven levels) as a random effect. For each carabid and spider data subset (i.e. activity density, species richness and category in a given functional trait), we fitted a maximal model that included all explanatory variables and meaningful interactions. Following this, non-significant interactions and variables were gradually removed (model simplification) to reach the most parsimonious model. ANOVA with a χ 2 test was used to test the significance of explanatory variables (the Anova function from the 'car' package, Fox and Weisberg 2010). To select the most parsimonious model within each subset, we calculated the weight of each model based on the AICc using the model.sel function of the 'MuMIn' package (Bartoń 2016). Model selection showed that the most parsimonious model was also the best one in all analyzed subsets. All analyses given above were conducted in R 3.4.1 (R Core Team 2017).

Assemblage characteristics for both predatory arthropod groups
In total, we collected 3878 carabids of 39 species (ESM -Tables 2 and 3

Estimated species richness and assemblage diversity
The species richness estimation revealed that it did not change substantially after the 20th carabid sample in setaside fields (27.25 species ± 0.7 SD, ESM- Fig. 1a), but no saturation was evident for carabids in winter wheat fields (ESM- Fig. 1b) and for spiders in both habitat types (ESM- Fig. 1c, ESM-Fig. 1d). The Rényi diversity profiles at the within-field scale revealed a negative effect of fertilizer use on carabids (Fig. 2a), but for spiders, this effect was not strong (Fig. 2b). At the between-field scale, the diversities of both groups were similar, with the diversity profiles intersecting: for rare species, the diversity was higher in the winter wheat fields than in the set-aside ones. Considering the more abundant species of both carabids and spiders, their diversity was higher in set-aside fields than in winter wheat ( Fig. 2c and d). At the landscape scale, carabid diversity was at the same level for rare species in both landscape types, however for the most abundant species, it was higher in simple landscapes (Fig. 2e). Spider diversity was slightly higher in the complex landscapes than the simple ones (Fig. 2f).

Responses to within-field, between-field and landscape factors
At the within-field scale, soil pH had a positive ( Fig. 3a and  b) while the SOM a slight negative effect ( Fig. 3c and d) Although the inorganic fertilizer reduced species richness of carabids by 35.3% (Fig. 4a), their activity density remained unaffected (Table 1). At between-field scale, the activity density of all carabids, and that of small, macropterous and open-habitat species was significantly higher in set-aside fields than in winter wheat fields, regardless of fertilization. The activity density of the assemblage was 49.5% lower in winter wheat than in set-asides fields (Fig. 4b). At the landscape scale, carabid species richness was 18.7% lower in complex landscapes than in the simple ones, while their activity density was reduced by 27.8% (Fig. 4c and d). Similarly, small, medium-sized and open habitat species were significantly more abundant in simple landscapes than in the complex ones (Table 1).
Spider species richness as well as activity density was significantly affected only at the between-field scale ( Table 1). Species richness was 25.1% lower and activity density was reduced by 44.1% in winter wheat fields than in adjacent setaside fields ( Fig. 5a and b). Small-sized species and stalkers also had low activity densities in winter wheat fields (Table 1). Medium-sized, running and generalist spiders did not significantly respond to any variables. Fertilizer treatment, soil properties or landscape complexity had no impact on the examined characteristics of spider assemblages.

Discussion
In this cultivated landscape, assemblages of both carabids and spiders were characterized by highly mobile, small, open-habitat-preferring species which can indicate a highly dynamic environment. The dominance of few but very abundant agrobiont species in both groups suggested the influence of strong environmental filters acting at different spatial scales. While spiders did not respond at the landscape scale, simple landscapes supported more abundant but functionally less diverse carabid assemblages. Highly mobile (macropterous) carabids can colonize arable crops from greater distances, while ground-dwelling spiders depend on the proximity of set-aside fields. Overall, winter wheat fields (where pest control should be delivered) were utilized mostly by

Within-field scale
Only carabids were affected at the within-field scale. Here, the species richness of carabids was negatively affected by inorganic fertilization, conforming to the prediction by our first hypothesis (i). The impact of fertilizer use on carabids is primarily indirect, acting through soil-and crop-related features, such as vegetation (crop) structure, diversity of weeds, and overall plant productivity (Siemann 1998;Kromp 1999). This can lead to changes in the availability of below-and above-ground food resources for carabids and consequently affects their spatial distribution (Hole et al. 2005;Pakeman and Stockan 2014;Mader et al. 2018). Individual carabid species such as Amara aenea (De Geer, 1774) and Chlaenius festivus (Panzer, 1796) might react to these microhabitat changes by avoiding fertilized patches, whilst other species showed slight preferences for either the fertilized or the control subplot (ESM -Table 3). In addition, lower species richness and diversity in the fertilized subplots suggest the existence of environmental filtering and a potential decrease in biological pest control within arable fields . Our findings are similar to those of Sądej et al. (2012), who observed fewer carabid species in plots fertilized by inorganic nitrogen but non-significant drop in activity density. Also, Honczarenko (1975) reported that some carabid species totally avoided plots with high doses of inorganic fertilizer and were concentrated in non-fertilized control plots.
In general, higher pH had a positive, while SOM a negative effect on carabids, providing partial support for our first hypothesis (i). Soils with neutral pH are characterized by a higher diversity of various ground-dwelling organisms including carabids (Kuperman 1996;Scheu 2001;Sądej et al. 2008;Tamburini et al. 2016). Although the fauna living in soils with a high content of SOM can serve as a good food source for ground-dwelling predators, thus resulting in higher numbers of carabids (Pfiffiner and Luka 2000; Sądej et al. 2012), we found the opposite pattern. It is likely that this is due to shifts in carabid assemblage composition mediated by open habitat species. Decreasing pH and increasing SOM content in mucky, peat-muck and peat soils negatively affect the activity density of mesophilous carabid species that were less tolerant to changing moisture conditions (Nietupski et al. 2010). In humus-rich soils, those mesophilous species are replaced by hygrophilous ones that were more associated with open habitats. We presume that the predominance of open habitat species in this study drove the overall response of the carabid assemblages via their negative response to increasing SOM content.

Between-field scale
At the between-field scale, in accordance with our second hypothesis (ii), we found that newly-established set-aside fields locally increased the activity density of carabids with the exception of medium-sized species (see below). The majority of carabid species used both habitats, but had lower activity densities in winter wheat than in set-aside fields. This might suggest that the spring colonization of winter wheat (at least partially) depended on adjacent set-aside fields. Set-aside fields, similarly to other non-crop habitats,  (Pfiffner and Luka 2000;Fournier and Loreau 2001;Bianchi et al. 2006). Concerning medium-sized carabids, their activity density did not significantly change between habitats. This might be due to the predominance of Poecilus cupreus (L. 1758) (ESM -Table 3), a common species occurring primarily in arable fields (Mader et al. 2018). This agrobiont species can be found at boundaries with set-aside fields and other non-crop habitats, as well as in centers of arable fields, in high abundances, suggesting insensitivity to the filtering effect at between-field scale (Anjum-Zubair et al. 2015;. For spiders, only our second hypothesis (ii) was supported, because spiders were affected only at the betweenfield scale. It seems that a major driver for spider distribution in this study was the presence of set-aside fields, while inorganic fertilization, soil conditions in winter wheat fields, and landscape complexity were not important. Whereas some species of carabids can hibernate in arable fields and do not need adjacent non-crop habitats for completing their life cycles (Holland et al. 2007), in-crop-only persistence is very rare in spiders (Clough et al. 2005) and has been observed only in few species (Mestre et al. 2018). The spillover of spiders into arable crops therefore heavily depends on the proximity of adjacent set-aside fields and/or other types of non-crop habitats (Clough et al. 2005;Bianchi et al. 2006;Gallé et al. 2018). However, a significant drop in both species richness and activity density in winter wheat fields suggest the existence of an environmental filter between habitats because only some species, mostly belonging to running spiders, were found in winter wheat (ESM -Table 4); the other species stayed in the set-aside fields. Running spiders that moved into the wheat fields did so regardless of inorganic fertilization or SOM. Active hunters, such as P. agrestis, one of the most dominant running spiders in agroecosystems in Central Europe (Kiss and Samu 2000), can easily disperse from adjoining overwintering sites and colonize even large fields (Samu and Szinetár 2002). This can be also the reason why running and open habitat spiders did not react at the between-field scale. Spiders of other hunting strategies, such as stalkers, could be positively affected by the complex vegetation structure of set-aside fields resulting in higher abundances in this habitat type (Luff and Rushton 1989;Cole et al. 2005;Schellhorn et al. 2014). These results might be indirect evidence that the structured vegetation and high plant diversity in our set-aside fields (the vegetation cover was a legume-grass mixture with several wild plants from the seed bank) and extensive management without any chemical input favored stalking and ambush-hunting species, such as Ozyptila simplex (O. Pickard-Cambridge, 1862), by providing better hunting conditions via a higher availability of shelters than in winter wheat fields.

Landscape scale
At the landscape scale, contrary to the third hypothesis (iii), we found no effect of landscape complexity on spiders, while carabid assemblages were more species-rich and abundant in simple landscapes than in complex ones. Carabid assemblages in agroecosystems might be unsaturated due to the depletion of specialist species that have a moderate dispersal capacity (Aviron et al. 2005;Duflot et al. 2014). Structurally complex landscapes mitigate this depletion due to the presence of semi-natural and other non-crop habitats, supporting recolonization of temporarily suitable crop fields and counteracting biotic homogenization (Hendrickx et al. 2007;Fahrig et al. 2011;Gáméz-Virués et al. 2015). However, our results showed the opposite pattern. For small, open habitat-preferring species, complex landscapes might act as a barrier and reduce their dispersal due to unfavorable Response of species richness (a) and activity density (b) of spiders to habitat type. Whiskers represent 95% confidence interval and capital letters at the bottom of bars indicate significant differences conditions in semi-natural habitats (Mauremooto et al. 1995;Al Hassan et al. 2013). In simple landscapes, a very high species-turnover of common agrobionts compensates this species depletion, leading to the species-rich and abundant but functionally less diverse assemblage (Aviron et al. 2005;Martin et al. 2019). In spiders, the lack of response at the landscape scale might be due to the distance considered in this study (1 km around winter wheat fields). It is possible that spider assemblages were affected at a higher or lower scale. In future studies, measuring the landscape heterogeneity at multiple distances (500 m, 1 km, 3 km from the study area) could help to identify the scale of landscape effects on ground-dwelling arthropods in agroecosystems.

Conclusions
Carabids and spiders are important providers of ecosystem services to promote ecological intensification. Particular AES programs should be planned in a broader context including multiple spatial scales to sustain these services. Set-aside fields may support natural enemies locally but the existence of scale-dependent environmental filters may diminish their efficiency. Although our results did not provide direct evidence, we suggest that dispersal ability and sensitivity to fertilization play a crucial role in trait-specific responses of carabids and spiders. It is also important to consider the cost/benefit ratio of the proposed AES habitat. Establishing and maintaining set-aside fields is relatively costly and not very popular among farmers due to reduction of the area under production and consequent lower yield (see Kleijn et al. 2019). Since winter wheat fields were utilized mostly by common agrobiont carabids and spiders as a result of scale-dependent environmental filters, management adjustments should be focused on arable fields rather than their boundaries. Increasing crop diversity and vegetation heterogeneity may support the taxonomic and functional diversity of carabids and spiders in crops where pest control should be achieved. Growing multiple crops simultaneously in the same field with reduced fertilizer use could be a potential option for farmers that want to support natural pest control (Kleijn et al. 2019). Nevertheless, these suggestions may only be appropriate for this study region. In addition, short-term sampling scheme provides only a snapshot in the complex dynamics of ground-dwelling arthropods in agroecosystems and thus the generalization of our results might be limited. The effect of environmental filtering on natural enemies may vary in relation to different parts of the season, crop type, AES habitat and/or region of Europe. Various AESs can be effective tools for ecosystem intensification but they have to be carefully designed and targeted (Batáry et al. 2015).