NeoBiota 11: 33–49, doi: 10.3897/neobiota.11.2060
Modelling the distribution of the invasive Roesel’s bush-cricket (Metrioptera roeselii) in a fragmented landscape
Sonja Preuss 1, Anna Cassel-Lundhagen 1, Åsa Berggren 1
1 Department of Ecology, PO Box 7044, Swedish University of Agricultural Sciences SE-75007 Uppsala, Sweden

Corresponding author: Sonja Preuss (sonja.preuss@slu.se)

Academic editor: I. Kühn 

received 11 September 2011 | accepted 21 November 2011 | Published 30 November 2011


(C) 2011 Sonja Preuss. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.


For reference, use of the paginated PDF or printed version of this article is recommended.

Abstract

The development of conservation strategies to mitigate the impact of invasive species requires knowledge of the species ecology and distribution. This is, however, often lacking as collecting biological data may be both time-consuming and resource intensive. Species distribution models can offer a solution to this dilemma by analysing the species-environment relationship with help of Geographic information systems (GIS). In this study, we model the distribution of the non-native bush-cricket Metrioptera roeselii in the agricultural landscape in mid-Swedenwhere the species has been rapidly expanding in its range since the 1990s. We extract ecologically relevant landscape variables from Swedish CORINE land-cover maps and use species presence-absence data from large-scale surveys to construct a species distribution model (SDM). The aim of the study is to increase the knowledge of the species range expansion pattern by examining how its distribution is affected by landscape composition and structure, and to evaluate SDM performance at two different spatial scales. We found that models including data on a scale of 1 × 1 km were able to explain more of the variation in species distribution than those on the local scale (10 m buffer on each side of surveyed road). The amount of grassland in the landscape, estimated from the area of arable land, pasture and rural settlements, was a good predictor of the presence of the species on both scales. The measurements of landscape structure – linear elements and fragmentation - gave ambivalent results which differed from previous small scaled studies on species dispersal behaviour and occupancy patterns. The models had good predictive ability and showed that areas dominated by agricultural fields and their associated grassland edges have a high probability being colonised by the species. Our study identified important landscape variables that explain the distribution of Metrioptera roeselii in Mid-Sweden that may also be important to other range expanding orthopteran species. This work will serve as a foundation for future analyses of species spread and ecological processes during range expansion.

Keywords

Orthoptera, presence-absence data, spatial scale, landscape structure, land use

Introduction

The development of effective strategies to manage the spread of invasive organisms requires data on species habitat preferences and knowledge of how landscape characteristics influence species dispersal and establishment (Cote and Reynolds 2002; Rosin et al. 2011). However, the collection of fine-detailed distribution data over large scales is time consuming and logistically challenging, hence data is missing for many species (Jimenez-Valverde et al. 2008). Management decisions have often to be taken swiftly (Morueta-Holme et al. 2010) and species distribution modelling becomes a handy tool when dealing with limited observation data and large spatial and temporal extents (Guisan and Thuiller 2005). By modelling species distribution as a function of ecologically relevant data on climate conditions and/ or landscape characteristics, it is possible to describe occupancy patterns and predict species range expansions (Hein et al. 2007; Early et al. 2008; De Groot et al. 2009; Bonter et al. 2010). Estimates of current and future species distributions rely on: (1) the strength of the relationship between environmental variables and the organism in question (Cote and Reynolds 2002), and (2) the availability of ecological relevant environmental data that can be applied at a range of geographic scales (Scott et al. 2002). It is also important to consider the impact of scale on the performance of the models (Scott et al. 2002), i.e. we need to know which environmental predictors give the best estimates for species presence at a given spatial scale.

Some species of orthopterans (grasshoppers and bush-crickets) have recently shown a rapid response to changed environmental conditions and are invading new areas outside their common range (Sword et al. 2008; Bazazi et al. 2011). Orthopterans are well suited for studying distribution patterns across a range of spatial and temporal scales, because they are relatively easy to survey and their ecology is well studied (Ingrisch and Köhler 1998; Gwynne 2001; Hein et al. 2003; Holzhauer et al. 2006). Metrioptera roeselii is an example of a range expanding species in northern Europe (Simmons and Thomas 2004; Gardiner 2009; Hochkirch and Damerau 2009; Species Gateway 2010). Detailed studies on the species’ ecology (e.g. Ingrisch 1984; Berggren et al. 2001; Poniatowski and Fartmann 2005; Holzhauer et al. 2006) and movement behaviour (Berggren et al. 2002; Berggren 2004, 2005) have increased the understanding of how Metrioptera roeselii responds to local biotic and abiotic factors. However it is currently unknown which of the factors are shaping the regional occupancy pattern of Metrioptera roeselii, and to what extent readily-available landscape data can be used to predict the regional distribution of the species

The aim of this study is to model the distribution of Metrioptera roeselii at a large scale (>2000 km2) using species presence-absence data from field surveys and digital landscape data available from the national cartographic agency. Since the predictive ability of occupancy models is known to be scale sensitive (Scott et al. 2002) we model the distribution of Metrioptera roeselii at two different spatial scales (‘landscape’ and ‘local’ scale) and compare model performance. At the ‘landscape’ scale, we measure the landscape composition and structure, factors that affect colonisation and establishment of populations (Werling and Gratton 2008). At the local scale we use land cover type as a predictor for species occurrence as it is thought to reflect closely species habitat requirements (Hirzel and Le Lay 2008).

The questions we sought to answer in this study were: (1) is there any difference in predictive ability of models which use landscape composition and structure versus those that only include local land cover type to explain the distribution pattern of Metrioptera roeselii and (2) which landscape variables explain best the occurrence of Metrioptera roeselii and are these variables consistent between the landscape and local scale?

Material and methods Study species

Metrioptera roeselii (Orthoptera: Tettigoniidae) (Hagenbach 1822) is a small (12–18 mm) predominantly short-winged and flightless bush-cricket commonly found in grasslands of central and northern Europe (Bellmann 2006). In Sweden Metrioptera roeselii occurs mainly in the Lake Mälaren region and the position of the population core area suggests that the species has been introduced via sea cargo (de Jong and Kindvall 1991). There are indications that the expansion of Metrioptera roeselii may cause the displacement of a native orthopteran species (Berggren and Low 2004), but its impact on the insect community as whole is largely unknown. Metrioptera roeselii is an omnivorous generalist that prefers tall grassland habitats. In the agricultural landscape the species is found in extensively grazed pastures, leys, grassy field margins, ditches, and road verges (Marshall and Haes 1988; Berggren et al. 2001). Forests, arable crop fields and intensively grazed pastures are considered to be unsuitable habitat for the species and urban areas are usually avoided (Ingrisch and Köhler 1998; de Jong and Kindvall 1991; Wissmann et al. 2009).

The reproductive season of Metrioptera roeselii in Scandinavia is between July and September. Males stridulate to attract females and the species-specific call makes the species easy to census (Marshall and Haes 1988). Metrioptera roeselii is a wing polymorphic species; extremely favourable weather conditions (mild springs and hot summers) and high population densities trigger the development of long winged morphs (macropters) (Poinatowski and Fartmann 2010). However, in normal years and at range margins the proportion of macropters in Metrioptera roeselii populations rarely exceeds two percent and the vast majority of individuals disperse by walking and jumping (Vickery 1965; Wissmann et al. 2009; pers. obs.).

Data collection

During 2008 and 2009 we surveyed an area of 2554 km2 in the Lake Mälaren region (mid-point 59°44'N, 16°52'E) for the presence of Metrioptera roeselii (Fig. 1). The landscape in this region consists of a mosaic of agricultural land (46%), forest (43%), scattered settlements and small towns (5%), lakes and waterways (3%) and a small proportion of other land use types (3%). In our surveys we sampled the land cover types proportionally to their occurrence in the landscape. We used known locations of Metrioptera roeselii (de Jong and Kindvall 1991; Berggren et al. 2001; Species Gateway 2010) as starting points for our surveys and surveyed the wider surroundings to map the current distribution of the species. We conducted auditory surveys by car (de Jong and Kindvall 1991; Berggren et al. 2001) on sunny days, between 10 am – 5 pm, from mid July until the end of August. Since the species’ call is strong and can be heard over distances of approximately 10 m (Fischer et al. 1997; Bellman 2006), it is possible to listen for stridulating males from the car window while driving slowly (~30 km/ h) along countryside roads (Berggren et al. 2001). We recorded our survey routes and observations of Metrioptera roeselii using a GPS (Garmin 60XL).

Figure 1.

Survey area (mid-point 59°44'N, 16°52'E) in south-central Sweden covering 2554 km2. The 1 × 1 km grid squares (n = 874) show presence (black) and absence (grey) of Metrioptera roeselii.

Variable selection

We used ArcGIS 9.2 (ESRI 2006) to plot and analyse the survey and landscape data. Information on landscape structure and landscape composition was extracted from a topographic map (Geographic Sweden Data (GSD) 1:50 000) and a Swedish CORINE (Coordination of Information on the Environment) land cover map (resolution 30 × 30 m) both available from the Swedish mapping, cadastral and land registration authority. We analysed the effect of landscape variables on the species’ distribution at two spatial scales: the landscape and the local scale. We placed a 1 × 1 km grid across the study area to create presence-absence squares from the species survey data and to design units in which we measured the predictor variables for the landscape scale analysis (Fig. 1). For the analysis at the local scale we use the same 1 × 1 km grid for the species data but extracted the land use data from a 10 m wide buffer strip running parallel to each side of the surveyed roads (i.e. the search area). We compared the models from the search area with the models at the landscape scale to test if we find similar effects of land use on species occurrence at a larger spatial scale.

The distribution of Metrioptera roeselii was treated as presence-absence data within the 1 × 1 km squares for both spatial scales of the analysis (n total = 874 with 318 absence and 556 presence squares). Squares where Metrioptera roeselii was absent were only included in the analysis if they were adjacent to a presence square. Based on our knowledge of the species dispersal behaviour (Berggren et al. 2001, 2002) we excluded distant and isolated absence squares from the analysis because we considered those squares to lie outside the species immediate colonisable area. We chose this conservative approach in order to minimise the number of false absences in the data which otherwise inflates the omission error, lowering the accuracy of the models (Guisan and Thuiller 2005). Because we were primarily interested in modelling the distribution of populations rather than dispersing in individuals, we only included squares in the analysis that contained at least two observations of male Metrioptera roeselii. Previous studies have shown that the species has a good colonising ability and propagules consisting of two males and two females can found sustainable populations (Berggren 2001). Because survey length affects detection probability of the species, we used survey length as a covariate in all models, and only included squares in the analyses in which more than 100 m of road was surveyed.

We used GIS to extract landscape variables that are of ecological relevance for Metrioptera roeselii (Berggren et al. 2001; Berggren et al. 2002; Berggren 2004) and which represent predefined categories in the maps that we used. The land cover categories were generic and consisted of sub-categories of land-use types that resembled each other in terms of vegetation- and management type: (1) arable land (under crop rotation; includes cultivation of cereals, fodder - and root crops, fallow land), (2) forest (includes broadleaved, coniferous and mixed forest, clear-cuts and young plantations), (3) pasture (includes dense herbaceous vegetation dominated by grasses under different grazing regimes), (4) urban areas (includes land with buildings and other man-made structures, small towns and villages), (5) rural settlements (includes solitary houses and farm buildings surrounded by grasslands and gardens), (6) linear elements (combined lengths of streams and roads), and (7) number of fragments of arable land (see Table 1).

Table 1.

Descriptive statistics for the major landscape features and predictor variables used in the regression analyses to explain the distribution of Metrioptera roeselii in south-central Sweden.

Presence squares (1 × 1 km) Absence squares (1 × 1 km)
Variable Min Mean Max SE Min Mean Max SE
Survey length [km] 0.06 1.02 2.62 0.02 0.01 0.69 2.06 0.02
Landscape scale
Arable land [ha] 0.00 53.09 100.00 1.10 0.00 39.96 100.00 1.68
Forest [ha] 0.00 36.91 99.12 1.10 0.00 45.94 100.00 1.75
Pasture [ha] 0.00 4.37 45.12 0.23 0.00 5.19 34.94 0.37
Urban [ha] 0.00 1.67 62.38 0.20 0.00 2.93 59.82 0.45
Rural settlements [ha] 0.00 1.65 17.06 0.10 0.00 0.58 8.88 0.08
Fragments† [count] 0.00 1.43 5.00 0.03 0.00 1.42 6.00 0.05
Linear Elements‡ [km] 0.31 3.45 13.59 0.06 0.20 3.37 10.21 0.09
Stream length [km] 0.00 0.95 3.94 0.04 0.00 0.98 3.36 0.05
Road Length [km] 0.11 2.50 11.91 0.05 0.14 2.39 10.21 0.08
Local scale
Arable land [ha] 0.00 1.21 4.11 0.03 0.00 0.56 3.89 0.04
Forest [ha] 0.00 0.46 3.46 0.03 0.00 0.54 3.76 0.04
Pasture [ha] 0.00 0.14 2.28 0.01 0.00 0.10 1.26 0.01
Rural settlements [m2] 0.00 0.20 3.15 0.02 0.00 0.12 2.20 0.02
Min = Minimum, Max = Maximum, SE is the standard error of the mean.† = Number of fragments of arable land, ‡ = the sum of the length of streams and roads.

We used Pearson’s product-moment correlations to test for the relationships between landscape variables using JMP version 8.0.1 (SAS Institute Inc. 2009). Arable land and forest were highly negatively correlated (r = - 0.86, p < 0.0001), suggesting they are mutually exclusive in the landscape. Thus, we choose to exclude forest and include arable land in the analyses as previous studies have shown that Metrioptera roeselii does not occur in forest areas and arable land under intensive cultivation but occurs and spreads along grassy field margins (Ingrisch and Köhler 1998; Berggren et al. 2001). Linear elements were positively correlated with urban areas (r = 0.56, p < 0.0001) as road length increases with urban development. We excluded urban areas from the analyses since we know from personal observations and records in the national species base (Species Gateway 2010) that Metrioptera roeselii is rarely found in urban areas due to the lack of suitable habitat. All other landscape variables showed low to moderate r-values (r ≤ 0.3) and were included in analyses. Moran’s I values indicated that the response variable was spatially structured which would cause our estimates of variable significance in the models to be exaggerated (Legendre 1993). However, our primary aim was not to elicit precise species-habitat relationships but rather to produce a general applicable model to predict the species distribution over a large spatial extent. We therefore chose a non-spatial modeling approach over explicitly accounting for spatial dependency in the species distribution model.

Statistical analyses

We used logistic regression models to investigate the relationship between the landscape variables and Metrioptera roeselii occurrence at two scales: the landscape scale (1 × 1 km units) and the local scale (10 m area either side of surveyed roads). For both analyses a balanced set of candidate models were considered (i.e. all possible combinations of the variables of interest) and these were ranked according to the relative strength of support for each model using Akaike’s information criterion (AIC). We used AIC weights (ωi) to generate weighted model-averaged parameter estimates when there was no clear best model by including all models within 5 AIC (Σ ωi = 0.95) from the highest-ranked model (Burnham and Anderson 2002). We also estimated the relative importance of the predictor variables by summing the AIC weights over all the models in which the variable was contained (Burnham and Anderson 2002). Parameter estimates and AIC for all models were calculated using the ‘glm’ function in the R 2.8.1 software (R Core Development Team 2008).

We used v-fold cross-validation (Witten and Frank 2000), to evaluate the prediction accuracy of the highest-ranked models from our analyses (i.e. survey scale and landscape scale). Of the 874 survey squares, 80% were randomly sub-sampled as the training set and used to parameterise the model. The coefficients of this model were then used to derive probabilities of occurrence for the remaining 20% of the survey squares. Among the number of data partitioning methods in model evaluation (Fielding and Bell 1997) this ratio of 80% training and 20% test data has been previously found useful (Dormann et al. 2008). The square-specific probabilities were used to calculate a random draw from a Bernoulli probability distribution for each square to produce a prediction (0 or 1) and these were compared to the observed data in the validation set (0 or 1) for each square. Differences in observation versus prediction were then recorded as a proportion of mismatches for the training data set. This was repeated 1000 times, with the proportion of mismatches being modeled as a distribution of errors; i.e. the proportional deviation of the predicted versus the observed – similar to a probability density curve. The median and 95% confidence intervals of these errors were then calculated using the cumulative distribution function (ecdf) in R 2.13.1 (R Development Core Team 2009).

Results

Models at the landscape scale had lower AIC values when compared to equivalent models at the local scale (Table 2), suggesting that variables measured at the landscape scale were better predictors of Metrioptera roeselii presence than those measured in the immediate survey area (local scale). There was strong support for arable land as an important positive predictor for this species, as it was the only variable present in all models with AIC support (Table 2). By comparing different scales in the analyses (landscape versus local) we show that the habitat variables were differently associated with the species presence depending on the spatial scale at which they were measured (Tables 2 and 3).

Table 2.

Model selection results for the effect of landscape variables on the occurrence of Metrioptera roeselii. The model selection statistics are number of parameters (K), Akaike’s information criterion (AIC), difference between model and minimum AIC values (∆AIC), and AIC weights (ωi). Only models with ∆AIC < 10 are shown.

Rank Model       K AIC ∆AIC ωi
Landscape scale
1 Sur + Ara + Rural + Pas + Lin + Frag 7 969.08 0 0.662
2 Sur + Ara + Rural + Pas + Lin 6 971.39 2.31 0.209
3 Sur + Ara + Rural + Pas + Frag 6 973.62 4.54 0.068
4 Sur + Ara + Rural + Pas 5 974.65 5.57 0.041
5 Sur + Ara + Rural + Lin + Frag 6 977.51 8.43 0.010
6 Sur + Ara + Rural + Lin 5 977.86 8.78 0.008
Local scale
1 Sur + Ara + Rural 4 983.48 0 0.581
2 Sur + Ara + Rural + Pas 5 984.65 1.17 0.324
3 Sur + Ara 3 987.82 4.34 0.066
4 Sur + Ara + Pas     4 989.49 6.01 0.029
Abbreviations used for the explanatory variables in the models: Sur = Survey length, Ara = Arable land, Rural = Rural settlements, Pas = Pasture, Lin = Linear Elements, Frag = Number of fragments of arable land.

At the landscape scale, Metrioptera roeselii presence was best explained by the full model, containing arable land, rural settlements, pasture, number of arable land fragments and linear elements (Table 2). The second- and third-ranked models differed in either number of fragments or linear elements, suggesting that structural landscape variables had weaker support in explaining Metrioptera roeselii occurrence. Contrary to expectation, occurrence of Metrioptera roeselii was negatively correlated with the amount of pasture and linear elements, and positively correlated with the number of fragments of arable land (Table 3). The three land-use variables (arable land, rural settlements and pasture) had the highest relative-importance weights (1.0), followed by linear elements (0.927) and number of fragments (0.778).

Table 3.

AIC-weighted model-averaged parameter estimates generated from the top three models (Σ ωi = 0.95) presented in Table 2.

Averaged Model Variable Coeff SE
Landscape Scale
(1 × 1 km) (Intercept) -1.478 0.282
Survey length 1.791 0.199
Arable land 0.014 0.003
Rural settlements 0.259 0.051
Pasture -0.042 0.013
Linear Elements‡ -0.124 0.050
Fragments† 0.159 0.078
Local Scale
(10 m buffer) (Intercept) -0.996 0.169
Survey length 0.643 0.228
Arable land 1.071 0.142
  Rural settlements 0.559 0.228
Pasture 0.100 0.111
† = Number of fragments of arable land, ‡ = the sum of the length of streams and roads.Coeff = coefficient; SE = Standard error of coefficient.

At the local scale, the two highest-ranked models contained arable land and rural settlements. This, in combination with their relative-importance weights (1.0 and 0.905 respectively), demonstrates the strong support for them as positive predictors (Tables 2 & 3). Although pasture was included in the second-highest-ranked model, an examination of Table 2 shows that its inclusion in models generally results in a lower ranking than models without it – suggesting very weak support for it as a predictor of Metrioptera roeselii presence (relative-importance of pasture = 0.353).

Cross-validation showed that the models were generally accurate in their predictions of species occurrence across the spatial scales for the environmental gradients examined in the study. The landscape-level model prediction for the probability of Metrioptera roeselii being detected in a square had an error which ranged from -0.091 to +0.080 (95% CI; Fig. 2a). At the survey scale, model prediction error for the probability of detection ranged between -0.075 to +0.097 (95% CI; Fig. 2b).

Figure 2.

Cross-validation accuracy of 1000 models using randomly selected training and validation sets (80% and 20% respectively). The curves show the relative deviation of prediction accuracy when comparing estimated to observed occurrence of Metrioptera roeselii being detected in a square at a the survey scale (vertical bars show the 95% CI for model prediction error: -0.075 to +0.097), and b the landscape scale (-0.091 to +0.080).

Discussion

In our study, the distribution of Metrioptera roeselii was best explained by models at the landscape scale. This indicates that measuring the landscape characteristics within 1 × 1 km units captures both the availability of habitat for the species and incorporates ecological functions of the landscape features (Crawford and Hoagland 2010). The weaker relationship between land use and species occurrence at the local scale could be attributable to the coarse grain size of the land-cover data failing to capture local aspects of habitat quality, i.e. vegetation heterogeneity, microclimate (Gardiner and Dover 2008) and its temporal variability (Gardiner et al. 2008; Poniatowski and Fartmann 2008) as well as important biotic interactions (Huston 2002) that are influencing the distribution of the species. Our study shows that landscape data extracted from digital map sources can be used to explain the regional distribution pattern of this expanding species. Determining biologically important variables and the optimal spatial scale is a prerequisite to predict the likelihood of occurrence of a species in non-surveyed sites with a resolution of 1 km 2 and form the base for monitoring species spread, serving conservation planning and future research on spatial processes shaping species distributions. The models can also be further developed and used for region-wide predictions in areas similar to the study area, assisting in devising management actions and possible control of undesired species expansion (Hutto and Young 2002; Scott et al. 2002). However, extrapolation of model results should be treated with caution. Abiotic factors such as land cover can generally be applied only within a limited spatial extent and time frame because the same variables can differ in habitat suitability moreover the same species may respond to different sets of variables in different parts of its distributional range (Guisan and Zimmerman 2000).

When modelling species distributions in fragmented landscapes it is important to incorporate the landscape structure into the analyses (Umetsu et al. 2008). The number of fragments of arable land was a positive predictor for the occurrence of Metrioptera roeselii, indicating that the field margins offer important edge habitat and serve as dispersal paths in the agricultural landscape (Berggren et al. 2001). Similar dispersal behaviour has been observed in the wood cricket Nemobius sylvestris that moves along habitat edges (Brouwers et al. 2011). Contrary to expectations, linear landscape elements (roads and streams) had a negative effect on species occurrence at the landscape scale. One possible explanation is that although linear elements have been associated with increased dispersal opportunities in small-scaled studies, at larger scales linear landscape features such as major roads and streams act as a barrier for the species dispersal if they separate suitable habitat areas (de Jong and Kindvall 1991). Due to the large spatial extent of our study it was not possible for us to explicitly incorporate spatial configuration and orientation of linear landscape features in the model.

At both spatial scales that we analysed, arable land and rural settlements turned out to be strong predictors for the presence of Metrioptera roeselii suggesting that these land use types can be used as surrogate measure for grassland habitat in the region. The positive effect of arable land on the occurrence of Metrioptera roeselii might be surprising at first since it is known that Metrioptera roeselii avoids crop fields because of the lack of shelter, food and egg laying places (Ingrisch and Köhler 1998). However, arable land is a generic land use description and vegetation cover varies with the type of crop cultivated. In Sweden, crop rotation is commonly practised (Söderberg 2006) and arable land becomes temporally a suitable habitat for orthopterans and other grassland living insects when crop fields are shifted into fallows or leys (Duelli et al. 1999).The ability to track resources is particularly important for species in dynamic landscapes. In areas with intensive agricultural production the grassy field margins and hedgerows often have high species richness and function as dispersal corridors and source habitats for colonisers of crop fields (Marshall and Moonen 2002; Meek et al. 2002). The present findings support our assumption that grassland insects like Metrioptera roeselii benefit from habitat heterogeneity in arable landscapes. Braschler et al. (2009) found that cricket (Ensifera) density was higher in fragmented plots, as uncut patches of grassy vegetation play an important role in maintaining insect diversity in the agricultural landscape by offering shelter from predators and serving as mating and egg laying sites. A previous study by Bieringer and Zulka (2003) showed that orthopteran species richness increases with distance to forest edge. We believe that the positive effect of arable land in our study was not simply because the bush-crickets avoided forest, but rather that agricultural areas contain a larger amount of suitable grassland vegetation than forests.

In cultivated landscapes, generalist species that are able to occupy a broad range of habitat types are less sensitive to local habitat loss (Marini et al. 2008, 2009a). Metrioptera roeselii is an example of a grassland generalist (Ingrisch and Köhler 1998) colonising a range of grassland types (Gardiner et al. 2008; Poniatowski and Fartmann 2005). Like its close relative Metrioptera bicolor (Kindvall 1996) it is able to sustain populations in small patches of habitat. Rural settlements, despite covering only a small area of the landscape, have been shown to provide important habitat for a range of species (Belfrage et al. 2005; Rosin et al. 2011) and may function as source patches for Metrioptera roeselii enabling the species to colonise surrounding areas. Extensive farming practices and small field sizes are positively correlated with habitat heterogeneity, which in turn has a positive effect on the local diversity of species with limited movement ability like pollinators and grassland living insects (Benton et al. 2003; Marini et al. 2009b; Steck et al. 2007).

We expected that the amount of pasture and the presence of Metrioptera roeselii would be positively correlated on both scales since Metrioptera roeselii has been found to colonise extensively grazed pastures (Poniatowski and Fartmann 2005). The negative correlation of pasture on Metrioptera roeselii occurrence at the landscape scale is difficult to interpret. A possible explanation could be that the overall proportion of pastures in the landscape is small and its distribution scattered which makes it more difficult for the species to colonize.

Species ecology, range size and rarity have an influence on model performance (Franklin et al. 2009; Syphard and Franklin 2009). Results from other studies (Heikkinen et al. 2006; Segurado and Araújo 2004) have shown that specialist species and species with a limited range are generally more accurately modeled than generalist species and species with a wide geographic range, Metrioptera roeselii is an example of the latter. The natural dynamics of the study species makes it more difficult to model its distribution because the assumption of the species being in equilibrium with the environment is violated and dispersal contributes to spatial autocorrelation in the data (Franklin et al. 2009). With these limitations in mind we thoroughly surveyed the range of environmental conditions present in the distribution area from the core of the study area to the margin aiming to obtain a large sample size as possible. Despite our surveys were conducted by car, we sampled all important habitat types (arable land, forests, pastures and human settlements) proportionally to their occurrence in the landscape. Aware of the trade-off between model generality, reality and precision (Guisan and Zimmerman 2000), we prioritized the former as our primary aim in this study was to develop predictive model for Metrioptera roeselii within the study region. The model can be further developed and applied to other grassland insects with similar traits.

Conclusions

Type of land use and structural landscape elements describing the amount of available habitat are important predictors for species occurrences (Hein et al. 2007; Kemp et al. 1990; Crawford and Hoagland 2010). The possibility to model Metrioptera roeselii distribution using survey data and available land-cover data on a scale that is easy to extract and utilise for managers is promising in that it will enable us to predict the direction and possible extent of future range expansion of the species. As many Orthopterans disperse and interact with the environment in a similar way (Hjermann and Ims 1996; Diekötter et al. 2007; Brouwers et al. 2011), the results from this study may also be valid for other related species that are now expanding their distribution areas. This is very useful, as many studies on grassland living insects face a similar dilemma: a limited availability of distribution data for species that are living in highly dynamic landscapes (Marini et al. 2009b). The possibility to utilise available distribution data in combination with land-cover data enables us to improve our understanding of the species ecology, to highlight areas of conservation concern and to predict species occurrences in a time of environmental change (Bonter et al. 2010).

Acknowledgements

We are grateful to Emilia Broberg and Lina Ahlbäck for their assistance in the field and with the GIS work. We thank Matthew Low for valuable comments on the manuscript. The research was funded by a research grant from FORMAS and The Magnus Bergvalls foundation to ÅB.

References
Bazazi S, Romanczuk P, Thomas S, Schimansky-Geier L, Hale JJ, Miler GA., Sword GA, Simpson SJ, Couzin ID (2011) Nutritional state and collective motion: from individuals to mass migration. Proceedings of the Royal Society of London B 27: 356-363. doi: 10.1098/rspb.2010.1447
Belfrage K, Björklund J, Salomonsson L (2005) The effects of farm size and organic farming on diversity of birds, pollinators, and plants in a Swedish landscape. Ambio 34: 582-588.
Bellmann H (2006) Der Kosmos Heuschreckenführer - Die Arten Mitteleuropas sicher bestimmen. Franckh-Kosmos Verlags-Gmbh & Co.KG, Stuttgart.
Benton TG, Vickery JA, Wilson JD (2003) Farmland biodiversity: is habitat heterogeneity the key? Trends in Ecology and Evolution 18: 182-188. doi: 10.1016/S0169-5347(03)00011-9
Berggren Å (2001) Colonization success in Roesel’s bush-cricket Metrioptera roeseli: the effects of propagule size. Ecology 82: 274-280.
Berggren Å (2004) Impact of grazing on individual male movement in Roesel’s bush-cricket Metrioptera roeseli: One possible clue to species range expansion. Journal of Insect Behavior 17: 419-429. doi: 10.1023/B:JOIR.0000042531.27859.ac
Berggren Å (2005) The effect of conspecifics on individual male movement in Roesel’s bush cricket, Metrioptera roeseli. Ecological Entomology 30: 480-483. doi: 10.1111/j.0307-6946.2005.00709.x
Berggren Å, Low M (2004) Exclusion of the native bog bush-cricket Metrioptera brachyptera by the currently invading Roesel’s bush-cricket Metrioptera roeseli. Entomologisk Tidskrift 125: 125-32.
Berggren Å, Birath B, Kindvall O (2002) Effect of corridors and habitat edges on dispersal behavior, movement rates, and movement angles in Roesel’s bush-cricket (Metrioptera roeseli). Conservation Biology 16: 1562-1569. doi: 10.1046/j.1523-1739.2002.01203.x
Berggren Å, Carlson A, Kindvall O (2001) The effect of landscape composition on colonization success, growth rate and dispersal in introduced bush-crickets Metrioptera roeseli. Journal of Animal Ecology 70: 663-670. doi: 10.1046/j.1365-2656.2001.00525.x
Bieringer G, Zulka KP (2003). Shading out species richness: edge effect of a pine plantation on the Orthoptera (Tettigoniidae and Acrididae) assemblage of an adjacent dry grassland. Biodiversity and Conservation 12 (7): 1481-1495. doi: 10.1023/A:1023633911828
Bonter DN, Zuckerberg B, Dickinson JL (2010) Invasive birds in a novel landscape: habitat associations and effects on established species. Ecography 33: 494-502.
Brouwers NC, Newton AC, Bailey S (2011) The dispersal ability of wood cricket (Nemobius sylvestris) (Orthoptera: Gryllidae) in a wooded landscape. European Journal of Entomology 108: 117-125.
Burnham KP, Anderson DR (2002) Model selection and multimodel inference: a practical information-theoretic approach. Springer-Verlag, New York, USA.
Cote IM, Reynolds JD (2002) Predictive ecology to the rescue? Science 298: 1181-1182. doi: 10.1126/science.1079074
Crawford PHC, Hoagland BW (2010) Using species distribution models to guide conservation at the state level: the endangered American burying beetle (Nicrophorus americanus) in Oklahoma. Journal of Insect Conservation 14: 511-521. doi: 10.1007/s10841-010-9280-8
De Groot M, Rebeusek F, Grobelnik V, Govedic M, Salamun A, Verovnik R (2009) Distribution modelling as an approach to the conservation of a threatened alpine endemic butterfly (Lepidoptera: Satyridae). European Journal of Entomology 106: 77-84.
de Jong J, Kindvall O (1991) Cikadavårtbitaren Metrioptera roeseli– nykomling eller hotad relikt? [The Roesel’s bush-cricket Metrioptera roeseli - new in Sweden or a threatened relict species?] Fauna och Flora 86: 215-221.
Diekötter T, Speelmans M, Dusoulier F, Van Wingerden WKRE, Malfait JP, Crist TO, Edwards PJ., Dietz H (2007) Effects of landscape structure on movement patterns of the flightless bush cricket Pholidoptera griseoaptera. Environmental Entomology 36: 90-98. doi: 10.1603/0046-225X(2007)36[90:EOLSOM]2.0.CO;2
Dormann CF, Purschke O, Marquez JR, Lautenbach G, Schroder SB (2008) Components of uncertainty in species distribution analysis: a case study of the great grey shrike. Ecology 89 (12): 3371-3386. doi: 10.1890/07-1772.1
Duelli P, Obrist MK, Schmatz DR (1999) Biodiversity evaluation in agricultural landscapes: above-ground insects. Agriculture, Ecosystems & Environment 74: 33-64. doi: 10.1016/S0167-8809(99)00029-8
Early R, Anderson B, Thomas CD (2008) Using habitat distribution models to evaluate large-scale landscape priorities for spatially dynamic species. Journal of Applied Ecology 45: 228-238. doi: 10.1111/j.1365-2664.2007.01424.x
Fielding AH, Bell JF (1997). A review of methods for the assessment of prediction errors in conservation presence/absence models. Environmental Conservation 24 (01): 38-49. doi: 10.1017/S0376892997000088
Franklin J, Wejnert KE, Hathaway SA, Rochester CJ, Fisher RN (2009). Effect of species rarity on the accuracy of species distribution models for reptiles and amphibians in southern California. Diversity and Distributions 15 (1): 167-177. doi: 10.1111/j.1472-4642.2008.00536.x
Fischer FP, Schulz U, Schubert H, Knapp P, Schmoger M (1997) Quantitative assessment of grassland quality: acoustic determination of population sizes of orthopteran indicator species. Ecological Applications 7: 909-920. doi: 10.1890/1051-0761(1997)007[0909:QAOGQA]2.0.CO;2
Gardiner T (2009) Macropterism of Roesel’s Bushcricket Metrioptera roeselii in Relation to Climate Change and Landscape Structure in Eastern England. Journal of Orthoptera Research 18: 95-102. doi: 10.1665/034.018.0110
Gardiner T, Dover J (2008) Is microclimate important for Orthoptera in open landscapes? Journal of Insect Conservation 12: 705-709. doi: 10.1007/s10841-007-9104-7
Gardiner T, Hill J, Marshall EJP (2008) Grass field margins and Orthoptera in eastern England. Entomologist’s Gazette 59: 251-257.
Guisan A, Thuiller W (2005) Predicting species distribution: offering more than simple habitat models. Ecological Letters 8: 993-1009. doi: 10.1111/j.1461-0248.2005.00792.x
Guisan A, Zimmermann NE (2000). Predictive habitat distribution models in ecology. Ecological Modelling 135(2–3): 147-186. doi: 10.1016/S0304-3800(00)00354-9
Gwynne DT (2001) Katydids and bush-crickets: reproductive behaviour and evolution of the Tettigoniidae. Cornell University Press, USA.
Heikkinen RK, Luoto M, Araujo MB, Virkkala R, Thuiller W, Sykes MT (2006). Methods and uncertainties in bioclimatic envelope modelling under climate change. Progress in Physical Geography 30 (6): 751-777. doi: 10.1177/0309133306071957
Hein S, Gombert J, Hovestadt T, Poethke H.J (2003) Movement patterns of the bush cricket Platycleis albopunctata in different types of habitat: matrix is not always matrix. Ecological Entomology 28: 432-438. doi: 10.1046/j.1365-2311.2003.00531.x
Hein S, Binzenhofer B, Poethke HJ, Biedermann R, Settele J, Schroder B (2007) The generality of habitat suitability models: A practical test with two insect groups. Basic and Applied Ecology 8: 310-320. doi: 10.1016/j.baae.2006.09.012
Hirzel AH, Le Lay G (2008) Habitat suitability modelling and niche theory. Journal of Applied Ecology 45: 1372-1381. doi: 10.1111/j.1365-2664.2008.01524.x
Hjermann DO, Ims RA (1996) Landscape ecology of the wart-biter Decticus verrucivorus in a patchy landscape. Journal of Animal Ecology 65: 768-780. doi: 10.2307/5675
Hochkirch A, Damerau M (2009) Rapid range expansion of a wing-dimorphic bush-cricket after the 2003 climatic anomaly. Biological Journal of the Linnean Society 97: 118-127. doi: 10.1111/j.1095-8312.2008.01199.x
Holzhauer S, Ekschmitt K, Sander A-C, Dauber J, Wolters V (2006) Effect of historic landscape change on the genetic structure of the bush-cricket Metrioptera roeseli. Landscape Ecology 21: 891-899. doi: 10.1007/s10980-005-0438-9
Huston MA (2002) Introductory essay: critical issues for improving predictions. In: Scott JM, Heglund PJ, Morrison ML, Haufler JB, Raphael MG, Wall WA, Samson FB (Eds) Predicting Species Occurrences: Issues of Accuracy and Scale. Island Press, Covelo, CA, pp. 7–21.
Hutto RL, Young JS (2002) Regional land bird monitoring: perspectives from the Northern Rocky Mountains. Wildlife Society B. 30: 738-750.
Ingrisch S (1984) The influence of environmental factors on dormancy and duration of egg development in Metrioptera roeseli (Orthoptera: Tettigoniidae). Oecologia 61: 254-258. doi: 10.1007/BF00396769
Ingrisch S, Köhler G (1998) Die Heuschrecken Mitteleuropas. Westarp-Wiss., Magdeburg.
Jimenez-Valverde A, Gomez JF, Lobo JM, Baselga A, Hortal J (2008) Challenging species distribution models: the case of Maculinea nausithous in the Iberian Peninsula. Annales Zoologici Fennici 45: 200-210.
Kemp WP, Harvey SJ, O‘Neill KM (1990) Patterns of vegetation and grasshopper community composition. Oecologia 83: 299-308.
Kindvall O (1996) Habitat heterogeneity and survival in a bush cricket metapopulation. Ecology 77: 207-214. doi: 10.2307/2265670
Legendre P (1993) Spatial Autocorrelation - Trouble or New Paradigm. Ecology 74: 1659-1673. doi: 10.2307/1939924
Marini L, Fontana P, Scotton M, Klimek S (2008) Vascular plant and Orthoptera diversity in relation to grassland management and landscape composition in the European Alps. Journal of Applied Ecology 45: 361-370. doi: 10.1111/j.1365-2664.2007.01402.x
Marini L, Fontana P, Klimek S, Battisti A, Gaston KJ (2009a) Impact of farm size and topography on plant and insect diversity of managed grasslands in the Alps. Biological Conservation 142: 394-403. doi: 10.1016/j.biocon.2008.10.034
Marini L, Fontana P, Battisti A, Gaston KJ (2009b) Agricultural management, vegetation traits and landscape drive orthopteran and butterfly diversity in a grassland-forest mosaic: a multi-scale approach. Insect Conservation and Diversity 2: 213-220. doi: 10.1111/j.1752-4598.2009.00053.x
Marshall EJR, Moonen AC (2002) Field margins in northern Europe: their functions and interactions with agriculture. Agriculture, Ecosystems & Environment 89: 5-21. doi: 10.1016/S0167-8809(01)00315-2
Marshall JA, Haes ECM (1988) Grasshoppers and Allied Insects of Great Britain and Ireland. Harley Books, Martins, Great Horkesley, Colchester, Essex, UK.
Meek B, Loxton D, Sparks T, Pywell R, Pickett H, Nowakowski M (2002) The effect of arable field margin composition on invertebrate biodiversity. Biological Conservation 106: 259-271. doi: 10.1016/S0006-3207(01)00252-X
Morueta-Holme N, Flojgaard C, Svenning JC (2010) Climate Change Risks and Conservation Implications for a Threatened Small-Range Mammal Species. Plos One 5.
Poniatowski D, Fartmann T (2005) Die Ökologie von Roesels Beißschrecke (Metrioptera roeselii) im Feuchtgrünland der Medebacher Bucht (Südwestfalen). Articulata 20: 85-111.
Poniatowski D, Fartmann T (2008) The classification of insect communities: Lessons from orthopteran assemblages of semi-dry calcareous grasslands in central Germany. European Journal of Entomology 105: 659-671.
Poniatowski D, Fartmann T (in press). Weather-driven changes in population density determine wing dimorphism in a bush-cricket species. Agriculture, Ecosystems and Environment. doi: 10.1016/j.agee.2010.10.006
R Development Core Team (2008) & 2009 R: a language and environment for statistical computing. - R Foundation for Statistical Computing. (http://www.R-project.org)
Rosin ZM, Skórka P, Lenda M, Moron D, Sparks TH, Tryjanowski P (2011) Increasing patch area, proximity of human settlement and larval food plants positively affect the occurrence and local population size of the habitat specialist butterfly Polyommatus coridon (Lepidoptera: Lycaenidae) in fragmented calcareous grasslands. European Journal of Entomology 108: 99-106.
SAS Institute Inc. (2009) JMP version 8.0.1. SAS Institute, Cary.
Scott MJ, Heglund PJ, Morrison ML, Haufler JB, Raphael MG, Wall WA, Samson FB (2002) Predicting species occurrences: issues of accuracy and scale. Island press, Washington, USA
Segurado P, Araujo MB (2004). An evaluation of methods for modelling species distributions. Journal of Biogeography 31 (10): 1555-1568. doi: 10.1111/j.1365-2699.2004.01076.x
Simmons AD, Thomas CD (2004) Changes in dispersal during species’ range expansions. American Naturalist 164: 378-395. doi: 10.1086/423430
Species Gateway (2010) http://www.artportalen.se/
Steck CE, Burgi M, Coch T, Duelli P (2007). Hotspots and richness pattern of grasshopper species in cultural landscapes. Biodiversity and Conservation 16 (7): 2075-2086. doi: 10.1007/s10531-006-9089-7
Sword GA, Lorch PD, Gwynne DT (2008) Radiotelemetric analysis of the effects of prevailing wind direction on Mormon cricket migratory band movement. Environmental Entomology 37: 889-896. doi: 10.1603/0046-225X(2008)37[889:RAOTEO]2.0.CO;2
Söderberg T (2006) Miljöeffekter av träda och olika växtföljder – rapport från projektets CAP:s mijlöeffekter. Jordbruksverkets rapport 2006: 4, Jönköping, Sweden.
Syphard AD, Franklin J (2009). Differences in spatial predictions among species distribution modeling methods vary with species traits and environmental predictors. Ecography 32 (6): 907-918. doi: 10.1111/j.1600-0587.2009.05883.x
Umetsu F, Metzger JP, Pardini R (2008) Importance of estimating matrix quality for modeling species distribution in complex tropical landscapes: a test with Atlantic forest small mammals. Ecography 31: 359-370. doi: 10.1111/j.0906-7590.2008.05302.x
Vickery VR (1965) Factors governing the distribution and dispersal of the recently introduced grasshopper, Metrioptera roeseli (Hgb.) (Orthoptera: Ensifera). Annales de la Societe Entomologique du Quebec 10: 165-171.
Werling BP, Gratton C (2008) Influence of field margins and landscape context on ground beetle diversity in Wisconsin (USA) potato fields. Agriculture, Ecosystems & Environment 128: 104-108. doi: 10.1016/j.agee.2008.05.007
Wissmann J, Schielzeth H, Fartmann T (2009) Landscape-scale Expansion of Roesel’s bush-cricket Metrioptera roeselii at the North-western Range Limit in Central Europe (Orthoptera: Tettigoniidae). Entomologia Generalis 31: 317-326.
Witten IH, Frank E (2000) Data Mining: Practical Machine Learning Tools and Techniques with Java Implementations. Academic Press: San Diego.