Biodiversity and Conservation 5, 493-502 (1996)

The SLOSS dilemma: a case study

ARTURO BAZ* and ANTONIO GARCIA-BOYERO Departamento de Biologia Animal, Universidad de Alcald, E-28871 Alcalti de Henares, Madrid, Spain.

Received 5 December 1994; revised and accepted 6 April 1995

Butterfly is examined on simulated archipelagoes of 2, 3, 4 and 5 holm oak forest fragments in the Guadalajara Province (central Spain). It is shown that there are more species on several small 'islands' than on a single island. Also, species number increases with the number of fragments that form the archipelago, and with the average distance between islands within the archipelago. Thus, we conclude, at least for in a system of fragmented holm oak forests in central Iberia, that the best strategy in order to maximize the conservation of species richness is the creation of a net of some small and scattered reserves. Keywords: SLOSS; butterflies: simulated archipelagoes; average distance; central Spain

Introduction SLOSS is the acronym used to debate whether the best strategy for species survival is to have a single large or several small remnant refuge patches as nature reserves (Diamond, 1975; Wilson and Willis, 1975). The debate has generated considerable controversy (see Wilcox and Murphy, 1985; Murphy and Wilcox, 1986; Lahti and Ranta, 1985, 1986), and among other approaches has been formulated in terms of maximizing species richness or in terms of minimizing extinction rates (Burkey, 1989). Diamond (1976), Jarvinen (1982) and Hubbell and Wright (1983) have questioned the appropriateness of species richness as a measure of conservation success. For some authors, the crucial question hinges on whether the function of a should be to support more species, or whether it is more important to weight species so that the reserve contains more species that would become extinct in the absence of the reserve. Another factor in the SLOSS argument is that large reserves may have minimal extinctions compared with small reserves. Chance environmental impacts are less likely to cause extinctions in large reserves, yet many small reserves, in contrast, may spread the risk. However, Spellerberg (1991) points out that large reserves are not always the most appropriate and there has been no need to involve theories from island biogeography to make that claim. Immediate conservation objectives and current circumstances will apply, but of greatest importance is the need to conserve, for as long as possible and as much as possible. These considerations have developed into the SLOPP (single large or plentiful patchy) debate which argues that it may be better to divide an area into very small patches each with many individuals, rather than a single large area (see Gilpin, 1988). This paper contains data about butterfly diversity in a forest fragmented system in order to clarify the best strategy for butterfly conservation in biotopes of central Spain from the point of view of maximizing species richness.

*To whom correspondence should be addressed. 0960-3115 © 1996 Chapman & Hall 494 Baz and Garcia-Boyero Methods Thirteen holm oak (Quercus rotundifolia Lam.) forest fragments in the Guadalajara Province of central Spain (Fig. la,b) were selected for study on the basis of variety in area, shape and complexity (Fig. lc). Each site was visited once every 15 days from 3 May to 15 September 1991, totalling 9 sample days. Samples were taken on sunny days between 10.00 and 16.00 h. The sampling scheme was based on sampling subunits (40 min) of collecting effort per site per sample (see Viejo, 1984; Baz, 1985). With this procedure, only the individuals that were caught are counted. Butterfly species that were identified while flying but were not caught are counted as one individual, and represent 1% of the total number of individuals. In transect count methods (Pollard, 1977) butterflies were identified in flight. However, the Mediterranean butterfly faunas are diverse, which renders these methods inefficient due to the practical impossibility of recognizing in flight many congeneric species at the species level. Consequently, in this work, we have selected a method more suitable for the peculiarities of the Mediterranean local faunas (see Baz and Garcfa-Boyero, 1995). The validity of the sampling method can be observed in Fig. 2 in which species accumulation curves among patches are presented. Fig. 2 shows that from the sixth day of sampling the number of species remain approximately constant. For test hypotheses about SLOSS we have followed the procedure of Simberloff and Gotelli (1984). Archipelagos of 2, 3, 4 and 5 small 'islands' have been simulated by randomly drawing pairs, trios, quartets and quintets of forest fragments, constraining the total area of each archipelago to less than that of the largest fragment (2115 ha) which has been used as the control. Thus assemblages of 13 pairs, 13 trios, etc of fragments have been established to match the 13 single fragments. Areas of the archipelagos are constrained to insure that, for each number of islands there was a uniform distribution of areas between the smallest possible for that number of islands and the largest possible (see Simberloff, 1986). For each randomly assembled archipelago, we have constructed a species list of all species found anywhere in the archipelago. Moreover, all distances were measured between pairs of islands. Then for every randomly assembled archipelago, average distance between pairs of islands within the archipelago was calculated (Simberloff, 1986). For an archipelago of only one island, the average distance was set to zero. Data about shape, isolation, tree density and structure of vegetation of these forest fragments can be found in Baz and Garcfa-Boyero (1995).

Results We obtained a grand total of 2282 individuals of 81 different species of butterflies (see Appendix) whose distribution among the holm oak fragments studied is shown in Table 1. Species number varies from 26 to 43 and no correlation was demonstrated between the number of butterfly species and the area (Ln species number = 3.49 + 0.023 Ln Area: r 2 = 11.02%;p > 0.05) (see Baz and Garc/a-Boyero, 1995 for a more comprehensive discussion of these results in relation with dynamics). Fig. 3a shows the relationship between number of species in an archipelago and total area of the archipelago, for archipelagos consisting of from one to five islands. As Wilson and Willis (1975) pointed out, when area and all other factors are held constant, a single SLOSS and butterflies 495 a

BRgHLNEGA

C

[]

1Km,

Figure 1. Map of the study site. (a) location of the Guadalajara Province. (b) location of the forest patches. Main roads and villages are also included. (c) Shape and relative size of the forest patches. island will have more species than will an archipelago comprising several islands. However, our data do not show such a pattern as can be observed in Fig. 3a. For all 65 archipelagos the regression equation is: Ln species number = 3.32 + 0.106 Ln Area where r 2 = 46.61%;p < 0.05 The same procedure has been used with the number of islands in an archipelago (see Fig. 496 Baz and Garcia-Bovero 90-

• @ -~---e 80-~

~m i /

/ ~ ¢48

E

Z m 30-

2

20-

10-

I T ..... ]F~ - I -- --T ...... ~ .... --T- ...... Y r 1 2 3 4 5 6 7 8 9 Sample day

Figure 2. Species accumulation curves among patches. A curve for all localities arc also included.

3b). The correlation between number of species and number of islands is much more than that for area and the number of species. The regression equation is

Ln species number = 3.63 + 0.316 Ln islands number where r 2 = 78.25%: p < 0.05

As the regression equation shows, the coefficient of Ln islands number is high and positive when the Wilson and Willis (1975) hypothesis predicts that it wilt be negative. Lastly, the relationship between species number and average distance among islands in an archipelago is defined by the following equation (see Fig, 3c)

Ln species number = 3.63 + 0.142 Ln distance where r 2 = 66.59%, p < 0.05 Thus, as the distance among islands increases, the butterfly fauna in an archipelago becomes richer. Since the 65 archipelagos are not all independent of one another (some of them contain SLOSS and butterflies 497 "A) 4.2 °%o oo ! 4

3.8 O0~_ [] 0 [] [] ~ 0 Lj [] 3.6 / [] E] [] 3.4 o 3.2

Log Area

4.4 pm 4.2

.| 4

3.8

3.6 B 3-4 [] 3.2 0 0 3 0.6 0 9 1.2 5 1.8 Log N [number of islands}

4.4 ~"Pl 4-2 0 0

I 4 J 3.8

3.6 B 3.4 []

3-2 i i , 3 4 Log(n+1) Distance

Figure 3. Relationships between numbers of butterfly species on 13 forest patches and 52 simulated archipelagos comprising 2-5 islands each and (a) area; (b) number of islands and (c) average distance. The following symbols are used: [S], 1 island; ~, 2 islands; ~, 3 islands; e, 4 islands; and O, 5 islands. 498 Baz and Garcia-Bovero some of the same islands) the degrees of freedom for the regression are less than 63. However, even with as few as 7, 3 and 4 degrees of freedom respectively (too few for this situation) the coefficient of determination would be significant with p < (I.05 (sec Simberloff and Gotelli, 1984; Simberloff, 1986. By means of stepwise multiple regression we have elaborated a model that allows the prediction of variation in the number of species in the simulated archipelagos, using the pool of 3 variables considered. The regression equation is:

Ln species number = 3.47 + 0.03 Ln Area + 0.19 Ln islands number + 0.05 Ln distance.

This relationship is significant (F = 112.66: p < 0.05 for the full regression) with an increase in variation accounted (r 2 adj. = 84%). A graphic representation of the validity of this model for predictions can be seen in Fig. 4.

Discussion In the SLOSS debate it can be observed that those contributors who focus on minimizing extinction rates (Wilson and Willis, 1975: Diamond, 1975, 1976: Terborgh, 1976: Fahrig and Merriam, 1985) generally come out in favour of large unfragmented reserves. On the contrary, those who focus on maximizing species richness (Simberloff and Abele, 1976.

-- Fitted [] Observed b t ...... T ..... 4.4

4.2b

[

"0 0 > Q 8

3.6~- I

a.4~- II I il J~ ' k D 3.2 ...... L ...... ~.~ 3.5 3.7 3.9 4.1 4.3 Prodicted l~gure 4. Observed species richness vs species richness predicted by the overall multiple regression equation. Bars represent 95% intervals for predictions. SLOSS and butterflies 499 Table 1. Area (ha) and species number for each sampling site

Sampling site 1 2 3 4 5 6 7 8 9 10 11 12 13 Area 3.6 3.8 7.2 26.1 39.1 71 90 153 181 189 439 670 2115 Number of species 38 26 34 43 38 39 33 43 37 33 41 41 34

1982; Jarvinen, 1982; Margules et al., 1982) claim that it is sometimes possible to pack more species into a set of smaller reserves. Obviously, the disagreement in the SLOSS debate emerges because the participants adopt different conservation goals. If our main goal is to maximize species diversity in a reserve, we must consider the advantage of capturing more species in the system containing highest habitat diversity. That groups of small sites have more species than a single large one can be explained by the habitat diversity hypothesis suggested by Game and Peterken (1984) for results of a study showing that several small woods in England typically contain more herb species than does a single large wood of equal area. If the habitat diversity hypothesis correctly explains the observation that groups of small sites often have more species than single large ones, one might expect that the further apart were the sites comprising an archipelago, the greater the habitat diversity encompassed and therefore the more species should be supported (Simberloff, 1986). Certainly, our data support this idea for butterflies on holm oak fragments, average distance explains more than 66% in the variation in species number. Consequently, for butterflies in the holm oak forests in Central Spain there is no evidence that one large refuge would be better than several small ones, if the goal is just to preserve the greatest number of species. Average distance and the number of islands are more important variables than area. In fact, while scattered small fragments totalling 1874 ha, contain almost all known butterfly species in the zone, a single area in order to contain the same number of species would need to be immensely larger (i.e., the fragment used as a control contains 34 species in 2115 ha for a total of 81, see Table 1). Despite the limitations of these results (communities are so idiosyncratic that results cannot automatically be generalized) (Simberloff and Abele, 1976,1982, 1984) our opinion is that in a fragmented landscape, historically managed as in central Spain, the best strategy in order to preserve the maximum number of species as possible is the creation of a network of small reserves. The Mediterranean is one of the richest areas of Europe in terms of with three-quarters of the total European insect fauna living in this area (Balleto and Casale, 1991). Butterflies are no exception, because Mediterranean countries contain the largest number of both total and endemic species (Munguira, 1995). Thus, the maintenance of butterfly diversity in the Mediterranean is of priority because as Dennis and Williams (1995) have pointed out, the survival of species in southern Europe is fundamental to the long-term maintenance (104-107 years) of species in northern Europe.

Acknowledgements We thank J.F. Sanchez-Rodriguez for his cooperation with field work. This work was funded by the University of Alcal~i, Research Project n ° 91 B/2 to A. Baz. 500 Baz and Garcia-Boyero References Balleto, E. and Casale, A. (1991) Mediterranean insect conservation. In The Conservation of Insects and their (N.M. Collins and J.A. Thomas, eds.) pp. 121-42. London: Academic Press. Baz, A. (1985) Ecologia y faunistica de las Mariposas de la comarca madrileha del rio Henares. Alcal~i de Henares, Spain: M.Sc. thesis, Universidad de Alcal~i. Baz, A. and Garcia-Boyero, A. (1995) The effects of forest fragmentation on butterfly communities in central Spain. J. Biogeogr. (in press). Burkey, T.V. (1989) Extinction in nature reserves. The effect of fragmentation and the importance of migration between reserve fragments. Oikos 55, 75-81. Dennis, R.L.H. and Williams, W.R, (1995) Implications of biogeographical structures for the conservation of European butterflies. In and Conservation of Butterflies (A.S. Pullin, ed.) pp. 213-29. London: Chapman & Hal|. Diamond, J.M. (1975) The island dilemma: Lessons of modern biogeographic studies for the design of natural reserves. Biol. Conserv. 7, 12945. Diamond, J.M. (1976) Island biogeography and conservation: Strategy and limitations. Science 193, 1027-9. Fahrig, J. and Merriam, G. (1985) Habitat patch connectivity and population survival. Ecology 66, 1762-8. Game, M. and Peterken, G.F. (1984) Nature reserve selection strategies in the woodland of central Lincolnshire, England. Biol. Conserv. 29, 137-81. Gilpin, M.E. (1988) A comment on Quinn and Hastings: extinction in subdivided habitats. Conserv. Biol. 2, 290-2. Hubbell, S.P. and Wright, S.J. (1983) Stochastic extinction and reserve size: a focal species approach. Oikos 41, 466-76. Kudrna, O. (1986) Aspects of the conservation of Butterflies in Europe. In: Butterflies ~fEurope. vol. 8. (O. Kudrna, ed.). Aula-Verlag. Wiesbaden. Jarvinen, O. (1982) Conservation of endangered plant populations: single large or several small reserves. Oikos 38, 301-7. Lahti, T. and Ranta, E. (1985) The SLOSS principle and conservation practice: An example. Oikos 44, 369-70. Lahti, T. and Ranta, E. (1986) Island biogeography and conservation: a reply to Murphy and Wilcox. Oikos 47, 388-9. Margules, C.R., Higgs, A.J. and Rafe, R.E. (1982) Modern biogeographic theory: are there any lessons for nature ? Biol. Conserv. 24, 115-28. Munguira, M.L. (1995) Conservation of butterfly habitats and diversity in European Mediterranean countries. In Ecology and Conservation of Butterflies (A.S. Pullin, ed,) pp, 277-89. Lond~m: Chapman & Hall. Murphy, D.D. and Wilcox, B.A. (1986) On island biogeography and conservation. Oikos 47, 385-7. Pollard, E. (1977) A method for assessing changes in the abundance of butterflies. Biol. Conserv. 12, 115-34. Simberloff, D. (1986) Design of nature reserves. In Wildlife Conservation Evaluation (M.B. Usher, ed.) pp. 315-37. London: Chapman & Hall. Simberloff, D. and Abele, L.G. (1976) Island biogeography theory and conservation practice. Science 191, 285-6. Simberloff, D. and Abele, L.G. (1982) Refuge design and island biogeography theory: Effects of fragmentation. Am. Nat. 120, 41-50. Simberloff, D. and Abele, L.G. (1984) Conservation and obfuscation: subdivision of reserves. Oikos 42, 399401. Simberloff, D. and Gotelli, N. (1984) Effects of insutarization on plant species richness in the prairie-forest ecotone. Biol. Conserv. 29, 27-46. Spellerberg, I.F. (1991) Biogeographical basis of conservation. In The Scientific Management o/ SLOSS and butterflies 501 Temperate Communities for Conservation (I.F. Spellerberg, F.B. Goldsmith and M.G. Morris, eds.) pp. 293-322. Oxford: Blackwell Scientific Publications. Terborgh, J.W. (1976) Island biogeography and conservation: strategy and limitations. Science 193, 1029-30. Viejo, J.L. (1984) Estudio faunfstico de los ropal6ceros del quejigar supramediterr~ineo de Madrid. SHILAP Rvta. Lep 12, 135--40. Wilcox, B,A. and Murphy, D.D. (1985) Conservation strategy: the effects of fragmentation on extinction. Am Nat. 125, 879-87. Wilson, E.O. and Willis, E.O. (1975) Applied biogeography. In Ecology and Evolution of Communities (M.L. Cody and J.M. Diamond, eds.) pp. 523-34. Cambridge, MS: Harvard University Press.

Appendix

List of butterflies recorded in Central Spain fragment forests during 1991. Also included are the total number of individuals. Scientific names follow the recommendations of Kudrna (1986).

1 2 3 4 5 6 7 8 9 10 11 12 13 Total Family Papilionidae Iphiclides feisthamelii (Duponchel, 1832) 1 1 2 1 1 6 Papilio machaon (Linnaeus, 1758) 1 1 2 Zerynthia rumina (Linnaeus, 1758) 1 1

Family Pieridae Anthocharis cardamines (Linnaeus, 1758) 3 31 232 1 15 Anthocharis euphenoides (Staudinger, 1869) 1 1 Aporia crataegi (Linnaeus, 1758) 1 1 1 6 1 9 19 Colias alfacariensis (Berger, 1948) 11315182 511262 57 Colias crocea (Geoffroy, 1785) 665111074453765 79 Euchloe ausonia (Hiibner, 1804) 1 1 1 2 5 Gonepteryx cleopatra (Linnaeus, 1767) 1 17 2 12 22 18 Leptidea sinapis (Linnaeus, 1758) 1 1 Pieris daplidice (Linnaeus, 1758) 10 18 5 11 11 4 4 16 20 15 12 4 9 139 Pieris rapae (Linnaeus, 1758) 7 5 11 1 6 14 12 5 6 15 2 84 Zegris eupheme (Esper, 1782) 1 1

Family Lycaenidae Aricia allous (Geyer, 1837) 2 31 4 13 2 16 Aricia cramera (Eschscholtz, 1821) 5 3 11 1 4 15 Callophrys rubi (Linnaeus, 1758) 1 111122 7 16 Celastrina argiolus (Linnaeus, 1758) 14 1 6 Cupido minimus (Fuessly, 1775) 1 1 Cyaniris semiargus (Rottemburg, 1775) 1 1 Glaucopsyche alexis (Poda, 1761) 11251 1141 11 4 41 Glaucopsyche melanops (Boisduval, 1828) 1 12 4 Lampides boeticus (Linnaeus, 1767) 1 1 1 3 L ycaena phlaeas (Linnaeus, 1761) 111 1 113 9 Nordmannia acaciae (Fabricius, 1787) 1 1 Nordmannia esculi (Hiibner, 1804) 346124 886 9 7 12 70 Nordmannia ilicis (Esper, 1779) 312624101663 5 7 9 74 Nordrnannia spini (Denis & Schiff., 1775) 1 211 3 1 9 Plebejus argus (Linnaeus, 1758) 1 2 3 Polyommatus albicans (Her.-Scha., 185l) 8 18 25 44 10 10 2 30 7 2 11 207 194 Polyommatus bellargus (Rott., 1775) 4 5 4 22 14 2 10 8 15 2 86 502 Baz and Garcia-Bo.vero 1 2 3 4 5 6 7 8 9 10 ll 12 13 Total

Polyommatus dorylas (Denis & Schiff., 1775) 1 l 2 Polyommatus escheri (Habner, 1822) 1 241 3 11 Polyommatus icarus (Rottemburg, 1775) 1331108416224151 7O Polyommatus thersites (Cantener, 1834) 11 111 5 Pseudophilotes panoptes (Htibner, 1813) 11564 91013 43 Quercusia quercus (Linnaeus, 1758) 9 21 5 1115 31 Syntarucus pirithous ( Linnae us, 1767) I

Family SatytMae Chazara briseis (Linnaeus, 1764) I 2 I 9 Coenonympha arcania (Linnaeus, 1761 ) I Coo, onympha dorus (Esper, 1782) 2 2 I 14 2 3 13 38 Coenonympha pamphilus (Linn., 17581 1 1 2 3 l0 2 1 4 1 6 2 33 Hipparchia fidia (Linnaeus. 1767) 1 1 Hipparchia hermione (Linnaeus, 1764) 1 1 2 Hipparchia semele (Linnaeus, 1758) 2 4 3 1 2 3 8 3 7 2 2 8 45 Hipparchia statilinus (Huffnagel, 17661 9 26 16 5 I1 14 23 19 21 33 4 12 16 209 Hyponephele lupina (Costa, 1836) 4 2 4 7 2 4 l 2 3 3 1 3 l 37 Hyponephele lycaon (Muschamps, 19151 2 3 I 3 3 4 I 4 I 9 31 Kanetisa circe (Fabricius, 17751 10 I 5 3 2 105 105 4 3 2 6 66 Lassiommara maera (Linnaeus, 1758) 1 1 2 Lassiommata megera (Linnaeus, 1767) 3 3 2 I 2 2 5 7 9 2 1 1 38 Maniola jurtina (Linnaeus, 1758) 4 6 I 5 2 3 4 3 4 3 1 36 Melanargia laehesis (Habner, 179(/) 24215 10 2(I 19108 t2197 2115 191 Pyronia bathseba (Fabricius, 1793) t2119 I 5 8 3 7 2610 1827 1t9 Pyroniu cecilia (Vallantin, 1894) 1 3 1 2 7 Pyronia tithonus (Linnaeus, 1771) 1 2 4 1 1 1 3 5 18 Satyrus actaea (Esper, 1781)) I 1

Family Nymphalidae A rgynnis pandora (Denis & Schiff., 1775) 1 2 I 21 7 Brenthis daphne (Denis & Schiff., 1775) 1 Brenthis hecate (Denis & Schiff., 1775) 1 1 Euphidryas aurinia (Rottemburg, 1775) 1 1 Euphidryas desfontainii (Oodart, 1819) 11015 2733 2 33 Limenitis reducta (Staudinger, 1901 ) 5 2 3 10 Melitaea cinxia (Linnaeus, 1758) 2 8 2 I 6 1 2O Melitaea didyma (Esper, 1779) 2 1 5 2 1 I 1 13 Melitaea phoebe (Denis & Schiff., t775) 3 1356311 252 32 Nymphalis polychloros (Linnaeus, 1758) 2 -1 Vanessa cardui (Linnaeus, 1758) 1 I 11 1 5

Family Hesperidae Carcharodus alceae ( Esper, 1780) 1 I Carcharodus boeticus (Rambur, I839) 2 Carcharodus lavatherae (Esper, 1780) 2 2 Erynnis tages (Linnaeus, 1758) 1 1 Hesperia comma (Linnaeus, 1758) 3 3 Pyrgus carthami (Habner, 1819) I I 2 Pyrgus cirsii (Rambur, 1839) l 1 2 I 2 2 9 Pyrgus malvae (Linnaeus, 17581 1 1 Spialia sertorius (Hoffmannsegg, 1804) 2 1 3 I 2 1 10 Syrichtus proto (Ochsenheimer, 1808) 3 6 4 4 2 5 5 14 2 10 6 60 Thymelicus acteon (Rottemburg, 1775) 1 3 4 Thymelicusflavus (Brunnich, 1763) 3 1 2 6 Thymelicus lineola (Ochsenheimer, 1808) 1 3 2 3 9