Time series analysis of lava flux

We have applied time series analytical techniques to the flux of lava from an extrusive eruption. Tilt data acting as a proxy for flux are used in a case study of the May–August 1997 period of the eruption at Soufrière Hills Volcano, Montserrat. We justify the use of such a proxy by simple calibratory arguments. Three techniques of time series analysis are employed: spectral, spectrogram and wavelet methods. In addition to the well-known ~9-hour periodicity shown by these data, a previously unknown periodic flux variability is revealed by the wavelet analysis as a 3-day cycle of frequency modulation during June–July 1997, though the physical mechanism responsible is not clear. Such time series analysis has potential for other lava flux proxies at other types of volcanoes.

A shift to randomness of brain oscillations in people with autism

BACKGROUND: Resting-state functional magnetic resonance imaging (fMRI) enables investigation of the intrinsic functional organization of the brain. Fractal parameters such as the Hurst exponent, H, describe the complexity of endogenous low-frequency fMRI time series on a continuum from random (H = .5) to ordered (H = 1). Shifts in fractal scaling of physiological time series have been associated with neurological and cardiac conditions. METHODS: Resting-state fMRI time series were recorded in 30 male adults with an autism spectrum condition (ASC) and 33 age- and IQ-matched male volunteers. The Hurst exponent was estimated in the wavelet domain and between-group differences were investigated at global and voxel level and in regions known to be involved in autism. RESULTS: Complex fractal scaling of fMRI time series was found in both groups but globally there was a significant shift to randomness in the ASC (mean H = .758, SD = .045) compared with neurotypical volunteers (mean H = .788, SD = .047). Between-group differences in H, which was always reduced in the ASC group, were seen in most regions previously reported to be involved in autism, including cortical midline structures, medial temporal structures, lateral temporal and parietal structures, insula, amygdala, basal ganglia, thalamus, and inferior frontal gyrus. Severity of autistic symptoms was negatively correlated with H in retrosplenial and right anterior insular cortex. CONCLUSIONS: Autism is associated with a small but significant shift to randomness of endogenous brain oscillations. Complexity measures may provide physiological indicators for autism as they have done for other medical conditions.

The role of biotic and abiotic factors in evolution of ant dispersal in the milkwort family (Polygalaceae)

A phylogenetic approach was taken to investigate the evolutionary history of seed appendages in the plant family Polygalaceae (Fabales) and determine which factors might be associated with evolution of elaiosomes through comparisons to abiotic (climate) and biotic (ant species number and abundance) timelines. Molecular datasets from three plastid regions representing 160 species were used to reconstruct a phylogenetic tree of the order Fabales, focusing on Polygalaceae. Bayesian dating methods were used to estimate the age of the appearance of ant-dispersed elaiosomes in Polygalaceae, shown by likelihood optimizations to have a single origin in the family. Topology-based tests indicated a diversification rate shift associated with appearance of caruncular elaiosomes. We show that evolution of the caruncular elaiosome type currently associated with ant dispersal occurred 54.0-50.5 million year ago. This is long after an estimated increase in ant lineages in the Late Cretaceous based on molecular studies, but broadly concomitant with increasing global temperatures culminating in the Late Paleocene-Early Eocene thermal maxima. These results suggest that although most major ant clades were present when elaiosomes appeared, the environmental significance of elaiosomes may have been an important factor in success of elaiosome-bearing lineages. Ecological abundance of ants is perhaps more important than lineage numbers in determining significance of ant dispersal. Thus, our observation that elaiosomes predate increased ecological abundance of ants inferred from amber deposits could be indicative of an initial abiotic environmental function.

Structural and equilibrium studies on mixed ligand complexes of Co(II), Ni(II), Cu(II) and Zn(II) with N-(2-benzimidazolyl)methyliminodiacetic acid and typical N,N donor ligands

Mixed ligand complexes: [Co(L)(bipy)] (.) 3H(2)O (1), [Ni(L)(phen)] (.) H2O (2), [Cu(L)(phen)] (.) 3H(2)O (3) and [Zn(L)(bipy)] (.) 3H(2)O (4), where L2- = two -COOH deprotonated dianion of N-(2-benzimidazolyl)methyliminodiacetic acid (H(2)bzimida, hereafter, H,L), bipy = 2,2' bipyridine and phen = 1,10-phenanthroline have been isolated and characterized by elemental analysis, spectral and magnetic measurements and thermal studies. Single crystal X-ray diffraction studies show octahedral geometry for 1, 2 and 4 and square pyramidal geometry for 3. Equilibrium studies in aqueous solution (ionic strength I = 10(-1) mol dm(-3) (NaNO3), at 25 +/- 1 degrees C) using different molar proportions of M(II):H2L:B, where M = Co, Ni, Cu and Zn and B = phen, bipy and en (ethylene diamine), however, provides evidence of formation of mononuclear and binuclear binary and mixed ligand complexes: M(L), M(H-1L)(-), M(B)(2+), M(L)(B), M(H-1L)(B)(-), M-2(H-1L)(OH), (B)M(H-1L)M(B)(+), where H-1L3- represents two -COOH and the benzimidazole NI-H deprotonated quadridentate (O-, N, O-, N), or, quinquedentate (O-, N, O-, N, N-) function of the coordinated ligand H,L. Binuclear mixed ligand complex formation equilibria: M(L)(B) + M(B)(2+) = (B)M(H-1L)M(B)(+) + H+ is favoured with higher pi-acidity of the B ligands. For Co(II), Ni(II) and Cu(II), these equilibria are accompanied by blue shift of the electronic absorption maxima of M(II) ions, as a negatively charged bridging benzimidazolate moiety provides stronger ligand field than a neutral one. Solution stability of the mixed ligand complexes are in the expected order: Co(II) < Ni(II) < Cu(II) > Zn(II). The Delta logK(M) values are less negetive than their statistical values, indicating favoured formation of the mixed ligand complexes over the binary ones. (c) 2005 Elsevier B.V. All rights reserved.

Mapping fibre orientation in complex-shaped biological systems with micrometre resolution by scanning X-ray microdiffraction

A fully automated procedure to extract and to image local fibre orientation in biological tissues from scanning X-ray diffraction is presented. The preferred chitin fibre orientation in the flow sensing system of crickets is determined with high spatial resolution by applying synchrotron radiation based X-ray microbeam diffraction in conjunction with advanced sample sectioning using a UV micro-laser. The data analysis is based on an automated detection of azimuthal diffraction maxima after 2D convolution filtering (smoothing) of the 2D diffraction patterns. Under the assumption of crystallographic fibre symmetry around the morphological fibre axis, the evaluation method allows mapping the three-dimensional orientation of the fibre axes in space. The resulting two-dimensional maps of the local fibre orientations - together with the complex shape of the flow sensing system - may be useful for a better understanding of the mechanical optimization of such tissues.

Isolation, fractionation and characterisation of proteins from Mucuna bean

The chemical composition and fractional distribution of protein isolates prepared from species of Mucuna bean were studied. Using six different extraction media, the yield of protein based on the Kjeldahl procedure varied from 8% to 34%, and the protein content varied from 75% to 95%. When the yields were high, the colour of the isolates generally tended to be dark and unsatisfactory. Hence, the use of chemical treatments and high pressure processing were explored. The solubility maxima for the protein isolates in water were found to occur at pH values of 2.0 and 11.0, while the pH corresponding to minimum solubility (i.e. isoelectric region) occurred at pH values of 4.0 and 5.0. The total essential amino acid in the isolates ranged from 495 to 557 mg g(-1) protein, which compares favourably with the recommended level for pre-school and school children. Methionine and cysteine were the limiting amino acids. A key nutritional attribute of the protein isolates was its high lysine content. The isolate can therefore complement cereal-based foods which are deficient in lysine. The proteins mainly consisted of albumins, glutelins and globulins. Prolamins were only present in trace concentration (< 0.3%). Gel filtration chromatograms of the isolates indicated the presence of major protein fractions with molecular weights of 40 and 15 kDa, while gel electrophoresis (SDS-PAGE) indicated a major broad zone with molecular weights of 36 +/- 7 and 17.3 +/- 13 kDa. (c) 2006 Elsevier Ltd. All rights reserved.

Unsupervised segmentation using Gabor wavelets and statistical features in LIDAR data analysis

In this paper, we address issues in segmentation Of remotely sensed LIDAR (LIght Detection And Ranging) data. The LIDAR data, which were captured by airborne laser scanner, contain 2.5 dimensional (2.5D) terrain surface height information, e.g. houses, vegetation, flat field, river, basin, etc. Our aim in this paper is to segment ground (flat field)from non-ground (houses and high vegetation) in hilly urban areas. By projecting the 2.5D data onto a surface, we obtain a texture map as a grey-level image. Based on the image, Gabor wavelet filters are applied to generate Gabor wavelet features. These features are then grouped into various windows. Among these windows, a combination of their first and second order of statistics is used as a measure to determine the surface properties. The test results have shown that ground areas can successfully be segmented from LIDAR data. Most buildings and high vegetation can be detected. In addition, Gabor wavelet transform can partially remove hill or slope effects in the original data by tuning Gabor parameters.

Face verification using Gabor wavelets and AdaBoost

This paper presents a new face verification algorithm based on Gabor wavelets and AdaBoost. In the algorithm, faces are represented by Gabor wavelet features generated by Gabor wavelet transform. Gabor wavelets with 5 scales and 8 orientations are chosen to form a family of Gabor wavelets. By convolving face images with these 40 Gabor wavelets, the original images are transformed into magnitude response images of Gabor wavelet features. The AdaBoost algorithm selects a small set of significant features from the pool of the Gabor wavelet features. Each feature is the basis for a weak classifier which is trained with face images taken from the XM2VTS database. The feature with the lowest classification error is selected in each iteration of the AdaBoost operation. We also address issues regarding computational costs in feature selection with AdaBoost. A support vector machine (SVM) is trained with examples of 20 features, and the results have shown a low false positive rate and a low classification error rate in face verification.

Facial feature extraction and selection by Gabor wavelets and boosting

In this paper, we present a feature selection approach based on Gabor wavelet feature and boosting for face verification. By convolution with a group of Gabor wavelets, the original images are transformed into vectors of Gabor wavelet features. Then for individual person, a small set of significant features are selected by the boosting algorithm from a large set of Gabor wavelet features. The experiment results have shown that the approach successfully selects meaningful and explainable features for face verification. The experiments also suggest that for the common characteristics such as eyes, noses, mouths may not be as important as some unique characteristic when training set is small. When training set is large, the unique characteristics and the common characteristics are both important.

The application of the Hilbert spectrum to the analysis of electromyographic signals

This paper investigates the application of the Hilbert spectrum (HS), which is a recent tool for the analysis of nonlinear and nonstationary time-series, to the study of electromyographic (EMG) signals. The HS allows for the visualization of the energy of signals through a joint time-frequency representation. In this work we illustrate the use of the HS in two distinct applications. The first is for feature extraction from EMG signals. Our results showed that the instantaneous mean frequency (IMNF) estimated from the HS is a relevant feature to clinical practice. We found that the median of the IMNF reduces when the force level of the muscle contraction increases. In the second application we investigated the use of the HS for detection of motor unit action potentials (MUAPs). The detection of MUAPs is a basic step in EMG decomposition tools, which provide relevant information about the neuromuscular system through the morphology and firing time of MUAPs. We compared, visually, how MUAP activity is perceived on the HS with visualizations provided by some traditional (e.g. scalogram, spectrogram, Wigner-Ville) time-frequency distributions. Furthermore, an alternative visualization to the HS, for detection of MUAPs, is proposed and compared to a similar approach based on the continuous wavelet transform (CWT). Our results showed that both the proposed technique and the CWT allowed for a clear visualization of MUAP activity on the time-frequency distributions, whereas results obtained with the HS were the most difficult to interpret as they were extremely affected by spurious energy activity. (c) 2008 Elsevier Inc. All rights reserved.

EMG signal filtering based on Empirical Mode Decomposition

This paper introduces a procedure for filtering electromyographic (EMG) signals. Its key element is the Empirical Mode Decomposition, a novel digital signal processing technique that can decompose my time-series into a set of functions designated as intrinsic mode functions. The procedure for EMG signal filtering is compared to a related approach based on the wavelet transform. Results obtained from the analysis of synthetic and experimental EMG signals show that Our method can be Successfully and easily applied in practice to attenuation of background activity in EMG signals. (c) 2006 Elsevier Ltd. All rights reserved.

Retrieval of global atmospheric electrical activity at a polluted urban site

The atmospheric electrical Potential Gradient (PG) arises from global thunderstorm activity, but surface measurements of the atmospheric Potential Gradient (PG) are influenced by global thunderstorms and local aerosol concentration changes. The local aerosol change can be monitored independently, and in some cases the concentration changes are closely related to PG changes. For these circumstances, a general theory to remove the local aerosol influence on PG measurements has been developed. Continuous measurements of PG and aerosol mass concentration were made during 24–31 Dec, 2005 within an urban environment at Reading, UK. The average diurnal variation of PG showed a double diurnal cycle, with maxima in the early morning and evening hours. The aerosol concentration has similar double maxima. Removing the aerosol using from the PG and aerosol correlation returns a single diurnal cycle, suggestive of the more global PG diurnal cycle.

A novel heteroditopic terpyridine-pincer ligand as building block for mono- and heterometallic Pd(II) and Ru(II) complexes

A palladium-catalyzed Stille coupling reaction was employed as a versatile method for the synthesis of a novel terpyridine-pincer (3, TPBr) bridging ligand, 4'-{4-BrC6H2(CH2NMe2)(2)-3,5}-2,2':6',2 ''-terpyridine. Mononuclear species [PdX(TP)] (X = Br, Cl), [Ru(TPBr)(tpy)](PF6)(2), and [Ru(TPBr)(2)](PF6)(2), synthesized by selective metalation of the NCNBr-pincer moiety or complexation of the terpyridine of the bifunctional ligand TPBr, were used as building blocks for the preparation of heterodi- and trimetallic complexes [Ru(TPPdCl)(tpy)](PF6)(2) (7) and [Ru(TPPdCl)(2)]-(PF6)(2) (8). The molecular structures in the solid state of [PdBr(TP)] (4a) and [Ru(TPBr)(2)](PF6)(2) (6) have been determined by single-crystal X-ray analysis. Electrochemical behavior and photophysical properties of the mono-and heterometallic complexes are described. All the above di- and trimetallic Ru complexes exhibit absorption bands attributable to (MLCT)-M-1 (Ru -> tpy) transitions. For the heteroleptic complexes, the transitions involving the unsubstituted tpy ligand are at a lower energy than the tpy moiety of the TPBr ligand. The absorption bands observed in the electronic spectra for TPBr and [PdCl(TP)] have been assigned with the aid of TD-DFT calculations. All complexes display weak emission both at room temperature and in a butyronitrile glass at 77 K. The considerable red shift of the emission maxima relative to the signal of the reference compound [Ru(tpy)(2)](2+) indicates stabilization of the luminescent (MLCT)-M-3 state. For the mono- and heterometallic complexes, electrochemical and spectroscopic studies (electronic absorption and emission spectra and luminescence lifetimes recorded at room temperature and 77 K in nitrile solvents), together with the information gained from IR spectroelectrochemical studies of the dimetallic complex [Ru(TPPdSCN)(tpy)](PF6)(2), are indicative of charge redistribution through the bridging ligand TPBr. The results are in line with a weak coupling between the {Ru(tpy)(2)} chromophoric unit and the (non)metalated NCN-pincer moiety.

Evaluation of peak-picking algorithms for protein mass spectrometry

Peak picking is an early key step in MS data analysis. We compare three commonly used approaches to peak picking and discuss their merits by means of statistical analysis. Methods investigated encompass signal-to-noise ratio, continuous wavelet transform, and a correlation-based approach using a Gaussian template. Functionality of the three methods is illustrated and discussed in a practical context using a mass spectral data set created with MALDI-TOF technology. Sensitivity and specificity are investigated using a manually defined reference set of peaks. As an additional criterion, the robustness of the three methods is assessed by a perturbation analysis and illustrated using ROC curves.

Aspects of the large scale tropical atmospheric circulation

An account is given of a number of recent studies with idealised models whose aim is to further understanding of the large-scale tropical atmospheric circulation. Initial-value integrations with a model with imposed heating are used to discuss aspects of the Asian summer monsoon, including constraints on cross-equatorial flow into the monsoon. The summer descent in the Mediterranean region and on the eastern sides of the summer subtropical anticyclones are seen to be associated with the monsoons to their east. An aqua-planet GCM is used to investigate the relationship between simple SST distributions and tropical convection and circulation. The existence of strong equatorial convection and Hadley cells is found to depend sensitively on the curvature of the meridional profile in SST. Zonally confined SST maxima produce convective maxima centred to the west and suppression of convection elsewhere. Strong equatorial zonal flow changes are found in some experiments and three mechanisms for producing these are investigated in a model with imposed heating. 1.