Estimation of survival rate and extinction probability for stage-structured populations with overlapping life stages

The development of methods providing reliable estimates of demographic parameters (e. g., survival rates, fecundity) for wild populations is essential to better understand the ecology and conservation requirements of individual species. A number of methods exist for estimating the demographics of stage-structured populations, but inherent mathematical complexity often limits their uptake by conservation practitioners. Estimating survival rates for pond-breeding amphibians is further complicated by their complex migratory and reproductive behaviours, often resulting in nonobservable states and successive cohorts of eggs and tadpoles. Here we used comprehensive data on 11 distinct breeding toad populations (Bufo calamita) to clarify and assess the suitability of a relatively simple method [the Kiritani-Nakasuji-Manly (KNM) method] to estimate the survival rates of stage-structured populations with overlapping life stages. The study shows that the KNM method is robust and provides realistic estimates of amphibian egg and larval survival rates for species in which breeding can occur as a single pulse or over a period of several weeks. The study also provides estimates of fecundity for seven distinct toad populations and indicates that it is essential to use reliable estimates of fecundity to limit the risk of under- or overestimating the survival rates when using the KNM method. Survival and fecundity rates for B. calamita populations were then used to define population matrices and make a limited exploration of their growth and viability. The findings of the study recently led to the implementation of practical conservation measures at the sites where populations were most vulnerable to extinction. © 2010 The Society of Population Ecology and Springer.

A metapopulation model of species boundaries

The majority, if not all, species have a limited geographic range bounded by a distribution edge. Violent ecotones such as sea coasts clearly produce edges for many species; however such ecotones, while sufficient for the formation of an edge, are not always necessary. We demonstrate this by simulation in discrete time of a spatially structured finite size metapopulation subjected to a spatial gradient in per-unit-time population extinction probability together with spatially structured dispersal and recolonisation. We find that relatively sharp edges separating a homeland or main geographical range from an outland or zone of relatively sparse and ephemeral colonisation can form in gradual environmental gradients. The form and placing of the edge is an emergent property of the metapopulation dynamics. The sharpness of the edge declines with increasing dispersal distance, and is dependent on the relative scales of dispersal distance and gradient length. The space over which the edge develops is short relative to the potential species range. The edge is robust against changes in both the shape of the environmental gradient and to a lesser extent to alterations in the kind of dispersal operating. Persistence times in the absence of environmental gradients are virtually independent of the shape of the dispersal function describing migration. The common finding of bell shaped population density distributions across geographic ranges may occur without the strict necessity of a niche mediated response to a spatially autocorrelated environment.

Extinction risk from climate change

Climate change over the past ,30 years has produced numerous shifts in the distributions and abundances of species1,2 and has been implicated in one species-level extinction3. Using projections of species’ distributions for future climate scenarios, we assess extinction risks for sample regions that cover some 20% of the Earth’s terrestrial surface. Exploring three approaches in which the estimated probability of extinction shows a powerlaw relationship with geographical range size, we predict, on the basis of mid-range climate-warming scenarios for 2050, that 15–37% of species in our sample of regions and taxa will be ‘committed to extinction’. When the average of the three methods and two dispersal scenarios is taken, minimal climate-warming scenarios produce lower projections of species committed to extinction (,18%) than mid-range (,24%) and maximum change (,35%) scenarios. These estimates show the importance of rapid implementation of technologies to decrease greenhouse

Fitness reduction and potential extinction of wild populations of Atlantic salmon, Salmo salar, as a result of interactions with escaped farm salmon.

The high level of escapes from Atlantic salmon farms, up to two million fishes per year in the North Atlantic, has raised concern about the potential impact on wild populations. We report on a twogeneration experiment examining the estimated lifetime successes, relative to wild natives, of farm, F1 and F2 hybrids and BC1 backcrosses to wild and farm salmon. Offspring of farm and hybrids (i.e. all F1 , F2 and BC1 groups) showed reduced survival compared with wild salmon but grew faster as juveniles and displaced wild parr, which as a group were significantly smaller. Where suitable habitat for these emigrant parr is absent, this competition would result in reduced wild smolt production. In the experimental conditions, where emigrants survived downstream, the relative estimated lifetime success ranged from 2% (farm) to 89% (BC1 wild) of that of wild salmon, indicating additive genetic variation for survival . Wild salmon primarily returned to fresh water after one sea winter (1SW) but farm and hybrids produced proportionately more 2SW salmon. However, lower overall survival means that this would result in reduced recruitment despite increased 2SW fecundity. We thus demonstrate that interaction of farm with wild salmon results in lowered fitness, with repeated escapes causing cumulative fitness depression and potentially an extinction vortex in vulnerable populations.

Escape times for rigid Brownian rotators in a bistable potential from the time evolution of the Green function and the characteristic time of the probability evolution

The greatest relaxation time for an assembly of three- dimensional rigid rotators in an axially symmetric bistable potential is obtained exactly in terms of continued fractions as a sum of the zero frequency decay functions (averages of the Legendre polynomials) of the system. This is accomplished by studying the entire time evolution of the Green function (transition probability) by expanding the time dependent distribution as a Fourier series and proceeding to the zero frequency limit of the Laplace transform of that distribution. The procedure is entirely analogous to the calculation of the characteristic time of the probability evolution (the integral of the configuration space probability density function with respect to the position co-ordinate) for a particle undergoing translational diffusion in a potential; a concept originally used by Malakhov and Pankratov (Physica A 229 (1996) 109). This procedure allowed them to obtain exact solutions of the Kramers one-dimensional translational escape rate problem for piecewise parabolic potentials. The solution was accomplished by posing the problem in terms of the appropriate Sturm-Liouville equation which could be solved in terms of the parabolic cylinder functions. The method (as applied to rotational problems and posed in terms of recurrence relations for the decay functions, i.e., the Brinkman approach c.f. Blomberg, Physica A 86 (1977) 49, as opposed to the Sturm-Liouville one) demonstrates clearly that the greatest relaxation time unlike the integral relaxation time which is governed by a single decay function (albeit coupled to all the others in non-linear fashion via the underlying recurrence relation) is governed by a sum of decay functions. The method is easily generalized to multidimensional state spaces by matrix continued fraction methods allowing one to treat non-axially symmetric potentials, where the distribution function is governed by two state variables. (C) 2001 Elsevier Science B.V. All rights reserved.

Comment on "correlated noise in a logistic growth model"

We argue the results published by Bao-Quan Ai et al [Phys. Rev E 67, 022903 (2003)] on "correlated noise in a logistic growth model " are not correct. Their conclusion that for larger values of the correlation parameter, lambda, the cell population is peaked at x=0, which denotes the high extinction rate is also incorrect. We find the reverse behaviour corresponding to their results, that increasing lambda, promotes the stable growth of tumour cells. In particular, their results for steady-state probability, as a function of cell number, at different correlation strengths, presented in figures 1 and 2 show different behaviour than one would expect from the simple mathematical expression for the steady-state probability. Additionally, their interpretation at small values of cell number that the steady state probability increases as they increase the correlation parameter is also questionable. Another striking feature in their figures (1 and 3) is that for the same values of the parameter lambda and alpha, their simulation produces two different curves both qualitatively and quantitatively.

A DFT study of the chain growth probability in Fischer-Tropsch synthesis

The chain growth probability (alpha value) is one of the most significant parameters in Fischer-Tropsch (FT) synthesis. To gain insight into the chain growth probability, we systematically studied the hydrogenation and C-C coupling reactions with different chain lengths on the stepped Co(0001) surface using density functional theory calculations. Our findings elucidate the relationship between the barriers of these elementary reactions and the chain length. Moreover, we derived a general expression of the chain growth probability and investigated the behavior of the alpha value observed experimentally. The high methane yield results from the lower chain growth rate for C-1 + C-1 coupling compared with the other coupling reactions. After C-1, the deviation of product distribution in FT synthesis from the Anderson-Schulz-Flory distribution is due to the chain length-dependent paraffin/olefin ratio. (C) 2008 Elsevier Inc. All rights reserved.

Linear quadratic and tumour control probability modelling in external beam radiotherapy

The standard linear-quadratic (LQ) survival model for external beam radiotherapy is reviewed with particular emphasis on studying how different schedules of radiation treatment planning may be affected by different tumour repopulation kinetics. The LQ model is further examined in the context of tumour control probability (TCP) models. The application of the Zaider and Minerbo non-Poissonian TCP model incorporating the effect of cellular repopulation is reviewed. In particular the recent development of a cell cycle model within the original Zaider and Minerbo TCP formalism is highlighted. Application of this TCP cell-cycle model in clinical treatment plans is explored and analysed.

Constraining the age of Lateglacial and early Holocene pollen zones and tephra horizons in southern Sweden with Bayesian probability methods

The sediment sequence from Hasseldala port in southeastern Sweden provides a unique Lateglacial/early Holocene record that contains five different tephra layers. Three of these have been geochemically identified as the Borrobol Tephra, the Hasseldalen Tephra and the 10-ka Askja Tephra. Twenty-eight high-resolution C-14 measurements have been obtained and three different age models based on Bayesian statistics are employed to provide age estimates for the five different tephra layers. The chrono- and pollen stratigraphic framework supports the stratigraphic position of the Borrobol Tephra as found in Sweden at the very end of the Older Dryas pollen zone and provides the first age estimates for the Askja and Hasseldalen tephras. Our results, however, highlight the limitations that arise in attempting to establish a robust, chronologically independent lacustrine sequence that can be correlated in great detail to ice core or marine records. Radiocarbon samples are prone to error and sedimentation rates in lake basins may vary considerably due to a number of factors. Any type of valid and 'realistic' age model, therefore, has to take these limitations into account and needs to include this information in its prior assumptions. As a result, the age ranges for the specific horizons at Hasseldala port are large and calendar year estimates differ according to the assumptions of the age-model. Not only do these results provide a cautionary note for overdependence on one age-model for the derivation of age estimates for specific horizons, but they also demonstrate that precise correlations to other palaeoarchives to detect leads or lags is problematic. Given the uncertainties associated with establishing age-depth models for sedimentary sequences spanning the Lateglacial period, however, this exercise employing Bayesian probability methods represents the best possible approach and provides the most statistically significant age estimates for the pollen zone boundaries and tephra horizons. Copyright (C) 2006 John Wiley & Sons, Ltd.

Laplace transforms of probability distributions and their inversions are easy on logarithmic scales

It is shown that, when expressing arguments in terms of their logarithms, the Laplace transform of a function is related to the antiderivative of this function by a simple convolution. This allows efficient numerical computations of moment generating functions of positive random variables and their inversion. The application of the method is straightforward, apart from the necessity to implement it using high-precision arithmetics. In numerical examples the approach is demonstrated to be particularly useful for distributions with heavy tails, Such as lognormal, Weibull, or Pareto distributions, which are otherwise difficult to handle. The computational efficiency compared to other methods is demonstrated for an M/G/1 queueing problem.

Anomalous extinction behaviour towards the type Ia SN 2003cg

We present optical and near-infrared photometry and spectroscopy of the Type Ia SN 2003cg, which exploded in the nearby galaxy NGC 3169. The observations cover a period between -8.5 and +414 d post-maximum. SN 2003cg is a normal but highly reddened Type Ia event. Its B magnitude at maximum B-max = 15.94 +/- 0.04 and Delta m(15)(B)(obs) = 1.12 +/- 0.04 [Delta m(15)(B)(intrinsic) = 1.25 +/- 0.05]. Allowing R-V to become a free parameter within the Cardelli et al. extinction law, simultaneous matches to a range of colour curves of normal SNe Ia yielded E(B - V) = 1.33 +/- 0.11, and RV = 1.80 +/- 0.19. While the value obtained for R-V is small, such values have been invoked in the past, and may imply a grain size which is small compared with the average value for the local interstellar medium.