Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2013; 8(4): e59227.
Published online 2013 April 1. doi:  10.1371/journal.pone.0059227
PMCID: PMC3613380

Quantifying High Resolution Transitional Breaks in Plant and Mammal Distributions at Regional Extent and Their Association with Climate, Topography and Geology

David L. Roberts, Editor



We quantify spatial turnover in communities of 1939 plant and 59 mammal species at 2.5 km resolution across a topographically heterogeneous region in south-eastern Australia to identify distributional breaks and low turnover zones where multiple species distributions overlap. Environmental turnover is measured to determine how climate, topography and geology influence biotic turnover differently across a variety of biogeographic breaks and overlaps. We identify the genera driving turnover and confirm the versatility of this approach across spatial scales and locations.


Directional moving window analyses, rotated through 360°, were used to measure spatial turnover variation in different directions between gridded cells containing georeferenced plant and mammal occurrences and environmental variables. Generalised linear models were used to compare taxic turnover results with equivalent analyses for geology, regolith weathering, elevation, slope, solar radiation, annual precipitation and annual mean temperature, both uniformly across the entire study area and by stratifying it into zones of high and low turnover. Identified breaks and transitions were compared to a conservation bioregionalisation framework widely used in Australia.


Detailed delineations of plant and mammal turnover zones with gradational boundaries denoted subtle variation in species assemblages. Turnover patterns often diverged from bioregion boundaries, though plant turnover adhered most closely. A prominent break zone contained either comparable or greater numbers of unique genera than adjacent overlaps, but these were concentrated in a small subsection relatively under-protected by conservation reserves. The environmental correlates of biotic turnover varied for different turnover zones in different subsections of the study area. Topography and temperature showed much stronger relationships with plant turnover in a topographically complex overlap, relative to a lowland overlap where weathering was most predictive. This method can quantify transitional turnover patterns from small to broad extents, at different resolutions for any location, and complements broad-scale bioregionalisation schemes in conservation planning.


Biogeographic breaks and transition zones have received considerable investigation at global and continental scales at coarse resolutions (e.g. [1][4]). Their detailed characterisation for distributions of multiple plant and animal species at smaller extents is less common. Some recent studies have adopted a narrower geographic focus, for instance Thiel-Egenter et al. [5] used a cluster analysis of plant species across the European Alps to assess the robustness of break zones, but this was at a comparatively coarse resolution.

It is increasingly pertinent to examine biogeographic phenomena at finer resolutions and at regional to local scales. This tests the generality of concepts derived by investigations over broad extents, as patterns and processes may not exhibit scale invariance [6], [7]. Practical concerns also advocate for translating insights achieved at broad scales to smaller extents. For instance, climatic effects vary over relatively small distances under topographic influence [8] with implications for how taxa respond to these changes. The effects of global processes such as climate change will ultimately need to be managed at regional and local scales [9], so it is important to understand how fine-scale variation in landscape features mediates climatic changes, particularly in areas of rich biodiversity.

Fine-scale quantification of biogeographic breaks also enables depiction of transitional zones in greater detail, facilitating their comparison with subtle variation in environmental gradients and potentially clarifying how different physical factors regulate species distributions across break transitions. This is useful for conservation planning because taxa occurring in transitional areas are often the most adaptive to a range of environmental conditions [10], [11]. Identifying continuous, fine-scale biogeographic distributional breaks may also promote design of efficient reserve networks. For instance, Shriner et al. [12] found that the numbers of species represented in equal area reserve networks increased as the grain size of distributional data decreased, such that best representation networks developed at coarser grains were almost 1000 times larger than equivalent networks using fine grained data. Fine resolution data may also be more informative to conservation planning in heterogeneous and fragmented landscapes [13] and in areas of high biological diversity [14], [15].

Quantifying the rate of species turnover across a landscape has been used to identify transitional turnover patterns and biogeographic breaks for bees [1] and birds [16]. Species turnover is the rate of compositional change in biotic communities across geographic space, i.e. the replacement of certain species with others among areas of any given size [17], [18]. High turnover areas experience greater rates of replacement between different species, and consequently show lower beta diversity and species richness, manifesting as a ‘break’ in distributions. Low turnover regions experience lower rates of community change and comprise more stable communities where the distributions of multiple different species overlap. Relating gridded maps of biotic turnover with corresponding environmental turnover can determine the relative influence of different physical factors on species turnover, and thus break formation. This was the general approach used by Di Virgilio et al. [19] to quantify breaks in plant, mammal and reptile distributions at 5 km resolution and their correspondence with geological and elevational turnover over a large portion of south-eastern Australia.

This research maps breaks and transitions in multiple plant and mammal species distributions at 2.5 km resolution across the south-east coast of New South Wales (NSW) in Australia, which is a smaller subsection of the same region previously examined by Di Virgilio et al. [19]. We provide a more detailed quantification of gradational breaks and overlaps in this region than was achieved by Di Virgilio et al. [19]. Break and overlap zones are also characterised more fully by identifying the specific flora and fauna exclusive to each zone and that may contribute to driving turnover. In changing the extent and resolution of these analyses, we also assess the versatility of this approach across multiple scales and resolutions and its potential for investigating transitional breaks at any location globally.

Secondly, Di Virgilio et al. [19] only examined inter-relationships between biotic and environmental factors uniformly across south-eastern NSW. Instead, we investigate whether the environmental influence on species turnover and break formation varies in different subsections of the study area, specifically in different geographic areas characterised by either high or low turnover. The environmental drivers considered are expanded to include climate (annual precipitation, annual mean temperature), topography (slope, solar radiation) and regolith weathering intensity, in addition to elevation and lithology. Some studies have identified climate as opposed to topographic variation as the principal driver of beta diversity in forest trees at a regional scale (e.g. [20]), whereas others have found a joint effect of climate and topography in regulating inter-annual vegetation variability [21]. We assess how the relative influence of environmental drivers varies for the different break and overlap zones.

Previous studies of trees in south-eastern NSW have shown the importance of variables such as climate, lithology and topographic position in predicting the richness of several species groups, e.g. [22], and how the environmental niche of eucalypt species varies along gradients of similar variables [23], [24]. The focus of our study is different in that it explicitly identifies transitional biogeographic breaks and overlaps for a greater variety of plant species and the role of spatial turnover of physical variables in driving these specific patterns.

Our results are interpreted in reference to the Interim Biogeographic Regionalisation for Australia version 6.1 (IBRA 6.1) [25], which is an existing bioregionalisation scheme widely used for conservation and reserve network planning in Australia [26], [27].

Materials and Methods

Study Area and Data Sources

The South East Corner (SEC) bioregion defined in IBRA 6.1 [25] straddles south-eastern New South Wales (NSW) and Victoria (VIC) in south-eastern Australia (Fig. 1). SEC in NSW comprises two parts: a montane region in the south-west and a larger area bordered by the Pacific Ocean, the South-Eastern Highlands (SEH) and Sydney Basin (SB) IBRA bioregions. We focused on SEC’s coastal portion in NSW (10,500 km2 in area) which is bordered by Cape Howe (149.974°, −37.494°) to the south and Ulladulla (150.467°, −35.358°) near its northern extent. SEC is further subdivided into three IBRA subregions.

Figure 1
Location of the study area in south-eastern New South Wales, south-eastern Australia.

SEC has a mesothermal climate, with long mild summers and little variation in annual rainfall. Average monthly temperatures range from minima of −3.5–8.4°C to maxima of 19.2–28.8°C. Average monthly rainfall minima and maxima are 29–102 mm and 58–155 mm respectively. SEC’s terrain is heterogeneous, comprising coastal lowlands extending ~15 km inland to the southern highlands to the west which rise to ~1250 m.

A 0.2° (~22 km) buffer was added around SEC (dashed white line in Fig. 1) to minimise edge effects, as data just outside the SEC boundary would influence the compositional dissimilarities between species and environmental values some distance inside the border. The buffer did not extend into Victoria because the biotic data had been sampled mainly within NSW. Biotic and environmental data contained within this buffer are included in moving window analyses, but excluded from generalised linear model (GLM) analyses. The buffer increased the study area to 17,240 km2.

The geographic locations of plants and mammals sampled within SEC plus the buffer were derived from the NSW Office of Environment and Heritage’s (OEH) Atlas of NSW Wildlife [28] (see This database records the geographic location of flora and fauna sightings in NSW from 1788 to 2009. Some records showed multiple observations of the same species at the same location and time; these were treated as a single point. Exotic species were excluded and records identified as potentially ambiguous e.g. with duplicate names, subspecies where a parent species was present etc., were omitted. The remaining plant genera were checked for taxonomic accuracy, with no errors found. Marine mammals (order Cetacea and families Otariidae and Phocidae) were excluded because they are restricted to the coastal zone.

The plant database contained 104,030 records, comprising 145 families, 595 genera and 1939 species and included a wide variety of different trees, shrubs, ferns, tussocks, herbs, grasses and other flowering plants. There were 41,392 mammal records (38 genera, 59 species). The mean positional accuracy of the unaltered Atlas mammal data is 1265 m with a standard deviation of 3391 m. Plant observations have a mean accuracy of 1137 m and standard deviation of 3451 m.

Elevation data was derived from the three arc-second Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) [29]. Slope and radiation data sets were derived from this DEM. Geological data were derived from the 1[ratio]1,000,000 scale surface geology of Australia GIS dataset [30]. Lithological classes were simplified using the same recoding scheme as Di Virgilio et al. [19]. Annual precipitation and annual mean temperature at 30 arc-seconds (~1 km) resolution were derived from the Worldclim data sets generated by Hijmans et al. [31]. Weathering intensity data were derived from the 100 m resolution weathering intensity index for Australia [32] which reflects changes to the physical characteristics of regolith as bedrock is altered into secondary minerals.

Statistical Analyses

The spatial analyses adapted those described by Di Virgilio et al. [19] and used the Biodiverse software, version 0.16 [33]. The overall approach was to quantify biotic and environmental turnover by progressively iterating directional moving windows throughout the study area to measure the dissimilarity between two sets of cells, where each cell recorded the number of unique occurrences of species or environmental values at its geographic location. Window dimensions were calibrated using correlograms of biotic turnover versus distance. Correlograms are widely used to investigate variation in natural phenomena across geographic space, for example genetic structure [34] and plant richness gradients [35].

All biotic and abiotic data were aggregated to 0.0225°×0.0225° cells (~2.5 km×2.5 km) by importing each data set into Biodiverse separately. All cells across the study area were used for the analyses, including those in which no species were recorded.

The Sørensen index was used to quantify the spatial turnover of biotic and lithological data:

equation image

where A is the number of unique occurrences of species or lithological classes found in both samples (i.e. both ‘neighbour sets’ of cells in the moving window); B is the number of species/classes unique to sample set one; and C is the number of species/classes unique to sample set two. This index produces a value between 0 (samples are identical) and 1 (samples share no species/classes). Alternative indices (Bray Curtis, Simpson’s Beta and Jaccard) were tested with the biotic data. The Jaccard index produced results similar to those of the Sørensen index, owing to the similarity in their calculation. Simpson’s Beta and Bray Curtis produced turnover patterns that shared some similarities with Sørensen-derived patterns, but were less detailed.

To investigate which genera are exclusive to discrete areas of high or low turnover, we identified the plant and mammal genera unique to each sample set and common to both. Once gradational areas of high and low plant turnover were identified via the moving window analyses, the outline of a discrete, high turnover ‘break’ zone was assigned to the sample one set and the perimeter of one of two discrete low turnover ‘overlap’ zones was assigned to the sample two set. The plant or mammal genera common to both sets and those exclusive to either sample set one or sample set two were then identified by comparing the genera in the break with those in one of the overlap areas.

Since the topographic, climatic and weathering variables are interval scaled data, the moving window analyses used a numeric dissimilarity index (eq. 2), which compared the set of numeric values in each sample set by calculating the mean absolute difference between all pairs of values in samples 1 and 2

equation image

where L1 and L2 are the sets of labels (values) in samples 1 and 2, respectively, and n1 and n2 are the overall abundances of labels in sample sets 1 and 2. w1i and w2j are the abundances of labels l1i and l2j in samples 1 and 2, respectively. This metric results in values ranging from zero (samples are identical) progressing theoretically to infinity as average differences increase.

We also calculated species richness and sample redundancy (eq. 3) for plants and mammals for each moving window. Richness is an important component of species turnover [36]. Sample redundancy [37] is a measure of sampling effort across SEC, and therefore enables reliability assessment (i.e. redundancy patterns can be reconciled against species turnover patterns in order to gauge the effects of poor sampling), and was calculated as:

equation image

where SR is the species richness and Ns is the number of observations of each species in the relevant sample set. This equation produces values from 1 in cases where there are many observations for each species and indicating areas that are well sampled, to zero, meaning there is no redundancy in the sampling because there is only one observation per species.

Moving Window Analyses

The mosaic of cells comprising a moving window were divided into one of two orthogonal, semi-elliptical neighbourhood sets (the sample one and two sets described above) and used the same geometric specification as the moving windows described by Di Virgilio et al. [19]. However, in this study the neighbourhood sets used different dimensions, calibrated via correlograms using the 0.0225° resolution biotic data.

The same window sizes were used for flora and mammals to aid comparison of turnover in each group. Since we aim to find the environmental correlates of species turnover, the same dimensions were used for moving window analyses of the environmental variables.

To gauge how biotic and abiotic turnover varies in different directions, moving window analyses were repeated in 25 increments of 15°; hence the window was rotated through 24 different angles from 0° to 360°.

Generalised Linear Model Regression

We investigated how environmental variables influenced species turnover across the whole study area and also how their relative influence on break and overlap formation varied once the study area was stratified into three different areas of high or low turnover (see Results for details of the zones/stratification). Generalized linear models (GLM; [38]) were used to regress untransformed plant or mammal turnover on the separate environmental turnover predictors generated by each of the 25 moving window orientations for the complete and stratified study area. A Gaussian error distribution for biotic turnover was assumed. We opted for univariate models because a multivariate approach may not capture important covariates.



The correlograms revealed a range of ~0.2° for flora, with a sill value of ~0.92, where the range is the distance on the x-axis at which the median turnover values begin to plateau and is the maximum distance to which there is spatially structured turnover (Fig. 2A). The mammal data had a range of ~0.25° and a sill value of ~0.8 (Fig. 2B). Therefore, the moving window outer ellipse had major and minor radii of 0.2° and 0.085° respectively. The inner ellipse major and minor radii were 0.085° and 0.0075°.

Figure 2
Species group correlograms showing the distribution of turnover values between pairs of cells over increasing geographic distance.

Species Richness

Plant richness was high across much of the SEC bioregion (Fig. 3, panels A–B), except for a large area of low richness in the north. Patches of higher plant richness tended not to extend across IBRA boundaries, e.g. the East Gippsland Lowlands (SEC 1). Most high richness areas were on low, uniform terrain, especially along the coast, except for a patch of high richness in the southwest.

Figure 3
Plant and mammal richness and sample redundancy patterns.

A greater proportion of SEC showed high mammal species richness (Fig. 3, panels E–F) with only small patches of low richness restricted to higher elevations. Nonetheless, some areas of high mammal richness corresponded with distinct regions of high plant richness e.g. the coastal lowlands from southern Bateman (SEC 3) to SEC’s southern border.

Richness patterns for both biotic groups changed slightly according to the orientation of the moving window because different combinations of cells are considered by different window orientations. For example, the shape of the northern low plant and mammal richness zones changed from the 90° to 0° windows, but their overall form was still consistent across orientations.

Sample Redundancy

Floristic sampling redundancy (Fig. 3, panels C–D) generally resembled plant richness patterns, with an area of medium to lower redundancy roughly corresponding to the northern low richness zone. High sampling redundancy extended from the southern Bateman bioregion to SEC’s southern border, interspersed with only a few patches of low redundancy. In comparison to plants, mammal redundancy patterns showed less variation with higher redundancy across a larger portion of SEC.

Biotic Turnover

Figure 4 shows the biotic turnover results for moving window analyses oriented at 90°. White areas on the maps are where either one or both of the moving window neighbourhoods contain no observations.

Figure 4
Plant and mammal species turnover maps showing location of turnover zones.


There was continuous variation in species turnover patterns with several breaks, overlaps and transition zones shown in greater detail than the 5 km resolution analyses of Di Virgilio et al. [19]. A prominent area of high plant turnover spanned northern SEC, with a strip of high turnover extending southward. The main body of this ‘northern break zone’ straddled a topographically complex area to the west, but flatter terrain along the coast. The southward extension corresponds to a low, flat area with undulating terrain along its southern boundary. Much of this break shows comparatively uniform turnover, although there is greater variation along its edges where there are transitions from high to low turnover in the south-east and south-west. Although Di Virgilio et al. [19] identified the location of this floristic break zone, our results show this break in much greater detail, with the presence of the southern extension and transitions along the zone boundaries more apparent.

A western low turnover area to the south-west of the break overlies high, topographically complex terrain separated from the break by a small valley (Fig. 1). In contrast, there is an eastern low turnover zone situated on the coastal plain. The breaks and overlaps varied in shape with moving window orientation, but again their overall form remained consistent (see Figures S1, S2, S3).

The outlines of the northern break zone, the hilly western overlap and the coastal eastern overlap are sketched in Fig. 5 and form the basis of the stratification of SEC into discrete turnover zones and the identification of genera exclusive to each zone. Although other discrete low turnover zones were present in SEC, these two overlaps were specifically selected because they flank the northern break on landscapes of varying topographic complexity.

Figure 5
Discrete high and low plant turnover zones in the study area.


Mammal turnover was lower across SEC and less sharply defined than plant turnover (Fig. 4), but there were some parallels with plant turnover patterns. Parts of the northern floristic break were evident, such as the western and north-eastern edges and southward extension, although all were fainter and smoother. Most of the coastline showed low mammal turnover, whereas high plant turnover extended from northern SEC to southern Bateman. Overall, mammal turnover adheres less to IBRA boundaries than plants, although some mammal turnover patterns do not extend across the boundaries such as the high turnover that aligns with the East Gippsland Lowlands subregion.

Biotic turnover was not as smooth as the corresponding richness and redundancy patterns. Breaks in mammal distributions appeared in both high and low redundancy areas and across extremes of richness. In contrast, the main portion of the floristic high turnover zone corresponds to areas of low to medium sampling redundancy and richness.

Genera driving high turnover

The northern break zone and the smaller eastern overlap contained similar numbers of plant species and genera, although eight times as many observations were sampled in the overlap (Table 1). 79 plant genera occurred exclusively within each zone, with 270 genera common to both regions (Table 2). Most of the 270 common genera are widespread throughout SEC. However, of the 79 genera exclusive to the break, 14 occur only in the break and nowhere else within SEC, while 45 genera are isolated from other populations elsewhere within SEC. The eastern overlap shows a similar trend: 11 genera occur nowhere else and 35 are isolated subpopulations.

Table 1
Total number of plant and mammal observations, genera and species sampled within the northern biogeographic break zone, eastern overlap and western overlap within the South East Corner (SEC) study area in New South Wales, Australia.
Table 2
The number of plant and mammal genera that occur exclusively either in the northern biogeographic break zone, or eastern overlap, and common to both areas, compared with the number of plant and mammal genera that occur exclusively either in the northern ...

A greater variety of plant genera and species were sampled in the northern break relative to the western overlap (Table 1). There were 135 genera unique to the break, 57 unique occurrences in the western overlap and 214 common genera (Table 2). The distribution of unique/common genera showed a similar pattern to that noted above. Plant genera identified by all comparisons are listed in Tables S1–S2 in File S1.

Similar numbers of mammal observations, genera and species were sampled in the break and eastern overlap (Table 1). There were two genera unique to the break, one unique genus in the overlap and 32 genera common to both areas (Table 2). Comparison of the break and western overlap revealed 8 genera unique to the former, whilst the overlap contained no exclusive genera and 26 genera which occurred in both areas (Table 2). In contrast to plants, mammal genera unique to the break were commonplace throughout SEC. See Tables S3–S4 in File S1 for the mammal genera identified.

Abiotic Turnover

Generally, environmental turnover did not correspond to the prominent biotic turnover zones (Fig. 6), though there were some fine-scale parallels. For instance, slope and solar radiation show a southern transition of low turnover that corresponds to the southward extension of the northern floristic break. Nonetheless, concordance between biotic and abiotic turnover patterns was inconsistent. The western plant overlap is on high elevation, topographically complex terrain showing higher elevation turnover, whilst the opposite applies to the eastern overlap. Elevation, slope, radiation and annual precipitation show substantial patches of higher turnover in north-western SEC, corresponding to part of the biotic northern break, whereas weathering turnover is low in this same area. With the exception of geology and weathering, abiotic turnover is low along the majority of the coastal lowlands.

Figure 6
Plant species turnover compared to environmental turnover.

Generalised Linear Models of Biotic and Abiotic Turnover

Turnover across the complete study area

Topographic and temperature turnover was substantially more predictive of variation in mammal turnover than for plants, as shown by the maximal R2 coefficients in Table 3. In contrast, weathering and geological turnover were poor predictors of both mammal and plant turnover, explaining only a small proportion of variation in biotic turnover. Biotic-environmental turnover relationships were positive across SEC, except for weathering turnover. All maximal R2 coefficients were significant at P<0.001.

Table 3
Maximum R2 coefficient of determination, regression coefficient (β1) values and t-scores produced by generalised linear model regressions between the corresponding biotic and environmental turnover results generated by 25 moving window orientations ...

Geographically stratified turnover

GLM regressions between biotic and environmental turnover in the stratified study area produced the maximal R2 coefficients and corresponding regression coefficients shown in Table 4. The two R2 coefficients in bold were not significant (flora-temperature P = 0.108; flora-elevation P = 0.026), whilst all other maximal coefficients were significant at P<0.01, with many having very small P values.

Table 4
Maximum R2 coefficient of determination and regression coefficient (β1) values and t-scores produced by generalised linear model regressions between the corresponding biotic and environmental turnover patterns generated by 25 moving window orientations, ...

Several differences between the biotic-environmental relationships emerged within the different turnover zones. Topographic and climatic turnover were more strongly predictive of mammal turnover in the northern break and eastern overlap compared to the western overlap, whilst geology and weathering turnover showed comparatively weaker associations with mammal turnover across all three zones. In comparison to mammals, environmental turnover had limited predictive capacity for variation in plant turnover, except for in the western overlap where weathering, elevation, slope and temperature turnover were substantially more predictive. Weathering turnover was also much more predictive of plant turnover than the other variables in the eastern overlap.

Most regression slopes were positive, but there were some negative relationships between biotic and environmental turnover. These occurred mainly in the two overlap zones, with only one negative slope in the northern break, and applied most consistently to relationships between plant, elevation, solar radiation and temperature turnover.


This research aimed to identify high resolution biogeographic breaks and transitions in multi-taxon distributions, the specific genera driving break formation and the environmental role in generating these patterns using a study area in south-eastern Australia as an example. Biotic and environmental turnover patterns in this region were depicted in substantially greater detail than was achieved by Di Virgilio et al. [19]. This is due the more localised extent of the present analyses, their finer resolution and the smaller dimensions of the moving window analyses. Gradational variation in turnover within and along the boundaries of break and overlap zones was shown more clearly and denoted the continuous change in species assemblages. Other studies such as [1], [2], [5], [16] have shown variation in turnover along the length of biogeographical boundaries, but over broad extents and at coarse resolutions. This more detailed characterisation of transitions in biotic and environmental breaks and overlaps could be informative to conservation planning, since priority is often assigned to areas of ecological transition where multiple taxa from adjacent assemblages overlap [39]. Typically conservation strategy relies on broad-scale ecoregions (e.g. IBRA) with solid boundaries, the limitations of which have been raised by others, e.g. [40].

The apparent breaks and overlaps do not appear closely linked to differing land-use patterns, except for an area of cleared/disturbed land which corresponds with the location of the southern extension of high plant turnover (Fig. 7). However, there are no further obvious parallels between cleared/disturbed land and the other main turnover patterns. National parks, reserves and NSW state forests are located in the same locations as the main portion of the northern break and each of the adjacent overlaps (Fig. 7). Since protected land and state forests exist in both high and low turnover areas, it seems unlikely that turnover patterns principally reflect land usage patterns. Moreover, the shape of the eastern overlap does not match the outline of the reserves and state forests.

Figure 7
Plant species turnover overlaid with the location of cleared/disturbed land and reserves, national parks and state forests.

We also quantified sample redundancy across SEC to reconcile the turnover patterns with the effects of poor sampling. Plant sampling redundancy ranges from low to medium levels in the north of the study area (Fig. 3, panels C and D), which is the approximate location of the large floral break zone, so it is unlikely that poor sampling alone explains the existence of this break. Redundancy variation is also smoother and more uniform than the turnover patterns, e.g. most of southern SEC shows uniformly high sample redundancy, but there are two disjunct areas of low turnover in this same area and not uniform, low turnover.

Despite some fine-scale similarities, overall the environmental turnover maps lacked consistent parallels with biotic turnover. For instance, high and low environmental turnover corresponded to separate areas of low biotic turnover. This appears to support the contention that greater environmental diversity may not always represent increased species richness [41].

Generalised linear model (GLM) regressions between biotic and abiotic turnover across SEC were much more predictive of mammal turnover than for plants. The positive relationships between most physical variables and biotic turnover suggest an intrinsic link between higher rates of environmental turnover and higher biotic turnover across SEC, although the opposite trend applied to weathering intensity. Relationships between plant and environmental turnover were weak across SEC, and varied only slightly. This contrasts with strong relationships between topography, temperature and mammal turnover, suggesting variation in these variables is influencing mammal turnover across SEC.

However, the R2 maxima derived from GLM regression fitting when SEC was stratified into discrete high and low turnover zones provided more insight into the relative contribution of different physical factors in affecting turnover. For instance, topography and climate were more predictive of mammal turnover in the northern break and eastern overlap, but markedly less so in the western overlap. Previous studies at continental scale have shown a joint role for topographic complexity and climate in influencing mammal turnover [42], but these results at a regional extent show how these associations vary in different areas with different geographic and environmental characteristics. Conversely, topography and climate lacked strong relationships with plant turnover in the low, uniform landscape of the eastern overlap and the northern break zone, but they were much more predictive of plant turnover in the hilly western overlap, which is broadly consistent with some previous findings. For instance, the response of inter-annual vegetation variability to climatic changes was closely linked to topographic variability [21]. Species distribution models of Eucalyptus fastigata in south-eastern NSW provided misleading range estimates when using only climatic variables and omitting topography [43]. The western overlap also contained 39 Eucalyptus species compared to 29 species in the eastern overlap, which accords with the finding of Austin et al. [22] that eucalypt species richness is lower on flats and in gullies compared to on ridges and slopes.

Negative relationships between environmental and biotic turnover were more common to plants and were exclusive to the two low turnover zones, with the exception of plant-temperature turnover in the northern break, which was only slightly negative (β1 = −0.021). This suggests that greater environmental variability in the break influences higher biotic turnover and potentially lower overall beta diversity and species richness. In contrast, increased variability in elevation, slope, solar radiation and temperature within the low turnover zones is associated with lower plant turnover and hence higher species richness because the distributions of multiple species overlap in these regions. Complex, heterogeneous environments have often been associated with higher species richness because of the greater variety of habitats and more opportunities for resource exploitation that they provide [44][46]. The negative biotic-environmental relationships in the two overlap zones are broadly consistent with this axiom. However, the large number of positive relationships between biotic and environmental turnover, especially within the break zone, suggest the opposite pattern with greater environmental heterogeneity engendering higher species turnover and consequently lower richness.

The northern break zone and the adjacent, smaller eastern overlap both contained similar genera and species numbers, despite many more observations being sampled in the latter. This overlap contains several cells of particularly low turnover, plus a band of low turnover extending ~25 km southwest and connecting with a separate low turnover zone at −37°. These low turnover cells indicate the presence of comparatively stable communities with high richness and diversity. They also coincide with areas of relatively low topographic and climatic variability. The connecting strip is 4 to 8 km in width and mostly comprises reserves and NSW state forests, though there is some incursion of cleared/disturbed land at its narrowest point. Nonetheless, this connecting strip of low turnover warrants investigation as a potential corridor linking the two overlaps, especially since corridors comprising natural landscape features may increase species movement in fragmented landscapes [47].

The eastern overlap may also be viable for additional conservation reserve siting. A key goal of conservation is to maximise biodiversity persistence [48], [49], which may be achieved by a variety of measures, such as maximising species abundances [50], minimising habitat fragmentation [51] and buffering against outside threats with compact reserve designs [52]. The eastern overlap appears to satisfy some of these criteria, especially since it contains a similar variety of genera and species to the much larger neighbouring break, but in a significantly smaller area, which may be more cost-effective. Currently this overlap hosts both dedicated nature reserves and NSW state forests, with the latter amenable to forestry operations regulated by the Forestry and National Park Estate Act, 1998 [53].

The northern break zone contained comparatively high proportions of unique plant and mammal genera and these same genera may play a significant role in driving higher rates of turnover. Figure 8 (panel A) shows the distribution density of plant genera common to the break and eastern overlap, with the majority occurring in the overlap. Genera exclusive to the eastern overlap are evenly distributed (Fig. 8, panel B). Both findings are expected because this is an area of comparatively stable, high diversity. That the high turnover zone hosts more than twice the number of unique plant genera than the western overlap and the same number as the eastern overlap is interesting, since break zones are characterised by higher rates of community change and lower overall diversity. However, the majority of these unique genera are clustered in a small portion of the break, in the coastal strip within the southern part of the Bateman bioregion (Fig. 8, panel D). This section of the Bateman bioregion has relatively low environmental turnover, but this increases further inland in the break zone where there are no exclusive genera. Plant genera sampled in the break are evenly distributed (Fig. 8, panel C), hence this pattern is not simply due to biased sampling.

Figure 8
Distribution densities of plant and mammal genera unique and common to break and overlap zones.

The distribution of mammal genera in the break and eastern overlap zones shows the same general pattern (Fig. 8 panel E) with more unique genera occurring exclusively in the break and also clustering in the southern part of the Bateman bioregion (Fig. 8, panel H). Although there are several state forests in south Bateman, ~11% of this region is protected by reserves (Fig. 7), hence this area may also warrant consideration in conservation planning.


We have mapped gradational biotic and environmental turnover patterns in significantly greater detail and for a smaller region than Di Virgilio et al. [19], demonstrating the applicability of this method to a range of scales and resolutions, potentially for any global location. There was a combination of concordance and divergence between species turnover patterns and the widely used IBRA bioregionalisation framework, which cautions against applying this framework too rigidly in natural resource planning. The environmental influence on turnover differed for mammals and plants in different subsections of SEC corresponding to separate break and overlap zones. For instance, topographic and temperature turnover showed stronger associations with plant turnover in a hilly, high altitude overlap zone compared to a homogenous coastal overlap. The capacity to quantify the gradational variation of discrete break and overlap zones, depict their transitional boundaries in high resolution, and identify genera exclusive to each zone could inform conservation planning, for instance by identifying candidate reserve areas.

Supporting Information

Figure S1

Plant species turnover maps in the South East Corner (SEC) bioregion of south-eastern New South Wales, Australia, for Sørensen moving window analyses rotated through 360° in 15° increments. Plant species turnover maps, moving window orientations 90° to 315°.


Figure S2

Plant species turnover maps in the South East Corner (SEC) bioregion of south-eastern New South Wales, Australia, for Sørensen moving window analyses rotated through 360° in 15° increments. Plant species turnover maps, moving window orientations 300° to 165°.


Figure S3

Plant species turnover maps in the South East Corner (SEC) bioregion of south-eastern New South Wales, Australia, for Sørensen moving window analyses rotated through 360° in 15° increments. Plant species turnover maps, moving window orientations 150° to 90°.


File S1

Comparison of the plant and mammal genera unique to northern break and eastern overlap zones, unique to the same break and a western overlap zone and common to the intersections of each pair. Table S1, Plant genera that occur exclusively in the northern biogeographic break zone (n = 79), exclusively in the eastern overlap zone (n = 79) and genera that occur in both areas (n = 270) within the South East Corner (SEC) study area in New South Wales, Australia. Table S2, Plant genera that occur exclusively in the northern biogeographic break zone (n = 135), exclusively in the western overlap zone (n = 57) and genera that occur in both areas (n = 214) within the South East Corner (SEC) study area in New South Wales, Australia. Table S3, Mammal genera that occur exclusively in the northern biogeographic break zone (n = 2), exclusively in the eastern overlap zone (n = 1) and genera that occur in both areas (n = 32) within the South East Corner (SEC) study area in New South Wales, Australia. Table S4, Mammal genera that occur exclusively in the northern biogeographic break zone (n = 8), exclusively in the western overlap zone (n = 0) and genera that occur in both areas (n = 26) within the South East Corner (SEC) study area in New South Wales, Australia.


Funding Statement

This research was supported by the Australian Research Council’s “Future Fellow” funding scheme (project number FT0992002, awarded to MCE) and a University of New South Wales Goldstar grant (RG114989 to SWL and MCE). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Williams PH (1996) Mapping variations in the strength and breadth of biogeographic transition zones using species turnover. Proceedings of the Royal Society of London Series B-Biological Sciences 263: 579–588
2. Ruggiero A, Lawton JH, Blackburn TM (1998) The geographic ranges of mammalian species in South America: spatial patterns in environmental resistance and anisotropy. Journal of Biogeography 25: 1093–1103
3. Morrone JJ (2004) The South American transition zone: Characterization and evolutionary relevance. Acta Entomológica Chilena 28: 41–50
4. Kreft H, Jetz W (2010) A framework for delineating biogeographical regions based on species distributions. Journal of Biogeography 37: 2029–2053
5. Thiel-Egenter C, Alvarez N, Holderegger R, Tribsch A, Englisch T, et al. (2011) Break zones in the distributions of alleles and species in alpine plants. Journal of Biogeography 38: 772–782
6. Rahbek C (2005) The role of spatial scale and the perception of large-scale species-richness patterns. Ecology Letters 8: 224–239
7. Burley HM, Laffan SW, Williams KJ (2012) Spatial non-stationarity and anisotropy of compositional turnover in eastern Australian Myrtaceae species. International Journal of Geographical Information Science 26: 2065–2081
8. O’Brien EM, Field R, Whittaker RJ (2000) Climatic gradients in woody plant (tree and shrub) diversity: water-energy dynamics, residual variation, and topography. Oikos 89: 588–600
9. Jenkins DG, Ricklefs RE (2011) Biogeography and ecology: two views of one world. Philosophical Transactions of the Royal Society B: Biological Sciences 366: 2331–2335 [PMC free article] [PubMed]
10. Smith TB, Kark S, Schneider CJ, Wayne RK, Moritz C (2001) Biodiversity hotspots and beyond: the need for preserving environmental transitions. Trends in Ecology & Evolution 16: 431–431
11. Araujo MB (2002) Biodiversity hotspots and zones of ecological transition. Conservation Biology 16: 1662–1663
12. Shriner SA, Wilson KR, Flather CH (2006) Reserve networks based on richness hotspots and representation vary with scale. Ecological Applications 16: 1660–1673 [PubMed]
13. Rouget M (2003) Measuring conservation value at fine and broad scales: implications for a diverse and fragmented region, the Agulhas Plain. Biological Conservation 112: 217–232
14. Knight AT, Driver A, Cowling RM, Maze K, Desmet PG, et al. (2006) Designing Systematic Conservation Assessments that Promote Effective Implementation: Best Practice from South Africa. Conservation Biology 20: 739–750 [PubMed]
15. Kremen C, Cameron A, Moilanen A, Phillips SJ, Thomas CD, et al. (2008) Aligning conservation priorities across taxa in Madagascar with high-resolution planning tools. Science 320: 222–226 [PubMed]
16. Williams PH, de Klerk HM, Crowe TM (1999) Interpreting biogeographical boundaries among Afrotropical birds: spatial patterns in richness gradients and species replacement. Journal of Biogeography 26: 459–474
17. Wilson MV, Shmida A (1984) Measuring beta diversity with presence absence data. Journal of Ecology 72: 1055–1064
18. Harrison S, Ross SJ, Lawton JH (1992) Beta-diversity on geographic gradients in Britain. Journal of Animal Ecology 61: 151–158
19. Di Virgilio G, Laffan SW, Ebach MC (2012) Fine scale quantification of floral and faunal breaks and their geographic correlates, with an example from south-eastern Australia. Journal of Biogeography 39: 1862–1876
20. Davidar P, Rajagopal B, Mohandass D, Puyravaud JP, Condit R, et al. (2007) The effect of climatic gradients, topographic variation and species traits on the beta diversity of rain forest trees. Global Ecology and Biogeography 16: 510–518
21. White AB, Kumar P, Tcheng D (2005) A data mining approach for understanding topographic control on climate-induced inter-annual vegetation variability over the United States. Remote Sensing of Environment 98: 1–20
22. Austin MP, Pausas JG, Nicholls AO (1996) Patterns of tree species richness in relation to environment in southeastern New South Wales, Australia. Australian Journal of Ecology 21: 154–164
23. Austin MP, Nicholls AO, Margules CR (1990) Measurement of the Realized Qualitative Niche: Environmental Niches of Five Eucalyptus Species. Ecological Monographs 60: 161–177
24. Austin MP, Smith TM, Van Niel KP, Wellington AB (2009) Physiological responses and statistical models of the environmental niche: a comparative study of two co-occurring Eucalyptus species. Journal of Ecology 97: 496–507
25. Department of Sustainability, Environment, Water, Population and Communities (2004) Interim Biogeographic Regionalisation of Australia (IBRA Version 6.1). Developed through the cooperative efforts of the Australia Government Department of the Environment and Heritage and State/Territory land management agencies. Canberra, ACT, Australia.
26. Thackway R, Cresswell ID (1997) A bioregional framework for planning the national system of protected areas in Australia. Natural Areas Journal 17: 241–247
27. National Reserve System Task Group (2009) Australia’s Strategy for the National Reserve System 2009-2030, Prepared by the National Reserve System Task Group Convened under the Natural Resource Management Ministerial Council. 1st ed. Canberra, ACT: Australian Government. 67.
28. OEH (2010) Atlas of NSW Wildlife. New South Wales (NSW) Office of Environment and Heritage (OEH) Available: Accessed 2010 Aug 12.
29. Jarvis A, Reuter HI, Nelson A, Guevara E (2008) Hole-filled SRTM for the globe Version 4, available from the CGIAR-CSI SRTM 90m Database. Available:
30. Raymond OL, Retter AJ, (eds) (2010) Surface geology of Australia 1:1,000,000 scale, 2010 edn [Digital Dataset] Geoscience Australia. Available:
31. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A (2005) Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 25: 1965–1978
32. Wilford J (2012) A weathering intensity index for the Australian continent using airborne gamma-ray spectrometry and digital terrain analysis. Geoderma 183: 124–142
33. Laffan SW, Lubarsky E, Rosauer DF (2010) Biodiverse, a tool for the spatial analysis of biological and related diversity. Ecography 33: 643–647
34. Pepper M, Doughty P, Arculus R, Keogh JS (2008) Landforms predict phylogenetic structure on one of the world’s most ancient surfaces. BMC Evolutionary Biology 8: 152. [PMC free article] [PubMed]
35. Camarero JJ, Gutierrez E, Fortin M (2006) Spatial patterns of plant richness across treeline ecotones in the Pyrenees reveal different locations for richness and tree cover boundaries. Global Ecology and Biogeography 15: 182–191
36. Lennon JJ, Koleff P, Greenwood JJD, Gaston KJ (2001) The geographical structure of British bird distributions: diversity, spatial turnover and scale. Journal of Animal Ecology 70: 966–979
37. Garcillan PP, Ezcurra E, Riemann H (2003) Distribution and species richness of woody dryland legumes in Baja California, Mexico. Journal of Vegetation Science 14: 475–486
38. Nelder JA, Wedderburn RWM (1972) Generalized Linear Models. Journal of the Royal Statistical Society Series A (General) 135: 370–384
39. Gaston KJ, Rodrigues ASL, van Rensburg BJ, Koleff P, Chown SL (2001) Complementary representation and zones of ecological transition. Ecology Letters 4: 4–9
40. McDonald R, McKnight M, Weiss D, Selig E, O’Connor M, et al. (2005) Species compositional similarity and ecoregions: Do ecoregion boundaries represent zones of high species turnover? Biological Conservation 126: 24–40
41. Araujo MB, Humphries CJ, Densham PJ, Lampinen R, Hagemeijer WJM, et al. (2001) Would environmental diversity be a good surrogate for species diversity? Ecography 24: 103–110
42. Qian H, Badgley C, Fox DL (2009) The latitudinal gradient of beta diversity in relation to climate and topography for mammals in North America. Global Ecology and Biogeography 18: 111–122
43. Austin MP, Van Niel KP (2011) Impact of landscape predictors on climate change modelling of species distributions: a case study with Eucalyptus fastigata in southern New South Wales, Australia. Journal of Biogeography 38: 9–19
44. MacArthur RH, Wilson RHMAEO (1967) The Theory of Island Biogeography: Princeton University Press.
45. Huston MA (1994) Biological Diversity: The Coexistence of Species: University Press.
46. Torras O, Gil-Tena A, Saura S (2008) How does forest landscape structure explain tree species richness in a Mediterranean context? Biodiversity and Conservation 17: 1227–1240
47. Gilbert-Norton L, Wilson R, Stevens JR, Beard KH (2010) A Meta-Analytic Review of Corridor Effectiveness. Conservation Biology 24: 660–668 [PubMed]
48. Pressey RL, Cabeza M, Watts ME, Cowling RM, Wilson KA (2007) Conservation planning in a changing world. Trends in Ecology & Evolution 22: 583–592 [PubMed]
49. Cabeza M, Moilanen A (2001) Design of reserve networks and the persistence of biodiversity. Trends in Ecology & Evolution 16: 242–248 [PubMed]
50. Rodrigues ASL, Gregory RD, Gaston KJ (2000) Robustness of reserve selection procedures under temporal species turnover. Proceedings of the Royal Society B-Biological Sciences 267: 49–55 [PMC free article] [PubMed]
51. Van Teeffelen AJA, Cabeza M, Moilanen A (2006) Connectivity, probabilities and persistence: Comparing reserve selection strategies. Biodiversity and Conservation 15: 899–919
52. Cabeza M, Araujo MB, Wilson RJ, Thomas CD, Cowley MJR, et al. (2004) Combining probabilities of occurrence with spatial reserve design. Journal of Applied Ecology 41: 252–262
53. Forestry and National Park Estate Act (1998) Australia: Government of New South Wales.

Articles from PLoS ONE are provided here courtesy of Public Library of Science