Dynamics of experimental populations of native and introduced blowflies (Diptera: Calliphoridae): mathematical modelling and the transition from asymptotic equilibrium to bounded oscillations

The equilibrium dynamics of native and introduced blowflies is modelled using a density-dependent model of population growth that takes into account important features of the life-history in these flies. A theoretical analysis indicates that the product of maximum fecundity and survival is the primary determinant of the dynamics. Cochliomyia macellaria, a blowfly native to the Americas and the introduced Chrysomya megacephala and Chrysomya putoria, differ in their dynamics in that the first species shows a damping oscillatory behavior leading to a one-point equilibrium, whereas in the last two species population numbers show a two-point limit cycle. Simulations showed that variation in fecundity has a marked effect on the dynamics and indicates the possibility of transitions from one-point equilibrium to bounded oscillations and aperiodic behavior. Variation in survival has much less influence on the dynamics.

A two-Dimensional finite volume method for variable density flow and solute transport through saturated-unsaturated media

Extensive groundwater withdrawal has resulted in a severe seawater intrusion problem in the Gooburrum aquifers at Bundaberg, Queensland, Australia. Better management strategies can be implemented by understanding the seawater intrusion processes in those aquifers. To study the seawater intrusion process in the region, a two-dimensional density-dependent, saturated and unsaturated flow and transport computational model is used. The model consists of a coupled system of two non-linear partial differential equations. The first equation describes the flow of a variable-density fluid, and the second equation describes the transport of dissolved salt. A two-dimensional control volume finite element model is developed for simulating the seawater intrusion into the heterogeneous aquifer system at Gooburrum. The simulation results provide a realistic mechanism by which to study the convoluted transport phenomena evolving in this complex heterogeneous coastal aquifer.

Studies on Plant Demography: Ranunculus Repens L., R. Bulbosus L. and R. Acris L.: III. A Mathematical Model Incorporating Multiple Modes of Reproduction

This paper presents a comparative population dynamics study of three closely related species of buttercups (Ranunculus repens, R. acris, and R. bulbosus). The study is based on an investigation of the behaviour of the seeds in soil under field conditions and a continuous monitoring of survival and reproduction of some 9000 individual plants over a period of 21/2 years in a coastal grassland in North Wales. The data were analysed with the help of an extension of Leslie's matrix method which makes possible an simultaneous treatment of vegetative and sexual reproduction. It was found that R. repens (a) depends more heavily on vegetative as compared with sexual reproduction, (b) shows indications of negatively density-dependent population regulation, and (c) exhibits little variation in population growth rates from site to site and from one year to the next. In contrast, R. bulbosus (a) depends exclusively on sexual reproduction, (b) shows indications of a positively density-dependent population behaviour, and (c) exhibits great variation in population growth rates from site to site and from one year to the next. R. acris exhibits an intermediate behaviour in all these respects. It is suggested that the attributes of R. repens are those expected of a species inhabiting a stable environment, while R. bulbosus exhibits some of the characteristics of a fugitive species.

A physics-based flexural phonon-dependent thermal conductivity model for single layer graphene

In this paper, we address a physics-based closed-form analytical model of flexural phonon-dependent diffusive thermal conductivity (kappa) of suspended rectangular single layer graphene sheet. A quadratic dependence of the out-of-plane phonon frequency, generally called flexural phonons, on the phonon wave vector has been taken into account to analyze the behavior of kappa at lower temperatures. Such a dependence has further been used for the determination of second-order three-phonon Umklapp and isotopic scatterings. We find that these behaviors in our model are best explained through the upper limit of Debye cut-off frequency in the second-order three-phonon Umklapp scattering of the long phonon waves that actually remove the thermal conductivity singularity by contributing a constant scattering rate at low frequencies and note that the out-of-plane Gruneisen parameter for these modes need not be too high. Using this, we clearly demonstrate that. follows a T-1.5 and T-2 law at lower and higher temperatures in the absence of isotopes, respectively. However in their presence, the behavior of kappa sharply deviates from the T-2 law at higher temperatures. The present geometry-dependent model of kappa is found to possess an excellent match with various experimental data over a wide range of temperatures which can be put forward for efficient electro-thermal analyses of encased/supported graphene.

Alpha decay half-lives of heavy nuclei within a generalized liquid drop model

Theoretical alpha-decay half-lives of the heaviest nuclei are calculated using the experimental Q value. The barriers in the quasi-molecular shape path is determined within a Generalized Liquid Drop Model (GLDM) and the WKB approximation is used. The results are compared with calculations using the Density-Dependent, M3Y (DDM3Y) effective interaction and the Viola-Seaborg-Sobiczewski (VSS) formulae. The calculations provide consistent estimates for the half-lives of the a decay chains of these superheavy elements. The experimental data stand between the GLDM calculations and VSS ones in the most time.

Probing the high density behavior of the symmetry energy by using proton elliptic flow

Based on the isospin- and momentum-dependent transport model IBUU04, we calculated the reaction of the Sn-132+Sn-124 systems in semi-central collisions at beam energies of 400/A MeV, 600/A MeV and 800/A MeV by adopting two different density dependent symmetry energies. It was found that the proton differential elliptic flow as a function of transverse momentum is quite sensitive to the density dependence of symmetry energy, especially for the considered beam energy range. Therefore the proton differential elliptic flow may be considered as a robust probe for investigating the high density behavior of symmetry energy in intermediate energy heavy ion collisions.

alpha decay half-lives of new superheavy nuclei within a generalized liquid drop model

The alpha decay half-lives of the recently produced isotopes of the 112, 114, 116 and 118 nuclei and decay products have been calculated in the quasi-molecular shape path using the experimental Q(alpha) value and a Generalized Liquid Drop Model including the proximity effects between nucleons in the neck or the gap between the nascent fragments. Reasonable estimates are obtained for the observed alpha decay half-lives. The results are compared with calculations using the Density-Dependent M3Y effective interaction and the Viola-Seaborg-Sobiczewski formulae. Generalized Liquid Drop Model predictions are provided for the alpha decay half-lives of other superheavy nuclei using the Finite Range Droplet Model Q(alpha) and compared with the values derived from the VSS formulae.

Functional responses of the intertidal amphipod Echinogammarus marinus: effects of prey supply, model selection and habitat complexity.

The influence of predation in structuring ecological communities can be informed by examining the shape and magnitude of the functional response of predators towards prey. We derived functional responses of the ubiquitous intertidal amphipod Echinogammarus marinus towards one of its preferred prey species, the isopod Jaera nordmanni. First, we examined the form of the functional response where prey were replaced following consumption, as compared to the usual experimental design where prey density in each replicate is allowed to deplete. E. marinus exhibited Type II functional responses, i.e. inversely density-dependent predation of J. nordmanni that increased linearly with prey availability at low densities, but decreased with further prey supply. In both prey replacement and non-replacement experiments, handling times and maximum feeding rates were similar. The non-replacement design underestimated attack rates compared to when prey were replaced. We then compared the use of Holling’s disc equation (assuming constant prey density) with the more appropriate Rogers’ random predator equation (accounting for prey depletion) using the prey non-replacement data. Rogers’ equation returned significantly greater attack rates but lower maximum feeding rates, indicating that model choice has significant implications for parameter estimates. We then manipulated habitat complexity and found significantly reduced predation by the amphipod in complex as opposed to simple habitat structure. Further, the functional response changed from a Type II in simple habitats to a sigmoidal, density-dependent Type III response in complex habitats, which may impart stability on the predator−prey interaction. Enhanced habitat complexity returned significantly lower attack rates, higher handling times and lower maximum feeding rates. These findings illustrate the sensitivity of the functional response to variations in prey supply, model selection and habitat complexity and, further, that E. marinus could potentially determine the local exclusion and persistence of prey through habitat-mediated changes in its predatory functional responses.

Aggregation of parasitism risk in an aphid-parasitoid system: Effects of plant patch size and aphid density

Few studies have linked density dependence of parasitism and the tritrophic environment within which a parasitoid forages. In the non-crop plant-aphid, Centaurea nigra-Uroleucon jaceae system, mixed patterns of density-dependent parasitism by the parasitoids Aphidius funebris and Trioxys centaureae were observed in a survey of a natural population. Breakdown of density-dependent parasitism revealed that density dependence was inverse in smaller colonies but direct in large colonies (>20 aphids), suggesting there is a threshold effect in parasitoid response to aphid density. The CV2 of searching parasitoids was estimated from parasitism data using a hierarchical generalized linear model, and CV2>1 for A. funebris between plant patches, while for T. centaureae CV2>1 within plant patches. In both cases, density independent heterogeneity was more important than density-dependent heterogeneity in parasitism. Parasitism by T. centaureae increased with increasing plant patch size. Manipulation of aphid colony size and plant patch size revealed that parasitism by A. funebris was directly density dependent at the range of colony sizes tested (50-200 initial aphids), and had a strong positive relationship with plant patch size. The effects of plant patch size detected for both species indicate that the tritrophic environment provides a source of host density independent heterogeneity in parasitism, and can modify density-dependent responses. (c) 2007 Gessellschaft fur Okologie. Published by Elsevier GmbH. All rights reserved.

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.

Stochastic dynamics in exotic and native blowflies: an analysis combining laboratory experiments and a two-patch metapopulation model

In this study we explored the stochastic population dynamics of three exotic blowfly species, Chrysomya albiceps, Chrysomya megacephala and Chrysomya putoria, and two native species, Cochliomyia macellaria and Lucilia eximia, by combining a density-dependent growth model with a two-patch metapopulation model. Stochastic fecundity, survival and migration were investigated by permitting random variations between predetermined demographic boundary values based on experimental data. Lucilia eximia and Chrysomya albiceps were the species most susceptible to the risk of local extinction. Cochliomyia macellaria, C. megacephala and C. putoria exhibited lower risks of extinction when compared to the other species. The simultaneous analysis of stochastic fecundity and survival revealed an increase in the extinction risk for all species. When stochastic fecundity, survival and migration were simulated together, the coupled populations were synchronized in the five species. These results are discussed, emphasizing biological invasion and interspecific interaction dynamics.

EVOLUTION OF DOMINANCE UNDER FREQUENCY-DEPENDENT INTRASPECIFIC COMPETITION IN AN ASSORTATIVELY MATING POPULATION

We study the evolution of higher levels of dominance as a response to negative frequency-dependent selection. In contrast to previous studies, we focus on the effect of assortative mating on the evolution of dominance under frequency-dependent intraspecific competition. We analyze a two-locus two-allele model, in which the primary locus has a major effect on a quantitative trait that is under a mixture of frequency-independent stabilizing selection, density-dependent selection, and frequency-dependent selection caused by intraspecific competition for a continuum of resources. The second (modifier) locus determines the degree of dominance at the trait level. Additionally, the population mates assortatively with respect to similarities in the ecological trait. Our analysis shows that the parameter region in which dominance can be established decreases if small levels of assortment are introduced. In addition, the degree of dominance that can be established also decreases. In contrast, if assortment is intermediate, sexual selection for extreme types can be established, which leads to evolution of higher levels of dominance than under random mating. For modifiers with large effects, intermediate levels of assortative mating are most favorable for the evolution of dominance. For large modifiers, the speed of fixation can even be higher for intermediate levels of assortative mating than for random mating.

Microevolutionary responses in experimental populations of plants to CO2-enriched environments: parallel results from two model systems.

Despite the critical role that terrestrial vegetation plays in the Earth's carbon cycle, very little is known about the potential evolutionary responses of plants to anthropogenically induced increases in concentrations of atmospheric CO2. We present experimental evidence that rising CO2 concentration may have a direct impact on the genetic composition and diversity of plant populations but is unlikely to result in selection favoring genotypes that exhibit increased productivity in a CO2-enriched atmosphere. Experimental populations of an annual plant (Abutilon theophrasti, velvetleaf) and a temperate forest tree (Betula alleghaniensis, yellow birch) displayed responses to increased CO2 that were both strongly density-dependent and genotype-specific. In competitive stands, a higher concentration of CO2 resulted in pronounced shifts in genetic composition, even though overall CO2-induced productivity enhancements were small. For the annual species, quantitative estimates of response to selection under competition were 3 times higher at the elevated CO2 level. However, genotypes that displayed the highest growth responses to CO2 when grown in the absence of competition did not have the highest fitness in competitive stands. We suggest that increased CO2 intensified interplant competition and that selection favored genotypes with a greater ability to compete for resources other than CO2. Thus, while increased CO2 may enhance rates of selection in populations of competing plants, it is unlikely to result in the evolution of increased CO2 responsiveness or to operate as an important feedback in the global carbon cycle. However, the increased intensity of selection and drift driven by rising CO2 levels may have an impact on the genetic diversity in plant populations.

Exposing extinction risk analysis to pathogens: Is disease just another form of density dependence?

In the United States and several other countries., the development of population viability analyses (PVA) is a legal requirement of any species survival plan developed for threatened and endangered species. Despite the importance of pathogens in natural populations, little attention has been given to host-pathogen dynamics in PVA. To study the effect of infectious pathogens on extinction risk estimates generated from PVA, we review and synthesize the relevance of host-pathogen dynamics in analyses of extinction risk. We then develop a stochastic, density-dependent host-parasite model to investigate the effects of disease on the persistence of endangered populations. We show that this model converges on a Ricker model of density dependence under a suite of limiting assumptions, including. a high probability that epidemics will arrive and occur. Using this modeling framework, we then quantify: (1) dynamic differences between time series generated by disease and Ricker processes with the same parameters; (2) observed probabilities of quasi-extinction for populations exposed to disease or self-limitation; and (3) bias in probabilities of quasi-extinction estimated by density-independent PVAs when populations experience either form of density dependence. Our results suggest two generalities about the relationships among disease, PVA, and the management of endangered species. First, disease more strongly increases variability in host abundance and, thus, the probability of quasi-extinction, than does self-limitation. This result stems from the fact that the effects and the probability of occurrence of disease are both density dependent. Second, estimates of quasi-extinction are more often overly optimistic for populations experiencing disease than for those subject to self-limitation. Thus, although the results of density-independent PVAs may be relatively robust to some particular assumptions about density dependence, they are less robust when endangered populations are known to be susceptible to disease. If potential management actions involve manipulating pathogens, then it may be useful to. model disease explicitly.

Unidirectional inversion of ibuprofen in Caco-2 cells: Developing a suitable model for presystemic chiral inversion study

Intestinal chiral inversion of ibuprofen is still lacking direct evidence. In a preliminary experiment, ibuprofen was found to undergo inversion in Caco-2 cells. This investigation was thus conducted to determine the characteristics and influence of some biochemical factors on the chiral inversion of ibuprofen in Caco-2 cells. The effects of substrate concentration (2.5-40 mu g/ml), cell density (0.5-2 x 10(6) cells/ well), content of serum (0-20%), coexistence of S ibuprofen (corresponding doses), sodium azide (10mm), exogenous Coenzyme A (CoA) (0.1 - 0.4 mm),. and palmitic acid (5-25 mu m) on inversion were examined. A stereoselective HPLC method based on the Chromasil-CHI-TBB column was developed for quantitative analysis of the drug in cell culture medium. The inversion ratio (F-i) and elimination rate constant were calculated as the indexes of inversion extent. Inversion of ibuprofen in Caeo-2 cells was found to be both dose and cell density dependent, indicating saturable characteristics. Addition of serum significantly inhibited the inversion, to an extent of 2.7 fold decrease at 20% content. Preexistence of S enantiomer exerted a significant inhibitory effect (p < 0.01 for all tests). Sodium azide decreased the inversion ratio from 0.43 to 0.32 (p < 0.01). Exogenous CoA and palmitic acid significantly promoted the inversion at all tested doses (p < 0.01 for all tests). This research provided strong evidence to the capacity and capability of intestinal chiral inversion. Although long incubation times up to 120 h were required, Caco-2 cells should be a suitable model for chiral inversion research of 2-APAs considering the human-resourced and well-defined characteristics from the present study.