Contributions of demography and dispersal parameters to the spatial spread of a stage-structured insect invasion

Stage-structured models that integrate demography and dispersal can be used to identify points in the life cycle with large effects on rates of population spatial spread, information that is vital in the development of containment strategies for invasive species. Current challenges in the application of these tools include: (1) accounting for large uncertainty in model parameters, which may violate assumptions of ‘‘local’’ perturbation metrics such as sensitivities and elasticities, and (2) forecasting not only asymptotic rates of spatial spread, as is usually done, but also transient spatial dynamics in the early stages of invasion. We developed an invasion model for the Diaprepes root weevil (DRW; Diaprepes abbreviatus [Coleoptera: Curculionidae]), a generalist herbivore that has invaded citrus-growing regions of the United States. We synthesized data on DRW demography and dispersal and generated predictions for asymptotic and transient peak invasion speeds, accounting for parameter uncertainty. We quantified the contributions of each parameter toward invasion speed using a ‘‘global’’ perturbation analysis, and we contrasted parameter contributions during the transient and asymptotic phases. We found that the asymptotic invasion speed was 0.02–0.028 km/week, although the transient peak invasion speed (0.03– 0.045 km/week) was significantly greater. Both asymptotic and transient invasions speeds were most responsive to weevil dispersal distances. However, demographic parameters that had large effects on asymptotic speed (e.g., survival of early-instar larvae) had little effect on transient speed. Comparison of the global analysis with lower-level elasticities indicated that local perturbation analysis would have generated unreliable predictions for the responsiveness of invasion speed to underlying parameters. Observed range expansion in southern Florida (1992–2006) was significantly lower than the invasion speed predicted by the model. Possible causes of this mismatch include overestimation of dispersal distances, demographic rates, and spatiotemporal variation in parameter values. This study demonstrates that, when parameter uncertainty is large, as is often the case, global perturbation analyses are needed to identify which points in the life cycle should be targets of management. Our results also suggest that effective strategies for reducing spread during the asymptotic phase may have little effect during the transient phase. Includes Appendix.

Demography and evolution of regional inequality

This paper addresses the importance of life cycle aspects in explaining the evolution of regional income inequality. The analysis of household microdata organized in age cohorts shows that Brazilian regional income inequality has different dynamics across generations, with income convergence being observed only for the older generations. The larger income share of younger generations produces a low speed of convergence in the country. When retirement payments, pensions, and other government transfers are excluded from income, convergence is not observed even for the older generations.

Demography in eusocial Hymenoptera

In my doctoral thesis I investigated the evolution of demographic traits within eusocial Hymenoptera. In the social bees, wasps and ants, eusociality has a unique effect on life span evolution as female larvae with the same genetic background can develop through phenotypic plasticity to a queen or a worker with vastly diverging life-history traits. Ant queens belong to the longest-lived insect species, while workers in most species live only a fraction of the queen’s life span. The average colony size of a species is positively correlated with social complexity, division of labor and diverging morphological female phenotypes all of which also affect life span. Therefore the demographic traits of interest in this thesis were life span and colony size. To understand the evolution of worker life span I applied a trade-off model that includes both hierarchical levels important in eusocial systems, namely the colony- and the individual-level. I showed that the evolution of worker life span may be an adaptive trait on the colony level to optimize resource allocation and therefore fitness in response to different levels of extrinsic mortality. A shorter worker life span as a result of reduced resource investments under high levels of extrinsic mortality increases colony fitness. In a further study I showed that Lasius niger colonies produce different aging phenotypes throughout colony development. Smaller colonies which apply a different foraging strategy than larger colonies produced smaller workers, which in turn have a longer life span as compared to larger workers produced in larger colonies. With the switch to cooperative foraging in growing colonies individual workers become less important for the colony caused by their increasing redundancy. Alternatively a trade of between growth and life span may lead to the results found in this study. A further comparative analysis to study the effect of colony size on life span showed a correlation between queen and worker life span when colony size is taken into account. While neither worker nor queen life span was associated with colony size, the differences between queen and worker life span increase with larger average colony sizes across all eusocial Hymenoptera. As colony size affects both queen and worker life span, I aimed to understand which factors lead to the small colony sizes displayed by some ant species. I therefore analyzed per-capita productivity at different colony sizes of eight cavity dwelling ant species. Most colonies of the study species grew larger than optimal productivity predicted. Larger colony size was shown to increase colony homeostasis, the predictability of future productivity and in turn the survival probability of the colony. I also showed that species that deploy an individual foraging mode may circumvent the density dependent decline in foraging success by splitting the colony to several nest sites.

Habitat-specific demography and source-sink dynamics in a population of Siberian jays

P>1. There are a number of models describing population structure, many of which have the capacity to incorporate spatial habitat effects. One such model is the source-sink model, that describes a system where some habitats have a natality that is higher than mortality (source) and others have a mortality that exceeds natality (sink). A source can be maintained in the absence of migration, whereas a sink will go extinct. 2. However, the interaction between population dynamics and habitat quality is complex, and concerns have been raised about the validity of published empirical studies addressing source-sink dynamics. In particular, some of these studies fail to provide data on survival, a significant component in disentangling a sink from a low quality source. Moreover, failing to account for a density-dependent increase in mortality, or decrease in fecundity, can result in a territory being falsely assigned as a sink, when in fact, this density-dependent suppression only decreases the population size to a lower level, hence indicating a 'pseudo-sink'. 3. In this study, we investigate a long-term data set for key components of territory-specific demography (mortality and reproduction) and their relationship to habitat characteristics in the territorial, group-living Siberian jay (Perisoreus infaustus). We also assess territory-specific population growth rates (r), to test whether spatial population dynamics are consistent with the ideas of source-sink dynamics. 4. Although average mortality did not differ between sexes, habitat-specific mortality did. Female mortality was higher in older forests, a pattern not observed in males. Male mortality only increased with an increasing amount of open areas. Moreover, reproductive success was higher further away from human settlement, indicating a strong effect of human-associated nest predators. 5. Averaged over all years, 76% of the territories were sources. These territories generally consisted of less open areas, and were located further away from human settlement. 6. The source-sink model provides a tool for modelling demography in distinct habitat patches of different quality, which can aid in identifying key habitats within the landscape, and thus, reduce the risk of implementing unsound management decisions.

Evolutionary demography of long-lived monocarpic perennials: a time-lagged integral projection model

1. The evolution of flowering strategies (when and at what size to flower) in monocarpic perennials is determined by balancing current reproduction with expected future reproduction, and these are largely determined by size-specific patterns of growth and survival. However, because of the difficulty in following long-lived individuals throughout their lives, this theory has largely been tested using short-lived species (< 5 years). 2. Here, we tested this theory using the long-lived monocarpic perennial Campanula thyrsoides which can live up to 16 years. We used a novel approach that combined permanent plot and herb chronology data from a 3-year field study to parameterize and validate integral projection models (IPMs). 3. Similar to other monocarpic species, the rosette leaves of C. thyrsoides wither over winter and so size cannot be measured in the year of flowering. We therefore extended the existing IPM framework to incorporate an additional time delay that arises because flowering demography must be predicted from rosette size in the year before flowering. 4. We found that all main demographic functions (growth, survival probability, flowering probability and fecundity) were strongly size-dependent and there was a pronounced threshold size of flowering. There was good agreement between the predicted distribution of flowering ages obtained from the IPMs and that estimated in the field. Mostly, there was good agreement between the IPM predictions and the direct quantitative field measurements regarding the demographic parameters lambda, R-0 and T. We therefore conclude that the model captures the main demographic features of the field populations. 5. Elasticity analysis indicated that changes in the survival and growth function had the largest effect (c. 80%) on lambda and this was considerably larger than in short-lived monocarps. We found only weak selection pressure operating on the observed flowering strategy which was close to the predicted evolutionary stable strategy. 6. Synthesis. The extended IPM accurately described the demography of a long-lived monocarpic perennial using data collected over a relatively short period. We could show that the evolution of flowering strategies in short- and long-lived monocarps seem to follow the same general rules but with a longevity-related emphasis on survival over fecundity.

Who Counts? Demography of Swidden Cultivators in Southeast Asia

Swidden cultivators are often found as a distinct category of farmers in the literature, but rarely appear in population censuses or other national and regional classifications. This has led to a worldwide confusion on how many people are dependent on this form of agriculture. The most often cited number of 200–300 million dates back to the early 1970s, but the source is obscure. We assess available, published data from nine countries in Southeast Asia and conclude that on this basis it is not possible to provide a firm estimate of the number of swidden cultivators in the region. A conservative range of 14–34 million people engaged in swidden cultivation in the region is suggested, however. We argue that along with improved knowledge of swidden livelihoods, there is an urgent need to develop techniques that will allow for better estimates of swidden populations in order to secure appropriate rural development and poverty reduction in swidden areas.