Predicting habitat distribution and frequency from plant species co-occurrence data

Aim Species frequency data have been widely used in nature conservation to aid management decisions. To determine species frequencies, information on habitat occurrence is important: a species with a low frequency is not necessarily rare if it occupies all suitable habitats. Often, information on habitat distribution is available for small geographic areas only. We aim to predict grid‐based habitat occurrence from grid‐based plant species distribution data in a meso‐scale analysis. Location The study was carried out over two spatial extents: Germany and Bavaria. Methods Two simple models were set up to examine the number of characteristic plant species needed per grid cell to predict the occurrence of four selected habitats (species data from FlorKart, http://www.floraweb.de ). Both models were calibrated in Bavaria using available information on habitat distribution, validated for other federal states, and applied to Germany. First, a spatially explicit regression model (generalized linear model (GLM) with assumed binomial error distribution of response variable) was obtained. Second, a spatially independent optimization model was derived that estimated species numbers without using spatial information on habitat distribution. Finally, an additional uncalibrated model was derived that calculated the frequencies of 24 habitats. It was validated using NATURA2000 habitat maps. Results Using the Bavarian models it was possible to predict habitat distribution and frequency from the co‐occurrence of habitat‐specific species per grid cell. As the model validations for other German federal states were successful, the models were applied to all of Germany, and habitat distribution and frequencies could be retrieved for the national scale on the basis of habitat‐specific species co‐occurrences per grid cell. Using the third, uncalibrated model, which includes species distribution data only, it was possible to predict the frequencies of 24 habitats based on the co‐occurrence of 24% of formation‐specific species per grid cell. Predicted habitat frequencies deduced from this third model were strongly related to frequencies of NATURA2000 habitat maps. Main conclusions It was concluded that it is possible to deduce habitat distributions and frequencies from the co‐occurrence of habitat‐specific species. For areas partly covered by habitat mappings, calibrated models can be developed and extrapolated to larger areas. If information on habitat distribution is completely lacking, uncalibrated models can still be applied, providing coarse information on habitat frequencies. Predicted habitat distributions and frequencies can be used as a tool in nature conservation, for example as correction factors for species frequencies, as long as the species of interest is not included in the model set‐up.     

Do McKinnon lists provide reliable data in bird species frequency? A comparison with transect-based data

Although occurrence-based listing methods could provide reliable lists of species composition for a site, the effective reliability of this method to provide more detailed information about species frequency (and abundance) has been rarely tested. In this paper, we compared the species frequencies obtained for the same set of species-rich sites (wetlands of central Italy) from two different methods: McKinnon lists and line transects. In all sites we observed: (i) rapid cumulating curves of line transect abundance frequencies toward the asymptote represented by the maximum value in McKinnon occurrence frequency; (ii) a large amount of species having a low frequency with line transect method showing a high range of variation in frequency obtained by McKinnon lists; (iii) a set of species having a subdominant (>0.02-<0.05) and dominant species (>0.05) frequency with line transect showed all the highest value in McKinnon frequency. McKinnon lists provides only a coarse-grained proxy of species frequency of individuals distinguishing only between common species (having the highest values of McKinnon frequency) and rare species (all the other species). Although McKinnon lists have some points of strength, this method does not discriminate the frequencies inside the subset of common species (sub-dominant and dominant species). Therefore, we suggest a cautionary approach when McKinnon frequencies should be used to obtain complex univariate metrics of diversity.     

Connections between ecology,biogeography, and paleobiology: Relationship between local abundance and geographic distribution in fossil and recent molluscs

Summary We used data on Contemporary and Pleistocene molluscs at one site in the Gulf of California to evaluate and extend earlier ideas about the relationship between local abundance and geographic distribution. For each species whose shells occurred in one Recent and two Pleistocene deposits, we measured its abundance in the sample and relative latitudinal position within its contemporary geographic range. Species near the edges of their ranges showed uniformly low abundances, whereas those near the centres exhibited a wide range of abundances. Species near the edges of their ranges also appear to have exhibited greater changes in abundance, including more colonization and extinction events, between the Pleistocene interglacial sample and the Recent one. The constraint of location in the geographic range on maximal local and regional abundance appears to offer an example of a connection between patterns and processes on local, regional, and geographical scales. Characteristics of community structure, such as relative abundance of individual species and frequency of local co-existence of multiple species, may be influenced by the location of the sample site with respect to the geographic ranges of the constituent species. These results demonstrate emergent, statistical features of population ecology and community organization that are manifest over geographic space and evolutionary time.     

Fire frequency and tree canopy structure influence plant species diversity in a forest-grassland ecotone

Disturbances and environmental heterogeneity are two factors thought to influence plant species diversity, but their effects are still poorly understood in many ecosystems. We surveyed understory vegetation and measured tree canopy cover on permanent plots spanning an experimental fire frequency gradient to test fire frequency and tree canopy effects on plant species richness and community heterogeneity within a mosaic of grassland, oak savanna, oak woodland, and forest communities. Species richness was assessed for all vascular plant species and for three plant functional groups: grasses, forbs, and woody plants. Understory species richness and community heterogeneity were maximized at biennial fire frequencies, consistent with predictions of the intermediate disturbance hypothesis. However, overstory tree species richness was highest in unburned units and declined with increasing fire frequency. Maximum species richness was observed in unburned units for woody species, with biennial fires for forbs, and with near-annual fires for grasses. Savannas and woodlands with intermediate and spatially variable tree canopy cover had greater species richness and community heterogeneity than old-field grasslands or closed-canopy forests. Functional group species richness was positively correlated with functional group cover. Our results suggest that annual to biennial fire frequencies prevent shrubs and trees from competitively excluding grasses and prairie forbs, while spatially variable shading from overstory trees reduces grass dominance and provides a wider range of habitat conditions. Hence, high species richness in savannas is due to both high sample point species richness and high community heterogeneity among sample points, which are maintained by intermediate fire frequencies and variable tree canopy cover.     

The anatomy, physiology, functional significance and evolution of specialized hearing organs of gerbilline rodents

Middle and inner ear anatomy correlates with neurophysiological responses to a wide range of sound frequencies for species of the Gerbillinae representing generalized, intermediate, and specialized anatomical conditions. Neurophysiological data were recorded from 81 specimens of 13 species representing six genera. Anatomical parameters involved in the process of hearing were correlated with the neurophysiological data to assess the effects of different degrees of anatomical specialization on hearing. The 13 species tested in this manner have graphic curves of auditory sensitivity of remarkably similar disposition over the frequencies tested and to those published for Kangaroo Rats. Ears with anatomical specializations show greater auditory sensitivity. The natural history of the Gerbillinae, particularly the kinds of predators, degree of predation, and habitat is reviewed and utilized to interpret the significance of the degree of auditory specialization in the forms studied and to evaluate the prevailing hypothesis that these specializations enhance the ability of these rodents to survive in open desert situations by detecting and evading predators. The middle ear anatomy of five additional genera and species was also studied. Thus, data on the entire spectrum of gerbilline middle ear morphology provide an evolutionary sequence. Certain anatomical parameters of the organ of Corti show a degree of specialization parallel to that of features of the middle ear. The morphological changes and possible functional roles of these features are considered. A very high correlation exists for degree of specialization and aridity of habitat, thus specialization increases with increasing aridity. This increased specialization may result from more effective predation in open xeric environments. Auditory acuity for a wide range of low frequency sounds augmented by auditory specialization is hence more advantageous here. There does not appear to be selection for hearing at particular frequencies in this range. The peaks of greatest auditory sensitivity appear to correspond to the resonant frequencies of the different components of the middle ear transformer and cavity.     

Data-driven habitat analysis of the Ponto-Caspian amphipod Dikerogammarus villosus in two invaded regions in Europe

When forecasting invasions, models built on a dataset from a certain region often have to be used for simulations in another geographic region. Results on the reliability and usefulness of such models are missing in literature. The present study compares habitat suitability models for the invasive amphipod species Dikerogammarus villosus developed based on data gathered in recently invaded rivers and channels in Flanders (Belgium), with similar models developed on the basis of long-term colonised systems in Croatia. The models were tested on their reliability in both regions. Two techniques, logistic regressions (LR) and classification trees (CT) were used to analyse the habitat preference of this species based on physical–chemical and morphological habitat characteristics. It was found that in Flanders, D. villosus prefers rivers with a non-natural bank structure, high oxygen saturation, low conductivity and good chemical water quality, which could be related to its distribution in large rivers and canals. In Croatian rivers, high oxygen saturation was the most important prerequisite for the species to be present. Despite the longer history of invasion in Croatia, the species seemed to have similar habitat preferences in both invaded regions. Both data-driven approaches yielded similar results, but CT performed somewhat better based on the used performance criteria (% Correctly Classified Instances, Kappa and Area Under Curve) and were easier to interpret compared to the LR. The CT models developed based on the data of Flanders performed moderately when applying on the data of Croatia, but had a lower performance when applied vice versa. The LR models did not perform well when applying on a dataset of another geographic area. Extrapolation of the logistic regression model seemed to be more difficult compared to classification tree models. Our results indicate that it is possible to determine the habitat preference of an invasive species and that these models could be applied to other regions in Europe in order to take preventive measures to control the further spread of invasive species. However, a major concern is that the models are developed based on a representative range of all relevant variables reflecting the stream conditions and that accurate data are important.     

Detecting adaptive evolution based on association with ecological gradients: Orientation matters!

Population genetic signatures of local adaptation are frequently investigated by identifying loci with allele frequencies that exhibit high correlation with ecological variables. One difficulty with this approach is that ecological associations might be confounded by geographic variation at selectively neutral loci. Here, we consider populations that underwent spatial expansion from their original range, and for which geographical variation of adaptive allele frequency coincides with habitat gradients. Using range expansion simulations, we asked whether our ability to detect genomic regions involved in adaptation could be impacted by the orientation of the ecological gradients. For three ecological association methods tested, we found, counter-intuitively, fewer false-positive associations when ecological gradients aligned along the main axis of expansion than when they aligned along any other direction. This result has important consequences for the analysis of genomic data under non-equilibrium population genetic models. Alignment of gradients with expansion axes is likely to be common in scenarios in which expanding species track their ecological niche during climate change while adapting to changing environments at their rear edge.     

An application of fuzzy set ordination to determine tree habitat suitability of sites from a regional data set

Abstract. Fuzzy set ordination is employed to evaluate sites on the basis of their suitability for particular tree species. The technique orders sites along an axis defined by the presences and absences of a given species of interest. A rationale is given in terms of noise reduction; in many situations the overall vegetation of a site will reflect habitat conditions better than the presence, absence, or quantitative performance of any single species. A data set of tree presence/absence covering a large part of the southeastern United States was analyzed and habitat suitability scores were calculated for each species. Monte-Carlo tests were used to measure the statistical power of the data set with regard to habitat preferences; 38 of the 49 species have cumulative frequency distributions showing significant departures from random expectation. Most statistically significant habitat preferences seem to be related to geographic range limits located within the study area, but some species found throughout the area also show significant departures from random expectation. The method may find applications in autecological studies of species, selection of representative site conditions for simulation modeling, and the solution of certain technical problems in ordination.     

Taxon-specific physico-chemical change points for stream benthic invertebrates

Freshwater organisms face numerous stressors, such as nutrient enrichment, contaminant pollution, sedimentation and alterations in stream hydrology and habitat structure. One of the most significant and widespread stressors in European freshwaters is expected to be water pollution from intensive land use. However, the information on critical threshold concentrations at which taxa decline or increase in frequency and abundance is missing for the large majority of river benthic invertebrate taxa. The main aim was to determine ecological change points for benthic invertebrate taxa at which rapid alterations in species frequency and abundance occur as a consequence of relatively small changes in the environmental gradient. These change points can be interpreted as critical threshold concentrations. A total of 468 river benthic invertebrate taxa and nine physico-chemical variables describing the daytime dissolved oxygen, chloride, nutrient concentrations and organic load were analyzed. We selected 751 river sites from a nationwide range of locations in Germany for this investigation. Depending on the physico-chemical variable, between 20.6% and 48.8% of the total number of tested taxa were assigned with a valid change point. All taxa were assigned to negative or positive response groups depending on the response direction. Except for daytime dissolved oxygen, negative responding taxa are referred to as decreasers and positive responding taxa as increasers, respectively. In total, 25.8–100% of the decreasers’ change points were below (and above in the case of daytime dissolved oxygen) the background values defined as quality criteria for German rivers by the water authorities. This indicates that stricter quality criteria may need to be set to reach the good ecological status according to the European Water Framework Directive. The calculated daytime dissolved oxygen change points were essentially in line with the species saprobic values and taxon-specific change points for physico-chemical variables fit well with the data provided in other international studies. We deliver valuable knowledge about the sensitivities and response schemes of river benthic invertebrate species. This information is especially useful for the development of efficient management and policy tools to predict the likelihood of occurrence of individual species under different levels of anthropogenic impact.     

Analysis of multilocus fingerprinting data sets containing missing data

Missing data are commonly encountered using multilocus, fragment‐based (dominant) fingerprinting methods, such as random amplified polymorphic DNA (RAPD) or amplified fragment length polymorphism (AFLP). Data sets containing missing data have been analysed by eliminating those bands or samples with missing data, assigning values to missing data or ignoring the problem. Here, we present a method that uses random assignments of band presence–absence to the missing data, implemented by the computer program famd (available from http://homepage.univie.ac.at/philipp.maria.schlueter/famd.html ), for analyses based on pairwise similarity and Shannon's index. When missing values group in a data set, sample or band elimination is likely to be the most appropriate action. However, when missing values are scattered across the data set, minimum, maximum and average similarity coefficients are a simple means of visualizing the effects of missing data on tree structure. Our approach indicates the range of values that a data set containing missing data points might generate, and forces the investigator to consider the effects of missing values on data interpretation.     

Functional Group Density as an index for assessing habitat quality in tallgrass prairie

We propose that patterns of plant functional group occurrences could be a reliable indicator of prairie vegetation quality. A method for assessing tallgrass prairie quality based on density and composition of plant functional groups was developed and tested by comparison with qualitative indices calculated from species data at 17 prairies in Illinois. Species sample data were recorded from quadrats while functional group data were recorded from segments of belt transects overlying the species sample transects. Prairies selected include remnants and restorations and represent a wide range of habitat quality including recognized natural areas, degraded remnants, and prairie plantings of varying age and success. For agglomerative clustering of prairie quality classes, a matrix of habitat indices and metrics was used based on species sample data from all sites. Three groups were identified in cluster analysis that were characteristic of high, medium, and low-quality prairies. Mean Functional Group Density (Mean FGD), an index developed based on the mean products of frequency and density among plant functional groups recorded from belt transects at each site, had the highest correlations to habitat quality indices among two other functional group indices tested. Mean FGD was highly correlated with the mean coefficient of conservatism and floristic quality index, indices calculated from species sample data that have been shown to be reliable indicators of habitat quality. In means comparison tests among prairie quality classes, Mean FGD differentiated high-quality from medium and low-quality prairies, but did not distinguish medium from low-quality sites (only two low-quality sites were identified), although rank order of Mean FGD was as predicted. There is a tradeoff in efficiency and precision between species-level and functional group sample data. Species-level data more precisely discriminate differences in low and medium-quality sites; however, functional group sampling is much more rapid requiring only 20-to-25% of the time required for collecting species-level data. Results from functional group sampling highlight differences in functional group composition among prairie quality classes. High-quality prairies are characterized by greater abundance of sedges and hemi-parasites while lower-quality prairies were affiliated more with non-native perennial forbs and annual/biennial species.     

Lineage Range Estimation Method Reveals Fine-Scale Endemism Linked to Pleistocene Stability in Australian Rainforest Herpetofauna

Areas of suitable habitat for species and communities have arisen, shifted, and disappeared with Pleistocene climate cycles, and through this shifting landscape, current biodiversity has found paths to the present. Evolutionary refugia, areas of relative habitat stability in this shifting landscape, support persistence of lineages through time, and are thus crucial to the accumulation and maintenance of biodiversity. Areas of endemism are indicative of refugial areas where diversity has persisted, and endemism of intraspecific lineages in particular is strongly associated with late-Pleistocene habitat stability. However, it remains a challenge to consistently estimate the geographic ranges of intraspecific lineages and thus infer phylogeographic endemism, because spatial sampling for genetic analyses is typically sparse relative to species records. We present a novel technique to model the geographic distribution of intraspecific lineages, which is informed by the ecological niche of a species and known locations of its constituent lineages. Our approach allows for the effects of isolation by unsuitable habitat, and captures uncertainty in the extent of lineage ranges. Applying this method to the arc of rainforest areas spanning 3500 km in eastern Australia, we estimated lineage endemism for 53 species of rainforest dependent herpetofauna with available phylogeographic data. We related endemism to the stability of rainforest habitat over the past 120,000 years and identified distinct concentrations of lineage endemism that can be considered putative refugia. These areas of lineage endemism are strongly related to historical stability of rainforest habitat, after controlling for the effects of current environment. In fact, a dynamic stability model that allows movement to track suitable habitat over time was the most important factor in explaining current patterns of endemism. The techniques presented here provide an objective, practical method for estimating geographic ranges below the species level, and including them in spatial analyses of biodiversity.     

Patterns in the geographic ranges of the world's woodpeckers

We used data on the world's woodpeckers to test for patterns in the geographic distributions of a single group of closely related species. The frequency distribution of woodpecker geographic range sizes is approximately lognormal. Most variation in range sizes is explained by differences between species within genera; that is, range size seems to be an evolutionarily labile trait. The largest woodpecker ranges are found in Eurasia, both when absolute differences are compared and when range size is measured as a proportion of estimated available habitat. Notably, there is a negative relationship between the mean range sizes attained by species in a genus or tribe in South America and the mean ranges attained by species in the same tribe or genus in North America. Large-bodied species tend to be more widely distributed and to live at higher latitudes, but both tendencies disappear if the taxonomic relatedness of species is controlled for. Species living at high latitudes also tend to be more widely distributed. This relationship seems largely due to the effect of North American woodpeckers, which show it even when the taxonomic relatedness of species is controlled. Small continents generally have more woodpecker species than do large ones. Woodpecker geographic range sizes are smaller the more woodpecker species inhabit an area. Species show less overlap in their geographic ranges with species of similar than with species of dissimilar body size. The implications of these results for our understanding of patterns in geographic range sizes are discussed.     

不同鸟类生境网络复合与优化——以苏锡常地区白鹭、鸳鸯、雉鸡为例

通过复合优化不同焦点物种的生境网络表征多数物种生境网络实现物种多样性保护,具有重要理论价值和实践意义。以苏锡常地区为研究区域,以白鹭、鸳鸯和雉鸡3种鸟类为焦点物种,通过采用生境斑块约束条件模型与最小成本路径法相结合的方法识别出生境网络;从满足集合覆盖问题视角出发,通过采用网络构成要素复合的方法进行网络复合,其中廊道与廊道的复合通过采用评价生态系统服务价值的方法实现;通过采用观测数据为生境斑块和迁移廊道增补依据的方法进行复合生境网络优化。结果表明:生境网络经过复合,廊道总面积、总长度及生态系统服务价值在分别减少了16%、68%和10%的情况下,保持了斑块间的连接度、兼顾到经济效益和生态效益最大化、覆盖了占总数86%以上的75处观测点;复合生境网络经过优化,覆盖了2010年及以后所有84处观测点,廊道总面积、总长度及生态系统服务价值分别增加了19%、21%和27%;说明生境网络复合优化定量分析方法和思路具有可行性,为网络复合研究提供了借鉴。     

Correlates of rediscovery and the detectability of extinction in mammals

Extinction is difficult to detect, even in well-known taxa such as mammals. Species with long gaps in their sighting records, which might be considered possibly extinct, are often rediscovered. We used data on rediscovery rates of missing mammals to test whether extinction from different causes is equally detectable and to find which traits affect the probability of rediscovery. We find that species affected by habitat loss were much more likely to be misclassified as extinct or to remain missing than those affected by introduced predators and diseases, or overkill, unless they had very restricted distributions. We conclude that extinctions owing to habitat loss are most difficult to detect; hence, impacts of habitat loss on extinction have probably been overestimated, especially relative to introduced species. It is most likely that the highest rates of rediscovery will come from searching for species that have gone missing during the 20th century and have relatively large ranges threatened by habitat loss, rather than from additional effort focused on charismatic missing species.     

Plant communities of New Brunswick in relation to environmental variation

Abstract. The objective of this study was to quantitatively describe vegetation-environment relationships at a regional scale within the Province of New Brunswick, Canada, using vegetation and environment data from 3947 provincial forestry sample plots. The major plant community types in the province were identified using cluster analysis. Relationships of these communities to climate, topography and soil variables were analyzed by Canonical Correspondence Analysis (CCA), using both a reduced data set consisting of cluster likelihood scores × sample plots and an unreduced species × sample plots data matrix. The vegetation types and major axes of environmental variation were mapped to examine the geographic distributions of these factors within the province. Eight communities were identified and described in terms of enhanced/reduced species (significantly higher or lower frequencies of occurrence in a specific community type relative to all plots) and common species (species in the community type with the highest frequencies of occurrence). The canonical axes explained 25 % of the variation in the vegetation cluster data. Vegetation composition was related to three major environmental gradients representing climate and elevation, soil moisture, and soil fertility. The geographic distributions of vegetation communities exhibited predictable but weak correspondence to the geographic distributions of individual environmental factors. Our findings emphasize the overriding importance of climate and topography and the secondary importance of soil factors in controlling vegetation pattern at the regional scale.     

Population Structure of the Giant Garter Snake, Thamnophis gigas

The giant garter snake, Thamnophis gigas, is a threatened species endemic to California’s Central Valley. We tested the hypothesis that current watershed boundaries have caused genetic differentiation among populations of T. gigas. We sampled 14 populations throughout the current geographic range of T. gigas and amplified 859 bp from the mitochondrial gene ND4 and one nuclear microsatellite locus. DNA sequence variation from the mitochondrial gene indicates there is some genetic structuring of the populations, with high F_ST values and unique haplotypes occurring at high frequency in several populations. We found that clustering populations by watershed boundary results in significant between-region genetic variance for mtDNA. However, analysis of allele frequencies at the microsatellite locus NSU3 reveals very low F_ST values and little between-region variation in allele frequencies. The discordance found between mitochondrial and microsatellite data may be explained by aspects of molecular evolution and/or T. gigas life history characteristics. Differences in effective population size between mitochondrial and nuclear DNA, or male-biased gene flow, result in a lower migration rate of mitochondrial haplotypes relative to nuclear alleles. However, we cannot exclude homoplasy as one explanation for homogeneity found for the single microsatellite locus. The mitochondrial nucleotide sequence data supports conservation practices that identify separate management units for T. gigas.     

Historical and recent processes shaping the geographic range of a rocky intertidal gastropod: phylogeography,ecology, and habitat availability

Factors shaping the geographic range of a species can be identified when phylogeographic patterns are combined with data on contemporary and historical geographic distribution, range‐wide abundance, habitat/food availability, and through comparisons with codistributed taxa. Here, we evaluate range dynamism and phylogeography of the rocky intertidal gastropod Mexacanthina lugubris lugubris across its geographic range – the Pacific coast of the Baja peninsula and southern California. We sequenced mitochondrial DNA (CO1) from ten populations and compliment these data with museum records, habitat availability and range‐wide field surveys of the distribution and abundance of M. l. lugubris and its primary prey (the barnacle Chthamalus fissus). The geographic range of M. l. lugubris can be characterized by three different events in its history: an old sundering in the mid‐peninsular region of Baja (~ 417,000 years ago) and more recent northern range expansion and southern range contraction. The mid‐peninsular break is shared with many terrestrial and marine species, although M. l. lugubris represents the first mollusc to show it. This common break is often attributed to a hypothesized ancient seaway bisecting the peninsula, but for M. l. lugubris it may result from large habitat gaps in the southern clade. Northern clade populations, particularly near the historical northern limit (prior to the 1970s), have high local abundances and reside in a region with plentiful food and habitat – which makes its northern range conducive to expansion. The observed southern range contraction may result from the opposite scenario, with little food or habitat nearby. Our study highlights the importance of taking an integrative approach to understanding the processes that shape the geographic range of a species via combining range‐wide phylogeography data with temporal geographic distributions and spatial patterns of habitat/food availability.     

皖南凹耳臭蛙地理分布格局与生境选择

凹耳臭蛙(Odorrana tormota)是一种小型山溪蛙类,由于叫声独特,种内具有超声通讯功能而引起国内外相关专家的广泛关注。首次针对皖南凹耳臭蛙地理分布、生境选择以及人为干扰强度对种群密度的影响进行了研究。共调查了9县1区32条不同的山溪,发现其中的16个有凹耳臭蛙分布,除黄山区的香溪和浮溪具文献记录外,其余14个分布点均为首次报道。基于现存凹耳臭蛙分布点的水系关系,推测该物种在皖南的分布呈现出以青弋江、新安江和水阳江三大水系为纽带、间断分布的格局。对其中10个分布点222只凹耳臭蛙生境选择的统计结果表明,在夜晚活动周期内,该物种对乔木、灌木丛、草本植物、岩石、沙滩和倒伏朽木等不同生境的选择具有显著差异(F5,54=25.75,P0.001),栖息在草本植物和灌丛中蛙的数量明显多于栖息在其它几种生境类型中的蛙类。暗示溪流沿岸的草本植物和灌丛对该物种的生存和繁衍具有重要意义,乔木对其非活动周期的隐蔽可能具有不可替代的作用。根据人类活动对其栖息地破坏程度,将人为干扰程度划分为强、中和弱3种不同强度,对不同干扰强度下种群相对密度进行单因素方差分析。结果表明,中等干扰强度下的种群密度显著高于其它两种强度下的种群密度(F2,13=8.155,P=0.005)。推测适度的人为干扰可能增加了栖息地内昆虫的种类和数量,对蛙类的生存和繁衍更为有利。但随着人为干扰强度的增强,对溪流沿岸植被的破坏更为剧烈,种群密度显著下降。