Otolith shape analysis as a tool for population discrimination of the white grunt (Haemulon plumieri) stock in the northern coast of the Yucatan Peninsula, Mexico

In this study, the population structure of the white grunt (Haemulon plumieri) from the northern coast of the Yucatan Peninsula was determined through an otolith shape analysis based on the samples collected in three locations: Celestún (N 20°49",W 90°25"), Dzilam (N 21°23", W 88°54") and Cancún (N 21°21",W 86°52"). The otolith outline was based on the elliptic Fourier descriptors, which indicated that the H. plumieri population in the northern coast of the Yucatan Peninsula is composed of three geographically delimited units (Celestún, Dzilam, and Cancún). Significant differences were observed in mean otolith shapes among all samples (PERMANOVA; F2, 99 = 11.20, P = 0.0002), and the subsequent pairwise comparisons showed that all samples were significantly differently from each other. Samples do not belong to a unique white grunt population, and results suggest that they might represent a structured population along the northern coast of the Yucatan Peninsula

Revisiting the compositional data. Some fundamental questions and new prospects in Archaeometry and Archaeology

In this paper we examine the problem of compositional data from a different startingpoint. Chemical compositional data, as used in provenance studies on archaeologicalmaterials, will be approached from the measurement theory. The results will show, in avery intuitive way that chemical data can only be treated by using the approachdeveloped for compositional data. It will be shown that compositional data analysis is aparticular case in projective geometry, when the projective coordinates are in thepositive orthant, and they have the properties of logarithmic interval metrics. Moreover,it will be shown that this approach can be extended to a very large number ofapplications, including shape analysis. This will be exemplified with a case study inarchitecture of Early Christian churches dated back to the 5th-7th centuries AD

Analyzing shapes as compositions of distances

We propose to analyze shapes as “compositions” of distances in Aitchison geometry asan alternate and complementary tool to classical shape analysis, especially when sizeis non-informative.Shapes are typically described by the location of user-chosen landmarks. Howeverthe shape – considered as invariant under scaling, translation, mirroring and rotation– does not uniquely define the location of landmarks. A simple approach is to usedistances of landmarks instead of the locations of landmarks them self. Distances arepositive numbers defined up to joint scaling, a mathematical structure quite similar tocompositions. The shape fixes only ratios of distances. Perturbations correspond torelative changes of the size of subshapes and of aspect ratios. The power transformincreases the expression of the shape by increasing distance ratios. In analogy to thesubcompositional consistency, results should not depend too much on the choice ofdistances, because different subsets of the pairwise distances of landmarks uniquelydefine the shape.Various compositional analysis tools can be applied to sets of distances directly or afterminor modifications concerning the singularity of the covariance matrix and yield resultswith direct interpretations in terms of shape changes. The remaining problem isthat not all sets of distances correspond to a valid shape. Nevertheless interpolated orpredicted shapes can be backtransformated by multidimensional scaling (when all pairwisedistances are used) or free geodetic adjustment (when sufficiently many distancesare used)

EEG Windowed statitical wavelet deviation for estimation of muscular artifacts

Electroencephalographic (EEG) recordings are, most of the times, corrupted by spurious artifacts, which should be rejected or cleaned by the practitioner. As human scalp EEG screening is error-prone, automatic artifact detection is an issue of capital importance, to ensure objective and reliable results. In this paper we propose a new approach for discrimination of muscular activity in the human scalp quantitative EEG (QEEG), based on the time-frequency shape analysis. The impact of the muscular activity on the EEG can be evaluated from this methodology. We present an application of this scoring as a preprocessing step for EEG signal analysis, in order to evaluate the amount of muscular activity for two set of EEG recordings for dementia patients with early stage of Alzheimer’s disease and control age-matched subjects.

Statistical analysis and modelling of weather radar beam propagation conditions in the Po Valley (Italy)

Ground clutter caused by anomalous propagation (anaprop) can affect seriously radar rain rate estimates, particularly in fully automatic radar processing systems, and, if not filtered, can produce frequent false alarms. A statistical study of anomalous propagation detected from two operational C-band radars in the northern Italian region of Emilia Romagna is discussed, paying particular attention to its diurnal and seasonal variability. The analysis shows a high incidence of anaprop in summer, mainly in the morning and evening, due to the humid and hot summer climate of the Po Valley, particularly in the coastal zone. Thereafter, a comparison between different techniques and datasets to retrieve the vertical profile of the refractive index gradient in the boundary layer is also presented. In particular, their capability to detect anomalous propagation conditions is compared. Furthermore, beam path trajectories are simulated using a multilayer ray-tracing model and the influence of the propagation conditions on the beam trajectory and shape is examined. High resolution radiosounding data are identified as the best available dataset to reproduce accurately the local propagation conditions, while lower resolution standard TEMP data suffers from interpolation degradation and Numerical Weather Prediction model data (Lokal Model) are able to retrieve a tendency to superrefraction but not to detect ducting conditions. Observing the ray tracing of the centre, lower and upper limits of the radar antenna 3-dB half-power main beam lobe it is concluded that ducting layers produce a change in the measured volume and in the power distribution that can lead to an additional error in the reflectivity estimate and, subsequently, in the estimated rainfall rate.

Particle shape and orientation in laser diffraction and static image analysis size distribution analysis of micrometer sized rectangular particles

Laser diffraction (LD) and static image analysis (SIA) of rectangular particles [United States Pharmacopeia, USP30-NF25, General Chapter <776>, Optical Miroscopy.] have been systematically studied. To rule out sample dispersion and particle orientation as the root cause of differences in size distribution profiles, we immobilize powder samples on a glass plate by means of a dry disperser. For a defined region of the glass plate, we measure the diffraction pattern as induced by the dispersed particles, and the 2D dimensions of the individual particles using LD and optical microscopy, respectively. We demonstrate a correlation between LD and SIA, with the scattering intensity of the individual particles as the dominant factor. In theory, the scattering intensity is related to the square of the projected area of both spherical and rectangular particles. In traditional LD the size distribution profile is dominated by the maximum projected area of the particles (A). The diffraction diameters of a rectangular particle with length L and breadth B as measured by the LD instrument approximately correspond to spheres of diameter ØL and ØB respectively. Differences in the scattering intensity between spherical and rectangular particles suggest that the contribution made to the overall LD volume probability distribution by each rectangular particle is proportional to A2/L and A2/B. Accordingly, for rectangular particles the scattering intensity weighted diffraction diameter (SIWDD) explains an overestimation of their shortest dimension and an underestimation of their longest dimension. This study analyzes various samples of particles whose length ranges from approximately 10 to 1000 μm. The correlation we demonstrate between LD and SIA can be used to improve validation of LD methods based on SIA data for a variety of pharmaceutical powders all with a different rectangular particle size and shape.

Caracteritació de models funcionals d'àmfores romanes mitjançant Finite Element Analysis

Projecte de recerca elaborat a partir d’una estada al Laboratory of Archaeometry del National Centre of Scientific Research “Demokritos” d’Atenes, Grècia, entre juny i setembre 2006. Aquest estudi s’emmarca dins d’un context més ampli d’estudi del canvi tecnològic que es documenta en la producció d’àmfores de tipologia romana durant els segles I aC i I dC en els territoris costaners de Catalunya. Una part d’aquest estudi contempla el càlcul de les propietats mecàniques d’aquestes àmfores i la seva avaluació en funció de la tipologia amforal, a partir de l’Anàlisi d’Elements Finits (AEF). L’AEF és una aproximació numèrica que té el seu origen en les ciències d’enginyeria i que ha estat emprada per estimar el comportament mecànic d’un model en termes, per exemple, de deformació i estrès. Així, un objecte, o millor dit el seu model, es dividit en sub-dominis anomenats elements finits, als quals se’ls atribueixen les propietats mecàniques del material en estudi. Aquests elements finits estan connectats formant una xarxa amb constriccions que pot ser definida. En el cas d’aplicar una força determinada a un model, el comportament de l’objecte pot ser estimat mitjançant el conjunt d’equacions lineals que defineixen el rendiment dels elements finits, proporcionant una bona aproximació per a la descripció de la deformació estructural. Així, aquesta simulació per ordinador suposa una important eina per entendre la funcionalitat de ceràmiques arqueològiques. Aquest procediment representa un model quantitatiu per predir el trencament de l’objecte ceràmic quan aquest és sotmès a diferents condicions de pressió. Aquest model ha estat aplicat a diferents tipologies amforals. Els resultats preliminars mostren diferències significatives entre la tipologia pre-romana i les tipologies romanes, així com entre els mateixos dissenys amforals romans, d’importants implicacions arqueològiques.

Shape complexity based on mutual information

Shape complexity has recently received attention from different fields, such as computer vision and psychology. In this paper, integral geometry and information theory tools are applied to quantify the shape complexity from two different perspectives: from the inside of the object, we evaluate its degree of structure or correlation between its surfaces (inner complexity), and from the outside, we compute its degree of interaction with the circumscribing sphere (outer complexity). Our shape complexity measures are based on the following two facts: uniformly distributed global lines crossing an object define a continuous information channel and the continuous mutual information of this channel is independent of the object discretisation and invariant to translations, rotations, and changes of scale. The measures introduced in this paper can be potentially used as shape descriptors for object recognition, image retrieval, object localisation, tumour analysis, and protein docking, among others

Role of Particle Shape on the Stress Propagation in Granular Packings

We present an experimental and numerical study on the influence that particle aspect ratio has on the mechanical and structural properties of granular packings. For grains with maximal symmetry (squares), the stress propagation in the packing localizes forming chainlike forces analogous to the ones observed for spherical grains. This scenario can be understood in terms of stochastic models of aggregation and random multiplicative processes. As the grains elongate, the stress propagation is strongly affected. The interparticle normal force distribution tends toward a Gaussian, and, correspondingly, the force chains spread leading to a more uniform stress distribution reminiscent of the hydrostatic profiles known for standard liquids

Diversity of the bacterial community in the surface soil of a pear orchard based on 16S rRNA gene analysis

A cultivation-independent approach based on polymerase chain reaction (PCR)-amplified partial small subunit rRNA genes was used to characterize bacterial populations in the surface soil of a commercial pear orchard consisting of different pear cultivars during two consecutive growing seasons. Pyrus communis L. cvs Blanquilla, Conference, and Williams are among the most widely cultivated cultivars in Europe and account for the majority of pear production in Northeastern Spain. To assess the heterogeneity of the community structure in response to environmental variables and tree phenology, bacterial populations were examined using PCR-denaturing gradient gel electrophoresis (DGGE) followed by cluster analysis of the 16S ribosomal DNA profiles by means of the unweighted pair group method with arithmetic means. Similarity analysis of the band patterns failed to identify characteristic fingerprints associated with the pear cultivars. Both environmentally and biologically based principal-component analyses showed that the microbial communities changed significantly throughout the year depending on temperature and, to a lesser extent, on tree phenology and rainfall. Prominent DGGE bands were excised and sequenced to gain insight into the identities of the predominant bacterial populations. Most DGGE band sequences were related to bacterial phyla, such as Bacteroidetes, Cyanobacteria, Acidobacteria, Proteobacteria, Nitrospirae, and Gemmatimonadetes, previously associated with typical agronomic crop environments

Calculation of Franck-Condon factors including anharmonicity: simulation of the C2 H4+ X̃2B 3u←C2 H4X̃1 Ag band in the photoelectron spectrum of ethylene

Our new simple method for calculating accurate Franck-Condon factors including nondiagonal (i.e., mode-mode) anharmonic coupling is used to simulate the C2H4+X2B 3u←C2H4X̃1 Ag band in the photoelectron spectrum. An improved vibrational basis set truncation algorithm, which permits very efficient computations, is employed. Because the torsional mode is highly anharmonic it is separated from the other modes and treated exactly. All other modes are treated through the second-order perturbation theory. The perturbation-theory corrections are significant and lead to a good agreement with experiment, although the separability assumption for torsion causes the C2 D4 results to be not as good as those for C2 H4. A variational formulation to overcome this circumstance, and deal with large anharmonicities in general, is suggested

Bilinear models for spatio-temporal point distribution analysis: application to extrapolation of left ventricular, biventricular and whole heart cardiac dynamics

In this work we describe the usage of bilinear statistical models as a means of factoring the shape variability into two components attributed to inter-subject variation and to the intrinsic dynamics of the human heart. We show that it is feasible to reconstruct the shape of the heart at discrete points in the cardiac cycle. Provided we are given a small number of shape instances representing the same heart atdifferent points in the same cycle, we can use the bilinearmodel to establish this. Using a temporal and a spatial alignment step in the preprocessing of the shapes, around half of the reconstruction errors were on the order of the axial image resolution of 2 mm, and over 90% was within 3.5 mm. From this, weconclude that the dynamics were indeed separated from theinter-subject variability in our dataset.

Correspondence analysis of raw data

Correspondence analysis has found extensive use in ecology, archeology, linguisticsand the social sciences as a method for visualizing the patterns of association in a table offrequencies or nonnegative ratio-scale data. Inherent to the method is the expression of the datain each row or each column relative to their respective totals, and it is these sets of relativevalues (called profiles) that are visualized. This relativization of the data makes perfect sensewhen the margins of the table represent samples from sub-populations of inherently differentsizes. But in some ecological applications sampling is performed on equal areas or equalvolumes so that the absolute levels of the observed occurrences may be of relevance, in whichcase relativization may not be required. In this paper we define the correspondence analysis ofthe raw unrelativized data and discuss its properties, comparing this new method to regularcorrespondence analysis and to a related variant of non-symmetric correspondence analysis.

Modifications in the Si valence band after ion-beam-induced oxidation

The changes undergone by the Si surface after oxygen bombardment have special interest for acquiring a good understanding of the Si+-ion emission during secondary ion mass spectrometry (SIMS) analysis. For this reason a detailed investigation on the stoichiometry of the builtup surface oxides has been carried out using in situ x-ray photoemission spectroscopy (XPS). The XPS analysis of the Si 2p core level indicates a strong presence of suboxide chemical states when bombarding at angles of incidence larger than 30°. In this work a special emphasis on the analysis and interpretation of the valence band region was made. Since the surface stoichiometry or degree of oxidation varies with the angle of incidence, the respective valence band structures also differ. A comparison with experimentally measured and theoretically derived Si valence band and SiO2 valence band suggests that the new valence bands are formed by a combination of these two. This arises from the fact that Si¿Si bonds are present on the Si¿suboxide molecules, and therefore the corresponding 3p-3p Si-like subband, which extends towards the Si Fermi level, forms the top of the respective new valence bands. Small variations in intensity and energy position for this subband have drastic implications on the intensity of the Si+-ion emission during sputtering in SIMS measurements. A model combining chemically enhanced emission and resonant tunneling effects is suggested for the variations observed in ion emission during O+2 bombardment for Si targets.

Analysis of the near-intrinsic and extrinsic photocapacitance due to the EL2 level in boron implanted GaAs

A detailed analysis of the photocapacitance signal at the near‐band and extrinsic energetic ranges in Schottky barriers obtained on horizontal Bridgman GaAs wafers, which were implanted with boron at different doses and annealed at several temperatures, has been carried out by using the optical isothermal transient spectroscopy, OITS. The optical cross sections have been determined as well as the quenching efficiency of the EL2 level which has been found to be independent of the annealing temperature. Moreover, the quenching relaxation presents two significant features: (i) a strong increase of the quenching efficiency from 1.35 eV on and (ii) a diminution of the quenching transient amplitude in relation with that shown by the fundamental EL2 level. In order to explain this behavior, different cases are discussed assuming the presence of several energy levels, the existence of an optical recuperation, or the association of the EL2 trap with two levels located, respectively, at Ev+0.45 eV and Ec−0.75 eV. The theoretical simulation, taking into account these two last cases, is in agreement with the experimental photocapacitance data at low temperature, as well as at room temperature where the EL2 filling phototransient shows an anomalous behavior. Moreover, unlike the previous data reported for the EL2 electron optical cross section, the values found using our experimental technique are in agreement with the behavior deduced from the theoretical calculation. The utilization of the OITS method has also allowed the determination of another level, whose faster optical contribution is often added to that of the EL2 level when the DLOS or standard photocapacitance is used.