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

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

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.

Arquitectures d'Aprenentatge Automàtic per al Reconeixement d'Estructures Lingüístiques en el Text

Estudi realitzat a partir d’una estada al Computer Science and Artificial Intelligence Lab, del Massachusetts Institute of Technology, entre 2006 i 2008. La recerca desenvolupada en aquest projecte se centra en mètodes d'aprenentatge automàtic per l'anàlisi sintàctica del llenguatge. Com a punt de partida, establim que la complexitat del llenguatge exigeix no només entendre els processos computacionals associats al llenguatge sinó també entendre com es pot aprendre automàticament el coneixement per a dur a terme aquests processos.

Interlinkages and forecasting made possible by the use of the DECOIN analytical tool kit

This technical report is a document prepared as a deliverable [D4.3 Report of the Interlinkages and forecasting prototype tool] of a EU project – DECOIN Project No. 044428 - FP6-2005-SSP-5A. The text is divided into 4 sections: (1) this short introductory section explains the purpose of the report; (2) the second section provides a general discussion of a systemic problem found in existing quantitative analysis of sustainability. It addresses the epistemological implications of complexity, which entails the need of dealing with the existence of Multiple-Scales and non-equivalent narratives (multiple dimensions/attributes) to be used to define sustainability issues. There is an unavoidable tension between a “steady-state view” (= the perception of what is going on now – reflecting a PAST --& PRESENT view of the reality) versus an “evolutionary view” (= the unknown transformation that we have to expect in the process of becoming of the observed reality and in the observer – reflecting a PRESENT --& FUTURE view of the reality). The section ends by listing the implications of these points on the choice of integrated packages of sustainability indicators; (3) the third section illustrates the potentiality of the DECOIN toolkit for the study of sustainability trade-offs and linkages across indicators using quantitative examples taken from cases study of another EU project (SMILE). In particular, this section starts by addressing the existence of internal constraints to sustainability (economic versus social aspects). The narrative chosen for this discussion focuses on the dark side of ageing and immigration on the economic viability of social systems. Then the section continues by exploring external constraints to sustainability (economic development vs the environment). The narrative chosen for this discussion focuses on the dark side of current strategy of economic development based on externalization and the “bubbles-disease”; (4) the last section presents a critical appraisal of the quality of energy data found in energy statistics. It starts with a discussion of the general goal of statistical accounting. Then it introduces the concept of multipurpose grammars. The second part uses the experience made in the activities of the DECOIN project to answer the question: how useful are EUROSTAT energy statistics? The answer starts with an analysis of basic epistemological problems associated with accounting of energy. This discussion leads to the acknowledgment of an important epistemological problem: the unavoidable bifurcations in the mechanism of accounting needed to generate energy statistics. By using numerical example the text deals with the following issues: (i) the pitfalls of the actual system of accounting in energy statistics; (ii) a critical appraisal of the actual system of accounting in BP statistics; (iii) a critical appraisal of the actual system of accounting in Eurostat statistics. The section ends by proposing an innovative method to represent energy statistics which can result more useful for those willing develop sustainability indicators.

Applications of the MuSIASEM approach to study changes in the metabolic pattern of Catalonia

This document presents an integrated analysis of the performance of Catalonia based on an analysis of how the energy consumption (measured at the societal level for the Catalan Society) is used within both the productive sectors of the economy and the household, to generate added value, jobs, and to guarantee a given level of material standard of living to the population. The trends found in Catalonia are compared to the trends of other European Countries to contextualize the performance of Catalonia with respect to other societies that have followed different paths of economic development. The first part of the document consists of the Multi-Scale Integrated Analysis of Societal and Ecosystem Metabolism (MuSIASEM) approach that has been used to provide this integrated analysis of Catalan Society across different scales (starting from an analysis of the specific sectors of the Catalan economy as an Autonomous Community and scaling up to an intra-regional (European Union 14) comparison) and across different dimensions of analyses of energy consumption coupled with added value generation. Within the scope of this study, we observe the various trajectories of changes in the metabolic pattern for Catalonia and the EU14 countries in the Paid Work Sectors composed of namely, the Agricultural Sector, the Productive Sector and the Services and Government Sector also in comparison with the changes in the household sector. The flow intensities of the exosomatic energy and the added value generated for each specific sector are defined per hour of human activity, thus characterized as exosomatic energy (MJ/hour) (or Exosomatic Metabolic Rate) and added value (€/hour) (Economic Labour Productivity) across multiple levels. Within the second part of the document, the possible usage of the MuSIASEM approach to land use analyses (using a multi-level matrix of categories of land use) has been conducted.

A first attempt of geographically-distributed Multi-scale integrated analysis of societal and ecosystem metabolism (MuSIASEM): mapping human time and energy throughput in metropolitan Barcelona

This study presents a first attempt to extend the “Multi-scale integrated analysis of societal and ecosystem metabolism (MuSIASEM)” approach to a spatial dimension using GIS techniques in the Metropolitan area of Barcelona. We use a combination of census and commercial databases along with a detailed land cover map to create a layer of Common Geographic Units that we populate with the local values of human time spent in different activities according to MuSIASEM hierarchical typology. In this way, we mapped the hours of available human time, in regards to the working hours spent in different locations, putting in evidence the gradients in spatial density between the residential location of workers (generating the work supply) and the places where the working hours are actually taking place. We found a strong three-modal pattern of clumps of areas with different combinations of values of time spent on household activities and on paid work. We also measured and mapped spatial segregation between these two activities and put forward the conjecture that this segregation increases with higher energy throughput, as the size of the functional units must be able to cope with the flow of exosomatic energy. Finally, we discuss the effectiveness of the approach by comparing our geographic representation of exosomatic throughput to the one issued from conventional methods.