Tectonic and metamorphic evolution of the Central Himalayan Domain in Southeast Zanskar (Kashmir, India)

La région du Zanskar, étudiée dans le cadre de ce travail, se situe au passage entre deux domaines himalayens fortement contrastés, la Séquence Cristalline du Haut Himalaya (HHCS), composée de roches métamorphiques et l'Himalaya Tethysien (TH), composé de séries sédimentaires. La transition entre ces deux domaines est marquée par une structure tectonique majeure, la Zone de Cisaillement du Zanskar (ZSZ), au sein de laquelle on observe une augmentation extrêmement rapide, mais néanmoins graduelle, du degré du métamorphisme entre le TH et le HHCS. Il a été établi que le HHCS n'est autre que l'équivalent métamorphique des séries sédimentaires de la base du TH. C'est principalement lors d'un épisode de mise en place de nappes à vergence sudouest, entre l'Eocène moyen et l'Oligocène, que les séries sédimentaires de la base du TH ont été entraînées en profondeur où elles ont subi un métamorphisme de type barrovien. Au début du Miocène, le HHCS à été exhumé en direction du sud-ouest sous forme d'une grande nappe, délimitée a sa base par le MCT (principal chevauchement central) et à son sommet par la Zone de Cisaillement du Zanskar. L'ensemble des zones barroviennes, de la zone à biotite jusqu'à la zone à disthène, a été cisaillée par les mouvements en faille normale au sommet du HHCS et se retrouve actuellement sur une épaisseur d'environ 1 kilomètre au sein de la ZSZ. La décompression associée à l'exhumation du HHCS a provoqué la fusion partielle d'une partie du HHCS et a donné naissance à des magmas de composition leucogranitiques. Grâce à la géothermobarometrie, et connaissant la géométrie de la ZSZ, il nous a été possible de déterminer que le rejet le long de cette structure d'extension est d'au moins 35?9 kilomètres. Une série d'arguments nous permet cependant de suggérer que ce rejet aurait pu être encore bien plus important (~100km). Les données géochronologiques nous permettent de contraindre la durée des mouvements d'extension le long de la ZSZ à 2.4?0.2 Ma entre 22.2?0.2 Ma et 19.8?0.1 Ma. Ce travail apporte de nouvelles données sur les processus métamorphiques, magmatiques et tectoniques liés aux phénomènes d'extension syn-orogeniques.<br/><br/>The southeastern part of Zanskar is located at the transition between two major Himalayan domains of contrasting metamorphic grade, the High Himalayan Crystalline Sequence (HHCS) and the Tethyan Himalaya (TH). The transition between the TH and the HHCS is marked by a very rapid, although perfectly gradual, decrease in metamorphic grade, which coincides with a major tectonic structure, the Zanskar Shear Zone (ZSZ). It is now an established fact that the relation between the HHCS and the TH is not one of basement-cover type, but that the metasedimentary series of the HHCS represent the metamorphic equivalent of the lowermost sedimentary series of the TH. This transformation of sedimentary series into metamorphic rocks, and hence the differentiation between the TH and the HHCS, is the consequence of crustal thickening associated to the formation of large scale southwest vergent nappes within the Tethyan Himalaya sedimentary series. This, Middle Eocene to Oligocene, episode of crustal thickening and associated Barrovian metamorphism is followed, shortly after, by the exhumation of the HHCS as a, large scale, south-west vergent, nappe. Foreword The exhumation of the HHCS nappe is marked by the activation of two contemporaneous structures, the Main Central Thrust at its base and the Zanskar Shear Zone at its top. Extensional movements along the ZSZ, caused the Barrovian biotite to the kyanite zones to be sheared and constricted within the ~1 km thick shear zone. Decompression associated with the exhumation of the HHCS induced the formation of leucogranitic magmas through vapour-absent partial melting of the highest-grade rocks. The combination of geothermobarometric data with a geometric model of the ZSZ allowed us to constrain the net slip at the top of the HHCS to be at least 35?9 kilometres. A set of arguments however suggests that these movements might have been much more important (~ 100 km). Geochronological data coupled with structural observations constrain the duration of ductile shearing along the ZSZ to 2.4?0.2 Ma between 22.2?0.2 Ma and 19.8?0.1 Ma. This study also addresses the consequences of synorogenic extension on the metamorphic, tectonic and magmatic evolution of the upper parts of the High Himalayan Crystalline Sequence.

Le pergélisol dans les terrains sédimentaires à forte déclivité: distribution, régime thermique et instabilités

Résumé: Dans le contexte d'un climat de plus en plus chaud, la localisation du pergélisol dans les terrains sédimentaires à forte déclivité et l'évaluation des mouvements de terrain qui y ont cours s'avèrent primordiales. S'insérant dans cette problématique, ce travail de thèse s'articule autour de deux axes de recherche différents. D'un point de vue statique, cette recherche propose une étude de la distribution et des caractéristiques du pergélisol dans les éboulis de la zone périglaciaire alpine. D'un point de vue dynamique, une analyse de l'influence des caractéristiques du pergélisol (teneur en glace, température du pergélisol, etc.) et des variations des températures de l'air et du sol sur les vitesses de fluage des corps sédimentaires gelés est effectuée. Afin de répondre à ce double objectif, l'approche "terrain" a été privilégiée. Pour déterminer la répartition et les caractéristiques du pergélisol, les méthodes traditionnelles de prospection du pergélisol ont été utilisées, à savoir la mesure de la température du sol à la base du manteau neigeux (BTS), la mesure de la température du sol en continu ainsi que la méthode géoélectrique. Les mouvements de terrain ont pour leur part été mesurés à l'aide d'un GPS différentiel. L'étude de la distribution du pergélisol a été effectuée dans une quinzaine d'éboulis situés dans les régions du Mont Gelé (Verbier-Nendaz) et d'Arolla principalement. Dans la plupart des cas, un pergélisol a pu être mis en évidence dans la partie inférieure des accumulations sédimentaires, alors que la partie médiane des éboulis n'est, le plus souvent, pas gelée. Si cette absence de pergélisol se prolonge parfois dans les portions sommitales des pentes, les mesures réalisées montrent que dans d'autres cas des sédiments gelés y sont à nouveau présents. Les résistivités électriques mesurées dans les portions gelées des éboulis étudiés sont dans la plupart des cas nettement inférieures à celles mesurées sur les glaciers rocheux. Des études préalables ont montré que des circulations d'air internes sont responsables de l'anomalie thermique négative et, lorsqu'il existe, du pergélisol que l'on trouve dans la partie inférieure d'éboulis situés plus de 1000 m plus bas que la limite inférieure régionale du pergélisol discontinu. L'étude de quatre sites de basse altitude (1400-1900 m), et notamment l'équipement du site de Dreveneuse (Préalpes Valaisannes) avec deux forages, des capteurs de température de surface et un anémomètre a permis de vérifier et de préciser le mécanisme de ventilation actif au sein des éboulis froids de basse altitude. Ce mécanisme fonctionne de la manière suivante: en hiver, l'air contenu dans l'éboulis, plus chaud et plus léger que l'air extérieur, monte à l'intérieur de l'accumulation sédimentaire et est expulsé dans ses parties sommitales. Cet effet de cheminée provoque une aspiration d'air froid à l'intérieur de la partie inférieure de l'éboulis, causant ainsi un sur-refroidissement marqué du terrain. En été, le mécanisme s'inverse, l'éboulis étant plus froid que l'air environnant. De l'air froid est alors expulsé au bas de la pente. Une ventilation ascendante hivernale a pu être mise en évidence dans certains des éboulis de haute altitude étudiés. Elle est probablement en grande partie responsable de la configuration particulière des zones gelées observées. Même si l'existence d'un effet de cheminée n'a pu être démontrée dans tous les cas, du fait notamment de la glace interstitielle qui entrave le cheminement de l'air, des indices laissant présager son possible fonctionnement existent dans la quasi totalité des éboulis étudiés. L'absence de pergélisol à des altitudes qui lui sont favorables pourrait en tous les cas s'expliquer par un réchauffement du terrain lié à des expulsions d'air relativement chaud. L'étude des mouvements de terrain a été effectuée sur une dizaine de sites, principalement sur des glaciers rocheux, mais également sur une moraine de poussée et - II - Résumé ? abstract quelques éboulis. Plusieurs glaciers rocheux présentent des formes de déstabilisation récente (niches d'arrachement, blocs basculés, apparition de la matrice fine à la surface, etc.), ce qui témoigne d'une récente accélération des vitesses de déplacement. Ce phénomène, qui semble général à l'échelle alpine, est probablement à mettre sur le compte du réchauffement du pergélisol depuis une vingtaine d'années. Les vitesses mesurées sur ces formations sont souvent plus élevées que les valeurs habituellement proposées dans la littérature. On note par ailleurs une forte variabilité inter-annuelle des vitesses, qui semblent dépendre de la variation de la température moyenne annuelle de surface. Abstract: In the context of a warmer climate, the localisation of permafrost in steep sedimentary terrain and the measurement of terrain movements that occur in these areas is of great importance. With respect to these problems, this PhD thesis follows two different research axes. From a static point of view, the research presents a study of the permafrost distribution and characteristics in the talus slopes of the alpine periglacial belt. From a dynamic point of view, an analysis of the influence of the permafrost characteristics (ice content, permafrost temperature, etc.) and air and soil temperature variations on the creep velocities of frozen sedimentary bodies is carried out. In order to attain this double objective, the "field" approach was favoured. To determine the distribution and the characteristics of permafrost, the traditional methods of permafrost prospecting were used, i.e. ground surface temperature measurements at the base of the snow cover (BTS), year-round ground temperature measurements and DC-resistivity prospecting. The terrain movements were measured using a differential GPS. The permafrost distribution study was carried out on 15 talus slopes located mainly in the Mont Gelé (Verbier-Nendaz) and Arolla areas (Swiss Alps). In most cases, permafrost was found in the lower part of the talus slope, whereas the medium part was free of ice. In some cases, the upper part of the talus is also free of permafrost, whereas in other cases permafrost is present. Electrical resistivities measured in the frozen parts of the studied talus are in most cases clearly lower than those measured on rock glaciers. Former studies have shown that internal air circulation is responsible for the negative thermal anomaly and, when it exists, the permafrost present in the lower part of talus slopes located more than 1000 m below the regional lower limit of discontinuous permafrost. The study of four low-altitude talus slopes (1400-1900 m), and notably the equipment of Dreveneuse field site (Valais Prealps) with two boreholes, surface temperature sensors and an anemometer permitted to verify and to detail the ventilation mechanism active in low altitude talus slopes. This mechanism works in the following way: in winter, the air contained in the block accumulation is warmer and lighter than the surrounding air and therefore moves upward in the talus and is expelled in its upper part. This chimney effect induces an aspiration of cold air in the interior of the lower part of talus, that causes a strong overcooling of the ground. In summer, the mechanism is reversed because the talus slope is colder than the surrounding air. Cold air is then expelled in the lower part of the slope. Evidence of ascending ventilation in wintertime could also be found in some of the studied high-altitude talus slopes. It is probably mainly responsible for the particular configuration of the observed frozen areas. Even if the existence of a chimney effect could not be demonstrated in all cases, notably because of interstitial ice that obstructs Résumé ? abstract - III - the air circulation, indices of its presence exist in nearly all the studied talus. The absence of permafrost at altitudes favourable to its presence could be explained, for example, by the terrain warming caused by expulsion of relatively warm air. Terrain movements were measured at about ten sites, mainly on rock glaciers, but also on a push moraine and some talus slopes. Field observations reveal that many rock glaciers display recent destabilization features (landslide scars, tilted blocks, presence of fine grained sediments at the surface, etc.) that indicate a probable recent acceleration of the creep velocities. This phenomenon, which seems to be widespread at the alpine scale, is probably linked to the permafrost warming during the last decades. The measured velocities are often higher than values usually proposed in the literature. In addition, strong inter-annual variations of the velocities were observed, which seems to depend on the mean annual ground temperature variations.

Scales and causes of heterogeneity in bars in a large multi-channel river: Rio Parana, Argentina

To date, published studies of alluvial bar architecture in large rivers have been restricted mostly to case studies of individual bars and single locations. Relatively little is known about how the depositional processes and sedimentary architecture of kilometre-scale bars vary within a multi-kilometre reach or over several hundreds of kilometres downstream. This study presents Ground Penetrating Radar and core data from 11, kilometre-scale bars from the Rio Parana, Argentina. The investigated bars are located between 30km upstream and 540km downstream of the Rio Parana - Rio Paraguay confluence, where a significant volume of fine-grained suspended sediment is introduced into the network. Bar-scale cross-stratified sets, with lengths and widths up to 600m and thicknesses up to 12m, enable the distinction of large river deposits from stacked deposits of smaller rivers, but are only present in half the surface area of the bars. Up to 90% of bar-scale sets are found on top of finer-grained ripple-laminated bar-trough deposits. Bar-scale sets make up as much as 58% of the volume of the deposits in small, incipient mid-channel bars, but this proportion decreases significantly with increasing age and size of the bars. Contrary to what might be expected, a significant proportion of the sedimentary structures found in the Rio Parana is similar in scale to those found in much smaller rivers. In other words, large river deposits are not always characterized by big structures that allow a simple interpretation of river scale. However, the large scale of the depositional units in big rivers causes small-scale structures, such as ripple sets, to be grouped into thicker cosets, which indicate river scale even when no obvious large-scale sets are present. The results also show that the composition of bars differs between the studied reaches upstream and downstream of the confluence with the Rio Paraguay. Relative to other controls on downstream fining, the tributary input of fine-grained suspended material from the Rio Paraguay causes a marked change in the composition of the bar deposits. Compared to the upstream reaches, the sedimentary architecture of the downstream reaches in the top ca 5m of mid-channel bars shows: (i) an increase in the abundance and thickness (up to metre-scale) of laterally extensive (hundreds of metres) fine-grained layers; (ii) an increase in the percentage of deposits comprised of ripple sets (to >40% in the upper bar deposits); and (iii) an increase in bar-trough deposits and a corresponding decrease in bar-scale cross-strata (<10%). The thalweg deposits of the Rio Parana are composed of dune sets, even directly downstream from the Rio Paraguay where the upper channel deposits are dominantly fine-grained. Thus, the change in sedimentary facies due to a tributary point-source of fine-grained sediment is primarily expressed in the composition of the upper bar deposits.

Conditioning of multiple-point statistics facies simulations to tomographic images

Geophysical tomography captures the spatial distribution of the underlying geophysical property at a relatively high resolution, but the tomographic images tend to be blurred representations of reality and generally fail to reproduce sharp interfaces. Such models may cause significant bias when taken as a basis for predictive flow and transport modeling and are unsuitable for uncertainty assessment. We present a methodology in which tomograms are used to condition multiple-point statistics (MPS) simulations. A large set of geologically reasonable facies realizations and their corresponding synthetically calculated cross-hole radar tomograms are used as a training image. The training image is scanned with a direct sampling algorithm for patterns in the conditioning tomogram, while accounting for the spatially varying resolution of the tomograms. In a post-processing step, only those conditional simulations that predicted the radar traveltimes within the expected data error levels are accepted. The methodology is demonstrated on a two-facies example featuring channels and an aquifer analog of alluvial sedimentary structures with five facies. For both cases, MPS simulations exhibit the sharp interfaces and the geological patterns found in the training image. Compared to unconditioned MPS simulations, the uncertainty in transport predictions is markedly decreased for simulations conditioned to tomograms. As an improvement to other approaches relying on classical smoothness-constrained geophysical tomography, the proposed method allows for: (1) reproduction of sharp interfaces, (2) incorporation of realistic geological constraints and (3) generation of multiple realizations that enables uncertainty assessment.

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.

Evolution of the Helvetic shelf (Switzerland) during the Barremian-early Aptian : paleoenvironmental, paleogeographic and paleoceanographic controlling factors

The Early Cretaceous has experienced the development of large shallow-water carbonate platform in tropical and subtropical regions, favoured by exceptionally warm climatic conditions, optimal trophic conditions and a suitable tectonic and paleogeographic context. This period was also characterized by shorter intervals, in which the widespread deposition of marine sediments enriched in organic matter occurred ("oceanic anoxic episodes": OAE). This study focuses on the Barremian- Aptian interval, during which the Urgonian platform developed throughout the northern Tethyan passive margin. Due to the Alpine orogeny, sediments belonging to this platform - named locally Schrattenkalk Formation, are presently outcropping in the Helvetic Alps. This study aims to reconstruct the paleogeographic evolution of the Helvetic platform, and to define the environmental and oceanographic factors, which influenced its development. Several key episodes in the life of this platform have been identified: - The installation of the platform, covering hemipelagic sediments of the Drusberg Member, near the limit between the early and late Barremian. - The temporary change of carbonate production type during the basal Aptian, with the deposition of the Rawil Member. - And finally the definitive interruption of photozoan carbonate platform sedimentation in the study area, during the early Aptian. The sedimentological, biostratigraphical and chemostratigraphic (8I3C) data lead to the sequential subdivision of eleven sections and one core, located throughout the different Helvetic nappes of Switzerland. The sequence stratigraphie framework, initially defined for the Urgonian carbonate platform of the Vercors area (SE France), is confirmed in the Helvetic nappes, where the same number of sequences was observed. Many similarities between these two areas are put forward in this work. The sequence stratigraphie framework helped to highlight the installation of a bioclastic body, included in the Schrattenkalk Formation, since the middle Early Barremian (sequence B2). The age of the installation of the rudist-rich limestone, which corresponds to the Urgonian facies sensu stricto, is attributed to the late Barremian (maximum flooding surface of the sequence B3). This age coincides with the one determined in other northern Tethyan areas for the installation of the Urgonian platform. The results of this study show a strong tectonic control of the platform architecture, with the presence of syn-sedimentary faults in a perpendicular position to the progradation direction of the platform. The presence of these faults was highlighted by the study of the evolution of the microfacies distribution and by thickness variations in different areas. Sea level fluctuations also played an important role in the various life phases of the platform. Three major falls in sea level have been identified. A significant emersion of the proximal domain has been observed, involving an important drop of the relative sea level, leading to the exposure of the Drusberg Member hemipelagic series. A second major drop in sea level is identified near the Barremian-Aptian boundary, and a third is registered on the top of the Upper Schrattenkalk Member on the whole platform; it is associated with a karst affecting the underlying limestones to a depth of over 20 meters. This observation sheds new light on the conditions linked to the demise of Urgonian platform, which was strongly influenced by this phase of emersion.