Parris 1999
Parris 1999
Abstract We investigated habitat variables influencing the species richness and composition of frog assemblages
at forest streams in southeast Queensland, Australia, to (i) determine whether these assemblage attributes can be
predicted from the nature of the local habitat, and (ii) assess the value of species richness as a measure of the
structure of frog assemblages. Stream size (as measured by the annual volume of precipitation in the upstream
catchment) was found to be a significant predictor of the species richness of an assemblage, while stream size,
understorey floristics in the riparian zone, the presence of palms and elevation were significantly correlated with
the composition of assemblages. No relationship was found between broad forest type and the species richness or
composition of frog assemblages. Species richness was not an informative description of the composition of frog
assemblages, as assemblages with identical species richness had a low average compositional similarity. Results of
this study indicate that frogs are responding to local, ground-level stream and vegetation variables, a finding of
both ecological interest and practical significance for the conservation management of forest-dwelling frogs.
Key words: amphibians, Australia, community ecology, conservation biology, forests, Mantel tests, multivariate
analysis, Queensland, species richness.
Habitat variable
recorded the presence or absence of each species at presence or absence of frog species at each site to
each site. One of three broad forest types, dry sclero- calculate the dissimilarity measure.
phyll (open) forest, wet sclerophyll (tall open) forest or These dissimilarity values were used to construct
subtropical rainforest (closed forest), was assigned to a site by site dissimilarity matrix for the 19 survey
each section of the transect on the basis of the forest sites. In order to determine the usefulness of species
overstorey. We recorded stream width and stream sub- richness as a measure of the composition of frog assem-
strate where the vegetation transect crossed the stream, blages, we calculated the average compositional dis-
and all signs of disturbance in the whole 1-ha site. similarity of assemblages with identical species richness.
Stream size was divided into two categories: small We then clustered the sites on the dissimilarity of their
(first or second order stream) and large (third order frog assemblages using the UPGMA (unweighted pair-
stream or larger). As stream order varies with the scale group method using arithmetic averages) algorithm,
of a map and the method used to calculate it, the loga- with b 5 0 (Sneath & Sokal 1973). A similar procedure
rithm of the average annual volume of precipitation in was followed to construct a dissimilarity matrix for
the upstream catchment was included as another mea- each of six habitat variables: the four used in the
sure of stream size. We calculated this by multiplying Poisson regression analyses plus stream size (small
the catchment area upstream from the survey site, mea- or large) and composition of the riparian understorey
sured from topographic maps, by the average yearly vegetation. The Bray–Curtis measure was used to
precipitation in the catchment determined using BIO- calculate the dissimilarity between sites for all the habi-
CLIM, a bioclimatic prediction system (Busby 1991; tat variables except elevation and logarithm of the
Nix & Switzer 1991). We then took the logarithm of average annual volume of precipitation in the upstream
the precipitation volume (in GL), as the values ranged catchment. A simple difference measure was used
over several orders of magnitude. for these variables (e.g. elevation at site A – elevation
at site B).
We investigated correlations between the matrix
Data analysis representing the frog assemblages and matrices repre-
senting habitat variables, using the Mantel test of
Species richness of assemblages matrix correspondence (Mantel 1967). This method
has been used in a number of previous ecological stud-
After identifying intercorrelated habitat variables and
ies (Sokal et al. 1980; Douglas & Endler 1982;
excluding them from further analysis, we investigated
Burgman 1987) and was recommended by Diniz-Filho
relationships between the species richness of frog
and Bini (1996). The Mantel test calculates the cor-
assemblages at our survey sites and habitat variables
relation between two dissimilarity matrices. The stat-
using Poisson regression. Four variables appropriate
istical significance of this correlation is calculated using
for univariate analysis were included: logarithm of
a Monte Carlo method, where the rows and columns
upstream catchment volume, broad forest type,
of one matrix are permuted randomly while those of
presence of palms and elevation. Both forward and
the other are held constant (Cliff & Ord 1981). After
backward step-wise selection of variables were used for
each permutation, the correlation between the matri-
the regression analyses.
ces is recalculated, and the resulting values form an
empirical distribution that is used for the significance
Composition of assemblages test. We used one million Monte Carlo permutations
We calculated the dissimilarity between survey sites in to test the significance of correlations. A positive cor-
terms of the composition of their frog assemblages relation indicates that sites that are similar in terms of
using the Bray–Curtis measure of compositional dis- a given habitat variable have similar frog assemblages.
similarity (Bray & Curtis 1957). This has been shown
to be a robust measure of the ecological distance
between sites (Faith et al. 1987). It is calculated as RESULTS
follows:
b1c
D5 Correlations between habitat variables
2a 1 b 1 c
where a is the number of matches (shared presences, Within the 17 habitat variables measured at each
or 1,1), and b and c are the number of (1,0) and (0,1) survey site were groups of intercorrelated variables.
mismatches, respectively. The Bray–Curtis measure Broad forest type in the riparian zone was correlated
gives a dissimilarity value (D) between 0 and 1 for each with broad forest types within 50 m of the stream,
pair of sites, with sites supporting identical assemblages composition of the overstorey vegetation and crown
having a dissimilarity of 0, and sites with no shared cover in the riparian zone. Broad forest type was also
species having a dissimilarity of 1. We used data on the correlated with a number of disturbance variables:
498 K . M . PARRIS AND M. A. McCA RT HY
evidence of fire, logging and grazing, and presence of compositional dissimilarity of frog assemblages with
the weed Lantana camara. Drier forests tended to show identical species richness were all > 0.5, and were
more evidence of these disturbances while rainforest clustered closely to the average dissimilarity of all pairs
showed less. Location (Blackall or Conondale Range) of assemblages. Although we can predict the species
was correlated with elevation, composition of the under- richness of a frog assemblage from the size of a stream,
storey vegetation, evidence of pigs and upstream catch- species richness is not an informative measure of the
ment volume. Sites in the Conondale Range were at a composition of an assemblage (Fig. 3).
higher elevation, had larger upstream catchment vol-
umes and more feral pigs than those in the Blackall
Range. There were also differences in understorey veg- Composition of frog assemblages
etation between the two locations. Stream size was cor-
related with stream width and substrate. Larger streams Cluster analysis of the 19 survey sites based on the
were wider and tended to have a rocky substrate, as compositional dissimilarity of their frog assemblages
higher flows of water wash away fine sediments from the produced four groups (A–D), although the lower two
stream bed. (groups C and D) were not well defined (Fig. 4). No
Fig. 2. Results of Poisson regression analysis of the rela- Fig. 3. A comparison between the average compositional
tionship between frog species richness at a survey site (l) and dissimilarity of frog assemblages with identical species rich-
the annual volume of precipitation in the catchment upstream ness (circles) and the average dissimilarity between all pairs
–
from the site (x): ln (l) 5 – 2.194 1 0.503 log(x), P < 0.01. of assemblages (D 5 0.75, n 5 171), represented by a line.
FROG A S S EM BLAG ES AT FOR ES T S T R EA M S 499
both the size of the stream and the adjacent forest were
influencing the suitability of a site for different frog
species. Of all the habitat variables tested, we found
stream size (small or large) was most strongly correl-
ated with the composition of frog assemblages at the
survey sites. Sites at streams of similar size tended to
have similar frog assemblages. The logarithm of the
average annual volume of precipitation in a stream’s
catchment, an alternative measure of stream size, was
also significantly correlated with the composition of
frog assemblages. No relationship was found between
broad forest type and the composition of frog assem-
blages (Table 3), indicating that ground layer vegetation
is more important for frogs than the nature of the forest
overstorey.
Elevation was found to be significantly correlated descriptions of the habitat of these frogs are expressed
with the composition of frog assemblages at a site in terms of broad forest type (Tyler 1992; Robinson
(Table 3). Our survey sites ranged in elevation from 100 1993). The assumption that diversity of forest types can
to 580 m, with those in the Conondale Range having act as a surrogate measure for faunal diversity (Austin
a higher average elevation than those in the Blackall & Meyers 1996; JANIS 1997), and that reservation of
Range (x– 5 475 m and 260 m, respectively). Elevation the range of broad forest types in a region or across
was therefore correlated with location. Three frog Australia will be sufficient to conserve the diversity of
species, Litoria tyleri, L. peronii and Pseudophryne forest-dwelling animals, is questionable. The agreed
raveni, were found only at sites in the Blackall Range, criteria for establishment of a representative reserve sys-
with L. tyleri and P. raveni at more than one site. Assa tem for Australian forests include use of information
darlingtoni, Litoria fallax, L. gracilenta and L. verreauxii on the distribution and genetic variation of individual
were found only in the Conondale Range, with Assa faunal species as well as the distribution of forest types
and L. fallax at more than one site (Table 2). The (JANIS 1997).
similarity of frog assemblages at similar elevations may Results of our study indicate that reservation of a
be a function of these interlocality differences rather range of stream sizes across the range of variation in
than an effect of elevation per se. Elevation has been understorey vegetation may be a better way to conserve
shown to have an effect on amphibian assemblages stream-breeding frogs in the forests of southeast
across broader geographical and elevational scales, such Queensland than reservation of the range of broad
as the whole of the American tropics (Duellman 1988). forest types. Data pertaining to stream size would be
Duellman’s (1988) study sites represented a 4000-m available in digital format and easy to incorporate into
range in elevation, compared to the 480 m range in our the GIS systems used by forest managers to assist with
study area. planning and reserve selection. However, mapping the
composition of the forest understorey would currently
require extensive field survey. This kind of information
Implications for management may be obtainable in the future with remote sensing
techniques such as airborne videography, which was
Conservation management of Australia’s frogs, as with recently used to assess and predict habitat complexity
other faunal groups, has largely focused on individual in eucalypt forests (Coops & Catling 1997).
species rather than on assemblages or communities.
Statistical models of habitat preferences have been
developed for some frog species in an effort to improve ACKNOWLEDGEMENTS
understanding of the environmental variables that influ-
ence their distribution (NSW National Parks & Wildlife We thank Peter West, Adrian Caneris, Ursula Grott,
Service 1994; Holloway 1997). While this approach is Denise Elias, Greg Czechura, Tony Norton, Scott
important, particularly for rare and endangered species, Osborn, Wyn Boon, Randall Donohue, Natasha West
it is not practical to manage every species individually. and Derek Cleland for their assistance in the field. Scott
Collating information from numerous single species Osborn assisted enormously with plant identifications.
models (e.g. 328 models for forest vertebrates in north- Andrew Claridge and David Lindenmayer provided
ern New South Wales, Australia; NSW National Parks valuable comments on the manuscript. This study was
& Wildlife Service 1994) into a manageable form could conducted under Q.D.P.I. permits number 788, 860,
be complicated, time-consuming and ultimately unsat- and 919, Q.N.P.W.S. permit number 2001, and
isfactory, because each species has an individual and Q.D.E.H. permit number HO/000139/95/SAA.
potentially conflicting set of habitat requirements.
Consequently, the aim of our study was to identify habi-
REFERENCES
tat variables influencing the structure of frog assem-
blages as a whole rather than the presence of their Adam P. (1994). Australian Rainforests. Oxford Biogeography
individual component species. Series no. 6. Oxford University Press, Oxford.
Species richness is not an informative measure of the Austin M. P. & Meyers J. A. (1996) Current approaches to mod-
structure of frog assemblages at forest streams in south- elling the environmental niche of eucalypts: implication for
east Queensland. Knowledge of the compositional management of forest biodiversity. For. Ecol. Manag. 85,
similarity of assemblages, and an understanding of the 95–106.
Barker J., Grigg G. C. & Tyler M. J. (1995). A Field Guide to
habitat variables that influence this, are likely to be
Australian Frogs. 2nd edn. Surrey Beatty and Sons, Chipping
more valuable for conservation management. Although Norton.
the composition of frog assemblages at forest streams Bray J. R. & Curtis J. T. (1957) An ordination of the upland forest
in southeast Queensland appears to be unrelated to the communities of southern Wisconsin. Ecol. Monog. 27,
character of the forest overstorey, many qualitative 325–49.
502 K . M . PARRIS AND M. A. McCA RT HY
Burgman M. A. (1987) An analysis of the distribution of plants MacArthur R. H. & Wilson E. O. (1967). The Theory of Island
on granite outcrops in southern Western Australia using Biogeography. Princeton University Press. Princeton.
Mantel tests. Vegetatio 71, 79–86. Mantel N. A. (1967) The detection of disease clustering and a
Busby J. R. (1991) BIOCLIM—A bioclimate analysis and pre- generalized regression approach. Cancer Res. 27, 209–20.
diction system. In: Nature Conservation: Cost Effective Martin A. A. (1967) Australian anuran life histories: some
Biological Surveys and Data Analysis (eds C. Margules & evolutionary and ecological aspects. In: Australian Inland
M. Austin), pp. 64–8. CSIRO, Canberra. Waters and Their Fauna: Eleven Studies (ed. A. H. Weatherley),
Cliff A. D. & Ord J. K. (1981). Spatial Processes: Models and pp. 175–91. Australian National University Press, Canberra.
Applications. Pion, London. McDonald K. R. & Davies M. (1990) Morphology and biology
Cogger H. G. (1994). Reptiles and Amphibians of Australia. 5th of the Australian tree frog Litoria pearsoniana (Copland)
edn. Reed, Sydney. (Anura: Hylidae). Trans. Roy. Soc. S. Aust. 114, 145–56.
Coops N. C. & Catling P. C. (1997) Utilising airborne multi- Morand A. & Joly P. (1995) Habitat variability and space util-
spectral videography to predict habitat complexity in euca- ization by the amphibian communities of the French upper-
lypt forests for wildlife management. Wildl. Res. 24, 691–703. Rhone floodplain. Hydrobiologia 301, 249–57.
Czechura G. V. (1991) The Blackall–Conondale Ranges: frogs, Nix H. A. & Switzer M. A. (1991). Rainforest Animals. Atlas of
reptiles and fauna conservation. In: The Rainforest Legacy, Vertebrates Endemic to Australia’s Wet Tropics. Kowari 1.
Vol. 2: Flora and Fauna of the Rainforests (eds G. Werren & Australian National Parks and Wildlife Service, Canberra.
P. Kershaw), pp. 311–24. Australian Government Printing NSW National Parks and Wildlife Service (1994) Fauna of north-
Service, Canberra. east NSW forests. North East Biodiversity Study Report no.
Czechura G. V. & Ingram G. J. (1990) Taudactylus diurnus and 3, New South Wales National Parks and Wildlife Service
the case of the disappearing amphibians. Mem. Queensl. Mus. (unpubl.).
29, 361–5. Oberdorff T., Hugueny B. & Guegan J. (1997) Is there an influ-
Diniz-Filho J. A. F. & Bini L. M. (1996) Assessing the relation- ence of historical events on contemporary fish species rich-
ship between multivariate community structure and environ- ness in rivers? Comparisons between Western Europe and
mental variables. Mar. Ecol. Prog. Ser. 143, 303–6. North America. J. Biogeog. 24, 461–7.
Douglas M. E. & Endler J. A. (1982) Quantitative matrix Parris K. M., Norton T. W. & Cunningham R. B. (1999) A
comparisons in ecological and evolutionary investigations. comparison of techniques for sampling amphibians in the
J. Theor. Biol. 99, 777–95. forests of southeast Queensland, Australia. Herpetologica 55,
Duellman W. E. (1988) Patterns of species diversity in anuran 271–83.
amphibians in the American tropics. Ann. Missouri Bot. Gard. Pearman P. B. (1997) Correlates of amphibian diversity in an
75, 79–104. altered landscape of Amazonian Ecuador. Cons. Biol. 11,
Faith D. P., Minchin P. R. & Belbin L. (1987) Compositional 1211–25.
dissimilarity as a robust measure of ecological distance. Petranka J. W. (1983) Fish predation: a factor affecting the spatial
Vegetatio 69, 57–68. distribution of a stream-breeding salamander. Copeia 1983,
Holloway S. (1997) Survey protocols for the stream-breeding 624–8.
frogs of Far East Gippsland: the application of habitat mod- Ricklefs R. E. & Schluter D. Eds. (1993). Species Diversity in
elling and an assessment of techniques. MSc thesis, The Ecological Communities. University of Chicago Press, Chicago.
University of Canberra, Canberra. Robinson M. (1993). A Field Guide to Frogs of Australia.
Hugueny B. & Paugy D. (1995) Unsaturated fish communities Australian Museum/Reed, Sydney.
in African rivers. Am. Nat. 146, 162–9. Sneath P. H. A. & Sokal R. R. (1973). Numerical Taxonomy.
Inger R. F. & Voris H. K. (1993) A comparison of amphibian Freeman, San Francisco.
communities through time and from place to place in Sokal R. R., Bird J. & Riska B. (1980) Geographic variation in
Bornean forests. J. Trop. Ecol. 9, 409–33. Pemphigus populicaulis (Insecta: aphididae) in eastern North
Ingram G. J. & McDonald K. R. (1993) An update on the decline America. Biol. J. Linn. Soc. 14, 163–200.
of Queensland’s amphibians. In: Herpetology in Australia: a Strahler A. N. (1964) Quantitative geomorphology of drainage
Diverse Discipline (eds D. Lunney & D. Ayres), pp. 297–303. basins and channel networks. In: Handbook of Applied
Royal Zoological Society of New South Wales, Mosman. Hydrology (ed. V. T. Chow), pp. 4–39. McGraw-Hill, New
JANIS (1997) Nationally agreed criteria for the establishment of a York.
comprehensive, adequate and representative reserve system for Tyler M. J. (1992). Encyclopedia of Australian Animals: Frogs.
forests in Australia. Commonwealth of Australia, Canberra. Angus and Robertson, Sydney.
Kats L. B., Petranka J. W. & Sih A. (1988) Antipredator defenses Zimmerman B. L. & Bierregaard R. O. (1986) Relevance of
and the persistence of amphibian larvae with fishes. Ecology the equilibrium theory of island biogeography and species-
69, 1865–70. area relations to conservation with a case from Amazonia.
Lack D. (1976). Island Birds. Blackwell Science, Oxford. J. Biogeog. 13, 133–43.