Oxidation of methane at the CH4/H2O-(CO2) transition zone in the external part of the Central Alps, Switzerland: Evidence from stable isotope investigations

With the aim of understanding the mechanisms that control the metamorphic transition from the CH4- to the H2O-(CO2)-dominated fluid zone in the Helvetic domain of the Central Alps of Switzerland, fluid inclusions in quartz, illite ``crystallinity'' index, vitrinite reflectance, and the stable isotope compositions of vein and whole rock minerals and fluids trapped in quartz were investigated along four cross-sections. Increasing temperature during prograde metamorphism led to the formation of dry gas by hydrocarbon cracking in the CH4-zone. Fluid immiscibility in the H2O-CH4-(CO2)-NaCl system resulted in cogenetic, CH4- and H2O-dominated fluid inclusions. In the CH4-zone, fluids were trapped at temperatures <= 270 +/- 5 degrees C. The end of the CH4-zone is markedby a sudden increase of CO2 content in the gas phase of fluid inclusions. At temperatures > 270 +/- 5 degrees C, in the H2O-zone, the total amount of volatiles within the fluid decreased below 1 mol% with no immiscibility. This resulted m total homogenization temperatures of H2O-(CO2-CH4)-NaCl inclusions below 180 degrees C. Hydrogen isotope compositions of methane in fluid inclusion have delta D values of less than -100 parts per thousand in the CH4-zone, typical for an origin through cracking of higher hydrocarbons, but where the methane has not equilibrated with the pore water. delta D values of fluid inclusion water are around -40 parts per thousand., in isotopic equilibrium with phyllosilicates of the whole rocks. Within the CH4 to H2O(CO2) transition zone, delta D(H2O) values in fluid inclusions decrease to -130 parts per thousand interpreted to reflect the contribution of deuterium depleted water from methane oxidation. In the H2O-zone, delta D(H2O) values increase again towards an average of -30 parts per thousand which is again consistent with isotopic equilibrium with host-rock phyllosilicates. delta C-13 values of methane in fluid inclusions from the CH4-zone are around -27 parts per thousand in isotopic equilibrium with calcite in veins and whole rocks. The delta C-13(CH4) values decrease to less than -35 parts per thousand at the transition to the H2O-zone and are no longer in equilibrium with the carbonates in the whole rocks. delta C-13 values of CO, are variable but too low to be in equilibrium with the wall rock fluids, compatible with a contribution of CO2 from closed system oxidation of methane. Differences in isotopic composition between host-rock and Alpine fissure carbonate are generally small, suggesting that the amount of CO2 produced by oxidation of methane was small compared to the C-budget in the rocks and local pore fluids were buffered by the wall rocks during precipitation of calcite within the fissures. (c) 2006 Elsevier B.V. All rights reserved.

Rift-related inheritance in orogens: a case study from the Austroalpine nappes in Central Alps (SE-Switzerland and N-Italy)

The Austroalpine nappe systems in SE-Switzerland and N-Italy preserve remnants of the Adriatic rifted margin. Based on new maps and cross-sections, we suggest that the complex structure of the Campo, Grosina/Languard, and Bernina nappes is inherited largely from Jurassic rifting. We propose a classification of the Austroalpine domain into Upper, Middle and Lower Austroalpine nappes that is new because it is based primarily on the rift-related Jurassic structure and paleogeography of these nappes. Based on the Alpine structures and pre-Alpine, rift-related geometry of the Lower (Bernina) and Middle (Campo, Grosina/Languard) Austroalpine nappes, we restore these nappes to their original positions along the former margin, as a means of understanding the formation and emplacement of the nappes during initial reactivation of the Alpine Tethyan margin. The Campo and Grosina/Languard nappes can be interpreted as remnants of a former necking zone that comprised pre-rift upper and middle crust. These nappes were juxtaposed with the Mesozoic cover of the Bernina nappe during Jurassic rifting. We find evidence for low-angle detachment faults and extensional allochthons in the Bernina nappe similar to those previously described in the Err nappe and explain their role during subsequent reactivation. Our observations reveal a strong control of rift-related structures during the subsequent Alpine reactivation on all scales of the former distal margin. Two zones of intense deformation, referred to as the Albula-Zebru and Lunghin-Mortirolo movement zones, have been reactivated during Alpine deformation and cannot be described as simple monophase faults or shear zones. We propose a tectonic model for the Austroalpine nappe systems that link inherited, rift-related structures with present-day Alpine structures. In conclusion, we believe that apart from the direct regional implications, the results of this paper are of general interest in understanding the control of rift structures during reactivation of distal-rifted margins.

Ecological-niche factor analysis: how to compute habitat-suitability maps without absence data?

We propose a multivariate approach to the study of geographic species distribution which does not require absence data. Building on Hutchinson's concept of the ecological niche, this factor analysis compares, in the multidimensional space of ecological variables, the distribution of the localities where the focal species was observed to a reference set describing the whole study area. The first factor extracted maximizes the marginality of the focal species, defined as the ecological distance between the species optimum and the mean habitat within the reference area. The other factors maximize the specialization of this focal species, defined as the ratio of the ecological variance in mean habitat to that observed for the focal species. Eigenvectors and eigenvalues are readily interpreted and can be used to build habitat-suitability maps. This approach is recommended in Situations where absence data are not available (many data banks), unreliable (most cryptic or rare species), or meaningless (invaders). We provide an illustration and validation of the method for the alpine ibex, a species reintroduced in Switzerland which presumably has not yet recolonized its entire range.

Multidimensional Deprivation in Contemporary Switzerland Across Social Groups and Time

We have investigated the phenomenon of deprivation in contemporary Switzerland through the adoption of a multidimensional, dynamic approach. By applying Self Organizing Maps (SOM) to a set of 33 non-monetary indicators from the 2009 wave of the Swiss Household Panel (SHP), we identified 13 prototypical forms (or clusters) of well-being, financial vulnerability, psycho-physiological fragility and deprivation within a topological dimensional space. Then new data from the previous waves (2003 to 2008) were classified by the SOM model, making it possible to estimate the weight of the different clusters in time and reconstruct the dynamics of stability and mobility of individuals within the map. Looking at the transition probabilities between year t and year t+1, we observed that the paths of mobility which catalyze the largest number of observations are those connecting clusters that are adjacent on the topological space.

Improved predictive mapping of indoor radon concentrations using ensemble regression trees based on automatic clustering of geological units.

PURPOSE: According to estimations around 230 people die as a result of radon exposure in Switzerland. This public health concern makes reliable indoor radon prediction and mapping methods necessary in order to improve risk communication to the public. The aim of this study was to develop an automated method to classify lithological units according to their radon characteristics and to develop mapping and predictive tools in order to improve local radon prediction. METHOD: About 240 000 indoor radon concentration (IRC) measurements in about 150 000 buildings were available for our analysis. The automated classification of lithological units was based on k-medoids clustering via pair-wise Kolmogorov distances between IRC distributions of lithological units. For IRC mapping and prediction we used random forests and Bayesian additive regression trees (BART). RESULTS: The automated classification groups lithological units well in terms of their IRC characteristics. Especially the IRC differences in metamorphic rocks like gneiss are well revealed by this method. The maps produced by random forests soundly represent the regional difference of IRCs in Switzerland and improve the spatial detail compared to existing approaches. We could explain 33% of the variations in IRC data with random forests. Additionally, the influence of a variable evaluated by random forests shows that building characteristics are less important predictors for IRCs than spatial/geological influences. BART could explain 29% of IRC variability and produced maps that indicate the prediction uncertainty. CONCLUSION: Ensemble regression trees are a powerful tool to model and understand the multidimensional influences on IRCs. Automatic clustering of lithological units complements this method by facilitating the interpretation of radon properties of rock types. This study provides an important element for radon risk communication. Future approaches should consider taking into account further variables like soil gas radon measurements as well as more detailed geological information.

Rockfall susceptibility assessment and remote geological mapping with LiDAR point clouds

Characterizing the geological features and structures in three dimensions over inaccessible rock cliffs is needed to assess natural hazards such as rockfalls and rockslides and also to perform investigations aimed at mapping geological contacts and building stratigraphy and fold models. Indeed, the detailed 3D data, such as LiDAR point clouds, allow to study accurately the hazard processes and the structure of geologic features, in particular in vertical and overhanging rock slopes. Thus, 3D geological models have a great potential of being applied to a wide range of geological investigations both in research and applied geology projects, such as mines, tunnels and reservoirs. Recent development of ground-based remote sensing techniques (LiDAR, photogrammetry and multispectral / hyperspectral images) are revolutionizing the acquisition of morphological and geological information. As a consequence, there is a great potential for improving the modeling of geological bodies as well as failure mechanisms and stability conditions by integrating detailed remote data. During the past ten years several large rockfall events occurred along important transportation corridors where millions of people travel every year (Switzerland: Gotthard motorway and railway; Canada: Sea to sky highway between Vancouver and Whistler). These events show that there is still a lack of knowledge concerning the detection of potential rockfalls, making mountain residential settlements and roads highly risky. It is necessary to understand the main factors that destabilize rocky outcrops even if inventories are lacking and if no clear morphological evidences of rockfall activity are observed. In order to increase the possibilities of forecasting potential future landslides, it is crucial to understand the evolution of rock slope stability. Defining the areas theoretically most prone to rockfalls can be particularly useful to simulate trajectory profiles and to generate hazard maps, which are the basis for land use planning in mountainous regions. The most important questions to address in order to assess rockfall hazard are: Where are the most probable sources for future rockfalls located? What are the frequencies of occurrence of these rockfalls? I characterized the fracturing patterns in the field and with LiDAR point clouds. Afterwards, I developed a model to compute the failure mechanisms on terrestrial point clouds in order to assess the susceptibility to rockfalls at the cliff scale. Similar procedures were already available to evaluate the susceptibility to rockfalls based on aerial digital elevation models. This new model gives the possibility to detect the most susceptible rockfall sources with unprecented detail in the vertical and overhanging areas. The results of the computation of the most probable rockfall source areas in granitic cliffs of Yosemite Valley and Mont-Blanc massif were then compared to the inventoried rockfall events to validate the calculation methods. Yosemite Valley was chosen as a test area because it has a particularly strong rockfall activity (about one rockfall every week) which leads to a high rockfall hazard. The west face of the Dru was also chosen for the relevant rockfall activity and especially because it was affected by some of the largest rockfalls that occurred in the Alps during the last 10 years. Moreover, both areas were suitable because of their huge vertical and overhanging cliffs that are difficult to study with classical methods. Limit equilibrium models have been applied to several case studies to evaluate the effects of different parameters on the stability of rockslope areas. The impact of the degradation of rockbridges on the stability of large compartments in the west face of the Dru was assessed using finite element modeling. In particular I conducted a back-analysis of the large rockfall event of 2005 (265'000 m3) by integrating field observations of joint conditions, characteristics of fracturing pattern and results of geomechanical tests on the intact rock. These analyses improved our understanding of the factors that influence the stability of rock compartments and were used to define the most probable future rockfall volumes at the Dru. Terrestrial laser scanning point clouds were also successfully employed to perform geological mapping in 3D, using the intensity of the backscattered signal. Another technique to obtain vertical geological maps is combining triangulated TLS mesh with 2D geological maps. At El Capitan (Yosemite Valley) we built a georeferenced vertical map of the main plutonio rocks that was used to investigate the reasons for preferential rockwall retreat rate. Additional efforts to characterize the erosion rate were made at Monte Generoso (Ticino, southern Switzerland) where I attempted to improve the estimation of long term erosion by taking into account also the volumes of the unstable rock compartments. Eventually, the following points summarize the main out puts of my research: The new model to compute the failure mechanisms and the rockfall susceptibility with 3D point clouds allows to define accurately the most probable rockfall source areas at the cliff scale. The analysis of the rockbridges at the Dru shows the potential of integrating detailed measurements of the fractures in geomechanical models of rockmass stability. The correction of the LiDAR intensity signal gives the possibility to classify a point cloud according to the rock type and then use this information to model complex geologic structures. The integration of these results, on rockmass fracturing and composition, with existing methods can improve rockfall hazard assessments and enhance the interpretation of the evolution of steep rockslopes. -- La caractérisation de la géologie en 3D pour des parois rocheuses inaccessibles est une étape nécessaire pour évaluer les dangers naturels tels que chutes de blocs et glissements rocheux, mais aussi pour réaliser des modèles stratigraphiques ou de structures plissées. Les modèles géologiques 3D ont un grand potentiel pour être appliqués dans une vaste gamme de travaux géologiques dans le domaine de la recherche, mais aussi dans des projets appliqués comme les mines, les tunnels ou les réservoirs. Les développements récents des outils de télédétection terrestre (LiDAR, photogrammétrie et imagerie multispectrale / hyperspectrale) sont en train de révolutionner l'acquisition d'informations géomorphologiques et géologiques. Par conséquence, il y a un grand potentiel d'amélioration pour la modélisation d'objets géologiques, ainsi que des mécanismes de rupture et des conditions de stabilité, en intégrant des données détaillées acquises à distance. Pour augmenter les possibilités de prévoir les éboulements futurs, il est fondamental de comprendre l'évolution actuelle de la stabilité des parois rocheuses. Définir les zones qui sont théoriquement plus propices aux chutes de blocs peut être très utile pour simuler les trajectoires de propagation des blocs et pour réaliser des cartes de danger, qui constituent la base de l'aménagement du territoire dans les régions de montagne. Les questions plus importantes à résoudre pour estimer le danger de chutes de blocs sont : Où se situent les sources plus probables pour les chutes de blocs et éboulement futurs ? Avec quelle fréquence vont se produire ces événements ? Donc, j'ai caractérisé les réseaux de fractures sur le terrain et avec des nuages de points LiDAR. Ensuite, j'ai développé un modèle pour calculer les mécanismes de rupture directement sur les nuages de points pour pouvoir évaluer la susceptibilité au déclenchement de chutes de blocs à l'échelle de la paroi. Les zones sources de chutes de blocs les plus probables dans les parois granitiques de la vallée de Yosemite et du massif du Mont-Blanc ont été calculées et ensuite comparés aux inventaires des événements pour vérifier les méthodes. Des modèles d'équilibre limite ont été appliqués à plusieurs cas d'études pour évaluer les effets de différents paramètres sur la stabilité des parois. L'impact de la dégradation des ponts rocheux sur la stabilité de grands compartiments de roche dans la paroi ouest du Petit Dru a été évalué en utilisant la modélisation par éléments finis. En particulier j'ai analysé le grand éboulement de 2005 (265'000 m3), qui a emporté l'entier du pilier sud-ouest. Dans le modèle j'ai intégré des observations des conditions des joints, les caractéristiques du réseau de fractures et les résultats de tests géoméchaniques sur la roche intacte. Ces analyses ont amélioré l'estimation des paramètres qui influencent la stabilité des compartiments rocheux et ont servi pour définir des volumes probables pour des éboulements futurs. Les nuages de points obtenus avec le scanner laser terrestre ont été utilisés avec succès aussi pour produire des cartes géologiques en 3D, en utilisant l'intensité du signal réfléchi. Une autre technique pour obtenir des cartes géologiques des zones verticales consiste à combiner un maillage LiDAR avec une carte géologique en 2D. A El Capitan (Yosemite Valley) nous avons pu géoréferencer une carte verticale des principales roches plutoniques que j'ai utilisé ensuite pour étudier les raisons d'une érosion préférentielle de certaines zones de la paroi. D'autres efforts pour quantifier le taux d'érosion ont été effectués au Monte Generoso (Ticino, Suisse) où j'ai essayé d'améliorer l'estimation de l'érosion au long terme en prenant en compte les volumes des compartiments rocheux instables. L'intégration de ces résultats, sur la fracturation et la composition de l'amas rocheux, avec les méthodes existantes permet d'améliorer la prise en compte de l'aléa chute de pierres et éboulements et augmente les possibilités d'interprétation de l'évolution des parois rocheuses.

The geologic and geomorphologic maps of Santa Coloma de Farners (Girona, Spain)

Notice about the geologic and geomorphologic maps of Santa Coloma de Farners, at scale 1:10,000 published by Unitat de Geologia of the Universitat de Girona

Geologic maps published by university de Girona's unity of geodynamics

In 1978 University de Girona's Unity of Geodynamics started publishing geological maps, in several themes and scales, of various areas of the Province of Girona (fig. 1). So far, up to 38 maps, gathered in several collections, have been published

Weighted Logratio Biplots, Correspondence Analysis and Spectral Maps

Starting with logratio biplots for compositional data, which are based on the principle of subcompositional coherence, and then adding weights, as in correspondence analysis, we rediscover Lewi's spectral map and many connections to analyses of two-way tables of non-negative data. Thanks to the weighting, the method also achieves the property of distributional equivalence

Using self organizing maps on compositional data

Self-organizing maps (Kohonen 1997) is a type of artificial neural network developed to explore patterns in high-dimensional multivariate data. The conventional version of the algorithm involves the use of Euclidean metric in the process of adaptation of the model vectors, thus rendering in theory a whole methodology incompatible with non-Euclidean geometries. In this contribution we explore the two main aspects of the problem: 1. Whether the conventional approach using Euclidean metric can shed valid results with compositional data. 2. If a modification of the conventional approach replacing vectorial sum and scalar multiplication by the canonical operators in the simplex (i.e. perturbation and powering) can converge to an adequate solution. Preliminary tests showed that both methodologies can be used on compositional data. However, the modified version of the algorithm performs poorer than the conventional version, in particular, when the data is pathological. Moreover, the conventional ap- proach converges faster to a solution, when data is \well-behaved". Key words: Self Organizing Map; Artificial Neural networks; Compositional data

PhotoLin: a program to identify and analyze linear structures in aerial photographs, satellite images and maps

A computer program, PhotoLin, written for an IBM-PC-compatible microcomputer is described which detects linear features in aerial photographs, satellite images and topographic maps. The program accepts images saved to PCX files as input and applies noise correction and smoothing filters and thinning routines. The output consists of a skeleton containing the median lines of linear features which can be represented on a map. The branches of the skeleton can be broken into sections of constant length for which the mean orientations are obtained for the preparation of rose diagrams. (C) 2001 Elsevier B.V. Ltd. All rights reserved.

Road maps towards an information society in Latin America and the Caribbean

Includes bibliography