Declining survival rates in a reintroduced population of the Mauritius kestrel: evidence for non-linear density dependence and environmental stochasticity

1. We studied a reintroduced population of the formerly critically endangered Mauritius kestrel Falco punctatus Temmink from its inception in 1987 until 2002, by which time the population had attained carrying capacity for the study area. Post-1994 the population received minimal management other than the provision of nestboxes. 2. We analysed data collected on survival (1987-2002) using program MARK to explore the influence of density-dependent and independent processes on survival over the course of the population's development. 3.We found evidence for non-linear, threshold density dependence in juvenile survival rates. Juvenile survival was also strongly influenced by climate, with the temporal distribution of rainfall during the cyclone season being the most influential climatic variable. Adult survival remained constant throughout. 4. Our most parsimonious capture-mark-recapture statistical model, which was constrained by density and climate, explained 75.4% of the temporal variation exhibited in juvenile survival rates over the course of the population's development. 5. This study is an example of how data collected as part of a threatened species recovery programme can be used to explore the role and functional form of natural population regulatory processes. With the improvements in conservation management techniques and the resulting success stories, formerly threatened species offer unique opportunities to further our understanding of the fundamental principles of population ecology.

The Evolution of Male-Biased Dispersal under the Joint Selective Forces of Inbreeding Load and Demographic and Environmental Stochasticity

Peer reviewed

Patchy populations in stochastic environments: Critical number of patches for persistence

We introduce a model for the dynamics of a patchy population in a stochastic environment and derive a criterion for its persistence. This criterion is based on the geometric mean (GM) through time of the spatial-arithmetic mean of growth rates. For the population to persist, the GM has to be greater than or equal to1. The GM increases with the number of patches (because the sampling error is reduced) and decreases with both the variance and the spatial covariance of growth rates. We derive analytical expressions for the minimum number of patches (and the maximum harvesting rate) required for the persistence of the population. As the magnitude of environmental fluctuations increases, the number of patches required for persistence increases, and the fraction of individuals that can be harvested decreases. The novelty of our approach is that we focus on Malthusian local population dynamics with high dispersal and strong environmental variability from year to year. Unlike previous models of patchy populations that assume an infinite number of patches, we focus specifically on the effect that the number of patches has on population persistence. Our work is therefore directly relevant to patchily distributed organisms that are restricted to a small number of habitat patches.

Population demography and the evolution of helping behaviors.

Limited dispersal may favor the evolution of helping behaviors between relatives as it increases their relatedness, and it may inhibit such evolution as it increases local competition between these relatives. Here, we explore one way out of this dilemma: if the helping behavior allows groups to expand in size, then the kin-competition pressure opposing its evolution can be greatly reduced. We explore the effects of two kinds of stochasticity allowing for such deme expansion. First, we study the evolution of helping under environmental stochasticity that may induce complete patch extinction. Helping evolves if it results in a decrease in the probability of extinction or if it enhances the rate of patch recolonization through propagules formed by fission of nonextinct groups. This mode of dispersal is indeed commonly found in social species. Second, we consider the evolution of helping in the presence of demographic stochasticity. When fecundity is below its value maximizing deme size (undersaturation), helping evolves, but under stringent conditions unless positive density dependence (Allee effect) interferes with demographic stochasticity. When fecundity is above its value maximizing deme size (oversaturation), helping may also evolve, but only if it reduces negative density-dependent competition.

The adaptive function of melanin-based coloration in the tawny owl (Strix aluco) ; the proximate role of the melanocortin system

Colour polymorphism is common in wild population. One of the main questioning of evolutionary biologists is to understand how different colour variants could have evolved and be maintained in fluctuating environments, a selective process that forces individuals to constantly adapt their strategies in order to survive. This issue is particularly true for traits that are genetically inherited. Natural selection erodes genotypes with lowest fitness (less adapted), reducing in turn global genetic variation within population. In this context, the study of the evolution and maintenance of melanin- based coloration is relevant since inter-individual variation in the deposition of these pigments is common in animal and plant kingdoms and is under strong genetic control. In this thesis, I focus on the specific case of the tawny owl (Strix aluco), a species displaying continuous variation in reddish pheomelanin-based coloration. Interestingly, empirical studies highlighted covariations between melanin-based coloration and important behavioural, physiological and life history traits. Recently, a genetic model pointed out the melanocortin system and their pleiotropic effects as a potential regulator of these covariations. Accordingly, this PhD thesis further investigates colour-specific behavioural, physiological, or life history strategies, while examining the proximate mechanisms underlying these reaction norms. We found that differently coloured tawny owls differently resolve fundamental trade-off between offspring number and quality (Chapter 1), light melanic individuals producing many low- quality offspring and dark, melanic ones producing few high-quality offspring. These reproductive strategies are likely to induce alternative physiological constraints. Indeed, we demonstrated that light melanic individuals produced higher levels of reactive oxygen species (ROS, Chapter 2), but also expressed higher levels of antioxidant (GSH, Chapters 2 & 3). Interestingly, we showed that light melanic breeding females could modulate their POMC prohormone levels according to the environmental conditions, while dark reddish ones produced constant levels of this prohormone {Chapter 4). Finally, we highlighted colour-specific patterns of prohormone convertase 1 (PCI) gene expression (Chapter 5), an enzyme responsible for POMC prohormone processing to ACTH and a- MSH, for instance. Altogether, these results provide strong evidence of colour-specific strategies, light and melanic tawny owls better coping with stressful and relaxed environments, respectively. Variation in melanin-based coloration is likely to be maintained by the heterogeneity of our study area and strong environmental stochasticity within and between years, these process favouring differently coloured tawny owls at different periods of time. From a proximate point of view, this PhD thesis supports the hypothesis that covariations between phenotypic traits and melanin-based coloration stems from the melanocortin system, especially the fundamental role of POMC gene expression and its processing to melanocortin peptides. - Le polymorphisme de couleur est une variation phénotypique très fréquente dans la nature. En biologie évolutive, une des problématiques clés est donc de comprendre comment différent morphes de couleur peuvent être apparus et maintenus au cours du temps dans des environnements aussi variables que les nôtres, surtout que ces fluctuations forcent ces morphes à s'adapter constamment pour assurer leur survie. Cette thématique est particulièrement réelle lorsque les variations phénotypiques sont héréditaires et donc sous forte influence génétique. La sélection naturelle a en effet le pouvoir d'éroder rapidement la variation génétique en éliminant les génotypes mal adaptés. Dans ce sens, l'étude de l'évolution, et de la maintenance de la coloration mélanique est donc tout à fait pertinente car la variation de coloration entre individus est très répandue à travers les règnes animal et végétal et sous forte influence génétique. Dans cette thèse, je me suis concentré sur le cas spécifique de la chouette hulotte (Strix aluco), une espèce présentant une variation continue dans la déposition de pigments pheomélaniques roux. De précédentes études ont déjà montré que cette variation de coloration était associée avec des variations de traits comportementaux, physiologiques ou d'histoire de vie. Récemment, une étude a souligné l'importance du système des mélanocortines et de leurs effets pléiotropes dans la régulation de ces covariations. En conséquence, cette thèse de doctoral a pour but d'étudier un peu plus les stratégies comportementales, physiologiques ou d'histoire de vie spécifiques à chaque morphe de couleur, tout en examinant un peu plus les mécanismes proximaux potentiellement à la base de ces normes de réactions. Nous constatons tout d'abord que les morphes de couleurs étaient associés à différentes stratégies dans la résolution de compromis telle que la production de beaucoup de jeunes ou des jeunes de qualité (Chapitre 1). Les morphes gris (dit peu mélaniques) ont tendance à produire beaucoup de jeunes mains de moindre qualité, alors que les morphes roux (dit fortement mélaniques) produisent moins de jeunes mais de meilleure qualité. Ces stratégies sont susceptibles alors d'induire certaines contraintes physiologiques. Par exemple, nous montrons que les morphes gris produisent plus de dérivés réactifs de l'oxygène (ROS, Chapitre 2), mais aussi plus d'antioxydants (GSH, Chapitres 2 & 3). Nous montrons ensuite que les femelles grises ont une plus grande capacité à moduler leur niveau de POMC prohormone dans le sang en fonction des conditions environnementales, alors que les femelles rousses gardent un niveau constant (Chapitre 4). Finalement, nous démontrons que les patterns d'expression du gène codant pour la prohormone convertase 1 varient chez des jeunes issus de parents gris ou roux (Chapitre 5). Ceci est particulièrement intéressant car cette enzyme permet de scinder la POMC prohormone en plusieurs peptides importants tels que l'ACTH ou l'a-MSH. En conclusion, ces résultats démontrent qu'il y a bel et bien des stratégies évolutives différentes entre les morphes de couleurs, les chouettes hulottes grises et rousses étant respectivement plus adaptés à des environnements stressants ou favorables. L'hétérogénéité de notre zone d'étude et la stochasticité environnementale qui caractérise ses habitats pourraient donc agir comme une source de sélection temporelle, laquelle favoriserait les différents morphes de couleurs à diverses périodes. D'un point de vue plus proximale maintenant, cette thèse de doctorat soutient l'hypothèse que les covariations observées entre la coloration mélanique et des traits phénotypiques importants sont modulées par les effets pléiotropes du système des mélanocortines, et met en avant le rôle prépondérant que pourrait jouer l'expression du gène POMC et sa post traduction en mélanocortines.

Linking landscape- and population ecology for large population management modelling: the case of Ibex (capra ibex) in Switzerland

La biologie de la conservation est communément associée à la protection de petites populations menacées d?extinction. Pourtant, il peut également être nécessaire de soumettre à gestion des populations surabondantes ou susceptibles d?une trop grande expansion, dans le but de prévenir les effets néfastes de la surpopulation. Du fait des différences tant quantitatives que qualitatives entre protection des petites populations et contrôle des grandes, il est nécessaire de disposer de modèles et de méthodes distinctes. L?objectif de ce travail a été de développer des modèles prédictifs de la dynamique des grandes populations, ainsi que des logiciels permettant de calculer les paramètres de ces modèles et de tester des scénarios de gestion. Le cas du Bouquetin des Alpes (Capra ibex ibex) - en forte expansion en Suisse depuis sa réintroduction au début du XXème siècle - servit d?exemple. Cette tâche fut accomplie en trois étapes : En premier lieu, un modèle de dynamique locale, spécifique au Bouquetin, fut développé : le modèle sous-jacent - structuré en classes d?âge et de sexe - est basé sur une matrice de Leslie à laquelle ont été ajoutées la densité-dépendance, la stochasticité environnementale et la chasse de régulation. Ce modèle fut implémenté dans un logiciel d?aide à la gestion - nommé SIM-Ibex - permettant la maintenance de données de recensements, l?estimation automatisée des paramètres, ainsi que l?ajustement et la simulation de stratégies de régulation. Mais la dynamique d?une population est influencée non seulement par des facteurs démographiques, mais aussi par la dispersion et la colonisation de nouveaux espaces. Il est donc nécessaire de pouvoir modéliser tant la qualité de l?habitat que les obstacles à la dispersion. Une collection de logiciels - nommée Biomapper - fut donc développée. Son module central est basé sur l?Analyse Factorielle de la Niche Ecologique (ENFA) dont le principe est de calculer des facteurs de marginalité et de spécialisation de la niche écologique à partir de prédicteurs environnementaux et de données d?observation de l?espèce. Tous les modules de Biomapper sont liés aux Systèmes d?Information Géographiques (SIG) ; ils couvrent toutes les opérations d?importation des données, préparation des prédicteurs, ENFA et calcul de la carte de qualité d?habitat, validation et traitement des résultats ; un module permet également de cartographier les barrières et les corridors de dispersion. Le domaine d?application de l?ENFA fut exploré par le biais d?une distribution d?espèce virtuelle. La comparaison à une méthode couramment utilisée pour construire des cartes de qualité d?habitat, le Modèle Linéaire Généralisé (GLM), montra qu?elle était particulièrement adaptée pour les espèces cryptiques ou en cours d?expansion. Les informations sur la démographie et le paysage furent finalement fusionnées en un modèle global. Une approche basée sur un automate cellulaire fut choisie, tant pour satisfaire aux contraintes du réalisme de la modélisation du paysage qu?à celles imposées par les grandes populations : la zone d?étude est modélisée par un pavage de cellules hexagonales, chacune caractérisée par des propriétés - une capacité de soutien et six taux d?imperméabilité quantifiant les échanges entre cellules adjacentes - et une variable, la densité de la population. Cette dernière varie en fonction de la reproduction et de la survie locale, ainsi que de la dispersion, sous l?influence de la densité-dépendance et de la stochasticité. Un logiciel - nommé HexaSpace - fut développé pour accomplir deux fonctions : 1° Calibrer l?automate sur la base de modèles de dynamique (par ex. calculés par SIM-Ibex) et d?une carte de qualité d?habitat (par ex. calculée par Biomapper). 2° Faire tourner des simulations. Il permet d?étudier l?expansion d?une espèce envahisseuse dans un paysage complexe composé de zones de qualité diverses et comportant des obstacles à la dispersion. Ce modèle fut appliqué à l?histoire de la réintroduction du Bouquetin dans les Alpes bernoises (Suisse). SIM-Ibex est actuellement utilisé par les gestionnaires de la faune et par les inspecteurs du gouvernement pour préparer et contrôler les plans de tir. Biomapper a été appliqué à plusieurs espèces (tant végétales qu?animales) à travers le Monde. De même, même si HexaSpace fut initialement conçu pour des espèces animales terrestres, il pourrait aisément être étndu à la propagation de plantes ou à la dispersion d?animaux volants. Ces logiciels étant conçus pour, à partir de données brutes, construire un modèle réaliste complexe, et du fait qu?ils sont dotés d?une interface d?utilisation intuitive, ils sont susceptibles de nombreuses applications en biologie de la conservation. En outre, ces approches peuvent également s?appliquer à des questions théoriques dans les domaines de l?écologie des populations et du paysage.<br/><br/>Conservation biology is commonly associated to small and endangered population protection. Nevertheless, large or potentially large populations may also need human management to prevent negative effects of overpopulation. As there are both qualitative and quantitative differences between small population protection and large population controlling, distinct methods and models are needed. The aim of this work was to develop theoretical models to predict large population dynamics, as well as computer tools to assess the parameters of these models and to test management scenarios. The alpine Ibex (Capra ibex ibex) - which experienced a spectacular increase since its reintroduction in Switzerland at the beginning of the 20th century - was used as paradigm species. This task was achieved in three steps: A local population dynamics model was first developed specifically for Ibex: the underlying age- and sex-structured model is based on a Leslie matrix approach with addition of density-dependence, environmental stochasticity and culling. This model was implemented into a management-support software - named SIM-Ibex - allowing census data maintenance, parameter automated assessment and culling strategies tuning and simulating. However population dynamics is driven not only by demographic factors, but also by dispersal and colonisation of new areas. Habitat suitability and obstacles modelling had therefore to be addressed. Thus, a software package - named Biomapper - was developed. Its central module is based on the Ecological Niche Factor Analysis (ENFA) whose principle is to compute niche marginality and specialisation factors from a set of environmental predictors and species presence data. All Biomapper modules are linked to Geographic Information Systems (GIS); they cover all operations of data importation, predictor preparation, ENFA and habitat suitability map computation, results validation and further processing; a module also allows mapping of dispersal barriers and corridors. ENFA application domain was then explored by means of a simulated species distribution. It was compared to a common habitat suitability assessing method, the Generalised Linear Model (GLM), and was proven better suited for spreading or cryptic species. Demography and landscape informations were finally merged into a global model. To cope with landscape realism and technical constraints of large population modelling, a cellular automaton approach was chosen: the study area is modelled by a lattice of hexagonal cells, each one characterised by a few fixed properties - a carrying capacity and six impermeability rates quantifying exchanges between adjacent cells - and one variable, population density. The later varies according to local reproduction/survival and dispersal dynamics, modified by density-dependence and stochasticity. A software - named HexaSpace - was developed, which achieves two functions: 1° Calibrating the automaton on the base of local population dynamics models (e.g., computed by SIM-Ibex) and a habitat suitability map (e.g. computed by Biomapper). 2° Running simulations. It allows studying the spreading of an invading species across a complex landscape made of variously suitable areas and dispersal barriers. This model was applied to the history of Ibex reintroduction in Bernese Alps (Switzerland). SIM-Ibex is now used by governmental wildlife managers to prepare and verify culling plans. Biomapper has been applied to several species (both plants and animals) all around the World. In the same way, whilst HexaSpace was originally designed for terrestrial animal species, it could be easily extended to model plant propagation or flying animals dispersal. As these softwares were designed to proceed from low-level data to build a complex realistic model and as they benefit from an intuitive user-interface, they may have many conservation applications. Moreover, theoretical questions in the fields of population and landscape ecology might also be addressed by these approaches.

Plurality of tree species responses to drought perturbation in Bornean tropical rain forest

Drought perturbation driven by the El Niño Southern Oscillation (ENSO) is a principal stochastic variable determining the dynamics of lowland rain forest in S.E. Asia. Mortality, recruitment and stem growth rates at Danum in Sabah (Malaysian Borneo) were recorded in two 4-ha plots (trees ≥ 10 cm gbh) for two periods, 1986–1996 and 1996–2001. Mortality and growth were also recorded in a sample of subplots for small trees (10 to <50 cm gbh) in two sub-periods, 1996–1999 and 1999–2001. Dynamics variables were employed to build indices of drought response for each of the 34 most abundant plot-level species (22 at the subplot level), these being interval-weighted percentage changes between periods and sub-periods. A significant yet complex effect of the strong 1997/1998 drought at the forest community level was shown by randomization procedures followed by multiple hypothesis testing. Despite a general resistance of the forest to drought, large and significant differences in short-term responses were apparent for several species. Using a diagrammatic form of stability analysis, different species showed immediate or lagged effects, high or low degrees of resilience or even oscillatory dynamics. In the context of the local topographic gradient, species’ responses define the newly termed perturbation response niche. The largest responses, particularly for recruitment and growth, were among the small trees, many of which are members of understorey taxa. The results bring with them a novel approach to understanding community dynamics: the kaleidoscopic complexity of idiosyncratic responses to stochastic perturbations suggests that plurality, rather than neutrality, of responses may be essential to understanding these tropical forests. The basis to the various responses lies with the mechanisms of tree-soil water relations which are physiologically predictable: the timing and intensity of the next drought, however, is not. To date, environmental stochasticity has been insufficiently incorporated into models of tropical forest dynamics, a step that might considerably improve the reality of theories about these globally important ecosystems.

Relaxation of species-specific neighborhood effects in Bornean rain forest under climatic perturbation

Evidence of negative conspecific density dependence (NDD) operating on seedling survival and sapling recruitment has accumulated recently. In contrast, evidence of NDD operating on growth of trees has been circumstantial at best. Whether or not local NDD at the level of individual trees leads to NDD at the level of the community is still an open question. Moreover, whether and how perturbations interfere with these processes have rarely been investigated. We applied neighborhood models to permanent plot data from a Bornean dipterocarp forest censused over two 10-11 year periods. Although the first period was only lightly perturbed, a moderately strong El Nino event causing severe drought occurred in the first half of the second period. Such events are an important component of the environmental stochasticity affecting the region. We show that local NDD on growth of small-to-medium-sized trees may indeed translate to NDD at the level of the community. This interpretation is based on increasingly negative effects of bigger conspecific neighbors on absolute growth rates of individual trees with increasing basal area across the 18 most abundant overstory species in the first period. However, this relationship was much weaker in the second period. We interpreted this relaxation of local and community-level NDD as a consequence of increased light levels at the forest floor due to temporary leaf and twig loss of large trees in response to the drought event. Mitigation of NDD under climatic perturbation acts to decrease species richness, especially in forest overstory and therefore has an important role in determining species relative abundances at the site.

Metapopulation extinction caused by mutation accumulation

Theory suggests that the risk of extinction by mutation accumulation can be comparable to that by environmental stochasticity for an isolated population smaller than a few thousand individuals. Here we show that metapopulation structure, habitat loss or fragmentation, and environmental stochasticity can be expected to greatly accelerate the accumulation of mildly deleterious mutations, lowering the genetic effective size to such a degree that even large metapopulations may be at risk of extinction. Because of mutation accumulation, viable metapopulations may need to be far larger and better connected than would be required under just stochastic demography.

Do harvest refuges buffer kangaroos against evolutionary responses to selective harvesting?

There is a wealth of literature documenting a directional change of body size in heavily harvested populations. Most of this work concentrates on aquatic systems, but terrestrial populations are equally at risk. This paper explores the capacity of harvest refuges to counteract potential effects of size-selective harvesting on the allele frequency,of populations. We constructed a stochastic, individual-based model parameterized with data on red kangaroos. Because we do not know which part of individual growth would change in the course of natural selection, we explored the effects of two alternative models of individual growth in which alleles affect either the growth rate or the maximum size. The model results show that size-selective harvesting can result in significantly smaller kangaroos for a given age when the entire population is subject to harvesting. In contrast, in scenarios that include dispersal from harvest refuges, the initial allele frequency remains virtually unchanged.

Habitat selection and population regulation in temporally fluctuating environments

Understanding and predicting the distribution of organisms in heterogeneous environments lies at the heart of ecology, and the theory of density-dependent habitat selection (DDHS) provides ecologists with an inferential framework linking evolution and population dynamics. Current theory does not allow for temporal variation in habitat quality, a serious limitation when confronted with real ecological systems. We develop both a stochastic equivalent of the ideal free distribution to study how spatial patterns of habitat use depend on the magnitude and spatial correlation of environmental stochasticity and also a stochastic habitat selection rule. The emerging patterns are confronted with deterministic predictions based on isodar analysis, an established empirical approach to the analysis of habitat selection patterns. Our simulations highlight some consistent patterns of habitat use, indicating that it is possible to make inferences about the habitat selection process based on observed patterns of habitat use. However, isodar analysis gives results that are contingent on the magnitude and spatial correlation of environmental stochasticity. Hence, DDHS is better revealed by a measure of habitat selectivity than by empirical isodars. The detection of DDHS is but a small component of isodar theory, which remains an important conceptual framework for linking evolutionary strategies in behavior and population dynamics.

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.

Thirty-year recovery trend in the once depleted Hawaiian green sea turtle stock

The green sea turtle is one of the long-lived species that comprise the charismatic marine megafauna. The green turtle has a long history of human exploitation with some stocks extinct. Here we report on a 30-year study of the nesting abundance of the green turtle stock endemic to the Hawaiian Archipelago. We show that there has been a substantial long-term increase in abundance of this once seriously depleted stock following cessation of harvesting since the 1970s. This population increase has occurred in a far shorter period of time than previously thought possible. There was also a distinct 3-4 year periodicity in annual nesting abundance that might be a function of regional environmental stochasticity that synchronises breeding behaviour throughout the Archipelago. This is one of the few reliable long-term population abundance time series for a large long-lived marine species, which are needed for gaining insights into the recovery process of long-lived marine species and long-term ecological processes. (C) 2003 Elsevier Ltd. All rights reserved.

Of sheep and rain: large-scale population dynamics of the red kangaroo

1. We analysed time-series data from populations of red kangaroos (Macropus rufus, Desmarest) inhabiting four areas in the pastoral zone of South Australia. We formulated a set of a priori models to disentangle the relative effects of the covariates: rainfall, harvesting, intraspecific competition, and domestic herbivores, on kangaroo population-growth rate. 2. The statistical framework allowed for spatial variation in the growth-rate parameters, response to covariates, and environmental variability, as well as spatially correlated error terms due to shared environment. 3. The most parsimonious model included all covariates but no area-specific parameter values, suggesting that kangaroo densities respond in the same way to the covariates across the areas. 4. The temporal dynamics were spatially correlated, even after taking into account the potentially synchronizing effect of rainfall, harvesting and domestic herbivores. 5. Counter-intuitively, we found a positive rather than negative effect of domestic herbivore density on the population-growth rate of kangaroos. We hypothesize that this effect is caused by sheep and cattle acting as a surrogate for resource availability beyond rainfall. 6. Even though our system is well studied, we must conclude that approximating resources by surrogates such as rainfall is more difficult than previously thought. This is an important message for studies of consumer-resource systems and highlights the need to be explicit about population processes when analysing population patterns.

Accounting for management costs in sensitivity analyses of matrix population models

Traditional sensitivity and elasticity analyses of matrix population models have been used to p inform management decisions, but they ignore the economic costs of manipulating vital rates. For exam le, the growth rate of a population is often most sensitive to changes in adult survival rate, but this does not mean that increasing that rate is the best option for managing the population because it may be much more expensive than other options. To explore how managers should optimize their manipulation of vital rates, we incorporated the cost of changing those rates into matrix population models. We derived analytic expressions for locations in parameter space where managers should shift between management of fecundity and survival, for the balance between fecundity and survival management at those boundaries, and for the allocation of management resources to sustain that optimal balance. For simple matrices, the optimal budget allocation can often be expressed as simple functions of vital rates and the relative costs of changing them. We applied our method to management of the Helmeted Honeyeater (Lichenostomus melanops cassidix; an endangered Australian bird) and the koala (Phascolarctos cinereus) as examples. Our method showed that cost-efficient management of the Helmeted Honeyeater should focus on increasing fecundity via nest protection, whereas optimal koala management should focus on manipulating both fecundity and survival simultaneously, These findings are contrary to the cost-negligent recommendations of elasticity analysis, which would suggest focusing on managing survival in both cases. A further investigation of Helmeted Honeyeater management options, based on an individual-based model incorporating density dependence, spatial structure, and environmental stochasticity, confirmed that fecundity management was the most cost-effective strategy. Our results demonstrate that decisions that ignore economic factors will reduce management efficiency.