Stability of a Seabird Population in the Presence of an Introduced Predator

We hypothesized that although large populations may appear able to withstand predation and disturbance, added stochasticity in population growth rate (λ) increases the risk of dramatic population declines. Approximately half of the Aleutian Islands' population of Least Auklets (Aethia pusilla) breed at one large colony at Kiska Island in the presence of introduced Norway rats (Rattus norvegicus) whose population erupts periodically. We evaluated two management plans, do nothing or eradicate rats, for this colony, and performed stochastic elasticity analysis to focus future research and management. Our results indicated that Least Auklets breeding at Kiska Island had the lowest absolute value of growth rate and more variable λ's (neither statistically significant) during 2001-2010, when compared with rat-free colonies at Buldir and Kasatochi islands. We found variability in the annual proportional change in population size among islands with Kiska Island having the fastest rate of decline, 78% over 20 years. Under the assumption that the eradication of rats would result in vital rates similar to those observed at rat-free Buldir and Kasatochi islands, we found the projected population decline decreased from 78% to 24% over 20 years. Overall, eradicating rats at Kiska Island is not likely to increase Least Auklet vital rates, but will decrease the amount of variation in λ, resulting in a significantly slower rate of population decline. We recommend the eradication of rats from Kiska Island to decrease the probability of dramatic population declines and ensure the future persistence of this important colony.

Stochastic modelling of rainfall from satellite data

Satellite-based rainfall monitoring is widely used for climatological studies because of its full global coverage but it is also of great importance for operational purposes especially in areas such as Africa where there is a lack of ground-based rainfall data. Satellite rainfall estimates have enormous potential benefits as input to hydrological and agricultural models because of their real time availability, low cost and full spatial coverage. One issue that needs to be addressed is the uncertainty on these estimates. This is particularly important in assessing the likely errors on the output from non-linear models (rainfall-runoff or crop yield) which make use of the rainfall estimates, aggregated over an area, as input. Correct assessment of the uncertainty on the rainfall is non-trivial as it must take account of • the difference in spatial support of the satellite information and independent data used for calibration • uncertainties on the independent calibration data • the non-Gaussian distribution of rainfall amount • the spatial intermittency of rainfall • the spatial correlation of the rainfall field This paper describes a method for estimating the uncertainty on satellite-based rainfall values taking account of these factors. The method involves firstly a stochastic calibration which completely describes the probability of rainfall occurrence and the pdf of rainfall amount for a given satellite value, and secondly the generation of ensemble of rainfall fields based on the stochastic calibration but with the correct spatial correlation structure within each ensemble member. This is achieved by the use of geostatistical sequential simulation. The ensemble generated in this way may be used to estimate uncertainty at larger spatial scales. A case study of daily rainfall monitoring in the Gambia, west Africa for the purpose of crop yield forecasting is presented to illustrate the method.

Introduction. Stochastic physics and climate modelling

Finite computing resources limit the spatial resolution of state-of-the-art global climate simulations to hundreds of kilometres. In neither the atmosphere nor the ocean are small-scale processes such as convection, clouds and ocean eddies properly represented. Climate simulations are known to depend, sometimes quite strongly, on the resulting bulk-formula representation of unresolved processes. Stochastic physics schemes within weather and climate models have the potential to represent the dynamical effects of unresolved scales in ways which conventional bulk-formula representations are incapable of so doing. The application of stochastic physics to climate modelling is a rapidly advancing, important and innovative topic. The latest research findings are gathered together in the Theme Issue for which this paper serves as the introduction.

Stochastic resonance in a nonlinear model of a rotating, stratified shear flow, with a simple stochastic inertia-gravity wave parameterization

We report on a numerical study of the impact of short, fast inertia-gravity waves on the large-scale, slowly-evolving flow with which they co-exist. A nonlinear quasi-geostrophic numerical model of a stratified shear flow is used to simulate, at reasonably high resolution, the evolution of a large-scale mode which grows due to baroclinic instability and equilibrates at finite amplitude. Ageostrophic inertia-gravity modes are filtered out of the model by construction, but their effects on the balanced flow are incorporated using a simple stochastic parameterization of the potential vorticity anomalies which they induce. The model simulates a rotating, two-layer annulus laboratory experiment, in which we recently observed systematic inertia-gravity wave generation by an evolving, large-scale flow. We find that the impact of the small-amplitude stochastic contribution to the potential vorticity tendency, on the model balanced flow, is generally small, as expected. In certain circumstances, however, the parameterized fast waves can exert a dominant influence. In a flow which is baroclinically-unstable to a range of zonal wavenumbers, and in which there is a close match between the growth rates of the multiple modes, the stochastic waves can strongly affect wavenumber selection. This is illustrated by a flow in which the parameterized fast modes dramatically re-partition the probability-density function for equilibrated large-scale zonal wavenumber. In a second case study, the stochastic perturbations are shown to force spontaneous wavenumber transitions in the large-scale flow, which do not occur in their absence. These phenomena are due to a stochastic resonance effect. They add to the evidence that deterministic parameterizations in general circulation models, of subgrid-scale processes such as gravity wave drag, cannot always adequately capture the full details of the nonlinear interaction.

QUAGMIRE v1.3: a quasi-geostrophic model for investigating rotating fluids experiments

QUAGMIRE is a quasi-geostrophic numerical model for performing fast, high-resolution simulations of multi-layer rotating annulus laboratory experiments on a desktop personal computer. The model uses a hybrid finite-difference/spectral approach to numerically integrate the coupled nonlinear partial differential equations of motion in cylindrical geometry in each layer. Version 1.3 implements the special case of two fluid layers of equal resting depths. The flow is forced either by a differentially rotating lid, or by relaxation to specified streamfunction or potential vorticity fields, or both. Dissipation is achieved through Ekman layer pumping and suction at the horizontal boundaries, including the internal interface. The effects of weak interfacial tension are included, as well as the linear topographic beta-effect and the quadratic centripetal beta-effect. Stochastic forcing may optionally be activated, to represent approximately the effects of random unresolved features. A leapfrog time stepping scheme is used, with a Robert filter. Flows simulated by the model agree well with those observed in the corresponding laboratory experiments.

Scale-invariant moving finite elements for nonlinear partial differential equations in two dimensions

A scale-invariant moving finite element method is proposed for the adaptive solution of nonlinear partial differential equations. The mesh movement is based on a finite element discretisation of a scale-invariant conservation principle incorporating a monitor function, while the time discretisation of the resulting system of ordinary differential equations is carried out using a scale-invariant time-stepping which yields uniform local accuracy in time. The accuracy and reliability of the algorithm are successfully tested against exact self-similar solutions where available, and otherwise against a state-of-the-art h-refinement scheme for solutions of a two-dimensional porous medium equation problem with a moving boundary. The monitor functions used are the dependent variable and a monitor related to the surface area of the solution manifold. (c) 2005 IMACS. Published by Elsevier B.V. All rights reserved.

Effects of chlorogenic acid and bovine serum albumin on the oxidative stability of low density lipoproteins in vitro

The ability of chlorogenic acid to inhibit oxidation of human low-density lipoprotein (LDL) was studied by in vitro copper-induced LDL oxidation. The effect of chlorogenic acid on the lag time before LDL oxidation increased in a dose dependent manner by up to 176% of the control value when added at concentrations of 0.25 -1.0 μM. Dose dependent increases in lag time of LDL oxidation were also observed, but at much higher concentrations, when chlorogenic acid was incubated with LDL (up to 29.7% increase in lag phase for 10 μM chlorogenic acid) or plasma (up to 16.6% increase in lag phase for 200 μM chlorogenic acid) prior to isolation of LDL, and this indicated that chlorogenic acid was able to bind, at least weakly, to LDL. Bovine serum albumin (BSA) increased the oxidative stability of LDL in the presence of chlorogenic acid. Fluorescence spectroscopy showed that chlorogenic acid binds to BSA with a binding constant of 3.88 x 104 M-1. BSA increased the antioxidant effect of chlorogenic acid, and this was attributed to copper ions binding to BSA, thereby reducing the amount of copper available for inducing lipid peroxidation.

Pilot study of the short-term physico-chemical stability of atenolol tablets stored in a multi-compartment compliance aid

Study objectives: There is a possibility that lower air, moisture and light protection could impact on physico-chemical stability of medicines inside multi-compartment compliance aids (MCCAs), although this has not yet been proved. The objectives of the study were to examine the physico-chemical stability of atenolol tablets stored in a compliance aid at room temperature, and at elevated temperature and humidity to simulate practice conditions. Methods: Atenolol 100 mg tablets in 28-chamber, plastic compliance aids with transparent lids were stored for four weeks at room temperature and at 40°C with 75% relative humidity. Tablets were also stored at room temperature in original packaging and Petri dishes. Physical tests were conducted to standards as laid down in the British Pharmacopoeia 2005, and dissolution to those of the United States Pharmacopoeia volume 24. Chemical stability was assessed by a validated high-performance liquid chromatography (HPLC) method. Results: Tablets at room temperature in original packaging, in compliance aids and Petri dishes remained the same in appearance and passed physico-chemical tests. Tablets exposed to 40°C with 75% relative humidity in compliance aids passed tests for uniformity of weight, friability and chemical stability but became pale and moist, softer (82 newtons ± 4; p< 0.0001) than tablets in the original packaging (118 newtons ± 6), more friable (0.14% loss of mass) compared with other tablets (0.005%), and failed the tests for disintegration (>15 minutes) and dissolution (only 15% atenolol released at 30 minutes). Conclusion: Although chemical stability was unaffected, storage in compliance aids at 40°C with 75% relative humidity softened atenolol tablets, prolonged disintegration time and hindered dissolution which could significantly reduce bioavailability. This formulation could be suitable for storage in compliance aids at 25°C, but not in hotter, humid weather.

Well-posed boundary value problems for integrable evolution equations on a finite interval

We consider boundary value problems posed on an interval [0,L] for an arbitrary linear evolution equation in one space dimension with spatial derivatives of order n. We characterize a class of such problems that admit a unique solution and are well posed in this sense. Such well-posed boundary value problems are obtained by prescribing N conditions at x=0 and n–N conditions at x=L, where N depends on n and on the sign of the highest-degree coefficient n in the dispersion relation of the equation. For the problems in this class, we give a spectrally decomposed integral representation of the solution; moreover, we show that these are the only problems that admit such a representation. These results can be used to establish the well-posedness, at least locally in time, of some physically relevant nonlinear evolution equations in one space dimension.

Tolerance traits and the stability of mutualism

Identifying factors which allow the evolution and persistence of cooperative interactions between species is a fundamental issue in evolutionary ecology. Various hypotheses have been suggested which generally focus on mechanisms that allow cooperative genotypes in different species to maintain interactions over space and time. Here, we emphasise the fact that even within mutualisms (interactions with net positive fitness effects for both partners), there may still be inherent costs, such as the occasional predation by ants upon aphids. Individuals engaged in mutualisms benefit from minimising these costs as long as it is not at the expense of breaking the interspecific interaction, which offers a net positive benefit. The most common and obvious defence traits to minimise interspecific interaction costs are resistance traits, which act to reduce encounter rate between two organisms. Tolerance traits, in contrast, minimise fitness costs to the actor, but without reducing encounter rate. Given that, by definition, it is beneficial to remain in mutualistic interactions, the only viable traits to minimise costs are tolerance-based 'defence' strategies. Thus, we propose that tolerance traits are an important factor promoting stability in mutualisms. Furthermore, because resistance traits tend to propagate coevolutionary arms races between antagonists, whilst tolerance traits do not, we also suggest that tolerance-based defence strategies may be important in facilitating the transition from antagonistic interactions into mutualisms. For example, the mutualism between ants and aphids has been suggested to have evolved from parasitism. We describe how phenotypic plasticity in honeydew production may be a tolerance trait that has prevented escalation into an antagonistic arms race and instead led to mutualistic coevolution.

Longevity, early emergence and body size in a pollinating fig wasp - implications for stability in a fig-pollinator mutualism

1. Fig trees (Ficus) are pollinated only by agaonid wasps, whose larvae also gall fig ovules. Each ovule develops into either a seed (when pollinated) or a wasp (when an egg is also laid inside) but not both. 2. Ovipositing wasps (foundresses) favour ovules near the centre of the enclosed inflorescence (syconium or 'fig'), leaving ovules near the outer wall to develop into seeds. This spatial stratification of wasps and seeds ensures reproduction in both partners, and thereby enables mutualism persistence. However, the mechanism(s) responsible remain(s) unknown. 3. Theory shows that foundresses will search for increasingly rare inner ovules and ignore outer ovules, as long as ovipositing in outer ovules is sufficiently slow and/or if inner ovules confer greater fitness to wasps. The fig-pollinator mutualism can therefore be stabilized by strong time constraints on foundresses and by offspring fitness gradients over variation in ovule position. 4. Female fig wasps cannot leave their galls without male assistance. We found that females in outer ovules were unlikely to be released. Inner ovules thus have added value to foundresses, because their female offspring are more likely to mate and disperse. 5. For those offspring that did emerge, gall position (inner/outer) and body size did not influence the order in which female pollinators exited syconia, nor did early emerging wasps enjoy increased life spans. 6. We also found that the life spans of female wasps nearly doubled when given access to moisture. We suggest that conflict resolution in the fig-pollinator mutualism may thus be influenced by tropical seasonality, because wasps may be less able to over-exploit ovules in dry periods due to time constraints.

Effect of minerals on casein micelle stability of cows' milk

The effects of minerals on casein micelle stability of individual cows' milk, throughout a complete lactation, were investigated. Calcium and calcium ions, magnesium, phosphorus, sodium, potassium and citrate contents were analysed, together with the following physical properties of milk; pH, ethanol stability, rennet clotting time and coagulum firmness. There was an inverse non-linear relationship between free calcium ion concentration and ethanol stability (ES; r=0.84). Rennet coagulation time showed a weaker relationship with free calcium ion concentration (r=0.44) but a stronger relationship with pH (r=0.66). In addition, samples containing higher amounts of free calcium ions produced a firmer gel. Citrate in natural samples acts as a stabilizing factor, as it slightly improves milk stability. Potassium, on the other hand, exhibited a negative correlation, but only with rennet clotting time (r=-0.52). Throughout lactation the average values were; free Ca 21 concentration 1.88 mm,, pH 6.63, ES 83.2% and clotting time 13.6 min. The equilibrium relationship between pH and free Ca2+ concentration was investigated by adjusting milk pH from 5.9 to 7.1, using acid and alkali. There was a good inverse linear relationship between pH and log (free Ca 21) for individual milk samples, with a gradient of -0.62 and a standard deviation of 0.042.

Association of n-3 polyunsaturated fatty acids with stability of atherosclerotic plaques: a randomised controlled trial

Background N-3 polyunsaturated fatty acids (PUFAs) from oily fish protect against death from cardiovascular disease. We aimed to assess the hypothesis that incorporation of n-3 and n-6 PUFAs into advanced atherosclerotic plaques increases and decreases plaque stability, respectively. Methods We did a randomised controlled trial of patients awaiting carotid endarterectomy. We randomly allocated patients control, sunflower oil (n-6), or fish-oil (n-3) capsules until surgery. Primary outcome was plaque morphology indicative of stability or instability, and outcome measures were concentrations of EPA, DHA, and linoleic acid in carotid plaques; plaque morphology; and presence of macrophages in plaques. Analysis was per protocol. Findings 188 patients were enrolled and randomised; 18 withdrew and eight were excluded. Duration of oil treatment was 7-189 days (median 42) and did not differ between groups. The proportions of EPA and DHA were higher in carotid plaque fractions in patients receiving fish oil compared with those receiving control (absolute difference 0.5 [95% CI 0.3-0.7], 0.4 [0.1-0.6], and 0.2 [0.1-0.4] g/100 g total fatty acids for EPA; and 0.3 [0.0-0.8], 0.4 [0.1-0.7], and 0.3 [0.1-0.6] g/100 g total fatty acids for DHA; in plaque phospholipids, cholesteryl esters, and triacylglycerols, respectively). Sunflower oil had little effect on the fatty acid composition of lipid fractions. Fewer plaques from patients being treated with fish oil had thin fibrous caps and signs of inflammation and more plaques had thick fibrous caps and no signs of inflammation, compared with plaques in patients in the control and sunflower oil groups (odds ratio 0.52 [95% CI 0.24-0.89] and 1.19 [1.02-1.57] vs control; 0.49 [0.23-0.90] and 1.16 [1.01-1.53] vs sunflower oil). The number of macrophages in plaques from patients receiving fish oil was lower than in the other two groups. Carotid plaque morphology and infiltration by macrophages did not differ between control and sunflower oil groups. Interpretation Atherosclerotic plaques readily incorporate n-3 PUFAs from fish-oil supplementation, inducing changes that can enhance stability of atherosclerotic plaques. By contrast, increased consumption of n-6 PUFAs does not affect carotid plaque fatty-acid composition or stability over the time course studied here. Stability of plaques could explain reductions in non-fatal and fatal cardiovascular events associated with increased n-3 PUFA intake.

Empirical mode decomposition: a novel technique for the study of tremor time series

Tremor is a clinical feature characterized by oscillations of a part of the body. The detection and study of tremor is an important step in investigations seeking to explain underlying control strategies of the central nervous system under natural (or physiological) and pathological conditions. It is well established that tremorous activity is composed of deterministic and stochastic components. For this reason, the use of digital signal processing techniques (DSP) which take into account the nonlinearity and nonstationarity of such signals may bring new information into the signal analysis which is often obscured by traditional linear techniques (e.g. Fourier analysis). In this context, this paper introduces the application of the empirical mode decomposition (EMD) and Hilbert spectrum (HS), which are relatively new DSP techniques for the analysis of nonlinear and nonstationary time-series, for the study of tremor. Our results, obtained from the analysis of experimental signals collected from 31 patients with different neurological conditions, showed that the EMD could automatically decompose acquired signals into basic components, called intrinsic mode functions (IMFs), representing tremorous and voluntary activity. The identification of a physical meaning for IMFs in the context of tremor analysis suggests an alternative and new way of detecting tremorous activity. These results may be relevant for those applications requiring automatic detection of tremor. Furthermore, the energy of IMFs was visualized as a function of time and frequency by means of the HS. This analysis showed that the variation of energy of tremorous and voluntary activity could be distinguished and characterized on the HS. Such results may be relevant for those applications aiming to identify neurological disorders. In general, both the HS and EMD demonstrated to be very useful to perform objective analysis of any kind of tremor and can therefore be potentially used to perform functional assessment.