Magmatic-dominated fluid evolution in the Jurassic Nambija gold skarn deposits (southeastern Ecuador)

The Jurassic (approximately 145 Ma) Nambija oxidized gold skarns are hosted by the Triassic volcanosedimentary Piuntza unit in the sub-Andean zone of southeastern Ecuador. The skarns consist dominantly of granditic garnet (Ad(20-98)) with subordinate pyroxene (Di(46-92)Hd(17-42)Jo(0-19)) and epidote and are spatially associated with porphyritic quartz-diorite to granodiorite intrusions. Endoskarn is developed at the intrusion margins and grades inwards into a potassic alteration zone. Exoskarn has an outer K- and Na-enriched zone in the volcanosedimentary unit. Gold mineralization is associated with the weakly developed retrograde alteration of the exoskarn and occurs mainly in sulfide-poor vugs and milky quartz veins and veinlets in association with hematite. Fluid inclusion data for the main part of the prograde stage indicate the coexistence of high-temperature (500A degrees C to > 600A degrees C), high-salinity (up to 65 wt.% eq. NaCl), and moderate- to low-salinity aqueous-carbonic fluids interpreted to have been trapped at pressures around 100-120 MPa, corresponding to about 4-km depth. Lower-temperature (510-300A degrees C) and moderate- to low-salinity (23-2 wt.% eq. NaCl) aqueous fluids are recorded in garnet and epidote of the end of the prograde stage. The microthermometric data (Th from 513A degrees C to 318A degrees C and salinity from 1.0 to 23 wt.% eq. NaCl) and delta(18)O values between 6.2aEuro degrees and 11.5aEuro degrees for gold-bearing milky quartz from the retrograde stage suggest that the ore-forming fluid was dominantly magmatic. Pressures during the early retrograde stage were in the range of 50-100 MPa, in line with the evidence for CO(2) effervescence and probable local boiling. The dominance of magmatic low-saline to moderately saline oxidizing fluids during the retrograde stage is consistent with the depth of the skarn system, which could have delayed the ingression of external fluids until relatively low temperatures were reached. The resulting low water-to-rock ratios explain the weak retrograde alteration and the compositional variability of chlorite, essentially controlled by host rock compositions. Gold was precipitated at this stage as a result of cooling and pH increase related to CO(2) effervescence, which both result in destabilization of gold-bearing chloride complexes. Significant ingression of external fluids took place after gold deposition only, as recorded by delta(18)O values of 0.4aEuro degrees to 6.2aEuro degrees for fluids depositing quartz (below 350A degrees C) in sulfide-rich barren veins. Low-temperature (< 300A degrees C) meteoric fluids (delta(18)O(water) between -10.0aEuro degrees and -2.0aEuro degrees) are responsible for the precipitation of late comb quartz and calcite in cavities and veins and indicate mixing with cooler fluids of higher salinities (about 100A degrees C and 25 wt.% eq. NaCl). The latter are similar to low-temperature fluids (202-74.5A degrees C) with delta(18)O values of -0.5aEuro degrees to 3.1aEuro degrees and salinities in the range of 21.1 to 17.3 wt.% eq. CaCl(2), trapped in calcite of late veins and interpreted as basinal brines. Nambija represents a deep equivalent of the oxidized gold skarn class, the presence of CO(2) in the fluids being partly a consequence of the relatively deep setting at about 4-km depth. As in other Au-bearing skarn deposits, not only the prograde stage but also the gold-precipitating retrograde stage is dominated by fluids of magmatic origin.

Fluid-rock-strain feedbacks in extensional shear zones

La présence de fluide météorique synchrone à l'activité du détachement (Farmin, 2003 ; Mulch et al., 2007 ; Gébelin et al., 2011), implique que les zones de cisaillement sont des systèmes ouverts avec des cellules de convections à l'échelle crustale et un intense gradient géothermique au sein du détachement (Morrison et Anderson, 1998, Gottardi et al., 2011). De plus, les réactions métamorphiques liées à des infiltrations fluides dans les zones de cisaillement extensionnel peuvent influencer les paramètres rhéologiques du système (White and Knipe, 1978), et impliquer la localisation de la déformation dans la croûte. Dans ce manuscrit, deux zones de cisaillement infiltrées par des fluides météoriques sont étudiées, l'une étant largement quartzitique, et l'autre de nature granitique ; les relations entre déformation, fluides, et roches s'appuient sur des approches structurales, microstructurales, chimiques et isotopiques. L'étude du détachement du Columbia river (WA, USA) met en évidence que la déformation mylonitique se développe en un million d'années. La phase de cisaillement principal s'effectue à 365± 30°C d'après les compositions isotopiques en oxygène du quartz et de la muscovite. Ces minéraux atteignent l'équilibre isotopique lors de leur recristallisation dynamique contemporaine à la déformation. La zone de cisaillement enregistre une baisse de température, remplaçant le mécanisme de glissement par dislocation par celui de dissolution- précipitation dans les derniers stades de l'activité du détachement. La dynamique de circulation fluide bascule d'une circulation pervasive à chenalisée, ce qui engendre localement la rupture des équilibres d'échange isotopiques. La zone de cisaillement de Bitterroot (MT, USA) présente une zone mylonitique de 600m d'épaisseur, progressant des protomylonites aux ultramylonites. L'intensité de la localisation de la déformation se reflète directement sur l'hydratation des feldspaths, réaction métamorphique majeure dite de « rock softening ». Une étude sur roche totale indique des transferts de masse latéraux au sein des mylonites, et d'importantes pertes de volume dans les ultramylonites. La composition isotopique en hydrogène des phyllosilicates met en évidence la présence (1) d'une source magmatique/métamorphique originelle, caractérisée par les granodiorites ayant conservé leur foliation magmatique, jusqu'aux protomylonites, et (2) une source météorique qui tamponne les valeurs des phyllosilicates des fabriques mylonitiques jusqu'aux veines de quartz non-déformées. Les compositions isotopiques en oxygène des minéraux illustrent le tamponnement de la composition du fluide météorique par l'encaissant. Ce phénomène cesse lors du processus de chloritisation de la biotite, puisque les valeurs des chlorites sont extrêmement négatives (-10 per mil). La thermométrie isotopique indique une température d'équilibre isotopique de la granodiorite entre 600-500°C, entre 500-300°C dans les mylonites, et entre 300 et 200°C dans les fabriques cassantes (cataclasites et veines de quartz). Basé sur les résultats issus de ce travail, nous proposons un modèle général d'interactions fluide-roches-déformation dans les zones de détachements infiltrées par des fluides météoriques. Les zones de détachements évoluent rapidement (en quelques millions d'années) au travers de la transition fragile-ductile ; celle-ci étant partiellement contrôlée par l'effet thermique des circulations de fluide météoriques. Les systèmes de détachements sont des lieux où la déformation et les circulations fluides sont couplées ; évoluant rapidement vers une localisation de la déformation, et de ce fait, une exhumation efficace. - The presence of meteoric fluids synchronous with the activity of extensional detachment zones (Famin, 2004; Mulch et al., 2007; Gébelin et al., 2011) implies that extensional systems involve fluid convection at a crustal scale, which results in high geothermal gradients within active detachment zones (Morrison and Anderson, 1998, Gottardi et al., 2011). In addition, the metamorphic reactions related to fluid infiltration in extensional shear zones can influence the rheology of the system (White and Knipe, 1978) and ultimately how strain localizes in the crust. In this thesis, two shear zones that were permeated by meteoric fluids are studied, one quartzite-dominated, and the other of granitic composition; the relations between strain, fluid, and evolving rock composition are addressed using structural, microstructural, and chemical/isotopic measurements. The study of the Columbia River detachment that bounds the Kettle core complex (Washington, USA) demonstrates that the mylonitic fabrics in the 100 m thick quartzite- dominated detachment footwall developed within one million years. The main shearing stage occurred at 365 ± 30°C when oxygen isotopes of quartz and muscovite equilibrated owing to coeval deformation and dynamic recrystallization of these minerals. The detachment shear zone records a decrease in temperature, and dislocation creep during detachment shearing gave way to dissolution-precipitation and fracturing in the later stages of detachment activity. Fluid flow switched from pervasive to channelized, leading to isotopic disequilibrium between different minerals. The Bitterroot shear zone detachment (Montana, USA) developed a 600 m thick mylonite zone, with well-developed transitions from protomylonite to ultramylonite. The localization of deformation relates directly to the intensity of feldspar hydration, a major rock- softening metamorphic reaction. Bulk-rock analyses of the mylonitic series indicate lateral mass transfer in the mylonite (no volume change), and significant volume loss in ultramylonite. The hydrogen isotope composition of phyllosilicates shows (1) the presence of an initial magmatic/metamorphic source characterized by the granodiorite in which a magmatic, and gneissic (protomylonite) foliation developed, and (2) a meteoric source that buffers the values of phyllosilicates in mylonite, ultramylonite, cataclasite, and deformed and undeformed quartz veins. The mineral oxygen isotope compositions were buffered by the host-rock compositions until chloritization of biotite started; the chlorite oxygen isotope values are negative (-10 per mil). Isotope thermometry indicates a temperature of isotopic equilibrium of the granodiorite between 600-500°C, between 500-300°C in the mylonite, and between 300 and 200°C for brittle fabrics (cataclasite and quartz veins). Results from this work suggest a general model for fluid-rock-strain feedbacks in detachment systems that are permeated by meteoric fluids. Phyllosilicates have preserved in their hydrogen isotope values evidence for the interaction between rock and meteoric fluids during mylonite development. Fluid flow generates mass transfer along the tectonic anisotropy, and mylonites do not undergo significant volume change, except locally in ultramylonite zones. Hydration of detachment shear zones attends mechanical grain size reduction and enhances strain softening and localization. Self-exhuming detachment shear zones evolve rapidly (a few million years) through the transition from ductile to brittle, which is partly controlled by the thermal effect of circulating surface fluids. Detachment systems are zones in the crust where strain and fluid flow are coupled; these systems. evolve rapidly toward strain localization and therefore efficient exhumation.

Geochemical and H-O-Sr-Nd isotope evidence for magmatic processes and meteoric-water interactions in the basal complex of La Gomera, Canary Islands

The plutonic rocks of the Basal Complex of La Gomera, Canary Islands, Spain, were studied by means of major and trace element contents and by H-O-Sr-Nd isotope compositions in order to distinguish primary magmatic characteristics and late-stage alteration products. Deciphering the effects of alteration allowed us to determine primary, plume-related compositions that indicated D- and (18)O-depletion relative to normal upper mantle, supporting the conclusions of earlier studies on the plutonic rocks of Fuerteventura and La Palma. Late-stage alteration took place during the formation of the intrusive series induced by interaction with meteoric water. Inferred isotopic compositions of the meteoric water indicate that the water infiltrated into the rock edifice at a height of about 1500 m above sea level, suggesting the existence of a subaerial volcano which was active during the intrusive activity and that it has been either distroyed or remain buried by later volcanic and landslide events.

H2O-δD-Fe-III relations of dehydrogenation and dehydration processes in magmatic amphiboles

Stable isotope compositions of a suite of magmatic amphiboles from alkaline basalts and andesitic rocks were examined to constrain the effects of degassing processes on the hydrogen isotope compositions. The Fe3+ (as Fe3+/Fe-total) and H2O contents, as well as the H isotope compositions of the amphiboles, differ markedly (27-58%, 0.5-2.2 wt%, -107 to -15 parts per thousand, respectively) but indicate systematic variations. The observed trends can be explained either as dehydrogenation or dehydration processes, both of which are coupled to oxidation processes, the latter most probably related to O2- substitution within amphiboles. The dehydrogenation-dehydration models can be used to assess the primary compositions of the magmas. As an important example, delta D values of amphiboles of Martian meteorites are discussed in a similar context. Copyright (c) 2006 John Wiley & Sons, Ltd.

Evolution of a polymineralic mantle shear zone and the role of second phases in the localization of deformation

The influence of second phases (e.g., pyroxenes) on olivine grain size was studied by quantitative microfabric analyses of samples of the Hilti massif mantle shear zone (Semail ophiolite, Oman). The microstructures range from porphyroclastic tectonites to ultramylonites, from outside to the center of the shear zone. Starting at conditions of ridge-related flow, they formed under continuous cooling leading to progressive strain localization. The dependence of the average olivine grain size on the second-phase content can be split into a second-phase controlled and a dynamic recrystallization-controlled field. In the former, the olivine grain size is related to the ratio between the second-phase grain size and volume fraction (Zener parameter). In the latter, dynamic recrystallization manifested by a balance between grain growth and grain size reduction processes yields a stable olivine grain size. In both fields the average olivine and second-phase grain size decreases with decreasing temperature. Combining the microstructural information with deformation mechanism maps suggests that the porphyroclastic tectonites (similar to 1100 degrees C) and mylonites (similar to 800 degrees C) formed under the predominance of dislocation creep. Since olivine-rich layers are intercalated with layer parallel, polymineralic bands in the mylonites, nearly equiviscous conditions can be assumed. In the ultramylonites, diffusion creep represents the major deformation mechanism in the polymineralic layers. It is this switch in deformation mechanism from dislocation creep to diffusion creep that forces strain to localize in the fine-grained polymineralic domains at low temperatures (<similar to 700 degrees C), underlining the role of the second phases on strain localization in cooling mantle rocks.

Active deformation fields: dense deformation field estimation for atlas-based segmentation using the active contour framework.

This paper presents a new and original variational framework for atlas-based segmentation. The proposed framework integrates both the active contour framework, and the dense deformation fields of optical flow framework. This framework is quite general and encompasses many of the state-of-the-art atlas-based segmentation methods. It also allows to perform the registration of atlas and target images based on only selected structures of interest. The versatility and potentiality of the proposed framework are demonstrated by presenting three diverse applications: In the first application, we show how the proposed framework can be used to simulate the growth of inconsistent structures like a tumor in an atlas. In the second application, we estimate the position of nonvisible brain structures based on the surrounding structures and validate the results by comparing with other methods. In the final application, we present the segmentation of lymph nodes in the Head and Neck CT images, and demonstrate how multiple registration forces can be used in this framework in an hierarchical manner.

Quantification of magmatic and hydrothermal processes in a peralkaline syenite-alkali granite complex based on textures, phase equilibria, and stable and radiogenic isotopes

The Puklen complex of the Mid-Proterozoic Gardar Province, South Greenland, consists of various silica-saturated to quartz-bearing syenites, which are intruded by a peralkaline granite. The primary mafic minerals in the syenites are augite +/- olivine + Fe-Ti oxide + amphibole. Ternary feldspar thermometry and phase equilibria among mafic silicates yield T = 950-750degreesC, a(SiO2) = 0.7-1 and an f(O2) of 1-3 log units below the fayalite-magnetite-quartz (FMQ) buffer at 1 kbar. In the granites, the primary mafic minerals are ilmenite and Li-bearing arfvedsonite, which crystallized at temperatures below 750degreesC and at f(O2) values around the FMQ buffer. In both rock types, a secondary post-magmatic assemblage overprints the primary magmatic phases. In syenites, primary Ca-bearing minerals are replaced by Na-rich minerals such as aegirine-augite and albite, resulting in the release of Ca. Accordingly, secondary minerals include ferro-actinolite, (calcite-siderite)(ss), titanite and andradite in equilibrium with the Na-rich minerals. Phase equilibria indicate that formation of these minerals took place over a long temperature interval from near-magmatic temperatures down to similar to300degreesC. In the course of this cooling, oxygen fugacity rose in most samples. For example, late-stage aegirine in granites formed at the expense of arfvedsonite at temperatures below 300degreesC and at an oxygen fugacity above the haematite-magnetite (HM) buffer. The calculated delta(18)O(melt) value for the syenites (+5.9 to +6.3parts per thousand) implies a mantle origin, whereas the inferred delta(18)O(melt) value of <+5.1parts per thousand for the granitic melts is significantly lower. Thus, the granites require an additional low-delta(18)O contaminant, which was not involved in the genesis of the syenites. Rb/Sr data for minerals of both rock types indicate open-system behaviour for Rb and Sr during post-magmatic metasomatism. Neodymium isotope compositions (epsilonNd(1170 Ma) = -3.8 to -6.4) of primary minerals in syenites are highly variable, and suggest that assimilation of crustal rocks occurred to variable extents. Homogeneous epsilon(Nd) values of -5.9 and -6.0 for magmatic amphibole in the granites lie within the range of the syenites. Because of the very similar neodymium isotopic compositions of magmatic and late- to post-magmatic minerals from the same syenite samples a principally closed-system behaviour during cooling is implied. In contrast, for the granites an externally derived fluid phase is required to explain the extremely low epsilon(Nd) values of about -10 and low delta(18)O between +2.0 and +0.5parts per thousand for late-stage aegirine, indicating an open system in the late-stage history. In this study we show that the combination of phase equilibria constraints with stable and radiogenic isotope data on mineral separates can provide much better constraints on magma evolution during emplacement and crystallization than conventional whole-rock studies.

Magmatic fluids in the Breccia-hosted epithermal Au-Ag deposit of Rosia Montana, Romania

The breccia-hosted epithermal Au-Ag deposit of Rosia Montana is located 7 kin northeast of Abrud, in the northern part of the South Apuseni Mountains, Romania. Estimated total reserves of 214.91 million metric toils (Mt) of ore at 1.46 g/t An and 6.9 g/t Ag (10.1 Moz of An and 47.6 Moz of Ag) make Rosia Montana one of the largest gold deposits in Europe. At this location, Miocene calc-alkaline magmatic and hydrothermal activity was associated with local extensional tectonics within a strike-slip regime related to the indentation of the Adriatic microplate into the European plate during the Carpathian orogenesis. The host rocks of the magmatic complex consist of pre-Mesozoic metamorphosed continental crust covered by Cretaceous turbiditic sediment (flysch). Magmatic activity at Rosia Montana and its surroundings occurred in several pulses and lasted about 7 m.y, Rosia Montana is a breccia-hosted epithermal system related to strong phreatomagmatic activity due to the shallow emplacement of the Montana dacite. The Montana dacite intruded Miocene volcaniclastic material (volcaniclastic breccias) and crops out at Cetate and Carnic Hills. Current mining is focused primarily on the Cetate open pit, which was mapped in detail, leading to the recognition of three distinct breccia bodies: the dacite breccia with a dominantly hydrothermal matrix, the gray polymict breccia with a greater proportion of sand-sized matrix support, and the black polymict breccia, which reached to the surface, contains carbonized tree trunks and has a dominantly barren elastic matrix. The hydrothermal alteration is pervasive. Adularia alteration with a phyllic overprint is ubiquitous; silicification and argillic alteration occur locally. Mineralization consists of quartz, adularia, carbonates (commonly Mn-rich), pyrite, Fe-poor sphalerite, galena, chalcopyrite, tetrahedrite, and native gold and occurs as disseminations, as well as in veins and filling vugs within the Montana dacite and the different breccias. The age of mineralization (12.85 +/- 0.07 Ma) was determined by Ar-40- Ar-39 dating on hydrothermal adularia crystals from vugs in the dacite breccia in the Cetate open pit. Microthermometric measurements of fluid inclusions in quartz phenocrysts from the Montana dacite revealed two fluid types that are absent from the hydrothermal breccia and must have been trapped at depth prior to dacite dome emplacement: brine inclusions (32-55 -wt % NaCl equiv, homogenizing at T-h > 460 degrees C) and intermediate density fluids (4.9-15.6 wt % NaCl equiv, T, between 345 degrees-430 degrees C). Secondary aqueous fluid inclusion assemblages in the phenocrysts have salinities of 0.2 to 2.2 wt percent NaCl equiv and T-h of 200 degrees to 280 degrees C. Fluid inclusion assemblages in hydrothermal quartz from breccias and veins have salinities of 0.2 to 3.4 wt percent NaCl equiv and T-h, from 200 degrees to 270 degrees C. The oxygen isotope composition of several zones of an ore-related epithermal quartz crystal indicate a very constant delta O-18 of 4.5 to 5.0 per mil for the mineralizing fluid, despite significant salinity and temperature variation over time. Following microthermometry, selected fluid inclusion assemblages were analyzed by laser ablation-inductively coupled-plasma mass spectrometry (LA-ICMS). Despite systematic differences in salinity between phenocryst-hosted fluids trapped at depth and fluids from quartz in the epithermal breccias, all fluids have overlapping major and trace cation ratios, including identical Na/K/Rb/Sr/Cs/Ba. Consistent with the constant near-magmatic oxygen isotope composition of the hydrothermal fluids, these data strongly indicate a common magmatic component of these chemically conservative solutes in all fluids. Cu, Pb, Zn, and Mn show variations in concentration relative to the relatively non-reactive alkalis, reflecting the precipitation of sulfide minerals together with An in the epithermal breccia, and possibly of Cu in an inferred subjacent porphyry environment. The magmatic-hydrothermal processes responsible for epithermal Au-Ag mineralization at Rosia Montana are, however, not directly related to the formation of the spatially associated porphyry Cu-Au deposit of Rosia Poieni, which occurred lout 3 m.y. later.

Structure and kinematics of the Siviez-Mischabel nappe in the Mattertal (western Alps of Switzerland) = Structure et cinématique de la nappe de Siviez-Mischabel le long de la vallée de Zermatt (Alpes de Suisse occidentale)

Résumé Cette étude porte sur le flanc inverse de la nappe de Siviez-Mischabel et sur les unités tectoniques sous jacentes (zone de Stalden supérieur et zone Houillère) dans la vallée menant à Zermatt. L'étude structurale du granite permien de Randa (orthogneiss oeillé) permet de mieux comprendre les effets de la déformation alpine sur les roches de socle. La cartographie détaillée de l'orthogneiss et de son encaissant, ainsi que l'étude lithostratigraphique des terrains sédimentaires associés permettent de proposer un schéma structural et cinématique du flanc inverse de la nappe de Siviez-Mischabel et de mieux comprendre ses relations avec les unités tectoniques sous-jacentes. L'analyse structurale de l'orthogneiss de Randa et de son encaissant révèle la superposition de plusieurs phases de déformation ductile. Cet orthogneiss formé sous des conditions métamorphiques du faciès schiste vert possède une forte schistosité alpine avec au moins deux linéations d'extension. La première, L1, orientée NW-SE est associée à la mise en place de la nappe. La seconde, L2, orientée SW-NE, se corrèle au cisaillement ductile du Simplon. La quantification de la déformation au moyen de la méthode de Fry sur les faciès porphyriques donne des ellipses à rapports axiaux compris entre 1.9 et 5.3, en accord avec les valeurs obtenues par d'autres marqueurs {tourmalines étirées, fibres). Les valeurs mesurées parallèlement à L1 ou L2 sont très semblables. La méthode de Fry a nécessité une étude théorique préalable afin de vérifier son applicabilité aux orthogneiss oeillés. La méthode requiert une distribution spatiale homogène et isotrope des marqueurs utilisés. Les tests statistiques effectués ont révélé que les phénocristaux de feldspath alcalin satisfont à cette condition et qu'ils peuvent être utilisés comme marqueur de la déformation au moyen de la méthode de Fry. Les valeurs obtenues révèlent l'importance du cisaillement ductile du Simplon sur la géométrie de la nappe dans la région d'étude. Le levé cartographique a permis d'améliorer la lithostratigraphie de la base de la nappe de Siviez-Mischabel. Trois formations en position renversée peuvent être observées sous les gneiss formant le coeur de la nappe. Ces trois formations forment le coeur du synclinal de St-Niklaus qui connecte la nappe de Siviez-Mischabel à la zone de Stalden supérieur. La datation par U-Pb de zircons détritiques et magmatiques par LA-ICP-MS permet de contraindre l'âge des formations observées (probablement Carbonifère à Trias précoce). Ces données ont des répercussions importantes sur la structure de la nappe dans la région, prouvant l'existence de plusieurs plis avec des séries normales et renversées bien préservées. La définition et la datation de ces formations, ainsi que leur identification dans la-Zone- Houillère avoisinante permettent de mieux comprendre la géométrie initiale et les relations tectoniques des nappes du Pennique moyen dans la vallée de Zermatt. Summary This study investigates the overturned limb of the Siviez-Mischabel nappe and underlying tectonic units (Upper Stalden zone and Houillère zone) in the Mattertal area. Detailed structural analysis in the Permian Randa granite (augen orthogneiss) allows a better understanding of the Alpine deformation effects on basement rocks. Detailed mapping of this orthogneiss and surrounding rocks, and the study of the lithostratigraphy in the related sedimentary horizons allow the proposition of a structural and kinematic model for the overturned limb of the Siviez-Mischabel and to better understand the relations with the underlying tectonic units. The structural analysis of the Randa orthogneiss and surrounding rocks revealed the superposition of several phases of ductile deformation. This orthogneiss formed under greenschist facies metamorphic conditions displays a strong Alpine foliation with at least two stretching lineations. The first lineation, L1, is oriented NW-SE and is related to the nappe emplacement northward. The second one, L2, is related to the Simplon ductile shear zone. Strain estimation using the Fry method has been performed on porphyritic facies of the Randa orthogneiss. The obtained ellipses have axial ratios varying between 1.9 and 5.3, in agreement with strain estimation obtained from other markers (stretched turmalines, fringes). The strain values are very similar if measured parallel to L1 or to L2. A theoretical approach was necessary to verify the relevant application of the Fry method to augen orthogneiss. This method requires that the distribution of the used markers has to be homogeneous and isotropic. Statistical tests have been done and revealed that K-feldspar phenocrysts satisfy these conditions and can be used as strain markers with the Fry method. The obtained strain measurements revealed the importance of the Simplon ductile shear zone on the geometry of the nappe in the studied area. Mapping has improved the lithostratigraphy at the base of the Siviez-Mischabel nappe. Three overturned formations can be observed below the gneisses forming the core of the nappe. These three formations form the St-Niklaus syncline, which connects the Siviez-Mischabel nappe to the underlying Upper Stalden zone. U-Pb dating of detrital and magmatic zircons by LA-ICPMS allowed the age of the observed formations to be constrained (presumably Carboniferous to Early Triassic). This data has critical implications for nappe structure in the region, composed of few recumbent folds with well preserved normal and overturned limbs. The definition and dating of these formations, as well as their identification in the adjacent "Houillère Zone" improve the understanding of the geometry and tectonic relations of the Middle Penninic nappes in the Mattertal.

Chemical and isotopic equilibrium between CO2 and CH4 in fumarolic gas discharges: Generation of CH4 in arc magmatic-hydrothermal systems

The chemical and isotopic composition of fumarolic gases emitted from Nisyros Volcano, Greece, and of a single gas sample from Vesuvio, Italy, was investigated in order to determine the origin of methane (CH,) within two subduction-related magmatic-hydrothermal environments. Apparent temperatures derived from carbon isotope partitioning between CH4 and CO2 of around 340degreesC for Nisyros and 470degreesC for Vesuvio correlate well with aquifer temperatures as measured directly and/or inferred from compositional data using the H2O-H-2-CO2-CO-CH4 geothermometer. Thermodynamic modeling reveals chemical equilibrium between CH4, CO2 and H2O implying that carbon isotope partitioning between CO2 and CH, in both systems is controlled by aquifer temperature. N-2/(3) He and CH4/(3) He ratios of Nisyros fumarolic gases are unusually low for subduction zone gases and correspond to those of midoceanic ridge environments. Accordingly, CH4 may have been primarily generated through the reduction of CO, by H, in the absence of any organic matter following a Fischer-Tropsch-type reaction. However, primary occurrence of minor amounts of thermogenic CH4 and subsequent re-equilibration with co-existing CO2 cannot be ruled out entirely- CO2/He-3 ratios and delta(13)C(CO2) values imply that the evolved CO2 either derives from a metasomatized mantle or is a mixture between two components, one outgassing from an unaltered mantle and the other released by thermal breakdown of marine carbonates. The latter may contain traces of organic matter possibly decomposing to CH4 during thermometamorphism. Copyright (C) 2004 Elsevier Ltd.

The Paleozoic magmatic history of the Maggia and Sambuco nappes, Lower Penninic, Central Lepontine Alps

Résumé : Les corps magmatiques sont des indicateurs essentiels dans toute reconstitution paléogéographique et/ou géodynamique d'un cycle orogénique, en particulier en contexte polycyclique, où la plupart des autres indices ont été oblitérés. Ils sont aisément datables et leurs caractéristiques géochimiques permettent de contraindre leur contexte tectonique de mise en place. Cette approche a été appliquée aux socles pré-mésozoïques des nappes penniques inférieures de Sambuco et de la Maggia, dans les Alpes centrales lepontines. Plusieurs événements magmatiques ont été identifiés dans le socle de Sambuco et datés par la méthode U-Pb sur zircon couplée à la technique LA-ICPMS. La suite calco-alcaline mafique rubanée de Scheggia est datée du Cambrien inférieur à 540-530 Ma ; le métagranite alumineux oeillé de Sasso Nero a un âge de 480-470 Ma, tout comme bien d'autres «older orthogneisses» des socles alpins. Il contient des zircons hérités d'âge panafricain à 630-610 Ma, indicateur d'une affiliation gondwanienne de ces terrains. Le pluton calco-alcalin du Matorello est daté à environ 300-310 Ma, et les filons lamprophyriques qu'il abrite à 300 Ma. La granodiorite de Cocco et le leucogranite de Ruscada, tous deux intrudés dans le socle de la nappe adjacente de la Maggia, ont des âges similaires à celui du Matorello. Ceci ajouté aux similitudes magmatiques observées entre Cocco et Matorello suggère une proximité paléogéographique des deux nappes au Permien-Carbonifère. Or ces dernières sont actuellement considérées appartenir à deux domaines paléogéographiques mésozoïques distincts : helvétique pour Sambuco et briançonnais pour Maggia, séparés par un bassin océanique. Si tel fut le cas, aucun mouvement décrochant ne doit avoir décalé les marges continentales de l'océan, retrouvées en parfaite coïncidence lors de sa fermeture. Le Matorello est un pluton recristallisé en faciès amphibolite et plissé par cinq phases successives de déformation non-coaxiales, qui ont conduit à son renversement complet, attesté par des indicateurs de paléogravité. Il préserve de spectaculaires phénomènes de coexistence liquide de magmas (essaims d'enclaves et Bills composites). Ce pluton était originellement tabulaire, construit par l'accumulation de multiples injections de magma en feuillets d'épaisseur métrique à décamétrique. Suivant le rythme de mise en place, les injections successives ont rapidement cristallisé avec des contours nets et bien définis (Bills composites) ou se sont mélangées avec les précédentes pour former une couche non consolidée de plusieurs dizaines de mètres d'épaisseur (granodiorite principale). Les injections individuelles sont délimitées par de subtils contrastes en granulométrie, proportions modales ou ségrégation de minéraux (schlieren), ou par des phénomènes d'érosion le long des surfaces de contact. Deux couches métriques à contour sinueux consistent en une accumulation compacte d'enclaves mafiques arrondies dans une matrice granodioritique fine. Le granoclassement des enclaves, la présence de figures de charge et de phénomènes érosifs en base de couche, ainsi que des schlieren de biotite entrecroisés évoquent l'injection de coulées de magma chargé d'enclaves et de faible viscosité en régime hydrodynamique turbulent dans un encaissant granodioritique encore largement liquide. La nature hybride des roches implique une chambre magmatique sous-jacente, en cours de différenciation et périodiquement réalimentée. Les magmas sont des liquides mafiques dérivés du manteau et des liquides anatectiques d'origine crustale, comme l'indique la gamme mesurée des rapports isotopiques initiaux du Sr (0.704 à 0.709) et des valeurs epsilon Nd (-2.1 à -4.7). Ces données montrent également que la contribution crustale est dominante, en accord avec les isotopes du plomb. Les phénomènes d'hybridation ont vraisemblablement eu lieu en base de croûte et dans la chambre magmatique sous-jacente au laccolite du Matorello. Les indicateurs de paléogravité du Matorello contribuent accessoirement à la compréhension de l'architecture actuelle de la nappe de Sambuco. Des plis isoclinaux à surface axiale verticale peuvent être mis en évidence par le contact entre les faciès dioritique et granodioritique. L'antiforme dont le Matorello forme le coeur est un synclinal, ce qui le positionne dans le Flanc inverse du grand pli couché que forme la nappe de Sambuco. Par ailleurs, des blocs de gneiss retrouvés dans le wildflysch sommital de la couverture de la nappe d'Antigorio ont été affiliés dans cette étude au pluton du Matorello. Ceci implique que le front de la nappe de Sambuco chevauchait déjà la partie est du bassin d'Antigorio au moment de sa fermeture. Par conséquent, ce n'est qu'en position externe que la nappe du Lebendun chevauche directement la nappe d'Antigorio. Abstract Magmatic bodies are important markers in paleo-geographic or geodynamic reconstructions of orogenic cycles, even more so in the case of polycyclic events where many of the other markers have been overwritten or destroyed. Plutons are relatively easy to date and their geochemical properties help constrain the tectonic context in which they were emplaced. This study focuses on the pre-mesozoic basement in the Sambuco and Maggia lower Penninic nappes located in the central Lepontine domain of the Alps. A number of magmatic events have been identified in the Sambuco basement. These events were dated using LA-ICPMS U/Pb on zircon grains. The mafic calc-alkaline banded Scheggia suite is dated as lower Cambrian, 540-530 Ma. The Al-rich Sasso-Nero lenticular gneiss is 480-470 Ma old (similarly to many older orfhogneisses of the Alpine basement) and contains 630-610 Ma old pan-African inherited zircons that illustrate the Gondwanian origin of these terranes.The calc-alkaline Matorello pluton is dated as 310-300 Ma whereas the lamprophyric bodies it contains are of 300 Ma. The Cocco granodiorite and the Ruscada leucogranite both intrude the basement of the adjacent Maggia nappe and are of similar ages to the Matorello. The ages as well as the geochemical similarities between the Cocco, Rucada and Matorello plutons suggest their paleo-geographic proximity at the Permian-Carboniferous boundary. However, these nappes are currently considered as belonging to two different Mesozoic paleo-geographic domains. Indeed, the Sambuco is considered as Helvetic whereas the Maggia is said to be Briançonnais, both separated by an oceanic basin. If this is the case, then it is essential that nostrike-slip movement has misaligned both continental margins since these coincide perfectly now that the oceanic domain closed. The Matorello pluton was originally a tabular intrusion, built up by the accumulation of multiple, several meter-thick, subhorizontal sheet-like injections of magma. Depending on their emplacement rate, the successive magma injections either solidified rapidly with sharp and rather well-defined boundaries (like the composite sills) or mingled with previous injections generating a thick molten layer up to several tens to hundred meters thick, like in the main granodioritic facies. These coalesced injections are hardly distinguishable, however subtle contrasts in granulometry, mineral modal proportions or mineral sorting (cross-bedded biotite-rich schlieren), as well as erosional features and/or crystal entrapment along contact surfaces allow to distinguish between the different injections. Two exceptional meter-thick layers display sinuous boundaries with the host granodiorite and consist of a densely packed accumulation of mafic enclaves in a granodioritic matrix. Gravitational sorting of the enclaves with load cast features at the base of the layers and sinuous biotite schlieren point to injection of low viscosity turbulent composite magma flows in the still largely molten granodiorite host. The hybrid nature of these rocks implies the existence of á periodically replenished and differentiated underlying magma chamber. Magmas are mafic liquids derived from the mantle and anatectic liquids of crustal origin, as shown by the (87Sr/86Sr), and epsilon Nd values (0.704-0.709 and -2.1 to -4.7 respectively. These data show that the crustal contribution is important, as confirmed by the Pb isotopes. The hybridisation processes seem to have occurred in the lower crust in magma chambers underlying the Matorello laccolith. The paleo-gravity markers in the Matorello help understand the architecture of the Sambuco nappe. Isoclinal folds with a vertical axial plane can be seen at the contact between dioritic and granodioritic facies. The antiform structure of which the Matorello is the heart is in fact a syncline. This places it in the inverse flanc of the large recumbent fold that constitutes the Sambuco nappe. The gneiss blocs found in the summital wildflysh cover of the Antigorio nappe have been linked to the Matorello pluton. This means that the front of the Sambuco nappe already overlapped the Antigorio basin when it closed. This implies that the Lebendun nappe can only overlap the Antigorio nappe in it's external position. Résumé grand public La chaîne alpine est la conséquence de la collision tertiaire entre deux masses continentales, l'Europe au nord et la péninsule apulienne africaine au sud, originellement séparées par l'océan mésozoïque téthysien. Cette collision a fermé un espace large de plusieurs centaines de km avec pour résultat l'écaillage de la croûte terrestre en unités tectoniques de dimensions variables, qui se sont empilées, imbriquées, éventuellement replissées en nappes de géométrie complexe. Cet amoncellement de 40 km d'épaisseur a vu sa température et sa pression lithostatique internes augmenter jusqu'à des valeurs de l'ordre de 680 °C et 6000 bars, induisant une recristallisation métamorphique des roches. L'un des objectifs de la géologie alpine est de reconstituer la géographie de la région aux temps mésozoïques de l'océan téthysien, en d'autres termes, de replacer chacune des unités tectoniques identifiées au sein de l'empilement alpin dans sa position originelle. Le défi est de taille et peut être comparé à celui de la reconstitution d'un vaste puzzle, dont certaines pièces seraient endommagées au niveau de leur contour ou leurs couleurs (métamorphisme), dissimulées par d'autres (enfouissement), voire tombées de la table de jeu (subduction, échappement latéral). Plusieurs approches ont été mises en oeuvre au cours du siècle écoulé. On citera en particulier la stratigraphie, la tectonique et le paléomagnétisme. Dans ce travail, nous avons essentiellement utilisé des techniques de datation isotopique absolue des roches (U/Pb sur zircon) qui, sur la base des connaissances acquises par l'ensemble des autres disciplines géologiques, nous ont permis de mieux contraindre ta paléogéographie mésozoïque du domaine «pennique inférieur » des Alpes centrales lépontines. Et au-delà? Nous savons tous que la disposition des continents à la surface de la Terre évolue constamment. Il est donc tentant d'essayer de remonter plus loin encore dans le temps et de reconstituer la physionomie de la marge sud européenne, tout au moins certains éléments de son histoire, au cours de l'ère paléozoïque. Les traces de ces événements très anciens sont naturellement ténues et dans ce contexte, les techniques de datation mentionnées ci-dessus deviennent les outils les plus performants. Ainsi, des datations u/Pb sur zircon nous ont permis de recenser plusieurs intrusions magmatiques, attribuées à quatre événements orogéniques anté-alpins. Des âges néoprotérozoïques (630-610 millions d'années ou Ma), cambrien inférieur (540-530 Ma), ordovicien inférieur (480-470 Ma) et carbonifère supérieur-permien inférieur (310-285 Ma) ont été obtenus dans le socle de la nappe de Sambuco. Des âges similaires à 300 Ma ont été obtenus dans la nappe voisine de la Maggia, qui permettent de relier ces deux unités. Aujourd'hui côte à côte, ces deux nappes devaient également se trouver proches l'une de l'autre il y a 300 Ma, lors de l'extension post-varisque. Les structures magmatiques spectaculaires préservées dans le pluton du Matorello (300 Ma) contraignent la géométrie actuelle de la nappe de Sambuco dans laquelle l'intrusion s'est mise en place. La forme originelle du pluton, aujourd'hui retourné et replissé plusieurs fois, s'avère être tabulaire, faite d'intrusions de faible épaisseur (1-300 m) s'étalant en forme de disque (30m à 2 km de diamètre). Les injections successives de magma se sont accumulées sous un toit dioritique précoce; elles sont issues, par le refais de fractures, d'une chambre magmatique plus profonde, périodiquement réalimentée par des magmas calco-alcalins d'origine mantellique contaminés parla croûte continentale profonde (εNd = -2.1 à -4.7). Des accumulations d'enclaves magmatiques arrondies et granoclassées dans des paléo-chenaux à fond érosif témoignent de conditions de mise en place hydrodynamiques à haute énergie. Ces enclaves sont emmenées de la chambre magmatique sous-jacente à la faveur d'épisodes de fracturation hydraulique liés à l'injection de magmas matelliques chauds dans des liquides différenciés riches en eau. Cette hypothèse est étayée par l'existence de filons composites. Une paléohorizontale a pu être déduite au sein du pluton, indiquant que cette partie de la nappe de Sambuco est verticalisée et isoclinalement replissée par la déformation alpine. Finalement, des blocs érodés du socle Sambuco ont été retrouvés dans le wildflysch sommital de la couverture sédimentaire mésozoïque de la nappe d'Antigorio sous-jacente. Ceci suggère que les blocs ont été fournis parle front de la nappe de Sambuco en train de chevaucher sur la nappe d'Antigorio au moment de la fermeture du bassin sédimentaire de cette dernière.

Three-dimensional hydrodynamic numerical modelling of fold nappe formation, basement-cover deformation and slab detachment with applications to the helvetic nappe system (W Switzerland)

Many three-dimensional (3-D) structures in rock, which formed during the deformation of the Earth's crust and lithosphere, are controlled by a difference in mechanical strength between rock units and are often the result of a geometrical instability. Such structures are, for example, folds, pinch-and-swell structures (due to necking) or cuspate-lobate structures (mullions). These struc-tures occur from the centimeter to the kilometer scale and the related deformation processes con-trol the formation of, for example, fold-and-thrust belts and extensional sedimentary basins or the deformation of the basement-cover interface. The 2-D deformation processes causing these structures are relatively well studied, however, several processes during large-strain 3-D defor-mation are still incompletely understood. One of these 3-D processes is the lateral propagation of these structures, such as fold and cusp propagation in a direction orthogonal to the shortening direction or neck propagation in direction orthogonal to the extension direction. Especially, we are interested in fold nappes which are recumbent folds with amplitudes usually exceeding 10 km and they have been presumably formed by ductile shearing. They often exhibit a constant sense of shearing and a non-linear increase of shear strain towards their overturned limb. The fold axes of the Morcles fold nappe in western Switzerland plunges to the ENE whereas the fold axes in the more eastern Doldenhorn nappe plunges to the WSW. These opposite plunge direc-tions characterize the Rawil depression (Wildstrubel depression). The Morcles nappe is mainly the result of layer parallel contraction and shearing. During the compression the massive lime-stones were more competent than the surrounding marls and shales, which led to the buckling characteristics of the Morcles nappe, especially in the north-dipping normal limb. The Dolden-horn nappe exhibits only a minor overturned fold limb. There are still no 3-D numerical studies which investigate the fundamental dynamics of the formation of the large-scale 3-D structure including the Morcles and Doldenhorn nappes and the related Rawil depression. We study the 3-D evolution of geometrical instabilities and fold nappe formation with numerical simulations based on the finite element method (FEM). Simulating geometrical instabilities caused by sharp variations of mechanical strength between rock units requires a numerical algorithm that can accurately resolve material interfaces for large differences in material properties (e.g. between limestone and shale) and for large deformations. Therefore, our FE algorithm combines a nu-merical contour-line technique and a deformable Lagrangian mesh with re-meshing. With this combined method it is possible to accurately follow the initial material contours with the FE mesh and to accurately resolve the geometrical instabilities. The algorithm can simulate 3-D de-formation for a visco-elastic rheology. The viscous rheology is described by a power-law flow law. The code is used to study the 3-D fold nappe formation, the lateral propagation of folding and also the lateral propagation of cusps due to initial half graben geometry. Thereby, the small initial geometrical perturbations for folding and necking are exactly followed by the FE mesh, whereas the initial large perturbation describing a half graben is defined by a contour line inter-secting the finite elements. Further, the 3-D algorithm is applied to 3-D viscous nacking during slab detachment. The results from various simulations are compared with 2-D resulats and a 1-D analytical solution. -- On retrouve beaucoup de structures en 3 dimensions (3-D) dans les roches qui ont pour origines une déformation de la lithosphère terrestre. Ces structures sont par exemple des plis, des boudins (pinch-and-swell) ou des mullions (cuspate-lobate) et sont présentés de l'échelle centimétrique à kilométrique. Mécaniquement, ces structures peuvent être expliquées par une différence de résistance entre les différentes unités de roches et sont généralement le fruit d'une instabilité géométrique. Ces différences mécaniques entre les unités contrôlent non seulement les types de structures rencontrées, mais également le type de déformation (thick skin, thin skin) et le style tectonique (bassin d'avant pays, chaîne d'avant pays). Les processus de la déformation en deux dimensions (2-D) formant ces structures sont relativement bien compris. Cependant, lorsque l'on ajoute la troisiéme dimension, plusieurs processus ne sont pas complètement compris lors de la déformation à large échelle. L'un de ces processus est la propagation latérale des structures, par exemple la propagation de plis ou de mullions dans la direction perpendiculaire à l'axe de com-pression, ou la propagation des zones d'amincissement des boudins perpendiculairement à la direction d'extension. Nous sommes particulièrement intéressés les nappes de plis qui sont des nappes de charriage en forme de plis couché d'une amplitude plurikilométrique et étant formées par cisaillement ductile. La plupart du temps, elles exposent un sens de cisaillement constant et une augmentation non linéaire de la déformation vers la base du flanc inverse. Un exemple connu de nappes de plis est le domaine Helvétique dans les Alpes de l'ouest. Une de ces nap-pes est la Nappe de Morcles dont l'axe de pli plonge E-NE tandis que de l'autre côté de la dépression du Rawil (ou dépression du Wildstrubel), la nappe du Doldenhorn (équivalent de la nappe de Morcles) possède un axe de pli plongeant O-SO. La forme particulière de ces nappes est due à l'alternance de couches calcaires mécaniquement résistantes et de couches mécanique-ment faibles constituées de schistes et de marnes. Ces différences mécaniques dans les couches permettent d'expliquer les plissements internes à la nappe, particulièrement dans le flanc inver-se de la nappe de Morcles. Il faut également noter que le développement du flanc inverse des nappes n'est pas le même des deux côtés de la dépression de Rawil. Ainsi la nappe de Morcles possède un important flanc inverse alors que la nappe du Doldenhorn en est presque dépour-vue. A l'heure actuelle, aucune étude numérique en 3-D n'a été menée afin de comprendre la dynamique fondamentale de la formation des nappes de Morcles et du Doldenhorn ainsi que la formation de la dépression de Rawil. Ce travail propose la première analyse de l'évolution 3-D des instabilités géométriques et de la formation des nappes de plis en utilisant des simulations numériques. Notre modèle est basé sur la méthode des éléments finis (FEM) qui permet de ré-soudre avec précision les interfaces entre deux matériaux ayant des propriétés mécaniques très différentes (par exemple entre les couches calcaires et les couches marneuses). De plus nous utilisons un maillage lagrangien déformable avec une fonction de re-meshing (production d'un nouveau maillage). Grâce à cette méthode combinée il nous est possible de suivre avec précisi-on les interfaces matérielles et de résoudre avec précision les instabilités géométriques lors de la déformation de matériaux visco-élastiques décrit par une rhéologie non linéaire (n>1). Nous uti-lisons cet algorithme afin de comprendre la formation des nappes de plis, la propagation latérale du plissement ainsi que la propagation latérale des structures de type mullions causé par une va-riation latérale de la géométrie (p.ex graben). De plus l'algorithme est utilisé pour comprendre la dynamique 3-D de l'amincissement visqueux et de la rupture de la plaque descendante en zone de subduction. Les résultats obtenus sont comparés à des modèles 2-D et à la solution analytique 1-D. -- Viele drei dimensionale (3-D) Strukturen, die in Gesteinen vorkommen und durch die Verfor-mung der Erdkruste und Litosphäre entstanden sind werden von den unterschiedlichen mechani-schen Eigenschaften der Gesteinseinheiten kontrolliert und sind häufig das Resulat von geome-trischen Istabilitäten. Zu diesen strukturen zählen zum Beispiel Falten, Pich-and-swell Struktu-ren oder sogenannte Cusbate-Lobate Strukturen (auch Mullions). Diese Strukturen kommen in verschiedenen Grössenordungen vor und können Masse von einigen Zentimeter bis zu einigen Kilometer aufweisen. Die mit der Entstehung dieser Strukturen verbundenen Prozesse kontrol-lieren die Entstehung von Gerbirgen und Sediment-Becken sowie die Verformung des Kontaktes zwischen Grundgebirge und Stedimenten. Die zwei dimensionalen (2-D) Verformungs-Prozesse die zu den genannten Strukturen führen sind bereits sehr gut untersucht. Einige Prozesse wäh-rend starker 3-D Verformung sind hingegen noch unvollständig verstanden. Einer dieser 3-D Prozesse ist die seitliche Fortpflanzung der beschriebenen Strukturen, so wie die seitliche Fort-pflanzung von Falten und Cusbate-Lobate Strukturen senkrecht zur Verkürzungsrichtung und die seitliche Fortpflanzung von Pinch-and-Swell Strukturen othogonal zur Streckungsrichtung. Insbesondere interessieren wir uns für Faltendecken, liegende Falten mit Amplituden von mehr als 10 km. Faltendecken entstehen vermutlich durch duktile Verscherung. Sie zeigen oft einen konstanten Scherungssinn und eine nicht-lineare zunahme der Scherverformung am überkipp-ten Schenkel. Die Faltenachsen der Morcles Decke in der Westschweiz fallen Richtung ONO während die Faltenachsen der östicher gelegenen Doldenhorn Decke gegen WSW einfallen. Diese entgegengesetzten Einfallrichtungen charakterisieren die Rawil Depression (Wildstrubel Depression). Die Morcles Decke ist überwiegend das Resultat von Verkürzung und Scherung parallel zu den Sedimentlagen. Während der Verkürzung verhielt sich der massive Kalkstein kompetenter als der Umliegende Mergel und Schiefer, was zur Verfaltetung Morcles Decke führ-te, vorallem in gegen Norden eifallenden überkippten Schenkel. Die Doldenhorn Decke weist dagegen einen viel kleineren überkippten Schenkel und eine stärkere Lokalisierung der Verfor-mung auf. Bis heute gibt es keine 3-D numerischen Studien, die die fundamentale Dynamik der Entstehung von grossen stark verformten 3-D Strukturen wie den Morcles und Doldenhorn Decken sowie der damit verbudenen Rawil Depression untersuchen. Wir betrachten die 3-D Ent-wicklung von geometrischen Instabilitäten sowie die Entstehung fon Faltendecken mit Hilfe von numerischen Simulationen basiert auf der Finite Elemente Methode (FEM). Die Simulation von geometrischen Instabilitäten, die aufgrund von Änderungen der Materialeigenschaften zwischen verschiedenen Gesteinseinheiten entstehen, erfortert einen numerischen Algorithmus, der in der Lage ist die Materialgrenzen mit starkem Kontrast der Materialeigenschaften (zum Beispiel zwi-schen Kalksteineinheiten und Mergel) für starke Verfomung genau aufzulösen. Um dem gerecht zu werden kombiniert unser FE Algorithmus eine numerische Contour-Linien-Technik und ein deformierbares Lagranges Netz mit Re-meshing. Mit dieser kombinierten Methode ist es mög-lich den anfänglichen Materialgrenzen mit dem FE Netz genau zu folgen und die geometrischen Instabilitäten genügend aufzulösen. Der Algorithmus ist in der Lage visko-elastische 3-D Ver-formung zu rechnen, wobei die viskose Rheologie mit Hilfe eines power-law Fliessgesetzes beschrieben wird. Mit dem numerischen Algorithmus untersuchen wir die Entstehung von 3-D Faltendecken, die seitliche Fortpflanzung der Faltung sowie der Cusbate-Lobate Strukturen die sich durch die Verkürzung eines mit Sediment gefüllten Halbgraben bilden. Dabei werden die anfänglichen geometrischen Instabilitäten der Faltung exakt mit dem FE Netz aufgelöst wäh-rend die Materialgranzen des Halbgrabens die Finiten Elemente durchschneidet. Desweiteren wird der 3-D Algorithmus auf die Einschnürung während der 3-D viskosen Plattenablösung und Subduktion angewandt. Die 3-D Resultate werden mit 2-D Ergebnissen und einer 1-D analyti-schen Lösung verglichen.