A test of statistical techniques for detecting density dependence in sequential censuses of animal populations

Summary Principal and reduced major axes, and Bulmer's (1975) tests have been suggested as methods for detecting the presence of density dependence in a series of population censuses that are unsuitable for analysis by alternative means e.g. by k-factor analysis. These alternative methods are tested using census data, some of which are previously unpublished, from natural populations known from independent evidence to be subject to density dependent processes. All the methods fail to detect density dependence reliably, irrespective of sample size and the dynamics of the population. We conclude that none of the methods tested is sufficiently reliable to be useful as a test of density dependence in sequential censues of animal populations.     

Ontogenetic shifts within the selfish herd: predation risk and foraging trade-offs change with age in colonial web-building spiders

We tested for the existence of density dependence in annual adult rangeland grasshopper (Orthoptera: Acrididae) data from Montana, USA (1951–1991). Statistical density dependence was, in the sense of a stochastic equilibrium or return tendency, detected in all of the grasshopper mean density time-series from the three major physiographic regions of the state, Northern Glaciated Plains, Southern Unglaciated Plains, and Western Mountains. Parameters were estimated for a model that described the stochastic equilibrium characteristics of regional mean densities. The analyses showed that rangeland grasshopper regional densities fluctuate according to gamma distribution with a mean of 6.1–6.3 grasshoppers per m². Further, when regions exhibit outbreaks, the resulting infestation period (duration of outbreak) is short, spanning only a few generations.     

Early survival in roe deer: causes and consequences of cohort variation in two contrasted populations

Time- and sex-specific summer survival of roe deer fawns was estimated using capture-mark-recapture methods in two enclosed populations living in contrasting conditions. The population of Trois Fontaines (eastern France) was roughly constant in size throughout the study period, while in Chizé (western France), the population experienced frequent summer droughts and numbers decreased continuously during the study. Early survival of fawns was low and highly variable over the years at both Chizé and Trois Fontaines, and demonstrated marked variations between cohorts that need to be taken into account when modelling roe deer population dynamics. In Trois Fontaines, fawn survival was positively correlated with early body growth and total rainfall in May and June. In Chizé, fawn survival decreased with increasing density and tended to increase with increasing rainfall in May and June and adult female body mass. These factors explained more than 75% of the variability in early survival observed in both populations. Variation between cohorts had different consequences for the two populations. At Trois Fontaines, cohort variation was limited to a numerical effect on early survival. However at Chizé, cohort variation was long-lasting and affected the phenotypic quality of survivors at later ages, and thereby future survival and breeding abilities (both numerical and quality effects). Male and female fawns had similar survival over their first summer in both populations. This result contrasts with the lower survival of young males often observed in ungulates. Two ultimate causes can be proposed to account for the low and variable survival of roe deer fawns over the first summer: the high energy expenditures incurred by does during each breeding attempt and/or the low absolute body size of newborn roe deer fawns. Received: 28 April 1997 / Accepted: 14 July 1997     

Temporal and spatial variation of larval parasitism in non-outbreaking populations of a folivorous moth

In order to assess the role of parasitoids in the regulation of non-outbreaking populations of Epirrita autumnata, a geometrid lepidopteran with outbreaking populations in northern Europe, we examined the temporal and spatial variation of larval parasitism in southwestern Finland during 6 successive years. The study was carried out on two spatial scales, among trees within sites of about 1 ha and among sites separated by distances of 2–10 km, using experimental and observational approaches respectively. The overall percent parasitism was independent of host density on both spatial scales, while temporally it fluctuated only little. Of the two main parasitoids, the commoner one, Protapanteles immunis, showed a variable response to host density on the larger spatial scale and negative density dependence on the smaller scale. Temporally, parasitism caused by this species was independent of host density. Another parasitoid, Phobocampe bicingulata, showed positive density dependence on the smaller spatial scale and had a variable response on the larger scale, but exhibited negative density dependence over time. The results of this study caution against drawing conclusions concerning population regulation on the grounds of spatial density dependence alone. Larval parasitoids apparently do not maintain low densities in the E. autumnata populations studied. However, they may suppress E. autumnata densities to a level low enough for density-dependent mortality factor(s) to become regulating. Among other mortality factors of E. autumnata, pupal predation has been found to be temporally positively density-dependent. Received: 19 October 1999 / Accepted: 10 January 2000     

Spatial scale and the detection of density dependence in spruce budworm outbreaks in eastern North America

Using two tests for direct density dependence and standard techniques of time series analysis, we identified density dependence in defoliation time series of the spruce budworm across its outbreak range in eastern North America over the years 1945–1988. We carried out analyses for the entire region and for grid cells of defoliation maps at five spatial scales created by aggregating the smallest grid cells. The rate of detection of direct density dependence, as assessed by two previously published methods, decreased with increasing spatial scale. Using both methods, density dependence was detected more frequently at the periphery of the outbreak range, where defoliation rate was lower. This result suggested that density-dependent regulation may be stronger in those areas. The first order autoregressive process was the basic model for defoliation dynamics overall and the most common model across spatial scales. Second-order processes were encountered much less frequently, and those commonly identified as resulting from delayed density dependence generally occurred across spatial scales at a rate expected by chance alone. Our results were similar to those of other published studies, which have found the detection of density dependence to decrease at larger spatial scales. The results also reinforced the importance of considering spatial scale when diagnosing population processes using time series of abundance for single species. Received: 26 December 1999 / Accepted: 17 March 2000     

Adaptive management of sika deer populations in Hokkaido, Japan: theory and practice

We investigated the utility of adaptive management (AM) in wildlife management, reviewing our experiences in applying AM to overabundant sika deer (Cervus nippon) populations in Hokkaido, Japan. The management goals of our program were: (1) to maintain the population at moderate density levels preventing population irruption, (2) to reduce damage to crops and forests, and (3) to sustain a moderate yield of hunting without endangering the population. Because of significant uncertainty in biological and environmental parameters, we designed a “feedback” management program based on controlling hunting pressure. Three threshold levels of relative population size and four levels of hunting pressure were configured, with a choice of four corresponding management actions. Under this program, the Hokkaido Government has been promoting aggressive female culling to reduce the sika deer population since 1998. We devised a harvest-based estimation for population size using relative population size and the number of deer harvested, and found that the 1993 population size (originally estimated by extrapolation of aerial surveys) had been underestimated. To reduce observation errors, a harvest-based Bayesian estimation was developed and the 1993 population estimate was again revised. Analyses of population trends and harvest data demonstrate that hunting is an important large-scale experiment to obtain reliable estimation of population size. A serious side effect of hunting on sika deer was inadvertent lead poisoning of large birds of prey. The prohibition of the use of lead bullets by the Hokkaido Government was successful in reducing the lead poisoning, but the problem still remains. Two case studies on sika population irruption show that the densities set by maximum sustainable yield may be too high to prevent damage to agriculture, forestry, and/or ecosystems. Threshold management based on feedback control is better for ecosystem management. Since volunteer hunters favor higher hunting efficiency in resource management (e.g., venison), it is necessary to support the development of professional hunters for culling operations for ecosystem management, where lower densities of deer should be set for target areas. Hunting as resource management and culling for ecosystem management should be synergistically combined under AM.     

Population and prehistory II: space-limited human populations in constant environments

We present a population model to examine the forces that determined the quality and quantity of human life in early agricultural societies where cultivable area is limited. The model is driven by the non-linear and interdependent relationships between the age distribution of a population, its behavior and technology, and the nature of its environment. The common currency in the model is the production of food, on which age-specific rates of birth and death depend. There is a single non-trivial equilibrium population at which productivity balances caloric needs. One of the most powerful controls on equilibrium hunger level is fertility control. Gains against hunger are accompanied by decreases in population size. Increasing worker productivity does increase equilibrium population size but does not improve welfare at equilibrium. As a case study we apply the model to the population of a Polynesian valley before European contact.     

Individuals,populations and the balance of nature: the question of persistence in ecology

Explaining the persistence of populations is an important quest in ecology, and is a modern manifestation of the balance of nature metaphor. Increasingly, however, ecologists see populations (and ecological systems generally) as not being in equilibrium or balance. The portrayal of ecological systems as “non-equilibrium” is seen as a strong alternative to deterministic or equilibrium ecology, but this approach fails to provide much theoretical or practical guidance, and warrants formalisation at a more fundamental level. This is available in adaptation theory, which allows population persistence to be explained as an epiphenomenon stemming from the maintenance, survival, movement and reproduction of individual organisms. These processes take place within a physicochemical and biotic environment that persists through structured annual cycles, but which is also spatiotemporally dynamic and subject to stochastic variation. The focus is thus shifted from the overproduction of offspring and the consequent density dependent population pressure thought to follow, to the adaptations and ecological circumstances that support those relatively few individuals that do survive.     

Detection and analysis of genetic variation in Hordeum spontaneum populations from Israel using RAPD markers

Randomly amplified polymorphic DNA (RAPD) markers were used to analyse genetic diversity within and between Hordeum spontaneum populations sampled from Israel. Nei's index of genetic differentiation was used to partition diversity into within and between population components. Fifty-seven per cent of the variation detected was partitioned within 10 H. spontaneum populations. Using principal component and multiple regression analysis, part of the variation detected between populations was seen to be associated with certain ecogeographical factors. Fifty-eight per cent of the distribution of the phenotypic frequencies of three RAPD phenotypes detected using a single primer in 20 H. spontaneum populations could be accounted for by four ecogeographical variables, suggesting adaptive variation at certain RAPD loci.     

Diurnal variation in sperm DNA fragmentation: analysis of 11,382 semen samples from two populations and in vivo animal experiments

ABSTRACT

Circadian rhythms have been found in some reproductive functions phenotypes but remain unclear for sperm DNA fragmentation index (DFI). The present study aims to investigate the diurnal variation of DFI in mice model and men sperm. Adult male mice were sacrificed for sperm DFI with Sperm Chromatin Structure Assay (SCSA) in 24 hours at 6 evenly distributed time points. A cosinor pattern of DFI was observed with a nadir at zeitgeber time 10 AM. In a community population with 630 semen samples collected between 8 AM and 20 PM, the temporal variation of DFI also fit a cosinor pattern with a ? 343° acrophase and a nadir at 11 AM (P = .031). In a reproductive-medical-center dataset of 10752 semen samples collected between 7 AM and 11 AM, the decreasing trend of DFI was also confirmed. For the males with multiple samples, intra-individual comparison between different timepoints was performed, and each consecutive hour after 7 AM was also associated with 2.5 (95% CI: ?1.0, 5.9)% lower DFI by SCSA or 4.9 (1.9, 7.8)% lower DFI by SCD. Our study reveals a daily diurnal variance in sperm DFI which may suggest a practical approach to get more qualified sperms for natural or assisted reproduction.

Abbreviations: BMI, Body mass index; DFI, DNA fragmentation index; MARHCS, Male Reproductive Health in the Chongqing College Students; RMC, Reproductive Medical Center; SCD, Sperm Chromatin Dispersion; SCSA, Sperm Chromatin Structure Assay.     

Genetic experiments with model populations: Fallacies in genetic analysis performed from samples of recombinants selected for two markers

Abstract: Samples selected for two markers (one from each parent) from the progeny of a genetic cross are altered both in the genotype and in the frequency distribution as compared to the original progeny population. The consequences of these alterations were analyzed in selected samples obtained from model progeny populations of hypothetical genetic crosses. In complete progeny populations, distribution of genotypes (pattern of genotypes), sequence of markers, and frequency distribution of individual genotypes exhibit genetically intelligent relationship to each other only under a sole condition: at the correct sequence of markers. In selected samples, the singularity of the interrelationship is relaxed. Consequences of the relaxation in the genetic analysis are multiple, the most striking of which is an insinuation of circularity for the linear gene map.     

Size variation in chewing lice <Emphasis Type="Italic">Docophorulus coarctatus</Emphasis>: how host size and louse population density vary together

Chewing lice of the species Docophorulus coarctatus were extracted from museum specimens of their host, the great grey shrike Lanius excubitor, by combing feathers from 36 freshly shot birds (shot between 1962 and 1974), and samples of ten individual lice (five female, five male) were randomly collected for measurements from each bird. Female lice were bigger than males for all studied measurements (P < 0.001 in all cases), although the size of both sexes obtained from individual hosts was positively correlated. The overall size of lice (derived from a principal components analysis) was positively correlated with the overall size of the avian host, and also with the population density of lice on the individual host. We suggest that variation in louse morphology is due to differences in selection pressure exerted by each host and by intraspecific competition due to conspecifics. This is, to the best of our knowledge, the first evidence that Harrison’s rule (parasites on larger host species are often bigger than those on smaller hosts) not only works in a multispecies comparison but also within a single host–single parasite system as well.     

r‐ and K‐selection in fluctuating populations is determined by the evolutionary trade‐off between two fitness measures: Growth rate and lifetime reproductive success

In a stable environment, evolution maximizes growth rates in populations that are not density regulated and the carrying capacity in the case of density regulation. In a fluctuating environment, evolution maximizes a function of growth rate, carrying capacity and environmental variance, tending to r‐selection and K‐selection under large and small environmental noise, respectively. Here we analyze a model in which birth and death rates depend on density through the same function but with independent strength of density dependence. As a special case, both functions may be linear, corresponding to logistic dynamics. It is shown that evolution maximizes a function of the deterministic growth rate r₀ and the lifetime reproductive success (LRS) R₀, both defined at small densities, as well as the environmental variance. Under large noise this function is dominated by r₀ and average lifetimes are small, whereas R₀ dominates and lifetimes are larger under small noise. Thus, K‐selection is closely linked to selection for large R₀ so that evolution tends to maximize LRS in a stable environment. Consequently, different quantities (r₀ and R₀) tend to be maximized at low and high densities, respectively, favoring density‐dependent changes in the optimal life history.     

Genetic analysis of fish genomes and populations: allozyme variation within and among Atlantic salmon from Downeast rivers of Maine

Analysis of genetic variation within and among samples of naturally produced Atlantic salmon ( n  = 372) from 7 Maine (U.S.A.) and one Canadian river were conducted based on 54 allozyme loci. Eight of the 54 loci proved polymorphic, and estimated heterozygosities ( H _S) based on all loci ranged from 0·012 to 0·026 (mean = 0·021, s . e . = 0·002). Only one of 56 tests revealed genotypic proportions that deviated significantly from Hardy–Weinberg expectations. Genetic distances ( D ) between samples ranged from 0·002 to 0·022. No obvious association existed between genetic and geographic distances. Cluster analysis of genetic distances revealed the Dennys River sample as the most differentiated when all samples were included in the analysis, though bootstrap support of the cluster analysis was generally weak. G ‐tests revealed significant differences in allele frequencies among samples at five of the polymorphic loci, and the G ‐value summed over all loci also indicated significant differences among samples. F _ST values indicated that 3·4% of the total genetic diversity was due to variability among samples, while 96·6% was due to variability within samples. These results indicate that the Atlantic salmon analyzed in this study had levels of genetic variability and differentiation among samples comparable to native populations from other areas collected across a similar geographic range.     

Density, climate and varying return points: an analysis of long-term population fluctuations in the threatened European tree frog

Experimental research has identified many putative agents of amphibian decline, yet the population-level consequences of these agents remain unknown, owing to lack of information on compensatory density dependence in natural populations. Here, we investigate the relative importance of intrinsic (density-dependent) and extrinsic (climatic) factors impacting the dynamics of a tree frog (Hyla arborea) population over 22 years. A combination of log-linear density dependence and rainfall (with a 2-year time lag corresponding to development time) explain 75% of the variance in the rate of increase. Such fluctuations around a variable return point might be responsible for the seemingly erratic demography and disequilibrium dynamics of many amphibian populations.     

Solar radiation stress in climbing snails: behavioural and intrinsic features define the Hsp70 level in natural populations of Xeropicta derbentina (Pulmonata)

Ectotherms from sunny and hot environments need to cope with solar radiation. Mediterranean land snails of the superfamily Helicoidea feature a behavioural strategy to escape from solar radiation-induced excessive soil heating by climbing up vertical objects. The height of climbing, and also other parameters like shell colouration pattern, shell orientation, shell size, body mass, actual internal and shell surface temperature, and the interactions between those factors may be expected to modulate proteotoxic effects in snails exposed to solar radiation and, thus, their stress response. Focussing on natural populations of Xeropicta derbentina, we conducted a ‘snapshot’ field study using the individual Hsp70 level as a proxy for proteotoxic stress. In addition to correlation analyses, an IT-model selection approach based on Akaike’s Information Criterion was applied to evaluate a set of models with respect to their explanatory power and to assess the relevance of each of the above-mentioned parameters for individual stress, by model averaging and parameter estimation. The analysis revealed particular importance of the individuals’ shell size, height above ground, the shell colouration pattern and the interaction height × orientation. Our study showed that a distinct set of behavioural traits and intrinsic characters define the Hsp70 level and that environmental factors and individual features strongly interact.     

On the evolution of heavy-metal tolerant populations in Armeria maritima: evidence from allozyme variation and reproductive barriers

The process of ecological differentiation leading to the evolution of heavy-metal tolerant populations in Armeria maritima was studied by comparing population genetic structure and pattern of gene flow between populations growing on heavy-metal contaminated against non-contaminated sites using allozyme markers. In addition the evolution of reproductive isolation among populations was studied by measuring pollen fertility in interpopulational hybrids. The allozyme data suggested that in A. maritima multiple independent evolutionary origins of heavy-metal tolerant populations have occurred in the absence of strong genetic bottlenecks. The pattern of gene flow among populations was consistent with the model of isolation by distance with considerable gene flow between neighbor populations, and no reduction of gene flow between tolerant and non-tolerant populations. Hence it appears that substantial gene flow has not hampered genetic differentiation, probably because of the high selection pressure for heavy-metal tolerance. The pattern of reproductive isolation among populations suggests that evolution of heavy-metal tolerant populations has not triggered the development of reproductive barriers against non-tolerant populations. However, partial reproductive isolation has occurred under geographic separation.