Rainfall variability at decadal and longer time scales: signal or noise?

Rainfall variability occurs over a wide range of temporal scales. Knowledge and understanding of such variability can lead to improved risk management practices in agricultural and other industries. Analyses of temporal patterns in 100 yr of observed monthly global sea surface temperature and sea level pressure data show that the single most important cause of explainable, terrestrial rainfall variability resides within the El Nino-Southern Oscillation (ENSO) frequency domain (2.5-8.0 yr), followed by a slightly weaker but highly significant decadal signal (9-13 yr), with some evidence of lesser but significant rainfall variability at interclecadal time scales (15-18 yr). Most of the rainfall variability significantly linked to frequencies tower than ENSO occurs in the Australasian region, with smaller effects in North and South America, central and southern Africa, and western Europe. While low-frequency (LF) signals at a decadal frequency are dominant, the variability evident was ENSO-like in all the frequency domains considered. The extent to which such LF variability is (i) predictable and (ii) either part of the overall ENSO variability or caused by independent processes remains an as yet unanswered question. Further progress can only be made through mechanistic studies using a variety of models.

Variability in time delay between two models of pulse oximeters for deriving the photoplethysmographic signals

Pulse oximetry is commonly used as an arterial blood oxygen saturation (SaO(2)) measure. However, its other serial output, the photoplethysmography (PPG) signal, is not as well studied. Raw PPG signals can be used to estimate cardiovascular measures like pulse transit time (PTT) and possibly heart rate (HR). These timing-related measurements are heavily dependent on the minimal variability in phase delay of the PPG signals. Masimo SET (R) Rad-9 (TM) and Novametrix Oxypleth oximeters were investigated for their PPG phase characteristics on nine healthy adults. To facilitate comparison, PPG signals were acquired from fingers on the same hand in a random fashion. Results showed that mean PTT variations acquired from the Masimo oximeter (37.89 ms) were much greater than the Novametrix (5.66 ms). Documented evidence suggests that I ms variation in PTT is equivalent to I mmHg change in blood pressure. Moreover, the PTT trend derived from the Masimo oximeter can be mistaken as obstructive sleep apnoeas based on the known criteria. HR comparison was evaluated against estimates attained from an electrocardiogram (ECG). Novametrix differed from ECG by 0.71 +/- 0.58% (p < 0.05) while Masimo differed by 4.51 +/- 3.66% (p > 0.05). Modem oximeters can be attractive for their improved SaO(2) measurement. However, using raw PPG signals obtained directly from these oximeters for timing-related measurements warrants further investigations.

Detecting finite bandwidth periodic signals in stationary noise using the signal coherence spectrum

All signals that appear to be periodic have some sort of variability from period to period regardless of how stable they appear to be in a data plot. A true sinusoidal time series is a deterministic function of time that never changes and thus has zero bandwidth around the sinusoid's frequency. A zero bandwidth is impossible in nature since all signals have some intrinsic variability over time. Deterministic sinusoids are used to model cycles as a mathematical convenience. Hinich [IEEE J. Oceanic Eng. 25 (2) (2000) 256-261] introduced a parametric statistical model, called the randomly modulated periodicity (RMP) that allows one to capture the intrinsic variability of a cycle. As with a deterministic periodic signal the RMP can have a number of harmonics. The likelihood ratio test for this model when the amplitudes and phases are known is given in [M.J. Hinich, Signal Processing 83 (2003) 1349-13521. A method for detecting a RMP whose amplitudes and phases are unknown random process plus a stationary noise process is addressed in this paper. The only assumption on the additive noise is that it has finite dependence and finite moments. Using simulations based on a simple RMP model we show a case where the new method can detect the signal when the signal is not detectable in a standard waterfall spectrograrn display. (c) 2005 Elsevier B.V. All rights reserved.

Variability in heart rate estimates from different pulse oximeters

Since its introduction, pulse oximetry has become a conventional clinical measure. Besides being arterial blood oxygen saturation (SpO2) measure, pulse oximeters can be used for other cardiovascular measurements, like heart rate (HR) estimations, derived from its photo plethysmographic (PPG) signals. The temporal coherence of the PPG signals and thereby HR estimates are heavily dependent on its minimal phase variability. A Masimo SET Rad-9TM, Novametrix Oxypleth and a custom designed PPG system were investigated for their relative phase variation. R-R intervals from electro-cardiogram (ECG) were recorded concurrently as reference. PPG signals obtained from the 3 systems were evaluated by comparing their respective beat-to-beat (B-B) intervals with the corresponding R-R estimates during a static test. For their relative B-B comparison to the ECG, Novametrix system differed 0.680.52% (p

Heart rate variability characterization using a time-frequency based instantaneous frequency estimation technique

In this paper, a new method for characterizing the newborn heart rate variability (HRV) is proposed. The central of the method is the newly proposed technique for instantaneous frequency (IF) estimation specifically designed for nonstationary multicomponen signals such as HRV. The new method attempts to characterize the newborn HRV using features extracted from the time–frequency (TF) domain of the signal. These features comprise the IF, the instantaneous bandwidth (IB) and instantaneous energy (IE) of the different TF components of the HRV. Applied to the HRV of both normal and seizure suffering newborns, this method clearly reveals the locations of the spectral peaks and their time-varying nature. The total energy of HRV components, ET and ratio of energy concentrated in the low-frequency (LF) to that in high frequency (HF) components have been shown to be significant features in identifying the HRV of newborn with seizures.

Time-Frequency Distribution Moments of Heart Rate Variability for Neonatal Seizure Detection

In this paper, we propose features extracted from the heart rate variability (HRV) based on the first and second conditional moments of time-frequency distribution (TFD) as an additional guide for seizure detection in newborn. The features of HRV in the low frequency band (LF: 0-0.07 Hz), mid frequency band (MF: 0.07-0.15 Hz), and high frequency band (HF: 0.15-0.6 Hz) have been obtained by means of the time-frequency analysis using the modified-B distribution (MBD). Results of ongoing time-frequency research are presented. Based on our preliminary results, the first conditional moment of HRV which is also known as the mean/central frequency in the LF band and the second conditional moment of HRV which is also known as the variance/instantaneous bandwidth (IB) in the HF band can be used as a good feature to discriminate the newborn seizure from the non-seizure

A matching pursuit-based signal complexity measure for the analysis of newborn EEG

This paper presents a new relative measure of signal complexity, referred to here as relative structural complexity, which is based on the matching pursuit (MP) decomposition. By relative, we refer to the fact that this new measure is highly dependent on the decomposition dictionary used by MP. The structural part of the definition points to the fact that this new measure is related to the structure, or composition, of the signal under analysis. After a formal definition, the proposed relative structural complexity measure is used in the analysis of newborn EEG. To do this, firstly, a time-frequency (TF) decomposition dictionary is specifically designed to compactly represent the newborn EEG seizure state using MP. We then show, through the analysis of synthetic and real newborn EEG data, that the relative structural complexity measure can indicate changes in EEG structure as it transitions between the two EEG states; namely seizure and background (non-seizure).

Distribution, expression, and motif variability of ankyrin domain genes in Wolbachia pipientis

The endosymbiotic bacterium Wolbachia pipientis infects a wide range of arthropods, in which it induces a variety of reproductive phenotypes, including cytoplasmic incompatibility (CI), parthenogenesis, male killing, and reversal of genetic sex determination. The recent sequencing and annotation of the first Wolbachia genome revealed an unusually high number of genes encoding ankyrin domain (ANK) repeats. These ANK genes are likely to be important in mediating the Wolbachia-host interaction. In this work we determined the distribution and expression of the different ANK genes found in the sequenced Wolbachia wMel genome in nine Wolbachia strains that induce different phenotypic effects in their hosts. A comparison of the ANK genes of wMel and the non-CI-inducing wAu Wolbachia strain revealed significant differences between the strains. This was reflected in sequence variability in shared genes that could result in alterations in the encoded proteins, such as motif deletions, amino acid insertions, and in some cases disruptions due to insertion of transposable elements and premature stops. In addition, one wMel ANK gene, which is part of an operon, was absent in the wAu genome. These variations are likely to affect the affinity, function, and cellular location of the predicted proteins encoded by these genes.

Comparison of BR and QR Eigenvalue Algorithms for Power System Small Signal Stability Analysis

The BR algorithm is a novel and efficient method to find all eigenvalues of upper Hessenberg matrices and has never been applied to eigenanalysis for power system small signal stability. This paper analyzes differences between the BR and the QR algorithms with performance comparison in terms of CPU time based on stopping criteria and storage requirement. The BR algorithm utilizes accelerating strategies to improve its performance when computing eigenvalues of narrowly banded, nearly tridiagonal upper Hessenberg matrices. These strategies significantly reduce the computation time at a reasonable level of precision. Compared with the QR algorithm, the BR algorithm requires fewer iteration steps and less storage space without depriving of appropriate precision in solving eigenvalue problems of large-scale power systems. Numerical examples demonstrate the efficiency of the BR algorithm in pursuing eigenanalysis tasks of 39-, 68-, 115-, 300-, and 600-bus systems. Experiment results suggest that the BR algorithm is a more efficient algorithm for large-scale power system small signal stability eigenanalysis.

Genetics and race variability of the lucerne - Colletotrichum trifolii pathosystem in Australia

The Australian-bred lucerne cultivars, Trifecta and Sequel, were found to possess useful levels of resistance to both Colletotrichum trifolii races 1 and 2. Race 2 has only been previously observed in the United States and surveys did not reveal its presence in Australia. Multilocus fingerprinting using random amplified polymorphic DNA (RAPDs) analysis revealed low diversity (<10% dissimilarity) within Australian C. trifolii collections, and between the Australian race 1 isolates and a US race 2 isolate. Studies on the inheritance of resistance to C. trifolii race 1 in individual clones from Trifecta and Sequel revealed the presence of 2 different genetic mechanisms. One inheritance was for resistance as a recessive trait, and the other indicated that resistance was dominant. The recessive system has never been previously reported, whereas in the US, 2 completely dominant and independent tetrasomic genes Anl and Ant have been reported to condition C. trifolii resistance. It was not possible to fit the observed segregations from our studies to a single-gene model. In contrast to US studies, clones of cv. Sequel exhibiting the recessive resistance reacted differently to spray and stem injection with C. trifolii inoculum, being resistant to the former and susceptible to the latter, providing additional evidence for the presence of a different genetic mechanism conditioning resistance to those previously reported in the US. As C. trifolii is one of the most serious diseases of lucerne worldwide, the future development of molecular markers closely linked to the dominant and recessive resistances identified in these studies, and the relationships between these resistances and Anl and Ans as determined by genetic mapping, appear to be useful areas of future study.

Analysis of small signal stability margins using genetic optimization

Power system small signal stability analysis aims to explore different small signal stability conditions and controls, namely: (1) exploring the power system security domains and boundaries in the space of power system parameters of interest, including load flow feasibility, saddle node and Hopf bifurcation ones; (2) finding the maximum and minimum damping conditions; and (3) determining control actions to provide and increase small signal stability. These problems are presented in this paper as different modifications of a general optimization to a minimum/maximum, depending on the initial guesses of variables and numerical methods used. In the considered problems, all the extreme points are of interest. Additionally, there are difficulties with finding the derivatives of the objective functions with respect to parameters. Numerical computations of derivatives in traditional optimization procedures are time consuming. In this paper, we propose a new black-box genetic optimization technique for comprehensive small signal stability analysis, which can effectively cope with highly nonlinear objective functions with multiple minima and maxima, and derivatives that can not be expressed analytically. The optimization result can then be used to provide such important information such as system optimal control decision making, assessment of the maximum network's transmission capacity, etc. (C) 1998 Elsevier Science S.A. All rights reserved.

Autoinhibition by an internal nuclear localization signal revealed by the crystal structure of mammalian importin alpha

Importin alpha is the nuclear import receptor that recognizes classical monopartite and bipartite nuclear localization signals (NLSs). The structure of mouse importin alpha has been determined at 2.5 Angstrom resolution. The structure shows a large C-terminal domain containing armadillo repeats, and a less structured N-terminal importin beta-binding domain containing an internal NLS bound to the NLS-binding site. The structure explains the regulatory switch between the cytoplasmic, high-affinity form, and the nuclear, low-affinity form for NLS binding of the nuclear import receptor predicted by the current models of nuclear import. Importin beta conceivably converts the low- to high-affinity form by binding to a site overlapping the autoinhibitory sequence. The structure also has implications for understanding NLS recognition, and the structures of armadillo and HEAT repeats.

Relation of TNF-related apoptosis-inducing ligand (TRAIL) receptor and FLICE-inhibitory protein expression to TRAIL-induced apoptosis of melanoma

Past studies have shown that apoptosis mediated by TNF-related apoptosis-inducing ligand (TRAIL) is regulated by the expression of two death receptors [TRAIL receptor 1 (TRAIL-RI) and TRAIL-R2] and two decoy receptors (TRAIL-R3 and TRAIL-R4) that inhibit apoptosis, In previous studies, me have shown that TRAIL but not other members of the tumor necrosis factor family induce apoptosis in approximately two-thirds of melanoma cell lines. Here, we examined whether the expression of TRAIL-R at the mRNA and protein level in a panel of 28 melanoma cell lines and melanocytes correlated with their sensitivity to TRAIL-induced apoptosis, We report that at least three factors appear to underlie the variability in TRAIL-induced apoptosis. (a) Pour of nine cell lines that were insensitive to TRAIL-induced apoptosis failed to express death receptors, and in two instances, lines were devoid of all TRAIL-Rs. Southern analysis suggested this was due to loss of the genes for the death receptors, (b) Despite the presence of mRNA for the TRAIL-R, some of the lines failed to express TRAIL-R protein on their surface. This was evident for TRAIL-RI and more so for the TRAIL decoy receptors TRAIL-R3 and -R4, Studies on permeabilized cells revealed that the receptors were located within the cytoplasm and redistribution from the cytoplasm may represent a posttranslational control mechanism. (c) Surface expression of TRAIL-RI and -R2 (but not TRAIL-R3 and -R4) showed an overall correlation with TRAIL-induced apoptosis. However, certain melanoma cell lines and clones were relatively resistant to TRAIL-induced apoptosis despite the absence of decoy receptors and moderate levels of TRAIL-RI and -R2 expression. This may indicate the presence of inhibitors within the cells, but resistance to apoptosis could not be correlated with expression of the caspase inhibitor FLICE-inhibitory protein. mRNA for another TRAIL receptor, osteoprotegerin, was expressed in 22 of the melanoma lines but not on melanocytes. Its role in induction of apoptosis remains to be studied. These results appear to have important implications for future clinical studies on TRAIL.

Deconvolving the delta O-18 seawater component from subseasonal coral delta O-18 and Sr/Ca at Rarotonga in the southwestern subtropical Pacific for the period 1726 to 1997

To reconstruct oceanographic variations in the subtropical South Pacific, 271-year long subseasonal time series of Sr/Ca and delta(18)O were generated from a coral growing at Rarotonga (21.5degreesS, 159.5degreesW). In this case, coral Sr/Ca appears to be an excellent proxy for sea surface temperature (SST) and coral delta(18)O is a function of both SST and seawater delta(18)O composition (delta(18)O(sw)). Here, we focus on extracting the delta(18)O(sw) signal from these proxy records. A method is presented assuming that coral Sr/Ca is solely a function of SST and that coral delta(18)O is a function of both SST and delta(18)O(sw). This method separates the effects of delta(18)O(sw) from SST by breaking the instantaneous changes of coral delta(18)O into separate contributions by instantaneous SST and delta(18)O(sw) changes, respectively. The results show that on average delta(18)O(sw) at Rarotonga explains similar to39% of the variance in delta(18)O and that variations in SST explains the remaining similar to61% of delta(18)O variance. Reconstructed delta(18)O(sw) shows systematic increases in summer months (December-February) consistent with the regional pattern of variations in precipitation and evaporation. The delta(18)O(sw) also shows a positive linear correlation with satellite-derived estimated salinity for the period 1980 to 1997 (r = 0.72). This linear correlation between reconstructed delta(18)O(sw) and salinity makes it possible to use the reconstructed delta(18)O(sw) to estimate the past interannual and decadal salinity changes in this region. Comparisons of coral delta(18)O and delta(18)O(sw) at Rarotonga with the Pacific decadal oscillation index suggest that the decadal and interdecadal salinity and SST variability at Rarotonga appears to be related to basin-scale decadal variability in the Pacific. Copyright (C) 2002 Elsevier Science Ltd.