Habitat Matrix Effects on Pond Occupancy in Newts - ESA Journals

2 downloads 0 Views 1MB Size Report
ship of newt abundance to both pond and landscape variables demonstrated the ... ing the role of metapopulation functioning in newt occupancy of ponds.
Habitat Matrix Effects on Pond Occupancy in Newts PIERRE JOLY,*‡ CLAUDE MIAUD,*§ ANTHONY LEHMANN,† AND ODILE GROLET* *UMR CNRS 5023 Ecologie des Eaux Douces et des Grands Fleuves, Université Claude Bernard Lyon 1, 69622 Villeurbanne Cedex, France †Laboratoire d’Ecologie et de Biologie Aquatique, Université de Genève, 18 chemin des Clochettes, 1206 Geneva, Switzerland

Abstract: In farmlands, the population viability of many amphibians is suspected to depend on the resistance the matrix of crop fields presents to movements between ponds and terrestrial sites and movements among ponds. Over recent decades the increase in cereal growing at the expense of cattle breeding has caused a drastic change in habitat matrix in many European regions. We investigated the effect of such change on populations of three newt species ( Triturus helveticus, T. alpestris, and T. cristatus) by comparing their abundances in sites that varied in amount of cultivated ground. A multivariate regression analysis of the relationship of newt abundance to both pond and landscape variables demonstrated the negative influence of cultivated ground on abundance. The width of the uncultivated sector linking the pond to the forest was a good predictor of abundance after the influences of both pond area and fish presence were removed. Moreover, newt presence was positively related to the number of ponds within that 50-ha surrounding area, highlighting the role of metapopulation functioning in newt occupancy of ponds. The relationship between newt abundance and width of uncultivated sectors agrees with present knowledge of the orientation mechanisms that underlie migration movements in urodeles. Such a relationship between connectedness and sector width shows that narrow, linear corridors such as hedgerows may not be useful in newt conservation. Our study also highlights the need to incorporate a behavioral component of habitat use into models of connectivity in conservation biology. Efecto de la Matriz del Hábitat en la Ocupación de Estanques por Tritones Resumen: Se considera que en tierras de cultivo la disponibilidad poblacional de muchos anfibios depende de la resistencia que la matriz de las tierras de cultivo presentan a los movimientos entre estanques y los sitios terrestres, y entre estanques. Durante las décadas recientes el incremento en la producción de cereal a cambio de la cría de ganado ha ocasionado un cambio drástico en la matriz del hábitat en muchas regiones de Europa. Investigamos los efectos de estos cambios en las poblaciones de tres especies de tritones ( Triturus helveticus, T. alpestris y T. cristatus) comparando sus abundancias en sitios que variaron en la cantidad de tierra cultivada. Un análisis de regresión multivariado de las relaciones entre la abundancia de tritones y las variables de los estanques y del paisaje demostró las influencias negativas de las tierras de cultivo sobre la abundancia. La amplitud de sectores sin cultivar que conectó a los estanques con el bosque fue un buen pronosticador de la abundancia después de remover las influencias del área del estanque y la presencia de peces. Más aún, la presencia de tritones estuvo positivamente relacionada con el número de estanques dentro de un área adyacente de 50 ha, resaltando el papel del funcionamiento de la metapoblación en la ocupación de estanques. La relación entre la abundancia de tritones y la amplitud de sectores sin cultivar concuerda con el conocimiento actual sobre los mecanismos de orientación que delimitan los movimientos de migración en urodelos. Esta relación entre la conectividad y la amplitud del sector muestra que los corredores estrechos y lineales como lo son los cercos podrían no ser útiles para la conservación de tritones. Nuestro estudio también resalta la necesidad de incorporar un componente conductual del uso del hábitat cuando se modela la conectividad en estudios de conservación biológica.

‡ email [email protected] § Current address: Université de Savoie, F-73370 Le Bourget du Lac, France Paper submitted April 19, 1999; revised manuscript accepted May 3, 2000.

239 Conservation Biology, Pages 239–248 Volume 15, No. 1, February 2001

240

Habitat Matrix Effects on Newts

Introduction At a local scale, the metapopulation concept supposes that species occurrence depends on population dynamics at local and regional scales and on movements among populations (Levins 1970; Hanski & Simberloff 1997). Within a metapopulation, movements of individuals are expected to be related primarily to distance between suitable habitat patches (Marsh et al. 1999). These movements can also be influenced by the physical quality of the matrix to be crossed, particularly in animals that move on the ground (Wiens 1997; Moilanen & Hanski 1998). In this respect, a goal of conservation biology is to relate local dynamics to landscape features, including the permeability of matrix habitats to dispersion. Directly ascribing a permeability coefficient to each habitat type is not an easy task (Pither & Taylor 1998), but the influence of the matrix structure can be indirectly estimated by comparing the abundance of individuals in habitat patches differing in composition and the configuration of the matrix surrounding them. For ground-dispersing animals that have to cross open fields to reach a target habitat fragment, the extension of cultivated grounds at the expense of grasslands creates new ecological conditions that can drastically influence the movements and survival of many species (Macdonald & Smith 1990). Because of their ground-dwelling habits, amphibians are suitable candidates for an analysis of the influence of matrix quality on species occurrence and abundance (Hecnar & M’Closkey 1996; Semlitsch & Bodie 1998). Besides metapopulation functioning, the complex amphibian life cycle implies regular migrations between terrestrial and aquatic habitats, either by breeding adults or by postmetamorphic juveniles. Dispersal through the matrix can have energy costs and increased mortality risks because much of the matrix may be regarded as unsuitable habitat (Ims & Yoccoz 1997) or habitat outside the species niche (Holt 1996). The permeability these areas offers to migrating amphibians constitutes a major influence on population viability and falls within the connectivity concept (Merriam 1984; Baudry & Merriam 1988). When both low local abundance and high dispersion prevail throughout the matrix, reproductive populations experience high sensitivity to both demographic and environmental stochasticity (Shaffer 1981, 1987; Foley 1992). In such small populations, immigration is expected to play a key role in compensating for extinction risks at the local scale. If immigration is favored by the shortness of the distances that separate suitable habitat patches, the connectedness of the matrix habitat also influences its success. Many amphibian species are in this situation, and they have proved suitable for empirical studies (Gill 1978a, 1978b; Sinsch 1992; Sinsch & Seidel 1995; Halley et al. 1996; Semlitsch 1998).

Conservation Biology Volume 15, No. 1, February 2001

Joly et al.

In temperate regions of Europe, newts (genus Triturus) occur on land used for cattle breeding as long as farm ponds (natural or artificial ponds used as water sources for domestic cattle) provide suitable conditions for their reproduction. In such agricultural ecosystems, forests are usually used by the newts as aestivation and hibernation habitats (Miaud 1990; Joly & Miaud 1993). Spatial separation of woods and ponds implies that open fields have to be crossed by newts. In the past few centuries, these open fields were largely pastures. During the last few decades, however, the area devoted to grain crops has increased at the expense of pastures. Where such a change in land use has occurred, some cattle ponds have been destroyed and others have become more or less surrounded by cultivated grounds. Our goal was to evaluate the influence of such a matrix modification on the population viability of newts. We compared the abundance of three newt species in habitat sites differing in degree of enclosure by cultivated grounds. We focused particularly on the matrix environments separating each studied pond from forested grounds.

Methods Study Area and Ecological Variables In the lowlands of Bresse and Dombes (northeast of Lyon, southeastern France), agriculture results in a complex landscape matrix of pastures, cultivated grounds, forests, and fish ponds (large ponds for fish breeding). In these regions, we selected 79 farm ponds from pasture areas to provide a sample of ponds that varied in degree of enclosure by cultivated land. To focus on the problem of the replacement of pastures by cultivated fields, this selection excluded all sites where connectedness was altered by causes other than cultivation, such as major roads or urbanized areas. For each pond we collected information on intrinsic variables thought capable of influencing newt presence: depth, area, percentage of gently sloping bank (slope of ⬍20%), proportion of area occupied by floating vegetation, and fish presence. Fish present in the farm ponds were predominantly pumpkinseed (Lepomis gibbosus), black bullhead (Ictalurus melas), roach (Rutilus rutilus), and tench (Tinca tinca). We studied the influence of landscape structure within circular, 50-ha areas, with a 400-m radius centered on each pond; 400 m corresponds to the average migration distance of newts (Dolmen 1980; Griffiths 1984; Cooke 1986; R. Jehle, unpublished results). Within each 50-ha area, we measured the following variables: proportion of matrix occupied by forest (main terrestrial habitat of the newts), proportion of matrix occupied by cultivated fields, length of hedgerow (possible substitutive terres-

Joly et al.

Habitat Matrix Effects on Newts

241

tablished in the field by plotting of land use on enlarged maps from the French National Geographic Institute. We estimated areas by counting pixels. Estimation of Newt Species Richness and Abundance In the studied region, three newt species coexist as a size-structured guild. From the smallest to the largest species, this guild is composed of Triturus helveticus (palmate newt), T. alpestris (Alpine newt), and T. cristatus (crested newt). The latter is considered a threatened species and is registered in Appendix II of the Bern Convention. Sampling took place in 1991–1993 from 15 April to 15 May. This period corresponds to the peak of breeding in the Bresse-Dombes region (Miaud 1990). Because dip-netting efficiency in cattleponds is relatively high during this period, with a capture probability between 0.4 and 0.8 (unpublished results), sampling was performed only once. We caught newts with a 5-mm mesh hand-net. Sampling duration varied from 30 to 60 minutes, depending on pond area (even when no newts were found). The variable used in the analysis was a standardized estimation of abundance, namely the mean number of newts caught during a 15-minute sampling bout. Data Analyses

Figure 1. The (a) measurement of the width of uncultivated sectors (in degrees) in the study area. This variable is calculated as the sum of all the azimuths in which a straight line between the pond and a forest patch did not cross a cultivated field. Sector width is a ⫹ b ⫹ c. Pond enclosure by cultivated fields is d ⫹ e. The (b) distribution of uncultivated sectors among the sampled sites.

trial habitat), number of ponds (suspected to influence metapopulation dynamics), and angle of the uncultivated sector. This last variable was estimated as the angular sum of all the directions from which a straight line between the pond and a forest did not cross cultivated field (Fig. 1a). Although such a sector proved nonexistent in 25 ponds, the other sites constituted a balanced sample of all degrees of sector width (Fig. 1b). All these variables were estimated from maps of the stations es-

The relationships between environmental characteristics and newt abundance were assumed to result from causative processes. In this context, many authors have taken advantage of generalized linear models (GLM: McCullagh & Nelder 1989; Guisan et al. 1998) and generalized additive models (GAM: Yee & Mitchell 1991; Leathwick 1995; Lehmann 1998). We used a GAM in our study to relate occurrence and abundance in each newt species to the environmental characteristics of the ponds and the surrounding landscape. We preferred Poisson models to binomial models (presence-absence to take into account the additional information supported by abundance on the health of local populations [Winston & Angermeier 1995]). We used S-PLUS software (Mathsoft Inc., Seattle), which specifies a Poisson family with a log function as link. Significant explanatory variables were selected for each species based on a stepwise procedure (Hastie & Tibshirani 1990; Chambers & Hastie 1993). Uncultivated sectors (ncs) that were highly correlated (r ⫽ 0.79) to the percentage of wood (woa) were preselected to account for the importance of forest habitat. The recent progress in computer technologies enables easy fitting of such models (Hastie & Tibshirani 1990; Chambers & Hastie 1993). The GLM and GAM present several advantages over classical multiple regression approaches in that they integrate data from different statistical distributions with the appropriate modeling of sta-

Conservation Biology Volume 15, No. 1, February 2001

242

Habitat Matrix Effects on Newts

Joly et al.

tistical error (e.g., normal as in multiple regression, binomial for presence-absence data, ordinal for classes of abundance, Poisson or negative binomial for species or individual counts). These more general formulations of regression models allow continuous variables (e.g., depth) and/or categorical variables (e.g., presence of fish) to be combined. As a nonparametric extension of GLM, GAM builds a model directly on a smoothed function from the explanatory variables instead of pre-establishing a parametric model. This allows one to test whether a response curve is bell-shaped or not and directly estimates the additive effect of each explanatory variable from its response curves in the prediction.

Results There were no newts in 21.5% of the study sites. Occurrence of each species was negatively related to mean body size, with T. helveticus occurring in 72.2% of the ponds, T. alpestris in 63.3%, and T. cristatus in 41.7%. The three species coexisted in 36.7% of the sampled sites. In 24.1% of the sites, the assemblage was composed only of T. alpestris and T. helveticus. Some sites were occupied by only a single species (16.5%), with T. helveticus alone more frequently than any other species alone (Fig. 2a). In the same way, abundance was also negatively related to species body size, with T. helveticus reaching the highest densities and T. cristatus the lowest (Fig. 2b). Abundance did not vary significantly between T. helveticus and T. alpestris, and only T. cristatus abundance proved significantly lower than those of the other species ( p ⬍ 0.005 in each comparison; KolmogorovSmirnov test; Siegel & Castellan 1988). The overall predictive power of the proposed models (GAM) was assessed by plotting the fitted against the observed abundances (Fig. 3). For each species, the number of included variables depended on the goodness-offit of the model. The best models fitted the data with r2 values of 0.80 for T. helveticus, 0.82 for T. alpestris, and 0.88 for T. cristatus: Triturus helveticus, s(poa) ⫹ s(dep) ⫹ s(fvg) ⫹ s(hed) ⫹ s(npo) ⫹ s(aga) ⫹ s(ucs) ⫹ fish; Triturus alpestris, s(poa) ⫹ s(dep) ⫹ s(gsb) ⫹ s(fvg) ⫹ s(hed) ⫹ s(npo) ⫹ s(aga) ⫹ s(ucs) ⫹ fish; and Triturus cristatus, s(poa) ⫹ dep ⫹ s(gsb) ⫹ s(fvg) ⫹ hed ⫹ s(aga) ⫹ ucs ⫹ fish; where s(x) is smoothed variable, poa is pond area, dep is depth, gsb is gently sloping bank, fvg is floating vegetation, woa is wood area, hed is hedgerow length, npo is number of ponds, aga is arable grounds area, and ucs is uncultivated sectors. Response curves for the three species in relation to the environment variables varied from one species to another and were not always linear (Fig. 4). Relative sensitivity of abundance to each habitat variable was estimated by deviance increase when the variable was

Conservation Biology Volume 15, No. 1, February 2001

Figure 2. Pond occupancy by newts: (a) number of sites according to newt assemblages (Ta, Triturus alpestris; Tc, T. cristatus; and Th, T. helveticus) and (b) number of sites according to newt abundance.

Joly et al.

Habitat Matrix Effects on Newts

243

Figure 3. Results of the generalized additive models analysis, expressed as the relationship between fitted and observed abundance of each newt species. Abundance is n individuals per 15-minute sampling effort. dropped from the model (Table 1). Overall, the abundances of the three species were negatively related to pond area, fish presence, and hedgerow length. We predicted that abundance would be positively associated with floating vegetation and uncultivated sectors. Proportion of cultivated area was related to the abundance of the three species following a bell curve: abundance increased with cultivated area until a threshold, above which it drastically declined. This relationship appeared more drastic in T. cristatus than in the other species. Relative sensitivity of abundance to each variable varied among species (Table 1). No major variable strongly affected the abundances of the three species. Four variables were shown to heavily influence the abundances of two species: floating vegetation for T. helveticus and T. cristatus, hedgerow length for T. helveticus and T. alpestris, arable-grounds area for T. helveticus and T. cristatus, and uncultivated sector for T. helveticus and T. alpestris. Other variables appeared less influential. To distinguish the influences of uncultivated sectors from that of forested areas (these two variables were positively correlated), we investigated the relationship between newt abundance and uncultivated sectors in a subsample of 20 sites where forested area varied within the limited range from 25% to 40% (the uncultivated sector varied from 0 to 360⬚). Within this selection of sites, high abundances were observed only when the width of the uncultivated sector exceeded 180⬚ for T. helveticus and T. alpestris and 140⬚ for T. cristatus (Fig. 5), suggesting a threshold effect. Considering only newt occurrence (and not abundance as in GAM analysis), a relationship was detected with number of ponds within the 50 ha (Fig. 6a). Whereas the probability of occupancy by a population of each newt species was high when the number of ponds exceeded 5/per 50 ha, it decreased as the ponds became more rare. The shape of such a relationship differed among species. Whereas neither linear regression

nor threshold were perceptible in T. helveticus (r 2 ⫽ 0.29), the occurrence of T. alpestris drastically increased when pond density increased from 2 to 3/per 50 ha (nonlinear regression; r2 ⫽ 0.43). In contrast, a linear regression model fitted the relationship between T. cristatus occurrence and pond density ( p ⫽ 0.04; r2 ⫽ 0.59) Moreover, the relationship between pond number and newt abundance also differed among species. In the smaller species, T. helveticus and T. alpestris, observed abundance was not high in situations where the number of ponds was the highest, but rather in those with 2 or 3 ponds/per 50 ha (Fig. 6b). There was only one pond frequency (5/per 50 ha) at which T. cristatus abundance was occasionally greater.

Discussion The distributions of the three newt species conform globally with the abundance-occupancy theory, with the commonest species also being locally the most abundant (Hanski 1982; Brown 1984; Gaston 1996). Despite low local abundances, T. cristatus was present in one-third of the sites. The relative rarity of this large species is probably related to the relative smallness of the studied ponds that prevent populations from reaching high numbers. Consequently, T. cristatus is probably the species most sensitive to habitat fragmentation (highest local extinction risk) (Shaffer 1987). Generalized additive models made it possible to jointly analyze influences of both aquatic and terrestrial components of the landscape environments. Variables describing the physico-chemical quality of water were not included in this analysis because of the usual low relevance of these variables compared to more macroscopic ecological ones in explaining newt distribution (Cooke & Frazer 1976; Beebee 1985; Giacoma 1988; Pavignano et

Conservation Biology Volume 15, No. 1, February 2001

244

Habitat Matrix Effects on Newts

al. 1990; Ildos & Ancona 1994). The four pond variables contributed significantly to the regressive model of each species, although their respective contributions varied. Pond area, extent of floating vegetation, and fish presence were the most influential variables in the models. Whereas most landscape variables (except forest area) influenced the abundances of the three species, number of ponds did not influence those of T. cristatus abundance. Pond area, hedgerow length, high proportion of land occupied by cultivated fields, and fish presence were negatively related to the abundance of each species. The negative relationship between abundance and pond area may result from a habitat-island effect (MacArthur & Wilson 1967), the highest abundances being predicted in the smallest patches (50–80 m2 ), which is consistent with the data we collected. But site area may also have lead to biased sampling efficiency, and the possibility of hand-netting yielding false absence or low abundance was probably higher in the largest habitat patches

Joly et al.

(ponds). Nevertheless, the variability of this relationship between area and abundance increased with increasing area, suggesting a higher sensitivity of abundance to factors other than area in the largest sites. More intriguing was the negative relationship between newt abundance and hedgerow length, suggesting that hedgerows do not constitute substitutive terrestrial habitats for the species of interest or that their effects on population dynamics were masked by other population processes. As expected, newt abundance was negatively related to the proportion of land occupied by cultivated fields. In two species (T. helveticus and T. cristatus), this relationship was bell-shaped, suggesting that newt abundance within a pond is enhanced in landscapes where the proportion of area in cultivation remains moderate. Such an enhancement of abundance probably results from a crowding effect. The negative influence of fish presence suggested by the GAMs confirms findings from previous studies (Aronsson & Stenson 1995). T. cristatus may have a greater vulnerability to fish presence than T.

Figure 4. Response curves of the abundances of Triturus helveticus, T. alpestris, and T. cristatus to the environmental characteristics (pond and surrounding area) kept in the generalized additive models analysis (center lines). The y-axes are log-scale (link function), are based on partial residuals, and indicate the relative influence of each explanatory variable on the prediction. Distance between the upper and lower curves indicates two times the pointwise standard errors for each curve or factor level. N.S., not significant and for fish variable; 0, absent; 1, present.

Conservation Biology Volume 15, No. 1, February 2001

Joly et al.

Habitat Matrix Effects on Newts

245

Table 1. Null deviance, residual deviance, degree of freedom, and residual df of each model of three newt species after the stepwise selection of explanatory variables. Parameter* Null deviance Residual deviance Degrees of freedom Residual df DI: pond area DI: depth DI: gently sloping banks DI: floating vegetation DI: hedgerow length DI: number of ponds DI: arable-grounds area DI: non-cultivated sector DI: fish

Triturus helveticus

Triturus alpestris

Triturus cristatus

1220 331 78 49.4 41 (4.6) 41 (4.6) ns 81 (9.1) 72 (8.1) 22 (2.5) 83 (9.3) 73 (8.2) 28 (3.1)

1054 261 78 45.6 40 (5) 17 (2.1) 25 (3.1) 10 (1.3) 65 (8.2) 25 (3.1) 14 (1.8) 66 (8.3) 10 (1.3)

230 40 78 58 12 (6.3) 3 (1.6) 25 (13.2) 44 (23.1) 11 (5.8) ns 37 (19.5) 4 (2.1) 16 (8.4)

*Residual deviance, deviance not supported by the selected variables. Difference between degrees of freedom and residual df is due to smoothing functions. Deviance increase (DI), marginal increases in deviance associated with dropping each selected variable from the species models, expressed by percentage in parentheses (deviance increase/[null deviance ⫺ residual deviance] ⫻ 100) (ns, nonsignificant influence).

helveticus and T. alpestris because of the nectonic behavior of its larvae. Of the environmental variables considered, only floating vegetation and uncultivated sectors were positively associated with newt abundance. Floating vegetation is used by the newts for egg laying, particularly in T. helveticus and T. cristatus (Miaud 1995). The positive relationship between the width of the uncultivated sector and newt abundances validates the hypothesis that this sector is the main component of connectivity between terrestrial and aquatic habitats in newts. Because the relationship between forest area and the abundances of both T. alpestris and T. cristatus was not significant, connectedness was the main landscape determinant of the abundances of these two species. The relationship between the width of the uncultivated sector and newt abundance suggests that the broader the sector the more individuals use it. This is in accordance with our knowledge of the orientation mechanisms that underlie urodele migrations. When migrating between their terrestrial and aquatic sites, these animals take a straight path (Shoop 1968; Madison & Shoop 1970; Shoop & Doty 1972; Douglas & Monroe 1981; Kleeberger & Werner 1982, 1983; Madison 1997). Most newt individuals are faithful to a given aquatic site ( Joly & Miaud 1989) and probably also to a terrestrial site (R. Jehle, unpublished results). At present, pilotage toward chemical cues emanating from the target site and magnetic compass are the only two orientation mechanisms that have been demonstrated experimentally (Hershey & Forester 1980; Phillips 1986a, 1986b, 1987; McGregor & Teska 1989; Phillips & Borland 1992; Joly & Miaud 1993). Both mecha-

Figure 5. Relationship between newt abundance and the width of the uncultivated sector (in degrees) in sites where woods covered from 25% to 40% of the 50-ha area surrounding the pond (Fig.1).

nisms determine straight trajectories. From one year to the next, most individuals follow the same path (P.J. & C.M., unpublished data). No evidence of the use of a cognitive map, which would allow animals to skirt a cultivated field, is presently available. Consequently, there is no evidence that linear corridors such as hedgerows or ditches can be used by newts as corridors allowing an efficient substitutive solution to reduction of the extent of fields permeable to their movements. Favorable migration habitat between terrestrial and aquatic sites has to be as broad as possible. In this context, newts differ from small rodents, for which linear corridors such as fencerows can ensure a valuable level of connectivity (Fahrig & Merriam 1985; Bennett 1990; Merriam & Lanoue 1990; La Polla & Barrett 1993).

Conservation Biology Volume 15, No. 1, February 2001

246

Habitat Matrix Effects on Newts

Figure 6. Relationship (a) between newt occurrence and number of ponds (n) within the 50-ha area surrounding the pond (1, n ⫽ 13; 2, n ⫽ 25; 3, n ⫽ 16; 4, n ⫽ 8; 5, n ⫽ 11; ⬎6, n ⫽ 6; hel, T. helveticus; alp, T. alpestris; cri, T. cristatus) and (b) for each species, between newt abundance and number of ponds in the surrounding area. Size of circles is positively related to numbers of sites. Conservation Biology Volume 15, No. 1, February 2001

Joly et al.

The width of the uncultivated sector probably modulates the number of newts able to complete their life cycles, whereas pond frequency in the immediate landscape probably influences metapopulation dynamics. Dispersal and migration, favored by a low interpond distance, are expected to ensure the continued presence of newts at a site, irrespective of their local reproductive success (Levins 1970; Gill 1978b; Hanski 1991). Our results (Fig. 6) confirm this hypothesis by establishing a positive relationship between newt occurrence in a pond and the number of ponds in close vicinity, despite variation in other environmental factors. The modality of the dependence of newt occurrence on pond frequency varied among species. Although no clear model was established for T. helveticus, a threshold effect was detected in T. alpestris and a linear relationship in T. cristatus. At least for these two species, abundance reached its highest levels in landscape configurations with intermediate pond frequency, suggesting that crowding occurs in configurations where spread to other ponds is restricted by isolation. Such pond frequencies also correspond to a habitat availability below which the species locally disappears (Wiens 1997). Other mechanisms can also be assumed to explain high abundances, such as high local connectivity, conspecific attraction, or other factors limiting population increase despite pond availability. In species with complex life cycles, processes operating in both aquatic and terrestrial phases are expected to contribute to the regulation of population numbers. The contribution of each phase may be expected to depend on the spatial configuration of these aquatic and terrestrial habitat types (Semlitsch 1998). Aquatic habitats are expected to influence population dynamics more prominently when they are scantily distributed, and the influence of terrestrial habitats is expected to predominate in configurations where ponds are plentiful. High newt abundances in situations of intermediate pond frequency support this hypothesis. But the ponds we sampled were selected to avoid the influence of ecological succession on abundance (newly created ponds or ponds filled with alluvium, which can be another important determinant of population dynamics) (Skelly & Meir 1997). Our study confirms the wide applicability of the connectivity concept to conservation biology (Hess 1996; Drechsler & Wissel 1998). It also demonstrates that the use of such a concept cannot be restricted to corridor design because the efficiency of corridors depends on the relationship between space-use behavior and landscape configuration for particular taxa. The integration of such behavioral contingencies in conservation biology represents a challenge for the coming years (Lima & Zollner 1996; Tischendorf et al. 1998). On the other hand, in analyzing only two levels of habitat determinants—local patch and close landscape—our study does

Joly et al.

not account for habitat influence at a broader scale. Further investigations should incorporate such an influence by using geographic information systems and landscape analysis in a metapopulation dynamics approach (Moilanen & Hanski 1998).

Acknowledgments This work was supported by grant 8E47 from the French Environment Ministry. It could not have been performed successfully without the efficient cooperation of W. Maille. A first draft was critically reviewed by G. Barker, B. Danielson, L. Fahrig, M. Groom, E. Pattee, P. SjögrenGulve, and anonymous referees. Literature Cited Aronsson, S., and J. A. E. Stenson. 1995. Newt-fish interactions in a small forest lake. Amphibia-Reptilia 16:177–184. Baudry, J., and G. Merriam. 1988. Connectivity and connectedness: functional versus structural patterns in landscapes. Pages 23–28 in K. F. Schreiber, editor. Connectivity in landscape ecology. Münstersche Geographische Arbeiten 29. Schöningh, Paderborn, Germany. Beebee, T. J. C. 1985. Discriminant analysis of amphibian habitat determinants in south-east England. Amphibia-Reptilia 6:35–43. Bennett, A. F. 1990. Habitat corridors and the conservation of small mammals in a fragmented forest environment. Landscape Ecology 4:109–122. Brown, J. H. 1984. On the relationship between abundance and distribution of species. The American Naturalist 124:255–279. Chambers, J. M., and T. J. Hastie. 1993. Statisitical models in Computer sciences series. Chapman and Hall, London. Cooke, A. 1986. Studies of the crested newt at Shillow Hill, 1984– 1986. Herpetofauna News 6:4–5. Cooke, A. S., and J. F. D. Frazer. 1976. Characteristics of newt breeding sites. Journal of Zoology 178:223–236. Dolmen, D. 1980. Local migration, rheotaxis, and philopatry by Triturus vulgaris within a locality in Central Norway. British Journal of Herpetology 6:151–158. Douglas, M. E., and B. L. Monroe. 1981. A comparative study of topographical orientation in Ambystoma (Amphibia: Caudata). Copeia: 460–463. Drechsler, M., and C. Wissel. 1998. Trade-offs between local and regional scale management of metapopulations. Biological Conservation 83:31–41. Fahrig, L., and G. Merriam. 1985. Habitat patch connectivity and population survival. Ecology 66:1762–1768. Foley, P. 1992. Predicting extinction times from environmental stochasticity and carrying capacity. Conservation Biology 8:124–137. Gaston, K. J. 1996. The multiple forms of the interspecific abundancedistribution relationship. Oikos 76:211–220. Giacoma, C. 1988. The ecology and distrtibution of newts in Italy. Annali de Museo Zoologico degli Universita di Napoli 26:49–84. Gill, D. E. 1978a. The metapopulation ecology of the red-spotted newt, Notophthalmus viridescens ( Rafinesque). Ecological Monographs 48:145–166. Gill, D. E. 1978b. Effective population size and interdemic migration rates in a metapopulation of the red-spotted newt, Notophthalmus viridescens ( Rafinesque). Evolution 32:839–849. Griffiths, R. A. 1984. Seasonal behaviour and intrahabitat movements in an urban population of smooth newts, Triturus vulgaris (Amphibia Salamandridae). Journal of Zoology ( London) 203:241–251.

Habitat Matrix Effects on Newts

247

Guisan, A., J-P. Theurillat, and F. Kienast. 1998. Predicting the distribution of plant species in an alpine environment. Journal of Vegetation Science 9:65–74. Halley, J. M., R. S. Oldham, and J. W. Arntzen. 1996. Predicting the persistence of amphibian populations with help of a spatial model. Journal of Applied Ecology 33:455–470. Hanski, I. 1982. Dynamics of regional distribution: the core and satellite hypothesis. Oikos 38:210–221. Hanski, I. 1991. Single-species metapopulation dynamics: concepts, models and observations. Biological Journal of the Linnean Society 42:17–38. Hanski, I., and D. Simberloff. 1997. The metapopulation approach, its history, conceptual domain, and application to conservation. Pages 5–42 in I. Hanski and M. E. Gilpin, editors. Metapopulation dynamics: ecology, genetics and evolution. Academic Press, San Diego. Hastie, T. J., and R. J. Tibshirani. 1990. Generalized additive models. Chapman and Hall, London. Hecnar, S. J., and R. T. M’Closkey. 1996. Regional dynamics and the status of amphibians. Ecology 77:2091–2097. Hershey, J. L., and D. C. Forester. 1980. Sensory orientation in Notophthalmus viridescens (Amphibia: Salamandridae). Canadian Journal of Zoology 58:266–276. Hess, G. R. 1996. Linking extinction to connectivity and habitat destruction in metapopulation models. The American Naturalist 148: 226–236. Holt, R. D. 1996. Adaptive evolution in source-sink environments: direct and indirect effects of density-dependence on niche evolution. Oikos 75:182–192. Ildos, A. S., and N. Ancona. 1994. Analysis of amphibian habitat preferences in a farmland area ( Po plain, Northern Italy). Amphibia-Reptilia 15:307–316. Ims, R. A., and N. G. Yoccoz. 1997. Studying transfer processes in metapopulations: emigration, dispersal and colonization. Pages 247– 265 in I. Hanski and M. E. Gilpin, editors. Metapopulation dynamics: ecology, genetics and evolution. Academic Press, San Diego. Joly, P., and C. Miaud. 1989. Fidelity to the breeding site in the Alpine newt Triturus alpestris. Behavioural Processes 19:47–56. Joly, P., and C. Miaud. 1993. How does a newt find its pond? The role of chemical cues in migrating newts (Triturus alpestris). Ethology Ecology and Evolution 5:447–455. Kleeberger, S. R., and J. K. Werner. 1982. Home range and homing behaviour of Plethodon cinereus in northern Michigan. Copeia: 409–415. Kleeberger, S. R., and J. K. Werner. 1983. Post-breeding migration and summer movement of Ambystoma maculatum. Journal of Herpetology 17:176–177. La Polla, V. N., and G. W. Barrett. 1993. Effects of corridor width and presence on the population dynamics of the meadow vole (Microtus pennsylvanicus). Landscape Ecology 8:25–37. Leathwick, J. R. 1995. Climatic relationships of some New Zealand forest tree species. Journal of Vegetation Science 6:237–248. Lehmann, A. 1998. GIS modeling of submerged macrophytes distribution using generalized additive models. Plant ecology 139:113–124. Levins, R. 1970. Extinction. Pages 77–107 in M. M. Gesternhaber editor. Some mathematical problems in biology. American Mathematical Society, Providence, Rhode Island. Lima, S. L., and P. A. Zollner. 1996. Towards a behavioral ecology of ecological landscapes. Trends in Ecology and Evolution 11:131–134. MacArthur, R. H., and E. O. Wilson. 1967. The theory of island biogeography. Princeton University Press, Princeton, New Jersey. Macdonald, D. W., and H. Smith. 1990. Dispersal, dispersion, and conservation in the agricultural ecosystem. Pages 18–64 in R. G. H. Bunce and D. C. Howard editors, Species dispersal in agricultural habitats. Belhaven Press, London. Madison, D. M. 1997. The emigration of radio-implanted spotted salamanders, Ambystoma maculatum. Journal of Herpetology 31: 542–551. Madison, D. M., and C. R. Shoop. 1970. Homing behavior, orientation,

Conservation Biology Volume 15, No. 1, February 2001

248

Habitat Matrix Effects on Newts

and home range of salamanders tagged with tantalum-182. Science 168:1484–1487. Marsh, D. M., E. H. Fegraus, and S. Harrison. 1999. Effects of breeding pond isolation on the spatial and temporal dynamics of pond use by the tungara frog, Physalemus pustulosus. Journal of Animal Ecology 68:804–814. McCullagh, P., and J. A. Nelder. 1989. Generalized linear models. Chapman and Hall, London. McGregor, J. H., and W. R. Teska. 1989. Olfaction as an orientation mechanism in migrating Ambystoma maculatum. Copeia: 779–781. Merriam, G. 1984. Connectivity: a fundamental ecological characteristic of landscape pattern. Pages 5–15 in J. Brandt and P. Agger, editors. Proceedings of the first international seminar on methodology in landscape ecological research and planning. Volume 1. Roskilde Universitetsforlag GeoRuc, Roskilde, Denmark. Merriam, G., and A. Lanoue. 1990. Corridor use by small mammals: field measurement for three experimental types of Peromyscus leucopus. Landscape Ecology 4:123–131. Miaud, C., 1990. La dynamique des populations subdivisées: étude comparative chez trois amphibiens urodèles (Triturus alpestris, T. helveticus et T. cristatus). Ph.D. thesis. University of Lyon, France. Miaud, C. 1995. Oviposition site selection in three species of European newts (Salamandridae) genus Triturus. Amphibia-Reptilia 16: 265–272. Moilanen, A., and I. Hanski. 1998. Metapopulation dynamics: effects of habitat quality and landscape structure. Ecology 79:2503–2515. Pavignano, I., C. Giacoma, and S. Castellano. 1990. A multivariate analysis of amphibian habitat determinants in north western Italy. Amphibia-Reptilia 11:311–324. Phillips, J. B. 1986a. Magnetic compass orientation in the eastern redspotted newt (Notophthalmus viridescens). Journal of Comparative Physiology (A) 158:103–109. Phillips, J. B. 1986b. Two magnetoreception pathways in a migratory salamander. Science 233:765–767. Phillips, J. B. 1987. Laboratory studies of homing orientation in the eastern spotted newt, Notophthalmus viridescens. Journal of Experimental Biology 131:215–229. Phillips, J. B., and S. C. Borland. 1992. Wavelength specific effects of light on magnetic compass orientation of the eastern red-spotted

Conservation Biology Volume 15, No. 1, February 2001

Joly et al.

newt Notophthalmus viridescens. Ethology Ecology and Evolution 4:33–42 Pither, J., and P. D. Taylor. 1998. An experimental assessment of landscape connectivity. Oikos 83:166–174. Semlitsch, R. D. 1998. Biological delineation of terrestrial buffer zones for pond-breeding salamanders. Conservation Biology 12:1113–1119. Semlitsch, R. D., and J. R. Bodie. 1998. Are small, isolated wetlands expendable? Conservation Biology 12:1129–1133. Shaffer, M. L. 1981. Minimum population sizes for species conservation. Bioscience 31:131–134. Shaffer, M. L. 1987. Minimum viable populations: coping with uncertainity. Pages 69–86 in M. E. Soulé, editor. Viable populations for conservation. Cambridge University Press, Cambridge, United Kingdom. Shoop, C. R. 1968. Migratory orientation of Ambystoma maculatum: movements near breeding ponds and displacements of migrating individuals. Biological Bulletin 135:230–238. Shoop, C. R., and T. L. Doty. 1972. Migratory orientation by marbled salamanders (Ambystoma opacum) near breeding area. Behavioral Biology 7:131–136. Siegel, S., and M. J. J. Castellan. 1988. Statistics for the behavioural sciences. McGraw-Hill, New York. Sinsch, U. 1992. Structure and dynamic of a natterjack toad metapopulation (Bufo calamita). Oecologia 90:489–499. Sinsch, U., and D. Seidel. 1995. Dynamics of local and temporal breeding assemblages in a Bufo calamita metapopulation. Australian Journal of Ecology 20:351–361. Skelly, D. K., and E. Meir. 1997. Rule-based models for evaluating mechanisms of distributional change. Conservation Biology 11:531–538. Tischendorf, L., U. Irmler, and R. Hingst. 1998. A simulation experiment on the potential of hedgerows as movement corridors for forest carabids. Ecological Modeling 106:107–118. Wiens, J. A. 1997. Metapopulation dynamics and landscape ecology. Pages 43–68 in I. Hanski and M. E. Gilpin, editors. Metapopulation biology. Academic Press, San Diego. Winston, M. R., and P. L. Angermeier. 1995. Assessing conservation value using centers of population density. Conservation Biology 9: 1518–1527. Yee, T. W., and N. D. Mitchell. 1991. Generalized additive models in plant ecology. Journal of Vegetation Science 2:587–602.