A new approach for objective identification of turns and steps in organism movement data relevant to random walk modelling

1. A first step in the analysis of complex movement data often involves discretisation of the path into a series of step-lengths and turns, for example in the analysis of specialised random walks, such as Lévy flights. However, the identification of turning points, and therefore step-lengths, in a tortuous path is dependent on ad-hoc parameter choices. Consequently, studies testing for movement patterns in these data, such as Lévy flights, have generated debate. However, studies focusing on one-dimensional (1D) data, as in the vertical displacements of marine pelagic predators, where turning points can be identified unambiguously have provided strong support for Lévy flight movement patterns. 2. Here, we investigate how step-length distributions in 3D movement patterns would be interpreted by tags recording in 1D (i.e. depth) and demonstrate the dimensional symmetry previously shown mathematically for Lévy-flight movements. We test the veracity of this symmetry by simulating several measurement errors common in empirical datasets and find Lévy patterns and exponents to be robust to low-quality movement data. 3. We then consider exponential and composite Brownian random walks and show that these also project into 1D with sufficient symmetry to be clearly identifiable as such. 4. By extending the symmetry paradigm, we propose a new methodology for step-length identification in 2D or 3D movement data. The methodology is successfully demonstrated in a re-analysis of wandering albatross Global Positioning System (GPS) location data previously analysed using a complex methodology to determine bird-landing locations as turning points in a Lévy walk. For this high-resolution GPS data, we show that there is strong evidence for albatross foraging patterns approximated by truncated Lévy flights spanning over 3·5 orders of magnitude. 5. Our simple methodology and freely available software can be used with any 2D or 3D movement data at any scale or resolution and are robust to common empirical measurement errors. The method should find wide applicability in the field of movement ecology spanning the study of motile cells to humans.

The future of the northeast Atlantic benthic flora in a high CO2 world

Seaweed and seagrass communities in the northeast Atlantic have been profoundly impacted by humans, and the rate of change is accelerating rapidly due to runaway CO2 emissions and mounting pressures on coastlines associated with human population growth and increased consumption of finite resources. Here, we predict how rapid warming and acidification are likely to affect benthic flora and coastal ecosystems of the northeast Atlantic in this century, based on global evidence from the literature as interpreted by the collective knowledge of the authorship. We predict that warming will kill off kelp forests in the south and that ocean acidification will remove maerl habitat in the north. Seagrasses will proliferate, and associated epiphytes switch from calcified algae to diatoms and filamentous species. Invasive species will thrive in niches liberated by loss of native species and spread via exponential development of artificial marine structures. Combined impacts of seawater warming, ocean acidification, and increased storminess may replace structurally diverse seaweed canopies, with associated calcified and noncalcified flora, with simple habitats dominated by noncalcified, turf-forming seaweeds.

Interaction Effects of Light, Temperature and Nutrient Limitations (N, P and Si) on Growth, Stoichiometry and Photosynthetic Parameters of the Cold-Water Diatom Chaetoceros wighamii

Light (20-450 μmol photons m^-2 s^-1), temperature (3-11°C) and inorganic nutrient composition (nutrient replete and N, P and Si limitation) were manipulated to study their combined influence on growth, stoichiometry (C:N:P:Chl a) and primary production of the cold water diatom Chaetoceros wighamii. During exponential growth, the maximum growth rate (~0.8 d^-1) was observed at high temperture and light; at 3°C the growth rate was ~30% lower under similar light conditions. The interaction effect of light and temperature were clearly visible from growth and cellular stoichiometry. The average C:N:P molar ratio was 80:13:1 during exponential growth, but the range, due to different light acclimation, was widest at the lowest temperature, reaching very low C:P (~50) and N:P ratios (~8) at low light and temperature. The C:Chl a ratio had also a wider range at the lowest temperature during exponential growth, ranging 16-48 (weight ratio) at 3°C compared with 17-33 at 11°C. During exponential growth, there was no clear trend in the Chl a normalized, initial slope (α*) of the photosynthesis-irradiance (PE) curve, but the maximum photosynthetic production (P_m) was highest for cultures acclimated to the highest light and temperature. During the stationary growth phase, the stoichiometric relationship depended on the limiting nutrient, but with generally increasing C:N:P ratio. The average photosynthetic quotient (PQ) during exponential growth was 1.26 but decreased to <1 under nutrient and light limitation, probably due to photorespiration. The results clearly demonstrate that there are interaction effects between light, temperature and nutrient limitation, and the data suggests greater variability of key parameters at low temperature. Understanding these dynamics will be important for improving models of aquatic primary production and biogeochemical cycles in a warming climate.

Phylogeographic analysis of the red seaweed Palmaria palmata reveals a Pleistocene marine glacial refugium in the English Channel

Phylogeography has provided a new approach to the analysis of the postglacial history of a wide range of taxa but, to date, little is known about the effect of glacial periods on the marine biota of Europe. We have utilized a combination of nuclear, plastid and mitochondrial genetic markers to study the biogeographic history of the red seaweed Palmaria palmata in the North Atlantic. Analysis of the nuclear rDNA operon (ITS1-5.8S-ITS2), the plastid 16S-trnI-trnA-23S-5S, rbcL-rbcS and rpl12-rps31-rpl9 regions and the mitochondrial cox2–3 spacer has revealed the existence of a previously unidentified marine refugium in the English Channel, along with possible secondary refugia off the southwest coast of Ireland and in northeast North America and/or Iceland. Coalescent and mismatch analyses date the expansion of European populations from approximately 128 000 bp and suggest a continued period of exponential growth since then. Consequently, we postulate that the penultimate (Saale) glacial maximum was the main event in shaping the biogeographic history of European P. palmata populations which persisted throughout the last (Weichselian) glacial maximum (c. 20 000 bp) in the Hurd Deep, an enigmatic trench in the English Channel.

Fault Diagnosis in Internal Combustion Engines Using Nonlinear Multivariate Statistics

This paper presents a statistical-based fault diagnosis scheme for application to internal combustion engines. The scheme relies on an identified model that describes the relationships between a set of recorded engine variables using principal component analysis (PCA). Since combustion cycles are complex in nature and produce nonlinear relationships between the recorded engine variables, the paper proposes the use of nonlinear PCA (NLPCA). The paper further justifies the use of NLPCA by comparing the model accuracy of the NLPCA model with that of a linear PCA model. A new nonlinear variable reconstruction algorithm and bivariate scatter plots are proposed for fault isolation, following the application of NLPCA. The proposed technique allows the diagnosis of different fault types under steady-state operating conditions. More precisely, nonlinear variable reconstruction can remove the fault signature from the recorded engine data, which allows the identification and isolation of the root cause of abnormal engine behaviour. The paper shows that this can lead to (i) an enhanced identification of potential root causes of abnormal events and (ii) the masking of faulty sensor readings. The effectiveness of the enhanced NLPCA based monitoring scheme is illustrated by its application to a sensor fault and a process fault. The sensor fault relates to a drift in the fuel flow reading, whilst the process fault relates to a partial blockage of the intercooler. These faults are introduced to a Volkswagen TDI 1.9 Litre diesel engine mounted on an experimental engine test bench facility.

MENTAL HEALTH AND CRIMINALITY AMONG PSYCHIATRIC OFFENDERS IN SOUTHEASTERN ONTARIO

OBJECTIVES: (1) Describe the population of mentally ill offenders over whom Ontario Review Board (ORB) held jurisdiction. (2) Assess the influences of psychopathology and criminal factors on criminal career. METHOD: This study was a retrospective case series design that reviewed all offenders who were court ordered for psychiatric evaluation at Mental Health Services Site of Providence Care in Kingston, Ontario from 1993 to 2007 (N=347). Eighty five subjects were found not criminally responsible on the account of mental disorder and were included in statistical analysis (n=85). Bivariate associations between five key variables and two outcome variables, seriousness of crime and recidivism, were examined. Logistic regressions were conducted to test the role of the predictor variables on the outcome variables. RESULTS: Age and change in principal psychiatric diagnosis over time were shown to be associated with seriousness of crime. Timing of psychiatric onset, early signs of deviance and change in diagnosis were shown to be associated with recidivism. On the whole, study population did not markedly vary in their distribution of variables by the outcome variables. Regression model included timing of psychiatric onset; psychiatric history; existence of criminal associate; child abuse history; and early signs of deviance. Recidivism was shown to be predicted by early signs of deviance (OR=8.154, p<0.05). Existence of criminal associates was shown to have substantial values of odds ratio at marginal significance (OR=7.577, p=0.13). CONCLUSION: Seriousness of crime is a complex factor that could not be sufficiently predicted by any one or combinations of study variables. Recidivism is better predicted by criminality factors than psychopathology. In the future, an exploratory analysis that more broadly examines the psychopathology and criminal factors in Canadian forensic population is needed. Findings from this study have important clinical and legal implications.

Dynamics of Viral RNA Synthesis during Measles Virus Infection

We propose a reference model of the kinetics of a viral RNA-dependent RNA polymerase (vRdRp) activities and its regulation during infection of eucaryotic cells. After measles virus infects a cell, mRNAs from all genes immediately start to accumulate linearly over the first 5 to 6 h and then exponentially until approximately 24 h. The change from a linear to an exponential accumulation correlates with de novo synthesis of vRdRp from the incoming template. Expression of the virus nucleoprotein (N) prior to infection shifts the balance in favor of replication. Conversely, inhibition of protein synthesis by cycloheximide favors the latter. The in vivo elongation speed of the viral polymerase is approximately 3 nucleotides/s. A similar profile with fivefold-slower kinetics can be obtained using a recombinant virus expressing a structurally altered polymerase. Finally, virions contain only encapsidated genomic, antigenomic, and 5'-end abortive replication fragment RNAs.

Calcium currents in interstitial cells from the guinea-pig bladder.

OBJECTIVE: To determine whether there are inward currents in interstitial cells (IC) isolated from the guinea-pig detrusor and if so, to characterise them using the patch-clamp technique and pharmacological agents. MATERIALS AND METHODS: Using the whole-cell patch-clamp technique, inward currents were studied in IC enzymatically isolated from the detrusor of the guinea-pig bladder. Currents were evoked by stepping positively from a holding potential of - 80 mV. RESULTS: Outward K+ currents were blocked by Cs+ internal solution to reveal inward currents, which activated at voltages more positive than - 50 mV, peaked at 0 mV, reversed near + 50 mV and were half-maximally activated at - 27 mV. The inward currents showed voltage-dependent inactivation and were half-maximally inactivated at - 36 mV. Fitting the activation and inactivation data with a Boltzmann function revealed a window current between - 40 mV and + 20 mV. The decay of the current evoked at 0 mV could be fitted with a single exponential with a mean time-constant of 88 ms. Replacing external Ca2+ with Ba2+ significantly increased this to 344 ms. The current amplitude was augmented by Ba2+, and by Bay K 8644. Inward currents were significantly reduced by 1 microm nifedipine, across the voltage range, but the blockade was more effective on the current evoked at 0 mV than that evoked by a step to - 20 mV, perhaps indicating voltage-dependence of the action of nifedipine or another component of inward current. Increasing the concentration of the drug to 10 microm caused no further significant reduction either at 0 mV or at -20 mV. However, in the presence of 1 microm nifedipine the latter current was significantly reduced by 100 microm Ni2+. Both currents were significantly reduced in Ca2+-free solution. CONCLUSIONS: IC from the guinea-pig detrusor possess inward currents with typical characteristics of L-type Ca2+ current. They also have a component of inward Ca2+ current, which was resistant to nifedipine, but sensitive to Ni2+. Further work is needed to characterise the latter conductance. PMID: 16686735 [PubMed - indexed for MEDLINE]

Goal Recognition through Goal Graph Analysis

We present a novel approach to goal recognition based on a two-stage paradigm of graph construction and analysis. First, a graph structure called a Goal Graph is constructed to represent the observed actions, the state of the world, and the achieved goals as well as various connections between these nodes at consecutive time steps. Then, the Goal Graph is analysed at each time step to recognise those partially or fully achieved goals that are consistent with the actions observed so far. The Goal Graph analysis also reveals valid plans for the recognised goals or part of these goals. Our approach to goal recognition does not need a plan library. It does not suffer from the problems in the acquisition and hand-coding of large plan libraries, neither does it have the problems in searching the plan space of exponential size. We describe two algorithms for Goal Graph construction and analysis in this paradigm. These algorithms are both provably sound, polynomial-time, and polynomial-space. The number of goals recognised by our algorithms is usually very small after a sequence of observed actions has been processed. Thus the sequence of observed actions is well explained by the recognised goals with little ambiguity. We have evaluated these algorithms in the UNIX domain, in which excellent performance has been achieved in terms of accuracy, efficiency, and scalability.

An exactly solvable four transition point model

Heavy particle collisions, in particular low-energy ion-atom collisions, are amenable to semiclassical JWKB phase integral analysis in the complex plane of the internuclear separation. Analytic continuation in this plane requires due attention to the Stokes phenomenon which parametrizes the physical mechanisms of curve crossing, non-crossing, the hybrid Nikitin model, rotational coupling and predissociation. Complex transition points represent adiabatic degeneracies. In the case of two or more such points, the Stokes constants may only be completely determined by resort to the so-called comparison- equation method involving, in particular, parabolic cylinder functions or Whittaker functions and their strong-coupling asymptotics. In particular, the Nikitin model is a two transition-point one-double-pole problem in each half-plane corresponding to either ingoing or outgoing waves. When the four transition points are closely clustered, new techniques are required to determine Stokes constants. However, such investigations remain incomplete, A model problem is therefore solved exactly for scattering along a one-dimensional z-axis. The energy eigenvalue is b(2)-a(2) and the potential comprises -z(2)/2 (parabolic) and -a(2) + b(2)/2z(2) (centrifugal/centripetal) components. The square of the wavenumber has in the complex z-plane, four zeros each a transition point at z = +/-a +/- ib and has a double pole at z = 0. In cases (a) and (b), a and b are real and unitarity obtains. In case (a) the reflection and transition coefficients are parametrized by exponentials when a(2) + b(2) > 1/2. In case (b) they are parametrized by trigonometrics when a(2) + b(2) <1/2 and total reflection is achievable. In case (c) a and b are complex and in general unitarity is not achieved due to loss of flux to a continuum (O'Rourke and Crothers, 1992 Proc. R. Sec. 438 1). Nevertheless, case (c) coefficients reduce to (a) or (b) under appropriate limiting conditions. Setting z = ht, with h a real constant, an attempt is made to model a two-state collision problem modelled by a pair of coupled first-order impact parameter equations and an appropriate (T) over tilde-tau relation, where (T) over tilde is the Stueckelberg variable and tau is the reduced or scaled time. The attempt fails because (T) over tilde is an odd function of tau, which is unphysical in a real collision problem. However, it is pointed out that by applying the Kummer exponential model to each half-plane (O'Rourke and Crothers 1994 J. Phys. B: At. Mel. Opt. Phys. 27 2497) the current model is in effect extended to a collision problem with four transition points and a double pole in each half-plane. Moreover, the attempt in itself is not a complete failure since it is shown that the result is a perfect diabatic inelastic collision for a traceless Hamiltonian matrix, or at least when both diagonal elements are odd and the off-diagonal elements equal and even.

Equilibrium polymerization of cyclic carbonate oligomers. II. Role of multiple active sites

Ring opening polymerization of bisphenol A polycarbonate is studied by Monte Carlo simulations of a model comprising a fixed number of Lennard-Jones particles and harmonic bonds [J. Chem. Phys. 115, 3895 (2001)]. Bond interchanges produced by a low concentration (0.10%less than or equal toc(a)less than or equal to0.36%) of chemically active particles lead to equilibrium polymerization. There is a continuous transition in both 2D and 3D from unpolymerized cyclic oligomers at low density to a system of linear chains at high density, and the polymeric phase is much more stable in three dimensions than in two. The steepness of the polymerization transition increases rapidly as c(a) decreases, suggesting that it is discontinuous in the limit c(a)-->0. The transition is entropy driven, since the average potential energy increases systematically upon polymerization, and there is a steady decline in the degree of polymerization as the temperature is lowered. The mass distribution functions for open chains and for rings are unimodal, with exponentially decaying tails that can be fitted by Zimm-Schulz functions and simpler exponential forms. (C) 2002 American Institute of Physics.

Light emission from scanning tunnelling microscope on polycrystalline Au films - what is happening at the single-grain level?

When operated with a metallic tip and sample the scanning tunnelling microscope constitutes a nanoscale, plasmonic light source yielding broadband emission up to a photon energy determined by the applied bias. The emission is due to tunnelling electron excitation and subsequent radiative decay of localized plasmon modes, which can be on the lateral scale of a single metal grain (similar to 25 nm) or less. For a Au-tip/Au-polycrystalline sample under ambient conditions it is found that the intensity and spectral content of the emitted light are not dependent on the lateral grain dimension, but are predominantly determined by the tip geometry. However, the intensity increases strongly with increasing film thickness (grain depth) up to 20-25 nm or approximately the skin depth of the Au film. Photon maps can show less emissive grains and two classes of this occurrence are distinguished. The first is geometrical in origin - a double-tip structure in this case - while the second is due to a contamination-induced lowering of the local work function that causes the tunnel gap to increase. It is suggested that differences in work-function lowering between grains presenting different crystalline facets, combined with an exponential decay in emitted light intensity with tip - sample distance, leads to grain contrast. These results are relevant to tip-enhanced Raman scattering and the fabrication of micro/nano-scale planar, light-emitting tunnel devices.

Effect of plasma scale length on multi-MeV proton production by intense laser pulses

The influence of the plasma density scale length on the production of MeV protons from thin foil targets irradiated at I lambda (2) = 5 x 10(19) Wcm(-2) has been studied. With an unperturbed foil, protons with energy >20 MeV were formed in an exponential energy spectrum with a temperature of 2.5 +/- 0.3 MeV. When a plasma with a scale length of 100 mum was preformed on the back of the foil, the maximum proton energy was reduced to

Novel Porphyrin-Incorporated Hydrogels for Photoactive Intraocular Lens Biomaterials

Novel surface-modified hydrogel materials have been prepared by binding charged porphyrins TMPyP (tetrakis-(4-N-methylpyridyl)porphyrin) and TPPS (tetrakis(4-sulfonatophenyl)porphyrin) to copolymers of HEMA (2-hydroxyethyl methacrylate) with either MAA (methacrylic acid) or DEAEMA (2-(diethylamino)ethylmethacrylate). The charged hydrogels display strong electrostatic interactions with the appropriate cationic or anionic porphyrins to give materials which are intended to be used to generate cytotoxic singlet oxygen (1O2) on photoexcitation and can therefore be used to reduce postoperative infection of the intraocular hydrogel-based replacement lenses that are used in cataract surgery. The UV/vis spectra of TMPyP in MAA:HEMA copolymers showed a small shift in the Soret band and a change from single exponential (161 ���­s) triplet decay lifetime in solution to a decay that could be fitted to a biexponential fit with two approximately equal components with ���´ ) 350 and 1300 ���­s. O2 bubbling reduced the decay to a dominant (90%) component with a much reduced lifetime of 3 ���­s and a minor, longer lived (20 ���­s) component. With D2O solvent the 1O2 lifetime was measured by 1270 nm fluorescence as 35 ���­s in MAA:HEMA, compared to 67 ���­s in solution, although absorbance-matched samples showed similar yield of 1O2 in the polymers and in aqueous solution. In contrast to the minor perturbation in photophysical properties caused by binding TMPyP to MAA:HEMA, TPPS binding to DEAEMA:HEMA copolymers profoundly changed the 1O2 generating ability of the TPPS. In N2-bubbled samples, the polymer-bound TPPS behaved in a similar manner to TMPyP in its copolymer host; however, O2 bubbling had only a very small effect on the triplet lifetime and no 1O2 generation could be detected. The difference in behavior may be linked to differences in binding in the two systems. With TMPyP in MAA:HEMA, confocal fluorescence microscopy showed significant penetration of the porphyrin into the core of the polymer film samples (>150 ���­m). However, for TPPS in DEAEMA:HEMA copolymers, although the porphyrin bound much more readily to the polymer, it remained localized in the first 20 ���­m, even in heavily loaded samples. It is possible that the resulting high concentration of TPPS may have cross-linked the hydrogels to such an extent that it significantly reduced the solubility and/or diffusion rate of oxygen into the doped polymers. This effect is significant since it demonstrates that even simple electrostatic binding of charged porphyrins to hydrogels can have an unexpectedly large effect on the properties of the system as a whole. In this case it makes the apparently promising TPPS/DEAEMA:HEMA system a poor candidate for clinical application as a postoperative antibacterial treatment for intraocular lenses while the apparently equivalent cationic system TMPyP/MAA:HEMA displays all the required properties.

Damage to plasmid DNA induced by low energy carbon ions

The damage induced in supercoiled plasmid DNA molecules by 1-6 keV carbon ions has been investigated as a function of ion exposure, energy and charge state. The production of short linear fragments through multiple double strand breaks has been demonstrated and exponential exposure responses for each of the topoisomers have been found. The cross section for the loss of supercoiling was calculated to be (2.2 +/- 0.5) x 10(-14) cm(2) for 2 keVC(+) ions. For singly charged carbon ions, increased damage was observed with increasing ion energy. In the case of 2 keV doubly charged ions, the damage was greater than for singly charged ions of the same energy. These observations demonstrate that ion induced damage is a function of both the kinetic and potential energies of the ion.