Green turtle somatic growth dynamics in a spatially disjunct Great Barrier Reef metapopulation

The growth dynamics of green sea turtles resident in four separate foraging grounds of the southern Great Barrier Reef genetic stock were assessed using a nonparametric regression modeling approach. Juveniles recruit to these grounds at the same size, but grow at foraging-ground-dependent rates that result in significant differences in expected size- or age-at-maturity. Mean age-at-maturity was estimated to vary from 25-50 years depending on the ground. This stock comprises mainly the same mtDNA haplotype, so geographic variability might be due to local environmental conditions rather than genetic factors, although the variability was not a function of latitudinal variation in environmental conditions or whether the food stock was seagrass or algae. Temporal variability in growth rates was evident in response to local environmental stochasticity, so geographic variability might be due to local food stock dynamics. Despite such variability, the expected size-specific growth rate function at all grounds displayed a similar nonmonotonic growth pattern with a juvenile growth spurt at 60-70 cm curved carapace length, (CCL) or 15-20 years of age. Sex-specific growth differences were also evident with females tending to grow faster than similar-sized males after the Juvenile growth spurt. It is clear that slow sex-specific growth displaying both spatial and temporal variability and a juvenile growth spurt are distinct growth behaviors of green turtles from this stock.

A stochastic metapopulation model accounting for habitat dynamics

A stochastic metapopulation model accounting for habitat dynamics is presented. This is the stochastic SIS logistic model with the novel aspect that it incorporates varying carrying capacity. We present results of Kurtz and Barbour, that provide deterministic and diffusion approximations for a wide class of stochastic models, in a form that most easily allows their direct application to population models. These results are used to show that a suitably scaled version of the metapopulation model converges, uniformly in probability over finite time intervals, to a deterministic model previously studied in the ecological literature. Additionally, they allow us to establish a bivariate normal approximation to the quasi-stationary distribution of the process. This allows us to consider the effects of habitat dynamics on metapopulation modelling through a comparison with the stochastic SIS logistic model and provides an effective means for modelling metapopulations inhabiting dynamic landscapes.

The role of habitat disturbance and recovery in metapopulation persistence

Classical metapopulation theory assumes a static landscape. However, empirical evidence indicates many metapopulations are driven by habitat succession and disturbance. We develop a stochastic metapopulation model, incorporating habitat disturbance and recovery, coupled with patch colonization and extinction, to investigate the effect of habitat dynamics on persistence. We discover that habitat dynamics play a fundamental role in metapopulation dynamics. The mean number of suitable habitat patches is not adequate for characterizing the dynamics of the metapopulation. For a fixed mean number of suitable patches, we discover that the details of how disturbance affects patches and how patches recover influences metapopulation dynamics in a fundamental way. Moreover, metapopulation persistence is dependent not only oil the average lifetime of a patch, but also on the variance in patch lifetime and the synchrony in patch dynamics that results from disturbance. Finally, there is an interaction between the habitat and metapopulation dynamics, for instance declining metapopulations react differently to habitat dynamics than expanding metapopulations. We close, emphasizing the importance of using performance measures appropriate to stochastic systems when evaluating their behavior, such as the probability distribution of the state of the. metapopulation, conditional on it being extant (i.e., the quasistationary distribution).

Limited effect of anthropogenic habitat fragmentation on molecular diversity in a rain forest skink, Gnypetoscincus queenslandiae

To examine the effects of recent habitat fragmentation, we assayed genetic diversity in a rain forest endemic lizard, the prickly forest skink (Gnypetoscincus queenslandiae), from seven forest fragments and five sites in continuous forest on the Atherton tableland of northeastern Queensland, Australia. The rain forest in this region was fragmented by logging and clearing for dairy farms in the early 1900s and most forest fragments studied have been isolated for 50-80 years or nine to 12 skink generations. We genotyped 411 individuals at nine microsatellite DNA loci and found fewer alleles per locus in prickly forest skinks from small rain forest fragments and a lower ratio of allele number to allele size range in forest fragments than in continuous forest, indicative of a decrease in effective population size. In contrast, and as expected for populations with small neighbourhood sizes, neither heterozygosity nor variance in allele size differed between fragments and sites in continuous forests. Considering measures of among population differentiation, there was no increase in F-ST among fragments and a significant isolation by distance pattern was identified across all 12 sites. However, the relationship between genetic (F-ST) and geographical distance was significantly stronger for continuous forest sites than for fragments, consistent with disruption of gene flow among the latter. The observed changes in genetic diversity within and among populations are small, but in the direction predicted by the theory of genetic erosion in recently fragmented populations. The results also illustrate the inherent difficulty in detecting genetic consequences of recent habitat fragmentation, even in genetically variable species, and especially when effective population size and dispersal rates are low.

Theory for designing nature reserves for single species

We examine the question of the optimal number of reserves that should be established to maximize the persistence of a species. We assume that the mean time to extinction of a single population increases as a power of the habitat area, that there is a certain amount of habitat to be reserved, and that the aim is to determine how this habitat is most efficiently divided. The optimal configuration depends on whether the management objective is to maximize the mean time to extinction or minimize the risk of extinction. When maximizing the mean time to extinction, the optimal number of independent reserves does not depend on the amount of available habitat for the reserve system. In contrast, the risk of extinction is minimized when individual reserves are equal to the optimal patch size, making the optimal number of reserves linearly proportional to the amount of available habitat. A model that includes dispersal and correlation in the incidence of extinction demonstrates the importance of considering the relative rate at which these two factors decrease with distance between reserves. A small number of reserves is optimal when the mean time to extinction increases rapidly with habitat area or when risks of extinction are high.

Estimates of sex- and age-class-specific survival probabilities for a southern Great Barrier Reef green sea turtle population

Sex- and age-class-specific survival probabilities of a southern Great Barrier Reef green sea turtle population were estimated using a capture - mark - recapture (CMR) study and a Cormack - Jolly - Seber (CJS) modelling approach. The CMR history profiles for 954 individual turtles tagged over a 9-year period ( 1984 - 1992) were classified into three age classes ( adult, subadult, juvenile) based on somatic growth and reproductive traits. Reduced-parameter CJS models, accounting for constant survival and time-specific recapture, fitted best for all age classes. There were no significant sex-specific differences in either survival or recapture probabilities for any age class. Mean annual adult survival was estimated at 0.9482 (95% CI: 0.92 - 0.98) and was significantly higher than survival for either subadults or juveniles. Mean annual subadult survival was 0.8474 ( 95% CI: 0.79 - 0.91), which was not significantly different from mean annual juvenile survival estimated at 0.8804 ( 95% CI: 0.84 - 0.93). The time-specific adult recapture probabilities were a function of sampling effort but this was not the case for either juveniles or subadults. The sampling effort effect was accounted for explicitly in the estimation of adult survival and recapture probabilities. These are the first comprehensive sex- and age-class-specific survival and recapture probability estimates for a green sea turtle population derived from a long-term CMR program.

Decreased relatedness between male prickly forest skinks (Gnypetoscincus queenslandiae) in habitat fragments

In species with low levels of dispersal the chance of closely related individuals breeding may be a potential problem; sex-biased dispersal is a mechanism that may decrease the possibility of cosanguineous mating. Fragmentation of the habitat in which a species lives may affect mechanisms such as sex-biased dispersal, which may in turn exacerbate more direct effects of fragmentation such as decreasing population size that may lead to inbreeding depression. Relatedness statistics calculated using microsatellite DNA data showed that rainforest fragmentation has had an effect on the patterns of dispersal in the prickly forest skink (Gnypetoscincus queenslandiae), a rainforest endemic of the Wet Tropics of north eastern Australia. A lower level of relatedness was found in fragments compared to continuous forest sites due to a significantly lower level of pairwise relatedness between males in rainforest fragments. The pattern of genetic relatedness between sexes indicates the presence of male-biased dispersal in this species, with a stronger pattern detected in populations in rainforest fragments. Male prickly forest skinks may have to move further in fragmented habitat in order to find mates or suitable habitat logs.

A new method for conservation planning for the persistence of multiple species

Although the aim of conservation planning is the persistence of biodiversity, current methods trade-off ecological realism at a species level in favour of including multiple species and landscape features. For conservation planning to be relevant, the impact of landscape configuration on population processes and the viability of species needs to be considered. We present a novel method for selecting reserve systems that maximize persistence across multiple species, subject to a conservation budget. We use a spatially explicit metapopulation model to estimate extinction risk, a function of the ecology of the species and the amount, quality and configuration of habitat. We compare our new method with more traditional, area-based reserve selection methods, using a ten-species case study, and find that the expected loss of species is reduced 20-fold. Unlike previous methods, we avoid designating arbitrary weightings between reserve size and configuration; rather, our method is based on population processes and is grounded in ecological theory.

Accounting for uncertainty in marine reserve design

Ecosystems and the species and communities within them are highly complex systems that defy predictions with any degree of certainty. Managing and conserving these systems in the face of uncertainty remains a daunting challenge, particularly with respect to developing networks of marine reserves. Here we review several modelling frameworks that explicitly acknowledge and incorporate uncertainty, and then use these methods to evaluate reserve spacing rules given increasing levels of uncertainty about larval dispersal distances. Our approach finds similar spacing rules as have been proposed elsewhere - roughly 20-200 km - but highlights several advantages provided by uncertainty modelling over more traditional approaches to developing these estimates. In particular, we argue that uncertainty modelling can allow for (1) an evaluation of the risk associated with any decision based on the assumed uncertainty; (2) a method for quantifying the costs and benefits of reducing uncertainty; and (3) a useful tool for communicating to stakeholders the challenges in managing highly uncertain systems. We also argue that incorporating rather than avoiding uncertainty will increase the chances of successfully achieving conservation and management goals.

Does colonization asymmetry matter in metapopulations?

Despite the considerable evidence showing that dispersal between habitat patches is often asymmetric, most of the metapopulation models assume symmetric dispersal. In this paper, we develop a Monte Carlo simulation model to quantify the effect of asymmetric dispersal on metapopulation persistence. Our results suggest that metapopulation extinctions are more likely when dispersal is asymmetric. Metapopulation viability in systems with symmetric dispersal mirrors results from a mean field approximation, where the system persists if the expected per patch colonization probability exceeds the expected per patch local extinction rate. For asymmetric cases, the mean field approximation underestimates the number of patches necessary for maintaining population persistence. If we use a model assuming symmetric dispersal when dispersal is actually asymmetric, the estimation of metapopulation persistence is wrong in more than 50% of the cases. Metapopulation viability depends on patch connectivity in symmetric systems, whereas in the asymmetric case the number of patches is more important. These results have important implications for managing spatially structured populations, when asymmetric dispersal may occur. Future metapopulation models should account for asymmetric dispersal, while empirical work is needed to quantify the patterns and the consequences of asymmetric dispersal in natural metapopulations.

Larval dispersal reveals regional sources and sinks in the Great Barrier Reef

We analysed simulated connectivity patterns for reef fish larvae in the Cairns section of the Great Barrier Reef, and identified 2 key subregions that exhibit regional scale source–sink dynamics. The source and sink were separated latitudinally by a boundary at 16.1°S, with the source subregion lying to the north. Larval transport between the 2 subregions was predominantly unidirectional, from north to south. Only a few local populations, described here as ‘gateway reefs’, were able to transport larvae from the sink subregion to the source subregion and thus maintain the connectedness of the metapopulation. The northern subregion was able to persist without external larval supply, but when conditions were recruitment limited, the southern subregion depended on larval supply from the north to persist. The relative autonomy of the northern subregion, and its importance in sustaining the southern subregion, will influence the effectiveness of conservation efforts.

Stochastic models for mainland-island metapopulations in static and dynamic landscapes

This paper has three primary aims: to establish an effective means for modelling mainland-island metapopulations inhabiting a dynamic landscape: to investigate the effect of immigration and dynamic changes in habitat on metapopulation patch occupancy dynamics; and to illustrate the implications of our results for decision-making and population management. We first extend the mainland-island metapopulation model of Alonso and McKane [Bull. Math. Biol. 64:913-958,2002] to incorporate a dynamic landscape. It is shown, for both the static and the dynamic landscape models, that a suitably scaled version of the process converges to a unique deterministic model as the size of the system becomes large. We also establish that. under quite general conditions, the density of occupied patches, and the densities of suitable and occupied patches, for the respective models, have approximate normal distributions. Our results not only provide us with estimates for the means and variances that are valid at all stages in the evolution of the population, but also provide a tool for fitting the models to real metapopulations. We discuss the effect of immigration and habitat dynamics on metapopulations, showing that mainland-like patches heavily influence metapopulation persistence, and we argue for adopting measures to increase connectivity between this large patch and the other island-like patches. We illustrate our results with specific reference to examples of populations of butterfly and the grasshopper Bryodema tuberculata.

Testing alternative models for the conservation of koalas in fragmented rural-urban landscapes

Predicting the various responses of different species to changes in landscape structure is a formidable challenge to landscape ecology. Based on expert knowledge and landscape ecological theory, we develop five competing a priori models for predicting the presence/absence of the Koala (Phascolarctos cinereus) in Noosa Shire, south-east Queensland (Australia). A priori predictions were nested within three levels of ecological organization: in situ (site level) habitat (< 1 ha), patch level (100 ha) and landscape level (100-1000 ha). To test the models, Koala surveys and habitat surveys (n = 245) were conducted across the habitat mosaic. After taking into account tree species preferences, the patch and landscape context, and the neighbourhood effect of adjacent present sites, we applied logistic regression and hierarchical partitioning analyses to rank the alternative models and the explanatory variables. The strongest support was for a multilevel model, with Koala presence best predicted by the proportion of the landscape occupied by high quality habitat, the neighbourhood effect, the mean nearest neighbour distance between forest patches, the density of forest patches and the density of sealed roads. When tested against independent data (n = 105) using a receiver operator characteristic curve, the multilevel model performed moderately well. The study is consistent with recent assertions that habitat loss is the major driver of population decline, however, landscape configuration and roads have an important effect that needs to be incorporated into Koala conservation strategies.