Butterfly monitoring using systematically placed transects in contrasting climatic regions – exploring an established spatial design for sampling

Butterfly monitoring schemes are recording programs initiated to monitor nationwide butterfly abundance and distribution patterns, often with help from volunteers. The method generates high-resolution data, but may be associated with a degree of habitat sampling bias if volunteers prefer to survey areas perceived to be high-quality butterfly habitats. This can result in habitats becoming underrepresented in the data set, leading to less information about the butterfly populations there. In the present study, we investigate the possibility of applying a spatial design used by the Swedish Bird Survey for nationwide, gridbased sampling, with a goal to get butterfly monitoring data covering a representative sample of different habitats. We surveyed four 2×2 km sampling squares, split into 100 m segments, in the southernmost region of Sweden (Scania) and four in the northernmost region (Norrbotten). The grid-based transects were compared with volunteer-selected transects in a GIS analysis using a refined Swedish version of CORINE land cover data to see how well these two transect designs represent true habitat coverage. A total of 53 km transect was monitored, resulting in 490 individuals and 29 different species recorded. We found that transect cover correlated significantly with overall land cover using both monitoring methods, though standardised transects outperformed volunteer-selected transects in habitat representation in Scania, but not in Norrbotten. Butterflies were found to aggregate significantly in specific habitats, but with contrasting results for the two geographically different regions. Grasslands in both regions generated a high number of recorded butterflies, although so did clear-cut and residential areas in Norrbotten as well. The Nature Conservation 14: 41–62 (2016) doi: 10.3897/natureconservation.14.7497 http://natureconservation.pensoft.net Copyright Elin Videvall et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. RESEARCH ARTICLE Launched to accelerate biodiversity conservation A peer-reviewed open-access journal Elin Videvall et al. / Nature Conservation 14: 41–62 (2016) 42 highest number of individuals recorded per transect was found in bogs in Scania. This study emphasises the value of complementing free site selection monitoring schemes with spatially representative schemes such as the Swedish Bird Survey, and sheds some light on general habitat preferences for Swedish butterflies in two contrasting climatic regions.


Introduction
Butterflies are the most widely studied of all insect groups (Dennis et al. 2006, van Swaay et al. 2008), and their sensitivity to environmental change together with the availability of butterfly data makes this group very useful as indicators for biodiversity (New 1997, Thomas 2005, Dennis et al. 2006, van Swaay et al. 2015).Furthermore, since many butterfly species require a warm microclimate for optimal growth and development (Wallisdevries andVan Swaay 2006, Eilers et al. 2013), they can also serve as indicators for climatic change (Thomas 2005, van Swaay et al. 2008, Betzholtz et al. 2013).While birds, another popular group of biodiversity indicators, range more widely and over larger areas, butterflies provide important additional and complementary area-specific information since they are more likely to reflect environmental changes occurring on a more detailed scale (van Swaay et al. 2006).
Habitat loss and fragmentation is a major driver behind the decline of many butterfly species worldwide (Bergman et al. 2004, Ekroos et al. 2010).Due to changes in land-use and intensification of agriculture throughout the last century there has been a loss of many open and half-open natural and semi-natural habitats (Nilsson et al. 2013, Cousins et al. 2015).This has had dramatic effects on many insect groups, and there has been a decline in butterfly numbers all over Europe (van Swaay et al. 2006, Konvicka et al. 2008, Van Dyck et al. 2009, Dover et al. 2011).The recently developed "Euro pean Grassland Butterfly Indicator" (van Swaay et al. 2015) based on butterfly data from 22 European countries suggests that grassland butterfly populations may have declined as much as 30% since 1990.
In order to effectively monitor population trends, butterfly monitoring schemes now run in several European countries, e.g. the UK, Finland, Germany and the Netherlands (van Swaay et al. 2008, 2015, Kühn et al. 2012, Botham et al. 2013, Heliölä et al. 2013).Schemes vary in the way monitoring sites are selected (van Swaay et al. 2015), with some schemes using random site placement whereas others have sites selected by coordinators.However, the majority of schemes use the more flexible approach where butterfly recorders can place their monitoring sites freely (van Swaay andWarren 2012, van Swaay et al. 2015).Free site selection is generally appreciated by recorders as they can be involved in the site selection process by influencing choice of site characteristics, accessibility, and being able to relate more closely to the sites that they monitor (van Swaay and Warren 2012).The main disadvantage of free site selection is that the geographic coverage commonly becomes non-random and hence not representative of habitats and butterfly populations in general (Dover et al. 1997, van Swaay and Warren 2012, van Swaay et al. 2015).The concern that open areas and in particular semi-natural grasslands tend to become overrepresented in monitoring data and that forests and other parts of the wider countryside are less well covered (van Swaay and Warren 2012) have recently inspired strategies that substantially reduce sampling bias (e.g, Brereton et al. 2011, Lang andBuhler 2012).
In the UK, species that are widespread across the general countryside compromise half of the butterfly fauna (Asher et al. 2001) but are underrepresented in traditional monitoring; declines for these species went largely undetected by monitoring schemes during the 20 th century (Brereton et al. 2011).This led to the launch of the Wider Countryside Butterfly Survey, a scheme which uses the grid-based sampling design of the British Bird Survey (Greenwood et al. 1995) to get representative trends across the whole countryside (Brereton et al. 2011, Botham et al. 2013, Roy et al. 2015).Grid-based, spatially representative sampling is increasingly used in bird monitoring today (e.g, Greenwood et al. 1995, Davey et al. 2013, Lehikoinen 2013) and offers an attractive development for butterfly monitoring (Kéry and Plattner 2007, Brereton et al. 2011, Lang and Buhler 2012, van Swaay and Warren 2012).However, there is a great need to evaluate how well such standardised designs suit butterfly monitoring in different climatic regions.Factors such as the openness of the countryside, topography, and the time and effort necessary to walk transects in different regions can differ markedly even within one single country (Brereton et al. 2011) and is particularly relevant in countries that span a large range of latitudes.For instance, a coniferous forest in Sweden's southern, continental region can be dense, dark and not particularly well suited for butterflies.In contrast, corresponding forests in the northern, boreal region can be much more open, sunlit, and hence attractive to butterflies.
The Swedish Butterfly Monitoring Scheme is a nationwide program with free site selection (Pettersson et al. 2011).Volunteers appreciate that they can place monitoring transects themselves and this flexible design has been central to the growth of the scheme.Because free site selection tends to result in some habitats being underrepresented (van Swaay and Warren 2012), it would be valuable to complement the Swedish scheme with a grid-based approach (cf.Lindenmayer and Likens 2009).A very promising way of doing so would be to count butterflies along transects that already form part of a standardised, nationwide monitoring design: the Swedish Bird Survey (Green and Lindström 2015).This scheme consists of 716 routes in a 25×25 km grid, covering Sweden as a whole.The layout ensures that all major habitats are proportionally represented.At the centre of each grid cell, the bird fauna is censused once per year along eight 1 km transects arranged in a 2×2 km square (Figure 1, Green and Lindström 2015).The standard method of butterfly monitoring, using "Pollard walk" transects, (Pollard and Yates 1993) could potentially be used along the bird monitoring transects (Brereton et al. 2011, Roy et al. 2015).Hence, the grid used by the Swedish Bird Survey offers an attractive design for obtaining butterfly monitoring data that cover major habitats proportionately.
With this study we have a twofold goal.The first is to evaluate if the large-scale, grid-based method used by the Swedish Bird Survey can be adapted and applied to butterfly monitoring.The second aim is to use the collected data to quantify butterfly abundance in relation to habitat characteristics in two contrasting climatic regions of Sweden: Scania and Norrbotten.

GIS analysis of land cover
While transects based on standardised grid-based designs are generally assumed to outperform free site selection in terms of true habitat representation, this assumption is rarely evaluated.Here, we quantified land cover in the Swedish Bird Survey transects visited in the present experiment (N = 4 in each region, total length: 32 km per region) as well as land cover in volunteer-selected free transects in the same regions (N = 5 in each region, total length 12.92 km in Scania and 12.22 km in Norrbotten).The free transects that were analysed comprised all sites in Norrbotten and a random, corresponding subsample of the transects in Scania.The free transects that we analysed covered less distance that the standardised transects and hence had smaller areas.To allow direct comparison despite differing total areas, we recalculated land cover to the smallest area (free transects, Scania: 4.45 km 2 ; Appendix: Table S3).Land cover data was obtained from the SMD database, which is a refined, Swedish version of the CORINE land cover database with the smallest mapping unit 1-25 ha, map resolution 25×25 m, and 2000 as reference year (Swedish Environmental Protection Agency, 2014).To quantify land cover class coverage, we clipped the SMD using a 200 m buffer around each transect following Davey et al. (2013) using QGIS (v.2.12, QGIS Development Team 2016).

Field data collection
Butterfly monitoring took place during the summer of 2010 from June 17 to July 7 in Scania, the southernmost part of Sweden, and from July 20-28 in Norrbotten, the northernmost part of Sweden.Each site consisted of the four outer sections of 2×2 km squares used by the Swedish Bird Survey (Figure 1), equalling an 8 km long transect.The survey squares were selected with the intention to cover different habitats along the transects, in order to thoroughly test the sampling method.The four sites visited in Scania were called Hyby, Tjörnarp, Slätteberga and Kongaö, and the four sites in Norrbotten were Sundom, Rosfors, Långberget and Bergnäset (Appendix: Table S1, Figure 1).Two persons (EV together with a colleague) walked all transects.Recording only took place during sunny weather, and monitoring did not start earlier than 10 am and ended not later than 5 pm (Central European summer time, UTC +2), which is the time most suitable for butterfly activity (van Swaay et al. 2008).Average tempera- ture throughout the day ranged from 17 °C to 25 °C, and the wind varied in different habitats, but never exceeded 5 in the Beaufort scale, which is acceptable in terms of butterfly monitoring (van Swaay et al. 2008).
All transects monitored were divided into segments of approximate 100 m length with an accompanying description for the habitat surrounding the transects.The distance of the transect segments were estimated with the help of maps.The different habitat categories were: deciduous forest, coniferous forest, grassland, residential area, fen/ bog, and clear-cut area (examples of three habitats can be seen in Figure 2).Road was added as an additional unique category since much monitoring had to be performed on small paths and roads due to accessibility.The habitat description of each 100 m segment consisted of percentages of each habitat category that the transect crossed (e.g.deciduous forest 80%, grassland 20%), in order to reflect the surroundings as closely as possible.As each butterfly individual was attributed to a specific 100 m segment, this typically resulted in non-integer numbers for butterfly individual and species counts per habitat category (i.e.0.8 butterflies in deciduous forest and 0.2 butterflies in grassland, based on the example above).Areas such as cultivated fields, highly dense forests and open water were avoided for reasons of safety, land-owner privacy, and accessibility.
Recordings were made of all butterflies (Rhopalocera) and burnet moths (Zygaenidae), as seen within an 'invisible box' of 5 m in front of the recorders, 2.5 m to each side and 5 m above, according to the 'Pollard walk' method (Pollard andYates 1993, Pettersson et al. 2011).The monitoring was paused during the time for identification of the species, and subsequently resumed.Observations were documented using a butterfly net and a camera, and validated to species level using colleagues and literature (Eliasson et al. 2005, Söderström 2006).The pace of walking depended on habitat, accessibility, and butterfly density, but was in general approximately 3 km/h.Double counting of individuals cannot be completely ruled out, but was avoided as far as possible.Individuals not caught and identified to species level were counted and included in the (total) butterfly abundance analyses but not in any species or biodiversity measurements.In total, 42 out of 490 individuals observed (8.6%) could not be identified to species level (Appendix: Table S2).Recorder bias should not be of major concern since one of the authors (EV) was present at all sites and single-handedly documented each individual recorded.Species names (English and scientific) follow Tolman (2001) and Eliasson et al. (2005), respectively.

Data analyses
For each region, we evaluated similarity in habitat coverage between the region as a whole and the two transect approaches using Spearman Rank Correlation tests within each region followed by two-sample tests of correlation coefficients using Fisher ztransformed values (Zar 1999).
Butterfly data was compiled using values per 100 m segment as the basis for calculations.To quantify biodiversity, we used the Simpson's Diversity Index (Magurran 2004), where both the number of species as well as the abundance of the species is taken into account.The Simpson's Index (D) measures the probability that two individuals randomly selected from a sample will belong to the same species, where n i = the number of individuals in the ith species, N = the total number of individuals and S = the number of species in the sample.We represent this biodiversity measurement as 1/D, called the Simpson's Reciprocal Index.In this variant of D, 1 is the lowest possible value, representing a community containing only one species, and the maximum possible value is the number of species in the sample (Magurran 2004).Further, we used the Simpson's Evenness Index (1/DS) representing the species evenness in sites with values ranging from 0 (aggregation of species) to 1 (completely even; Magurran 2004).
Alpha (α) diversity generally measures species diversity of a defined area or habitat, whereas beta (β) diversity is used as a measure of the difference between two or more defined areas (Magurran 2004, Anderson et al. 2011).Following this, we use the term α diversity for the site-specific diversity measure (number of species per site), and β diversity for the difference between the number of species at each site and the number of species in the region).Each diversity measure results in four replicates per region.In order to assess whether α and β diversity differed between regions, we used a Wilcoxon two-sample test.
Chi-square tests of butterfly abundance were performed for each habitat category, to test if butterfly distributions differed from random expectations, i.e. if individuals distributed among habitat categories according to their relative coverage along the transects.All statistical analyses were performed in R (v. 2.15, R Core Team 2013).

GIS land cover analysis
The land cover representation of the Swedish Bird Survey transects and the volunteerselected transects correlated significantly with overall land cover within Scania (SMD land cover classes: SBS transects: r = 0.890, N = 14, p < 0.001; in free transects: r = 0.524, N = 16, p = 0.04, Appendix: Table S3) as well as in Norrbotten (SMD land cover classes: SBS transects: r = 0.733, N = 20, p < 0.001; free transects: r = 0.422, N = 23, p = 0.045, Appendix: Table S3).In Scania, the correlation between overall land cover and Swedish Bird Survey transects was significantly higher than that between overall land cover and free selection transects (Z = 2.054, p = 0.03).There was no difference in the relation between overall land cover and the two transect categories in Norrbotten (Z = 1.472, p = 0.14).

Use of bird monitoring squares
The systematically placed 2×2 km squares used by the Swedish Bird Survey correspond to 8 km butterfly transect per site (80 segments).In total we monitored 53 km transect (83%) out of the 64 km transect present in the eight squares combined.Using transects along the borders of Bird Survey squares proved to be slightly more than what was normally possible to cover within one day of butterfly recording.The landscape along the transects was sometimes difficult to traverse, and some parts of the transects were completely inaccessible.The overall distance monitored in the two regions was very similar, with 27 km in Scania and 26 km in Norrbotten.Even though some of the squares were not completely surveyed, we got a substantial amount of data with an average of 6.6 km transect monitored per site.

Habitat coverage in the two regions
Field estimates of total habitat coverage was divided fairly equal between deciduous forest (25%), coniferous forest (26%), and grassland (30%), among all transects monitored (Figure 3).The remaining three habitat categories had less coverage, with clearcut area at 4%, residential area had 9%, and fen/bog covered 6% of all transects.A large part (39%) of all monitoring was performed alongside smaller roads or paths.Scania had a higher percentage of deciduous forest while the dominant habitat in Norrbotten was coniferous forest (Figure 3).

Butterfly monitoring
A total of 490 butterfly individuals were recorded, with 250 counted in Scania (9.3 individuals per km transect) and 240 in Norrbotten (9.2 individuals per km transect).We recorded 29 different butterfly species (Appendix: Table S2), but no burnet moths.We found 22 species in Scania and 16 species in Norrbotten.Of the 22 species recorded in Scania, 13 were uniquely found in this region, and in Norrbotten, 6 out of 16 species were only found in this region.The site with the highest number of species recorded was Långberget in Norrbotten with 12 different species, whereas Kongaö in Scania was the site with most individual recordings (147 individuals) (Appendix: Table S1).

Butterfly abundance in different habitats
Butterfly abundance in different habitats was compared to the expected number of individuals relative to habitat coverage in each region (Table 1 and Figure 3).Within the habitat category deciduous forest, a significantly lower number of individuals was observed in the region of Scania than expected if the butterflies had distributed randomly over the transects (χ 2 = 24.2,p < 0.001) (Table 1), but they were not relatively fewer than expected in the deciduous forests of Norrbotten (χ 2 = 0.48, p = 0.488).The opposite was true for coniferous forests: in Norrbotten we observed significantly fewer individuals than expected (χ 2 = 45.8, p < 0.001), but in Scania the number of  individuals recorded in this habitat matched the expected number very well (χ 2 = 0.07, p = 0.787).The habitat category fen/bog inhabited significantly more individuals than expected in Scania (χ 2 = 19.2,p < 0.001), but significantly fewer individuals in Norrbotten (χ 2 = 20.2,p < 0.001) (Table 1 and Figure 3).In Norrbotten, we found a significantly higher number of individuals in clear-cut areas (χ 2 = 6.73, p = 0.009) and in residential areas (χ 2 = 6.94, p = 0.008) than expected, but not in the clear-cut areas of Scania (χ 2 = 1.14, p = 0.286), and in Scania's residential areas we even recorded significantly fewer individuals than expected (χ 2 = 5.71, p = 0.017) (Table 1).Not surprisingly, significantly more individuals were found in grasslands than expected, in both Scania (χ 2 = 8.35, p = 0.004) and Norrbotten (χ 2 = 24.2,p < 0.001) (Figure 3), though the aggregation of butterflies in grasslands were denser in Norrbotten.

Butterfly species richness in different habitats
The habitat with highest number of species recorded in total was grassland, with 16 different species recorded in Scania and 12 species in Norrbotten, however this might be partly due to relatively high coverage of grassland monitored, with 10 km grassland visited in Scania (37%), and 5.8 km in Norrbotten (22.3%) (Figure 3).Nonetheless, grassland displayed significantly higher number of species in Norrbotten than expected (N = 12, χ 2 = 4.58, p = 0.032; Table 2).Despite lower habitat coverage of coniferous forests and fens/bogs in Scania, these habitats still harboured several butterfly species (N = 10 and N = 4, respectively).The same was true for Norrbotten's clear-cut areas and deciduous forests (N = 9 and N = 5, respectively; Figure 3 and Table 2).The habitat with the lowest species diversity, despite moderate habitat coverage, was fens/bogs in Norrbotten (N = 2, Figure 3).

Biodiversity measures
The biodiversity in the two regions was measured using the Simpson Reciprocal Diversity Index (1/D) and the Simpson Evenness Index (1/DS).The Simpson Reciprocal Index for Scania (6.86) was slightly higher (although non-significantly) than the index for Norrbotten (5.10) (this difference was tested between the sites using a Wilcoxon rank sum test: W = 12, p = 0.31).The Simpson Evenness Index for Scania (0.31), was similar and not significantly different than the corresponding number for Norrbotten (0.32) (tested between the sites with a Wilcoxon rank sum test: W = 11, p = 0.47).
The mean species number within a region, the α diversity, did not differ significantly between Scania and Norrbotten (Wilcoxon rank sum test: W = 10, p = 0.661).The β diversity, defined here as the mean difference in species number between each site and the total species number of that region, was significantly higher in Scania than in Norrbotten (Wilcoxon rank sum test: W = 16, p = 0.028).

Butterfly recording using systematically placed transects
In this study we have tested the possibility of using systematically placed transects in butterfly monitoring schemes in order to get butterfly recordings with reduced volunteer habitat bias.Volunteer recorders are most often free to select the location for monitoring (Roy et al. 2007, van Swaay et al. 2008), which may result in taxonomically and geographically biased records.Dennis and Thomas (2000) and Dennis et al. (1999) demonstrated that species richness and occurrence are positively correlated with recording intensity.The volunteers' visits are biased by access (e.g. the distance from their home), the location of potential 'butterfly hot spots' (either diversity or rarity hot spots), and areas with a greater number of butterfly resources (such as seminatural grasslands).This type of bias allows for good coverage of environmentally protected sites but may not provide trends representative for species in other habitats (e.g.woodland species, Roy et al. 2007).
Our GIS analyses showed that standardised transects mirrored overall land cover better than free transects in Scania, but not significantly better in Norrbotten.The Norrbotten landscape is generally less urbanised than Scania (cf.Appendix: Table S3), and one likely explanation that free transects do not represent true habitat coverage as well in Scania is that such transects are more likely to be placed near where volunteers live (cf.Dennis and Thomas 2000;Pettersson et al. 2011) and human influence on adjacent habitats is more pronounced in Scania than in Norrbotten.This is true not only for the extent of urban structures but also for semi-natural grasslands and other habitats heavily influenced by human activities.
We found that the transects used by the Swedish Bird Survey provided good coverage of traditionally underrepresented butterfly habitats such as forests, clear-cuts and wetlands in both regions (Figure 3, cf. Brereton et al. 2011, Roy et al. 2015).Even though the habitat categorisation was broad, our results clearly indicate that there are several habitats other than grasslands harbouring butterflies.Similarly, Berg et al. (2011) demonstrated that several typically overlooked habitats in forest-dominated landscapes can have at least as high numbers of butterfly species as semi-natural grasslands.Thus, it is evident that systematically placed transects throughout the country can be an important consideration in order to reduce bias in habitat coverage.It is highly likely that some of the habitats monitored, such as deciduous and coniferous forests, are strongly underrepresented in traditional butterfly monitoring.
During the year that the present study was performed, there were still relatively few free transects in Scania and Norrbotten and a direct comparison of simultaneously collected butterfly recordings from multiple standardized and free transects was not possible.Now substantially more free transects are monitored and such a direct comparison of observations made throughout the season would be a logical next step.

Butterfly abundance and species richness in different habitats
We found major differences in butterfly abundance and species richness in the different habitats and between the two regions.In grasslands, we found more butterfly individuals than expected in both regions (Table 1), and the total number of species found was the highest in this habitat, for both regions (Table 2).Although grassland was also the second most monitored habitat in terms of total transect length for both regions (Figure 3), the number of species recorded per 100 m transect was still relatively high, and our results confirm the importance of grassland habitats for both butterfly abundance and species richness.
The clear-cut and residential areas in Norrbotten harboured more butterfly individuals than expected, but not in Scania; this region had instead significantly fewer individuals recorded in residential areas (Table 1).For the two forest habitat categories, we found fewer individuals in deciduous forests in Scania (but not in Norrbotten), and fewer individuals in coniferous forests in Norrbotten (but not in Scania) (Table 1).The contrasting differences between the regions in butterfly abundance in different habitats may be due to the drastic latitudinal and climatic differences of the two regions.Interestingly, we recorded significantly more individuals and species than expected in the habitat fen/bog in Scania, where a relatively small area harboured a substantial number of Cranberry Blues (Plebeius optilete), but also other species.These results stress the value of using representative, grid-based geographic sampling so that diversity hotspots and common habitats become neither over-nor underrepresented in monitoring schemes.
The β diversity was significantly different between the regions Scania and Norrbotten, indicating that the sites in Norrbotten harboured many of the same species, as opposed to Scania where the sites often had different species composition.This is most likely because Norrbotten has a smaller species pool compared to Scania, which harbours several rare and local butterfly species (Eliasson et al. 2005(Eliasson et al. , Öckinger et al. 2006)).Alternatively, the 2×2 km squares in Norrbotten could be more similar to each other in terms of habitat composition compared to the Scanian study squares.
Some butterfly species that we recorded in Scania were not seen in Norrbotten (N = 13), and in Norrbotten 6 out of 16 species were not found in Scania.This is likely due to the climate differences resulting in different species distributions, but it is also plausible that many species were not detected by chance.We monitored the two regions intentionally during different dates to take into account their differences in spring arrival.Recordings in Scania were performed between June 17 th and July 7 th , and recordings in Norrbotten in late July (20 th -28 th ).Because of the large latitudinal differences between the regions (Appendix: Table S1), spring can arrive approximately seven weeks later in Norrbotten than in Scania (Alexandersson 2002), so we believe the later monitoring in Norrbotten would take most of this into account.

Implications for future butterfly monitoring
A 'reduced effort' monitoring scheme is based on a higher number of transects and counted only a few times per year, as opposed to the more traditional scheme, which are to a greater extent based on more regular visits and free site selection (van Swaay et al. 2008).A reduced effort scheme with fewer visits per year but more sites monitored might make volunteers more willing to record in areas where they are less likely to record high numbers of individuals or many species (Roy et al. 2007).Stratified random sampling of sites or systematically placed sites (like in this study), will generate a lot of data for widespread butterfly species but fewer records of rare or localised species (van Swaay et al. 2008).This can, however, be mitigated by adding additional transects at sites where these species are known to occur (van Swaay et al. 2008, van Swaay andWarren 2012).
A Swedish butterfly monitoring scheme with systematically placed transects throughout the country with the intention to cover different habitats would produce a representative picture of the nation's butterfly population without introducing bias from habitat choice by the recorder.Free and systematic site selection should however not be seen as mutually exclusive.The Swedish Bird Survey started with free site selection in 1969 and added its grid-based network of geographically representative transects to the monitoring scheme in 1996.Similarly, the UK butterfly monitoring scheme has added a complementary, Wider Countryside monitoring scheme (Brereton et al. 2011, Roy et al. 2015).In other words, free selection transects and geographically representative ones can certainly exist side by side.In fact, grid networks encouraging spatially representative selection of free sites are used by some monitoring schemes (e.g.Åström et al. 2014).
Following this, we suggest that it would be valuable to complement monitoring schemes with free site selection such as the Swedish butterfly monitoring by adding standardised, grid-based sampling schemes.An exciting possibility resulting from joint monitoring of different organisms in a grid-based design is that more general biodiversity trends such as those indicated by Thomas et al. (2004) could be addressed at high spatial resolution.In general, we believe that transects for recording butterflies need to be smaller than the 2×2 km squares used by the Swedish Bird Survey.If national butterfly monitoring relies on recordings by volunteers, the transect cannot be too long as it might deter volunteers from participating.Monitoring squares the size of 1×1 km or 750×750 m, as used in study by Jonason et al. (2010), is a much more reasonable size for a butterfly monitoring scheme.and Martin Stjernman for providing the Swedish Bird Survey maps.Tobias Roth acted as reviewer for this manuscript and provided valuable comments on a previous version.LP is funded by the Swedish Environmental Protection Agency as part of the Swedish Butterfly Monitoring Scheme.
Table S3.Land cover information for sites in the two regions: a) area in km 2 for Scania and Norrbotten as a whole, for Swedish Butterfly Monitoring Scheme (free) transects (N = 5 in each region), and Swedish Bird Survey (standardised) transects (see Supplementary Material Table S1

Figure 1 .A
Figure 1.A Map over Sweden with systematically placed bird monitoring squares throughout the country, used by the Swedish Bird Survey.The two regions featured in this study are shown, with monitoring squares visited portrayed in red.The squares are 2×2 km wide, indicate location, and not shown to scale B Example map of a bird monitoring square: site Kongaö in Scania, Sweden, indicated as a red, 2×2 km square.Butterfly transects followed the periphery of bird monitoring squares (N = 8 sites) as closely as possible for as long as weather conditions and time of day permitted (Pettersson et al. 2011).Lantmäteriet, I2014/00579.

Figure 3 .
Figure 3.A Transect monitored (km) in different habitats for two geographically different regions in Sweden: the southernmost region, Scania, and the northernmost region, Norrbotten B Number of butterfly individuals recorded per 100 m transect in different habitats (mean values ± SE).Numbers in parentheses indicate total number of individuals monitored per habitat (non-integers due to proportional habitat coverage per 100 m transect) C Total number of butterfly species recorded in the different habitats.

Table 1 .
Butterfly abundance in different habitats for two geographically different Swedish regions, Scania and Norrbotten.

Table 2 .
Species richness in different habitats in two geographically different Swedish regions, Scania and Norrbotten.

,
N = 4 in each region), b) Land cover data following the SMD database which is a refined Swedish version of the CORINE land cover database (Swedish Environmental Protection Agency, 2014).Code denotes SMD raster code, CLC Code the corresponding CORINE code, SMD Code the refined Swedish Code, CLC/SMD Description explains the SMD code.Numbers presented are coverage in each habitat category for each subgroup, recalculated for ease of comparison to the smallest area of all six subgroups (Free Transects, Scania: 4.45 km 2 )