Application of compositional data analysis to geochemical data of marine sediments

In an earlier investigation (Burger et al., 2000) five sediment cores near the RodriguesTriple Junction in the Indian Ocean were studied applying classical statistical methods(fuzzy c-means clustering, linear mixing model, principal component analysis) for theextraction of endmembers and evaluating the spatial and temporal variation ofgeochemical signals. Three main factors of sedimentation were expected by the marinegeologists: a volcano-genetic, a hydro-hydrothermal and an ultra-basic factor. Thedisplay of fuzzy membership values and/or factor scores versus depth providedconsistent results for two factors only; the ultra-basic component could not beidentified. The reason for this may be that only traditional statistical methods wereapplied, i.e. the untransformed components were used and the cosine-theta coefficient assimilarity measure.During the last decade considerable progress in compositional data analysis was madeand many case studies were published using new tools for exploratory analysis of thesedata. Therefore it makes sense to check if the application of suitable data transformations,reduction of the D-part simplex to two or three factors and visualinterpretation of the factor scores would lead to a revision of earlier results and toanswers to open questions . In this paper we follow the lines of a paper of R. Tolosana-Delgado et al. (2005) starting with a problem-oriented interpretation of the biplotscattergram, extracting compositional factors, ilr-transformation of the components andvisualization of the factor scores in a spatial context: The compositional factors will beplotted versus depth (time) of the core samples in order to facilitate the identification ofthe expected sources of the sedimentary process.Kew words: compositional data analysis, biplot, deep sea sediments

Application of compositional data analysis to geochemical data of marine sediments

In an earlier investigation (Burger et al., 2000) five sediment cores near the Rodrigues Triple Junction in the Indian Ocean were studied applying classical statistical methods (fuzzy c-means clustering, linear mixing model, principal component analysis) for the extraction of endmembers and evaluating the spatial and temporal variation of geochemical signals. Three main factors of sedimentation were expected by the marine geologists: a volcano-genetic, a hydro-hydrothermal and an ultra-basic factor. The display of fuzzy membership values and/or factor scores versus depth provided consistent results for two factors only; the ultra-basic component could not be identified. The reason for this may be that only traditional statistical methods were applied, i.e. the untransformed components were used and the cosine-theta coefficient as similarity measure. During the last decade considerable progress in compositional data analysis was made and many case studies were published using new tools for exploratory analysis of these data. Therefore it makes sense to check if the application of suitable data transformations, reduction of the D-part simplex to two or three factors and visual interpretation of the factor scores would lead to a revision of earlier results and to answers to open questions . In this paper we follow the lines of a paper of R. Tolosana- Delgado et al. (2005) starting with a problem-oriented interpretation of the biplot scattergram, extracting compositional factors, ilr-transformation of the components and visualization of the factor scores in a spatial context: The compositional factors will be plotted versus depth (time) of the core samples in order to facilitate the identification of the expected sources of the sedimentary process. Kew words: compositional data analysis, biplot, deep sea sediments

A Log-linear model of grain size influence on the geochemistry of sediments

Sediment composition is mainly controlled by the nature of the source rock(s), and chemical (weathering) and physical processes (mechanical crushing, abrasion, hydrodynamic sorting) during alteration and transport. Although the factors controlling these processes are conceptually well understood, detailed quantification of compositional changes induced by a single process are rare, as are examples where the effects of several processes can be distinguished. The present study was designed to characterize the role of mechanical crushing and sorting in the absence of chemical weathering. Twenty sediment samples were taken from Alpine glaciers that erode almost pure granitoid lithologies. For each sample, 11 grain-size fractions from granules to clay (ø grades &-1 to &9) were separated, and each fraction was analysed for its chemical composition.The presence of clear steps in the box-plots of all parts (in adequate ilr and clr scales) against ø is assumed to be explained by typical crystal size ranges for the relevant mineral phases. These scatter plots and the biplot suggest a splitting of the full grain size range into three groups: coarser than ø=4 (comparatively rich in SiO2, Na2O, K2O, Al2O3, and dominated by “felsic” minerals like quartz and feldspar), finer than ø=8 (comparatively rich in TiO2, MnO, MgO, Fe2O3, mostly related to “mafic” sheet silicates like biotite and chlorite), and intermediate grains sizes (4≤ø &8; comparatively rich in P2O5 and CaO, related to apatite, some feldspar).To further test the absence of chemical weathering, the observed compositions were regressed against three explanatory variables: a trend on grain size in ø scale, a step function for ø≥4, and another for ø≥8. The original hypothesis was that the trend could be identified with weathering effects, whereas each step function would highlight those minerals with biggest characteristic size at its lower end. Results suggest that this assumption is reasonable for the step function, but that besides weathering some other factors (different mechanical behavior of minerals) have also an important contribution to the trend.Key words: sediment, geochemistry, grain size, regression, step function

A Log-linear model of grain size influence on the geochemistry of sediments

Sediment composition is mainly controlled by the nature of the source rock(s), and chemical (weathering) and physical processes (mechanical crushing, abrasion, hydrodynamic sorting) during alteration and transport. Although the factors controlling these processes are conceptually well understood, detailed quantification of compositional changes induced by a single process are rare, as are examples where the effects of several processes can be distinguished. The present study was designed to characterize the role of mechanical crushing and sorting in the absence of chemical weathering. Twenty sediment samples were taken from Alpine glaciers that erode almost pure granitoid lithologies. For each sample, 11 grain-size fractions from granules to clay (ø grades <-1 to >9) were separated, and each fraction was analysed for its chemical composition. The presence of clear steps in the box-plots of all parts (in adequate ilr and clr scales) against ø is assumed to be explained by typical crystal size ranges for the relevant mineral phases. These scatter plots and the biplot suggest a splitting of the full grain size range into three groups: coarser than ø=4 (comparatively rich in SiO2, Na2O, K2O, Al2O3, and dominated by “felsic” minerals like quartz and feldspar), finer than ø=8 (comparatively rich in TiO2, MnO, MgO, Fe2O3, mostly related to “mafic” sheet silicates like biotite and chlorite), and intermediate grains sizes (4≤ø <8; comparatively rich in P2O5 and CaO, related to apatite, some feldspar). To further test the absence of chemical weathering, the observed compositions were regressed against three explanatory variables: a trend on grain size in ø scale, a step function for ø≥4, and another for ø≥8. The original hypothesis was that the trend could be identified with weathering effects, whereas each step function would highlight those minerals with biggest characteristic size at its lower end. Results suggest that this assumption is reasonable for the step function, but that besides weathering some other factors (different mechanical behavior of minerals) have also an important contribution to the trend. Key words: sediment, geochemistry, grain size, regression, step function

Analisi de compostos orgànics persistents i estudi de la procedència d'hidrocarburs en sediments del mar Mediterrani català

Els sediments del fons marí són molt importants en els ecosistemes aquàtics. Aquests sediments són considerats absorbents naturals, ja que s’hi acumulen moltes substàncies nocives, com ara metalls pesants i contaminants orgànics persistents (COPs). En aquest estudi s’han analitzat mostres de sediments del mar Mediterrani amb la finalitat d’examinar els COPs i la possible procedència d’hidrocarburs derivats del petroli. Les mostres han estat extretes davant de la desembocadura del riu Besòs (se n’ha fet una extracció a 8 km i una altra a 54 km de distància respecte la costa) i davant del municipi català anomenat Vilanova i la Geltrú (se n’ha fet, també, una extracció a 8-10 km i una altra a 30-40 km de distància respecte la costa). D’aquestes mostres, obtingudes per perforació, se n’han analitzat el primer i segon centímetres. S’han quantificat els hidrocarburs aromàtics policíclics (PAHs), els bifenils policlorats (PCBs) i els compostos organoclorats (OCPs) utilitzats antigament com a pesticides. S’ha posat de manifest que en les zones més properes a la costa, les concertacions COPs són més elevades a la regió més meridional, tot i estar davant de la desembocadura del riu Besòs. Segurament això és degut a la influència de la deriva litoral (corrent costaner marí en direcció NE-SW). Finalment, l’índex Phytane/Pristane, l’índex de preferència de carboni o, en anglès Carbon Preference Index (CPI) i la barreja complexa sense resoldre o, en anglès, Unresolved Complex Mixture (UCM) apunten a la conclusió que l’origen dels hidrocarburs analitzats (n-alcans) podria ser petrogènic.

Analysis of Pleistocene paleodrainage evolution in the Po Basin (Italy) by multivariate statistical techniques

In order to obtain a high-resolution Pleistocene stratigraphy, eleven continuouslycored boreholes, 100 to 220m deep were drilled in the northern part of the PoPlain by Regione Lombardia in the last five years. Quantitative provenanceanalysis (QPA, Weltje and von Eynatten, 2004) of Pleistocene sands was carriedout by using multivariate statistical analysis (principal component analysis, PCA,and similarity analysis) on an integrated data set, including high-resolution bulkpetrography and heavy-mineral analyses on Pleistocene sands and of 250 majorand minor modern rivers draining the southern flank of the Alps from West toEast (Garzanti et al, 2004; 2006). Prior to the onset of major Alpine glaciations,metamorphic and quartzofeldspathic detritus from the Western and Central Alpswas carried from the axial belt to the Po basin longitudinally parallel to theSouthAlpine belt by a trunk river (Vezzoli and Garzanti, 2008). This scenariorapidly changed during the marine isotope stage 22 (0.87 Ma), with the onset ofthe first major Pleistocene glaciation in the Alps (Muttoni et al, 2003). PCA andsimilarity analysis from core samples show that the longitudinal trunk river at thistime was shifted southward by the rapid southward and westward progradation oftransverse alluvial river systems fed from the Central and Southern Alps.Sediments were transported southward by braided river systems as well as glacialsediments transported by Alpine valley glaciers invaded the alluvial plain.Kew words: Detrital modes; Modern sands; Provenance; Principal ComponentsAnalysis; Similarity, Canberra Distance; palaeodrainage

Oscil.lacions del nivell del mar i sediments associats durant l'Holocè a la Costa Brava Meridional (Girona)

Several marine deposits in the southern Costa Brava shoreline have been studied. They appear at different heights above and below the present sea level. Three groups are defined in relation to their origin: conglomerate levels at the bottom of the cliffs, not-cementedemerged beaches and cemented submerged beaches. The age of the emerged beaches has been accurately determined by means of the ceramic content and radiocarbon dating. Chronological succession of the deposits and their stratigraphic and paleontological characteristics allow to define a sequence of the sea level changes during the latest stages of the versiliane transgression in the studied area

Analysis of Pleistocene paleodrainage evolution in the Po Basin (Italy) by multivariate statistical techniques

In order to obtain a high-resolution Pleistocene stratigraphy, eleven continuously cored boreholes, 100 to 220m deep were drilled in the northern part of the Po Plain by Regione Lombardia in the last five years. Quantitative provenance analysis (QPA, Weltje and von Eynatten, 2004) of Pleistocene sands was carried out by using multivariate statistical analysis (principal component analysis, PCA, and similarity analysis) on an integrated data set, including high-resolution bulk petrography and heavy-mineral analyses on Pleistocene sands and of 250 major and minor modern rivers draining the southern flank of the Alps from West to East (Garzanti et al, 2004; 2006). Prior to the onset of major Alpine glaciations, metamorphic and quartzofeldspathic detritus from the Western and Central Alps was carried from the axial belt to the Po basin longitudinally parallel to the SouthAlpine belt by a trunk river (Vezzoli and Garzanti, 2008). This scenario rapidly changed during the marine isotope stage 22 (0.87 Ma), with the onset of the first major Pleistocene glaciation in the Alps (Muttoni et al, 2003). PCA and similarity analysis from core samples show that the longitudinal trunk river at this time was shifted southward by the rapid southward and westward progradation of transverse alluvial river systems fed from the Central and Southern Alps. Sediments were transported southward by braided river systems as well as glacial sediments transported by Alpine valley glaciers invaded the alluvial plain. Kew words: Detrital modes; Modern sands; Provenance; Principal Components Analysis; Similarity, Canberra Distance; palaeodrainage

Oscil.lacions del nivell del mar i sediments associats durant l'Holocè a la Costa Brava Meridional (Girona)

Several marine deposits in the southern Costa Brava shoreline have been studied. They appear at different heights above and below the present sea level. Three groups are defined in relation to their origin: conglomerate levels at the bottom of the cliffs, not-cemented emerged beaches and cemented submerged beaches. The age of the emerged beaches has been accurately determined by means of the ceramic content and radiocarbon dating. Chronological succession of the deposits and their stratigraphic and paleontological characteristics allow to define a sequence of the sea level changes during the latest stages of the versiliane transgression in the studied area

A Multi-state model for the analysis of functional decline and mortality of frail elderly patients

L’objectiu d’aquest estudi, que correspon a un projecte de recerca sobre la pèrdua funcional i la mortalitat de persones grans fràgils, és construir un procés de supervivència predictiu que tingui en compte l’evolució funcional i nutricional dels pacients al llarg del temps. En aquest estudi ens enfrontem a l’anàlisi de dades de supervivència i mesures repetides però els mètodes estadístics habituals per al tractament conjunt d’aquest tipus de dades no són apropiats en aquest cas. Com a alternativa utilitzem els models de supervivència multi-estats per avaluar l’associació entre mortalitat i recuperació, o no, dels nivells funcionals i nutricionals considerats normals. Després d’estimar el model i d’identificar els factors pronòstics de mortalitat és possible obtenir un procés predictiu que permet fer prediccions de la supervivència dels pacients en funció de la seva història concreta fins a un determinat moment. Això permet realitzar un pronòstic més precís de cada grup de pacients, la qual cosa pot ser molt útil per als professionals sanitaris a l’hora de prendre decisions clíniques.

Approximating subtree distances between Phylogenies

We give a 5-approximation algorithm to the rooted Subtree-Prune-and-Regraft (rSPR) distance between two phylogenies, which was recently shown to be NP-complete by Bordewich and Semple [5]. This paper presents the first approximation result for this important tree distance. The algorithm follows a standard format for tree distances such as Rodrigues et al. [24] and Hein et al. [13]. The novel ideas are in the analysis. In the analysis, the cost of the algorithm uses a \cascading" scheme that accounts for possible wrong moves. This accounting is missing from previous analysis of tree distance approximation algorithms. Further, we show how all algorithms of this type can be implemented in linear time and give experimental results.

Distribución de radionúclidos naturales (7Be y 210Pb) y antropogénicos (137Cs y 239,240Pu) en el Océano Ártico

La distribución de los radionúclidos (7Be, 137Cs, 210Pb y 239,240Pu) en los diferentes compartimentos (columna de agua, sedimentos del fondo y hielo) depende no sólo de las propiedades fisicoquímicas del propio radionúclido, sino también de los procesos que tienen lugar en el Océano Ártico. El hielo marino desempeña un papel relevante en la dispersión de los radionúclidos desde las áreas de formación del hielo a lo largo del Ártico debido a la liberación de éstos durante su tránsito o definitivamente en las áreas de ablación, especialmente en el Estrecho de Fram. Además, el hielo durante su deriva incorpora radionúclidos debido a la deposición atmosférica. Los radionúclidos conservativos (137Cs) permanecen disueltos en la columna de agua, mientras que los reactivos (210Pb e isótopos de Pu) se hallan predominantemente en los sedimentos de las plataformas continentales. En cambio, en las cuencas árticas centrales la mayor parte se encuentran en la columna de agua debido a la limitación en el arrastre por parte de las partículas (scavenging). En el Estrecho de Fram, como resultado de la liberación masiva de sedimentos se activa el proceso de scavenging y se observa como los inventarios de todos los radionúclidos estudiados en sedimentos aumentan en comparación con los valores de las cuencas centrales árticas. El inventario de radionúclidos en sedimentos del hielo (SIS) presenta valores inferiores a los otros compartimentos pero no es despreciable por la relevancia de la dispersión.

Reconstrucció de les temperatures superficials de l’aigua i de l’aire del Skagerrak (Escandinàvia) dels darrers 2000 anys

S’han utilitzat uns proxies nous per a estimar les temperatures superficials de l’aigua i de l’aire dels darrers 2000 anys de la zona de l’estret del Skagerrak per tal d’estudiar l’evolució del clima en aquesta zona i comparar-la amb la tendència del pal d’hoquei observada en diversos registres de l’hemisferi Nord. Els proxies utilitzats estan basats el els GDGTS i són el TEX86, que permet estimar la SST, el MAAT que permet estimar la temperatura mitjana anual de l’aire, el BIT que permet estimar la proporció de material terrestre en sediments marins o lacustres i el CBT que permet estimar el pH del sòl. Aquestes estimacions s’han comparat amb altres registres de la mateixa zona d’estudi. S’observa que el senyal de temperatura estimat pel TEX86 prové d’una zona en la columna d’aigua situada entre la zona fòtica i el fons marí. L’estimació del MAAT permet observar un òptim climàtic medieval, una petita edat del gel i un increment accentuat de les temperatures en els darrers 200 anys en la zona d’estudi. D’altra banda sembla que el MAAT estima la temperatura anual de l’aire a nivell regional, no local. L’estimació del BIT permet observar que la major part del testimoni sedimentari analitzat prové de material autòcton. S’observen dos pics al voltant de 445 DC i 1462 DC que podrien indicar episodis de pluges torrencials i/o canvis en la dinàmica de les corrents marines al Skagerrak. Finalment el pH estimat mitjançant el CBT sembla ser no representatiu dels sòls de Noruega, per tant aquesta forma d’estimar el pH podria ser no vàlida per als sediments marins.

Reconstrucció de les temperatures superficials marines dels últims dos milions d'anys al Corrent de Benguela utilitzant l'índex TEX86

Per tal de reconstruir i obtenir una evolució climàtica a partir de les temperatures superficials marines dels últims dos milions d’anys al Corrent de Benguela (costa oest sud-africana) s’han analitzat 60 mostres del testimoni amb ODP 175-1084. Per fer-ho, s’ha utilitzat la nova calibració de l’índex TEX86 (Kim et al., 2007) i s’han representat els resultats juntament amb els valors obtinguts amb altres índexs referents a les temperatures mitjanes anuals de l’aire (MAAT) i el grau d’aportació sedimentària d’origen continental als sediments marins (BIT). També s’han comparat amb els registres de temperatura d’altres estudis en la mateixa àrea obtinguts a partir de proxies diferents. Els resultats del TEX86 mostren certa concordança amb alguns dels valors obtinguts en altres estudis i proposen hipòtesis que relacionen de manera directa els tres índexs calculats.

Anàlisi de la composició en metalls pesants de sòls adjacents a carreteres per fluorescència de raigs X portàtil

L’objectiu del present estudi és la caracterització i avaluació del grau de contaminació per metalls pesants dels sòls adjacents a una carretera de l’àrea metropolitana de Barcelona. Aquesta caracterització inclou la comparació dels valors mesurats amb els nivells màxims legals i la determinació de l’origen antropogènic o litogènic dels metalls mitjançant mètodes estadístics. Per a assolir aquest objectiu, s’emprarà la tècnica de la fluorescència de raigs X portàtil (FPXRF), de manera que l’estudi també inclourà la validació dels resultats obtinguts mitjançant aquesta tècnica per comparació amb un mètode de referència, així com l’estudi de la principal limitació que presenta la tècnica, que és la presència d’humitat a les mostres.