The Influence of Landscape Heterogeneity - Ground Beetles (Coleoptera: Carabidae) in Fthiotida, Central Greece

Abstract Pitfall traps were used to sample Carabidae in agricultural land of the Spercheios valley, Fthiotida, Central Greece. Four pairs of cultivated fields were sampled. One field of each pair was located in a heterogeneous area and the other in a more homogeneous area. Heterogeneous areas were composed of small fields. They had high percentages of non-cropped habitats and a high diversity of land use types. Homogeneous areas were composed of larger fields. They had lower percentages of non-cropped habitats and a lower diversity of land use types. One pair of fields had been planted with cotton, one with maize, one with olives and one with wheat. Altogether 28 carabid species were recorded. This paper describes the study areas, the sampling methods used and presents the data collected during the study. Neither heterogeneous nor homogeneous areas had consistently higher abundance levels, activity density levels, species richness levels or diversity levels. However, significant differences were seen in some of the comparisons between heterogeneous and homogeneous areas.


Introduction
The level of heterogeneity in agricultural landscapes can influence farmland wildlife (Benton et al. 2003;Tscharntke et al. 2005;Fahrig et al. 2011). Heterogeneous farmland, forming a mosaic of different land use types, has been seen to benefit invertebrate groups such as the Aranea (Sunderland and Samu 2000), the Lepidoptera (Weibull et al. 2000) and the Carabidae (Östman et al. 2001), in terms of their abundance, diversity and condition.
Locally, soil moisture and soil type have a large impact on carabid distributions (Meissner 1984;Thiele 1977). However, some aspects of landscape heterogeneity will also influence the Carabidae. These aspects are field size (Kromp 1999), the presence of non-cropped habitat (Pollard 1968;Sotherton 1985) and land use diversity (Östman et al. 2001).
Small fields are thought to be easier for carabids to recolonize after disturbance, due to the shorter dispersal distances involved. Additionally, landscapes with small fields are likely to have high levels of land use diversity, which will create refuges for carabids in times of disturbance. This is because cultivation practices take place at different times in different crops, so a diverse landscape will always have some undisturbed habitat, while cultivation practices take place elsewhere (Kromp 1999).
Another aspect of landscape heterogeneity that is important for carabids, is the presence of non-cropped habitat. This may take the form of grassy field margins, hedgerows and areas of semi-natural habitat, such as fallow, woodland and wasteland. Although carabids live in the crops during the vegetation season, they are known to use field margins to hibernate in (Pollard 1968;Sotherton 1985). Where this occurs, carabids are less vulnerable to ploughing, harrowing and insecticide application, which usually take place in the autumn and winter (Holland and Luff 2000). In addition to this, field margins are known to be important as breeding sites for many carabid species (Desender and Turin 1989).
In Greece, agriculture is often extensive and small in scale, with large-scale, intensive farming occurring only on the flatter and more fertile land. This means that the country has relatively high levels of habitat richness (Eurostat 2012). Although this situation may be beneficial to wildlife, economies of scale mean that such areas are expensive to farm and have, for a long time, been prone to land abandonment and agricultural intensification (Kasimis et al. 2003;Louloudis and Maraveyas 1997). Today, these problems have been compounded with the effects of the economic crisis, meaning that small-scale farmers and the landscapes they have helped to create are now under even greater pressure (Eurostat 2012). So far, little work has been conducted on the carabid communities of farmland in Central Greece, in spite of the benefits associated with such diverse agricultural landscapes.
This study aims to identify closely matched and situated areas of heterogeneous and homogeneous farmland. These areas will then be compared using landscape analysis, so that the different aspects of their heterogeneity can be examined. Finally, matched fields of the same crop types within the heterogeneous and homogeneous areas will be compared, to see how they differ in terms of their carabid abundance, activity density, species richness and diversity. In this way, it may be determined whether high levels of landscape heterogeneity benefit Carabidae at a local level, within individual fields.

Study Areas and Sampled Fields
The study was conducted on agricultural land in the Spercheios valley, Fthiotida, Central Greece. Four pairs of fields were sampled for Carabidae between May and October in 2007. One pair of the sampled fields had been planted with cotton, one with maize, one with olives and one with wheat.
One field of each pair was located in a heterogeneous area, which had small field sizes, a large amount of non-cropped habitat and a high level of land use diversity. The other field of each pair was located in a more homogeneous area, which had larger field sizes, less non-cropped habitat and a lower level of land use diversity (Suppl. material 2).
The heterogeneous and homogeneous areas were matched, in order to be as similar as possible regarding all factors apart from their heterogeneity. They were matched according to their mean elevation, their distances to villages and roads, as well as to large expanses of woodlandand wasteland. Fig. 1 shows the relative positions of the different study areas. Here it can be seen that the sizes of the study areas varied and their shapes were irregular. This was because their outlines delineated the places that were considered most heterogeneous or homogeneous when compared to each other and to the surrounding countryside. When the landscape analysis was performed; however, comparisons took place between equally proportioned areas.
The sampled fields chosen from within the heterogeneous and homogeneous areas were also matched regarding their crop type, soil type, previous crop type, agrochemical treatment, harvesting time, elevation and whether or not they were irrigated. The data used to match the fields, along with the size and location of each field, are provided in Table 1   The land use maps in Figs 2, 3, 4 show the positions of the individual fields within the study areas. Each land use type is represented by a polygon of a different pattern and the sampled fields are marked with red circles. These maps were made using satellite photographs, taken from an eye altitude of 1.5 km. As homogeneous area 3b was very large, the map of this area was divided in two, so that two maps were made, one of area 3bi, containing the sampled wheat field and one of area 3bii, containing the sampled cotton field. The key to these maps is shown in Fig. 5. Land use maps of areas 1a, 1b, 2a and 2b. The sampled maize and olive fields within these areas are marked with red circle.
The Influence of Landscape Heterogeneity -Ground Beetles (Coleoptera: ...  Land use maps of areas 3bi and 3bii. The sampled wheat and cotton fields within these areas are marked with red circle. The program FRAGSTATS (McGarigal et al. 2002) was used to conduct landscape analysis on each of the study areas, using the data presented in Figs 2, 3, 4. This resulted in the calculation of six landscape metrics (Suppl. material 2). The data for the landscape analysis were obtained by placing equal sized rectangles over the land use maps as a separate layer. The map data within these rectangles were then converted into numerical data and input into FRAGSTATS. For the metrics "land use diversity" and "land use richness" the data from 10 ha rectangles were used. This was because the results of these metrics would vary depending on the percentage of background landscape included in the analysis. For the other metrics; however, where this was not the case, 50 ha rectangles were used, meaning that all of the collected data could be analyzed.
The metric "land use aggregation" which concerns land use configuration, is dimensionless, so therefore lacks units. The same is true for the metric "land use similarity", which indicates how similar the sampled fields were, regarding land use, to other fields in the surrounding areas (McGarigal et al. 2002).

Sampling Procedure
The locations of the pitfall traps within the sampled fields are listed as geographic coordinates in Table 1. These coordinates represent the mid points of the lines of traps. Most fields were large enough to accommodate a single line of 10 traps, situated 10 m apart, 5 m from one of the field margins. However, the olive grove and the maize field in the heterogeneous areas were too small to accommodate a single line of 10 traps. So in these fields two rows of 5 traps were used, positioned on opposite sides of the fields, 10 m apart and 5 m from each of the field margins.
Not all of the pairs of fields were sampled during every 15 day sampling period. This was because access to the fields depended on the cultivation practices taking place at the time. Irrigation, spraying, harvesting, ploughing, pruning and fertilizer application would all prevent access to the fields and the setting of traps. Initially, a pair of alfalfa fields was also sampled, but the frequency of harvesting meant that most of these traps were destroyed before they could be collected. So sampling in these fields was not continued. The Influence of Landscape Heterogeneity -Ground Beetles (Coleoptera: ...
If one field of a pair was inaccessible for a given 15 day period, then the other field of that pair was not sampled either. This insured that comparisons between heterogeneous and homogeneous areas were fair, as the same amount of data was obtained in both fields, at the same time of year. In all, 440 traps were set and 391 were recovered successfully. Successful traps were those that were not flooded by irrigation, destroyed by other farming practices or dug up by animals. The dates covered by the 15 day sampling periods, as well as the numbers of traps set and successfully recovered are shown in Table 2.  22 May to 6 June 10 10 9 23 June to 8 July 10 10 9 9 July to 24 July 10 10 9 25 July to 9 Aug From the 391 successful traps, 320 (40 traps from each field) were chosen for use in the data analysis (Table 2). The number 40 was chosen because the least sampled field, "cotton a" provided a total of 42 successful traps. To begin with, two traps, chosen at random, were removed from the "cotton a" data set. Then 40 equivalent traps were chosen from the "cotton b" data set. These traps were taken from the same sampling periods as those in the "cotton a" data set. This procedure was repeated for all of the field pairs. It meant that the comparisons between fields in heterogeneous and homogeneous areas were always fair, as they were made using data from the same number of traps, taken during the same time of year.
The traps themselves were made out of 250 ml plastic cups. These had a depth of 10 cm and a rim diameter of 7.3 cm. They were part filled with ethylene glycol, and covered with wooden lids to prevent flooding and the capture of larger, non-target species. This left a gap of 2 cm between the trap rims and their lids. The Carabidae were identified to species level using Arndt et al. (2011) and Trautner and Geigenmüller (1987). Then a subsample of the specimens, containing all of the sampled species was checked by a taxonomist from the University of Athens. A small number of specimens were also sent to the Natural History Museum in London for verification (Acknowledgements). All specimens were stored in alcohol and are now kept in the private collection of the author.

Data Analysis
For each field, the data from 40 traps were combined, then the number of carabid species found in each field were recorded, along with their relative abundance (n). The annual activity density (ADa) (Brandmayr et al. 2005) was also calculated for each species, using the following equation: The diversity of carabid species was calculated for each field using the Simpson's Diversity Index (D), which is presented here as the complement (1-D).
Carabid abundance and species richness were also calculated for each trap used in the data analysis. The resulting 40-trap data sets were tested using the Anderson-Darling test. This showed that the data sets were rarely normally distributed, even after transformation, meaning that nonparametric statistics were used for significance testing. Mann-Whitney U tests were used to determine the significance of differences between heterogeneous and homogeneous areas.
To compare between the different crop types, carabid abundance and species richness were again calculated for each trap. Then the data from both fields of each crop type were combined, resulting in four, 80-trap data sets, one for each crop type. Kruskal-Wallis tests were then used to determined the significance of differences between the cotton, maize, olive and wheat cultivations. tot tot Locality:

Results
Suppl. material 1 shows the abundance (n) and the annual activity density (ADa) of each of the carabid species sampled in the different fields. Table 3 summarizes these data and also presents the carabid species richness and diversity levels for all of the sampled fields. Suppl. material 2 presents the results of the six landscape metrics calculated for each of the study areas. Notes: It digs burrows under stones and is mostly phytophagous (Trautner and Geigenmüller 1987). In this study it was found rarely (n = 2) in the wheat field of the homogeneous area. th rd Table 3.
The total abundance (N), total annual activity density (ADa), species richness and diversity (1-D) of Carabidae in each field.
Notes: This species was found in the cotton field in the heterogeneous area (n = 1), the maize field in the heterogeneous area (n = 51), the olive grove in the homogeneous area (n = 4), the wheat field in the heterogeneous area (n = 1) and the wheat field in the homogeneous area (n = 5).

Notes:
In agricultural areas, it is often found on arable land, pastureland and alfalfa. It is usually absent in fields with abundant weed cover (Holland and Luff 2000). It likes open areas and is an autumn breeder (Thiele 1977). In this study, it was found in the maize field in the homogeneous area (n = 1), the olive grove in the homogeneous area (n = 3) and the wheat field in the homogeneous area (n = 1).

Notes:
In this study, this species was only found rarely (n = 3) in the olive grove in the homogeneous area.
Notes: In this study, it was found in the olive grove in the heterogeneous area (n = 3), the olive grove in the homogeneous area (n = 7), the wheat field in the heterogeneous area (n = 3) and the wheat field in the homogeneous area (n = 12).
Notes: It is a phytophagous and xerophilous species (Pavlícek et al. 2005). In this study, it was found in the maize field in the homogeneous area (n = 2), the olive grove in the heterogeneous area (n = 1) and the wheat field in the heterogeneous area (n = 5).

Discussion
Neither the heterogeneous nor the homogeneous areas had consistently higher abundance, activity density, species richness, or diversity levels. This suggests that the level of heterogeneity of the study areas did not have a great influence on the carabid communities of the sampled fields. Areas with small field sizes, large amounts of noncropped habitat and high land use diversity did not appear to benefit the Carabidae.
Additionally, there did not seem to be an association between any of the landscape metrics in Suppl. material 2 and the levels of carabid abundance, activity density, richness or diversity. Even the "percentage of non-cropped habitat", which is a factor thought to be particularly influential for the Carabidae and more broadly for farmland biodiversity (Holland and Luff 2000;Desender and Turin 1989;Fahrig et al. 2011) appeared not to have a consistent influence.
The results of this study do not agree with those reviewed by Benton et al. (2003), Tscharntke et al. (2005) and Fahrig et al. (2011). However, this may be because the relationship between heterogeneity and biodiversity is thought not to be linear (Fahrig et al. 2011), in which case, this study's design may have been inappropriate. Fahrig et al. (2011) suggest that when studying the influence of landscape heterogeneity, areas should be chosen that form a gradient in terms of their heterogeneity. Here though, a paired design was used and different crop types were sampled in the different area pairs. Fahrig et al. (2011) also recommend sampling a much larger number of areas (40-60). In this study, a small number of areas (only 6) were sampled, so replication occurred within the sampled fields, rather than at the area level. While this provided information at a local level, about the sampled fields themselves, it would not have provided information at a landscape level. It is likely therefore, that if a larger number of areas had been sampled, the overall influence of landscape heterogeneity would have been clearer.
A related issue is that too few land use types may have been sampled in this study. Heterogeneity is believed to enhance biodiversity, through different species being associated with different land use types, at different times in their lives (Dunning et al. 1992;Pollard 1968;Sotherton 1985;Desender and Turin 1989). This suggests that if more land use types had been sampled, a more consistent pattern may have been seen in the results.
Additionally, the need to compare closely situated areas, matched regarding other factors apart from their heterogeneity, meant that it was difficult to choose areas that differed greatly in all aspects of their heterogeneity. For example, although "land use diversity" was high and "land use similarity" was low in all of the heterogeneous areas, "land use richness" did not always follow this pattern. For the wheat comparison, "land use richness" was slightly higher in the homogeneous area, while for the cotton comparison "land use richness" was the same in both areas (Suppl. material 2). This implies that the heterogeneous and homogeneous areas may not have been different enough to influence the Carabidae present within them. Having said that though, highly homogeneous farmland was not present in the Spercheios Valley, so comparing widely different heterogeneity levels would have meant comparing areas from different regions, something which may have resulted in unfair comparisons.
It is also possible, as mentioned by Fahrig et al. (2011), that while appearing to be heterogeneous in satellite photographs, the heterogeneous areas in this study were not functionally heterogeneous from the point of view of the Carabidae. The different land use types in the heterogeneous areas may have actually provided similar resources for the Carabidae. Woodland and tree crops, for example, may have been functionally similar land use types, due to their low disturbance levels and the rich plant communities they supported.
Despite these issues, significant differences were seen in some of the comparisons between heterogeneous and homogeneous areas (Following Subsections). These results are interesting as they provide information about how the Carabidae within individual fields are affected by heterogeneity at the landscape level.

Relative Abundance
For the cotton fields, carabid abundance per trap was significantly higher in the homogeneous area (U = 1178, p = 0.0003). For the maize fields though, there was not a significant difference between the heterogeneous and homogeneous area (U = 831.5, p = 0.7642). The olive groves had significantly higher carabid abundance in the heterogeneous area (U = 586, p = 0.0404). The result of this comparison was one of the few that showed a positive influence of heterogeneity. For the wheat fields though, there was not a significant difference between the heterogeneous and homogeneous area (U = 693, p = 0.3077). Finally, when the data from all of the fields in each type of area were combined, there was no significant difference between heterogeneous and homogeneous areas (U = 12721 p = 0.9283).
On the whole, crop type appeared to have a greater influence on the Carabidae than did heterogeneity. The results of the Kruskal-Wallis test show that there was a highly significant difference in carabid abundance between the different crop types (H = 34, p = <0.0001). This was probably due to the specific microclimates created by each crop type and to the differences in husbandry that each crop required (Holland and Luff 2000). For instance, the maize fields had the highest carabid abundance levels of all crop types, probably because they were irrigated throughout the growing season and because they were not treated with insecticides or herbicides. The cotton fields, on the other hand, had the lowest overall abundance levels, probably due to the treatment of these fields with the insecticide Phosalone. Differences between crop types; however, should not have affected comparisons between heterogeneity levels, as crop type and husbandry were matched between heterogeneous and homogeneous areas.
When the abundances of individual species were considered, significant differences were seen between some heterogeneous and homogeneous areas, but these were not consistent for all of the crop type comparisons. Relative abundance patterns also varied depending on the carabid species, with some species showing higher abundances in a heterogeneous area, and others in a homogeneous area. This may have been due to differences in dispersal ability, causing variation in they way individual species experienced heterogeneity (Tscharntke et al. 2005).
For the cotton comparison, Pterostichus (Platysma) niger was found in significantly higher numbers in the homogeneous area (U = 1041, p = 0.0209). Then when the influence of crop type was examined, significantly greater numbers of this species were seen in the maize fields (H = 136.13, p = <0.0001). This may have been largely due to the maize fields being irrigated regularly throughout the summer, which would have provided favourable conditions for this species.
For the olive comparison, Microlestes luctuosus was found in significantly greater numbers in the heterogeneous area (U = 592, p = 0.0455). This species is known to prefer areas of tall vegetation (Fadda et al. 2008), which it would have found to a greater degree in the heterogeneous area, due to the presence of greater percentages of non-cropped habitat (Suppl. material 2). There was also a significant difference in the abundance of M. luctuosus between the different crop types (H = 73.6, p = <0.0001). It was trapped in the greatest numbers in the wheat fields. This may have been because these fields remained undisturbed as stubble from early June, when harvest took place, until October, when they were ploughed. Leaving stubbles intact for many months is something that has been seen to be beneficial for many groups of species, including carabids (Evans et al. 2002).
For the maize comparison, Poecilus cupreus was found in significantly greater numbers in the heterogeneous area (U = 521, p = 0.0074). This may have been due to a preference for areas of non-cropped habitat at certain times of life (Thiele 1977) and could indicate a requirement for high levels of heterogeneity when inhabiting agricultural land.
For the cotton comparison, Tapinopterus (Tapinopterus) taborskyi was found in significantly greater numbers in the homogeneous area (U = 1050, p = 0.0164). This may suggest a preference for low heterogeneity levels. However, as this species has only ever been found on Oiti mountain (Arndt et al. 2011), which is adjacent to the Spercheios Valley, this seems unlikely. Perhaps its high abundance level in the cotton field in the homogeneous area was due to there being an easy dispersal route from Oiti in that direction, possibly due to the prevailing wind direction.
For the maize comparison, Cylindera germanica was found in significantly higher numbers in the homogeneous area (U = 1199, p = 0.0001). This was probably because C. germanica is accustomed to living in agricultural areas (Varvara and Apostol 2008) so a low level of heterogeneity would not be detrimental to it. There was also a significant difference seen in the abundance of this species between the different crop types (H = 17.19, p = 0.0006), with far greater numbers being trapped in the maize fields. This was probably due to the frequent irrigation required by the maize crop. This would make sense, as C. germanica is also known to be a common inhabitant of damp areas, such as flood plains (Arndt et al. 2011).
The most common species in this study, Harpalus (Pseudoophonus) rufipes, did not show any significant differences between heterogeneous and homogeneous areas. This is probably due to it being a very common inhabitant of agricultural land (Luff 2007;Thiele 1977), and so being well adapted to living in all kinds of agricultural landscapes. However, significantly greater numbers of this species were seen in the maize fields (H = 136.13, p = <0.0001). Again this was probably due to the frequent irrigation of this crop type, which would have provided moisture during the dry summer period.

Species Richness
For the cotton fields, there were significantly higher numbers of carabid species per trap in the homogeneous area (U = 1190, p = 0.0002), something that does not indicate a positive influence of landscape heterogeneity. For the maize and wheat fields, no significant differences were seen between heterogeneous and homogeneous areas (maize U = 673.5, p = 0.2263, wheat U = 759.5, p = 0.7039). For the olive groves too, there was not a significant difference between the heterogeneous and homogeneous area (U = 644.5, p = 0.1362). Additionally, when the data from all of the fields in each type of area were combined, there was no significant difference between heterogeneous and homogeneous areas (U = 13255, p = 0.5823).
The results of the Kruskal-Wallis test; however, showed that there was a highly significant difference in species numbers per trap between the different crop types (H = 90, p = <0.0001). The olive groves and the wheat fields had the highest overall richness levels of the four crop types. This may have been because these fields were organically farmed (Eyre et al. 2009;Clark 1999;Pfiffner and Niggli 1996). Another important factor may have been that the olive groves were left relatively undisturbed for most of the year, allowing the development of rich plant communities. The only cultivation practices that took place during the sampling season were pruning and organic manure application. The wheat fields were also left undisturbed for a large part of the summer. After harvest in early June, their stubbles were left intact until they were ploughed in October, which may have conferred benefits to the Carabidae living within them (Evans et al. 2002). The cotton fields; however, had the lowest species richness levels of all the crop types, something that was probably due to the treatment of these fields with the insecticide Phosalone.

Diversity
Carabid diversity levels were not consistently higher in either the heterogeneous or the homogeneous areas. Nor was there a clear association between carabid diversity levels and the levels of any of the landscape metrics in Suppl. material 2.
Neither were diversity levels consistently higher in any one crop type. Although overall diversity levels were highest in the olives groves, the least disturbed of all of the different cultivations. The highest diversity levels for individual fields were seen in the wheat and olive cultivations in the homogeneous areas. However, the lowest diversity level was seen in the wheat field in the heterogeneous area.
Again these results suggest that the level of heterogeneity had little influence on the Carabidae in the sampled fields. Although, as previously discussed, these findings may have been due to limitations imposed by this study's design.