Analog Modelling of N-S anticlines of the External Sierras (Jaca basin, Southern Pyrenees)

Report for the scientific sojourn carried out at the Uppsala Universitet, Sweden, from April to July the 2007. Two series of analogue models are used to explore ductile-frictional contrasts of the basal décollement in the development of oblique and transverse structures simultaneously to thin-skinned shortening. These models simulate the evolution of the Central External Sierras (Southern Pyrenees, Spain), which constitute the frontal emerging part of the southernmost Pyrenean thrust sheet. They are characterized by the presence of transverse N-S to NW-SE anticlines, which are perpendicular to the Pyrenean structural trend and developed in the hangingwall of the Santo Domingo thrust system, detaching on an unevenly distributed Triassic materials (evaporitic-dolomitic interfingerings). Model setup performs a décollement made by three patches of silicone neighbouring pure brittle sand. Model series A test the thickness ratio between overburden and décollement. Model series B test the width of frictional detachment areas. Model results show how deformation reaches further in areas detached on ductile layer whereas frictional décollement areas assimilate the strain by means of an additional uplift. This replicates the structural style of Central External Sierras: higher structural relief of N-S anticlines with regard to orogen-parallel structures, absence of a representative ductile décollement in the core, tilting towards the orogen and foreland-side closure not thrusted by the frontal emerging South-Pyrenean thrust.

Complex landslide behaviour and structural controlled by the structures: A 3D conceptual model example of Åknes rockslide, Norway

Åknes is an active complex large rockslide of approximately 30?40 Mm3 located within the Proterozoic gneisses of western Norway. The observed surface displacements indicate that this rockslide is divided into several blocks moving in different directions at velocities of between 3 and 10 cm year?1. Because of regional safety issues and economic interests this rockslide has been extensively monitored since 2004. The understanding of the deformation mechanism is crucial for the implementation of a viable monitoring system. Detailed field investigations and the analysis of a digital elevation model (DEM) indicate that the movements and the block geometry are controlled by the main schistosity (S1) in gneisses, folds, joints and regional faults. Such complex slope deformations use pre-existing structures, but also result in new failure surfaces and deformation zones, like preferential rupture in fold-hinge zones. Our interpretation provides a consistent conceptual three-dimensional (3D) model for the movements measured by various methods that is crucial for numerical stability modelling. In addition, this reinterpretation of the morphology confirms that in the past several rockslides occurred from the Åknes slope. They may be related to scars propagating along the vertical foliation in folds hinges. Finally, a model of the evolution of the Åknes slope is presented.

Western Préalpes Médianes Romandes: Timing and structure. A review.

Between the original position and their present day location as klippen, the Prealpes Medianes underwent a complex history of paleotectonics and alpine tectonics. Due to the opening of the Piemont ocean the Brianconnais sedimentation realm of the Prealpes Medianes evolved as a rim basin of the northern passive margin during Jurassic to Eocene times. Different paleotectonic features (normal faults, synsedimentary growth structures, inversion structures) developed and were active above a basal detachment in evaporitic layers. The tectonic movements were a consequence of thermal events in the crust. Isolated from the Iberic continent at the end of the Late Cretaceous, the Brianconnais exotic terrain was incorporated into the accretionary prism of the closing Piemont ocean and the incipient alpine orogeny during the Lutetian-Bartonian. The Prealpes Medianes were detached from their homeland during the Bartonian-Priabonian and were transported onto the foreland. The tectonic style is one of a thin-skinned foreland fold and thrust belt. Fault associated fold development above a main decollement, together with internal deformation, represent the Prealpes Medianes main structural features. The very low-grade metamorphic conditions have their origin in the heat flux induced by tectonic burial by overriding nappes in the accretionary prism. After having been transported on top of the developing Helvetic nappes the Prealpes were emplaced in their present day position in front of the Alpine mountain belt during Oligocene times. Post-emplacement and out of sequence thrusting, possibly younger than Oligocene, is observed and can be related to thrusting in the sedimentary substratum and the basement.

Homotopy Batalin-Vilkovisky algebras

This paper provides an explicit cofibrant resolution of the operad encoding Batalin-Vilkovisky algebras. Thus it defines the notion of homotopy Batalin-Vilkovisky algebras with the required homotopy properties. To define this resolution we extend the theory of Koszul duality to operads and properads that are defined by quadratic and linear relations. The operad encoding Batalin-Vilkovisky algebras is shown to be Koszul in this sense. This allows us to prove a Poincaré-Birkhoff-Witt Theorem for such an operad and to give an explicit small quasi-free resolution for it. This particular resolution enables us to describe the deformation theory and homotopy theory of BV-algebras and of homotopy BV-algebras. We show that any topological conformal field theory carries a homotopy BV-algebra structure which lifts the BV-algebra structure on homology. The same result is proved for the singular chain complex of the double loop space of a topological space endowed with an action of the circle. We also prove the cyclic Deligne conjecture with this cofibrant resolution of the operad BV. We develop the general obstruction theory for algebras over the Koszul resolution of a properad and apply it to extend a conjecture of Lian-Zuckerman, showing that certain vertex algebras have an explicit homotopy BV-algebra structure.

Segmentation of brain structures in presence of a space-occupying lesion.

Brain deformations induced by space-occupying lesions may result in unpredictable position and shape of functionally important brain structures. The aim of this study is to propose a method for segmentation of brain structures by deformation of a segmented brain atlas in presence of a space-occupying lesion. Our approach is based on an a priori model of lesion growth (MLG) that assumes radial expansion from a seeding point and involves three steps: first, an affine registration bringing the atlas and the patient into global correspondence; then, the seeding of a synthetic tumor into the brain atlas providing a template for the lesion; finally, the deformation of the seeded atlas, combining a method derived from optical flow principles and a model of lesion growth. The method was applied on two meningiomas inducing a pure displacement of the underlying brain structures, and segmentation accuracy of ventricles and basal ganglia was assessed. Results show that the segmented structures were consistent with the patient's anatomy and that the deformation accuracy of surrounding brain structures was highly dependent on the accurate placement of the tumor seeding point. Further improvements of the method will optimize the segmentation accuracy. Visualization of brain structures provides useful information for therapeutic consideration of space-occupying lesions, including surgical, radiosurgical, and radiotherapeutic planning, in order to increase treatment efficiency and prevent neurological damage.

40Ar/39Ar dating of Penninic Front tectonic displacement (W Alps) during the Lower Oligocene (31-34 Ma)

Direct absolute dating of the Penninic Frontal Thrust tectonic motion is achieved using the Ar-40/Ar-39 technique in the Pelvoux Crystalline Massif (Western Alps). The dated phengites were formed syn-kinematically in shear zones. They underline the brittle-ductile stretching lineation, pressure-shadow fibres and slickensides consistent with underthrusting of the European continental slab below the propagating Penninic Thrust. Chlorite-phengite thermobarometry yields 10-15 km and T similar to 280 degrees C, while Ar-40/Ar-39 phengite ages mainly range between 34 and 30 Ma, with one younger age at 27 Ma. This Early Oligocene age range matches a major tectonic rearrangement of the Alpine chain. Preservation of prograde Ar-40/Ar-39 ages is ascribed to passive exhumation of the Pelvoux shear zone network, sandwiched between more external thrusts and the Penninic Front reactivated as an E-dipping detachment fault. Partial resetting in the Low Temperature part of argon spectra below 24 Ma is ascribed to brittle deformation and alteration of phengites.

Laser scanning-based recognition of rotational movements on a deep seated gravitational instability: The Cinque Torri case (North-Eastern Italian Alps)

The Cinque Torri group (Cortina d'Ampezzo, Italy) is an articulated system of unstable carbonatic rock monoliths located in a very important tourism area and therefore characterized by a significant risk. The instability phenomena involved represent an example of lateral spreading developed over a larger deep seated gravitational slope deformation (DSGSD) area. After the recent fall of a monolith of more than 10 000 m3, a scientific study was initiated to monitor the more unstable sectors and to characterize the past movements as a fundamental tool for predicting future movements and hazard assessment. To achieve greater insight on the ongoing lateral spreading process, a method for a quantitative analysis of rotational movements associated with the lateral spreading has been developed, applied and validated. The method is based on: i) detailed geometrical characterization of the area by means of laser scanner techniques; ii) recognition of the discontinuity sets and definition of a reference frame for each set, iii) correlation between the obtained reference frames related to a specific sector and a stable external reference frame, and iv) determination of the 3D rotations in terms of Euler angles to describe the present settlement of the Cinque Torri system with respect to the surrounding stable areas. In this way, significant information on the processes involved in the fragmentation and spreading of a former dolomitic plateau into different rock cliffs has been gained. The method is suitable to be applied to similar case studies.

Stable isotope geochemistry and formation mechanisms of quartz veins; Extreme paleoaltitudes of the Central Alps in the Neogene

Quartz veins ranging in size from less than 50 cm length and 5 cm width to greater than 10 m in length and 5 m in width are found throughout the Central Swiss Alps. In some cases, the veins are completely filled with milky quartz, while in others, sometimes spectacular void-filling quartz crystals are found. The style of vein filling and size is controlled by host rock composition and deformation history. Temperatures of vein formation, estimated using stable isotope thermometry and mineral equilibria, cover a range of 450 degrees C down to 150 degrees C. Vein formation started at 18 to 20 Ma and continued for over 10 My. The oxygen isotope values of quartz veins range from 10 to 20 permil, and in almost all cases are equal to those of the hosting lithology. The strongly rock-buffered veins imply a low fluid/rock ratio and minimal fluid flow. In order to explain massive, nearly morromineralic quartz formation without exceptionally large fluid fluxes, a mechanism of differential pressure and silica diffusion, combined with pressure solution, is proposed for early vein formation. Fluid inclusions and hydrous minerals in late-formed veins have extremely low delta D values, consistent with meteoric water infiltration. The change from rock-buffered, static fluid to infiltration from above can be explained in terms of changes in the large-scale deformation style occurring between 20 and 15 Ma. The rapid cooling of the Central Alps identified in previous studies may be explained in part, by infiltration of cold meteoric waters along fracture systems down to depths of 10 km or more. An average water flux of 0.15 cm 3 cm(-2)yr(-1) entering the rock and reemerging heated by 40 degrees C is sufficient to cool rock at 10 km depth by 100 degrees C in 5 million years. The very negative delta D values of < -130 permil for the late stage fluids are well below the annual average values measured in meteoric water in the region today. The low fossil delta D values indicate that the Central Alps were at a higher elevation in the Neogene. Such a conclusion is supported by an earlier work, where a paleoaltitude of 5000 meters was proposed on the basis of large erratic boulders found at low elevations far from their origin.

Sperm Ultrastructure of the Spider Crab Maja brachydactyla (Decapoda: Brachyura)

This study describes the morphology of the sperm cell of Maja brachydactyla, with emphasis on localizing actin and tubulin. The spermatozoon of M. brachydactyla is similar in appearance and organization to other brachyuran spermatozoa. The spermatozoon is a globular cell composed of a central acrosome, which is surrounded by a thin layer of cytoplasm and a cup-shaped nucleus with four radiating lateral arms. The acrosome is a subspheroidal vesicle composed of three concentric zones surrounded by a capsule. The acrosome is apically covered by an operculum. The perforatorium penetrates the center of the acrosome and has granular material partially composed of actin. The cytoplasm contains one centriole in the subacrosomal region. A cytoplasmic ring encircles the acrosome in the subapical region of the cell and contains the structures-organelles complex (SO-complex), which is composed of a membrane system, mitochondria with few cristae, and microtubules. In the nucleus, slightly condensed chromatin extends along the lateral arms, in which no microtubules have been observed. Chromatin fibers aggregate in certain areas and are often associated with the SO-complex. During the acrosomal reaction, the acrosome could provide support for the penetration of the sperm nucleus, the SO-complex could serve as an anchor point for chromatin, and the lateral arms could play an important role triggering the acrosomal reaction, while slightly decondensed chromatin may be necessary for the deformation of the nucleus.

Oligo-Miocene extensional tectonics and fluid flow across the Northern Snake Range detachment system, Nevada

The Northern Snake Range (Nevada) represents a spectacular example of a metamorphic core complex and exposes a complete section from the mylonitic footwall into the hanging wall of a fossil detachment system. Paired geochronological and stable isotopic data of mylonitic quartzite within the detachment footwall reveal that ductile deformation and infiltration of meteoric fluids occurred between 27 and 23 Ma. Ar-40/Ar-39 ages display complex recrystallization-cooling relationships but decrease systematically from 26.9 +/- 0.2 Ma at the top to 21.3 +/- 0.2 Ma at the bottom of footwall mylonite. Hydrogen isotope (delta D) values in white mica are very low (-150 to -145 %) within the top 80-90 m of detachment footwall, in contrast to values obtained from the deeper part of the section where values range from -77 to -64 %, suggesting that time-integrated interaction between rock and meteoric fluid was restricted to the uppermost part of the mylonitic footwall. Pervasive mica-water hydrogen isotope exchange is difficult to reconcile with models of Ar-40 loss during mylonitization solely by volume diffusion. Rather, we interpret the Ar-40/Ar-39 ages of white mica with low-delta D values to date syn-mylonitic hydrogen and argon isotope exchange, and we conclude that the hydrothermal system of the Northern Snake Range was active during late Oligocene (27-23 Ma) and has been exhumed by the combined effects of ductile strain, extensional detachment faulting, and erosion.

Sedimentary record of the northward flight of India and its collision with Eurasia (Ladakh Himalaya, India)

Stratigraphic and petrographic analysis of the Cretaceous to Eocene Tibetan sedimentary succession has allowed us to reinterpret in detail the sequence of events which led to closure of Neotethys and continental collision in the NW Himalaya. During the Early Cretaceous, the Indian passive margin recorded basaltic magmatic activity. Albian volcanic arenites, probably related to a major extensional tectonic event, are unconformably overlain by an Upper Cretaceous to Paleocene carbonate sequence, with a major quartzarenite episode triggered by the global eustatic sea-level fall at the Cretaceous/Tertiary boundary. At the same time, Neotethyan oceanic crust was being subducted beneath Asia, as testified by calc-alkalic volcanism and forearc basin sedimentation in the Transhimalayan belt. Onset of collision and obduction of the Asian accretionary wedge onto the Indian continental rise was recorded by shoaling of the outer shelf at the Paleocene/Eocene boundary, related to flexural uplift of the passive margin. A few My later, foreland basin volcanic arenites derived from the uplifted Asian subduction complex onlapped onto the Indian continental terrace. All along the Himalaya, marine facies were rapidly replaced by continental redbeds in collisional basins on both sides of the ophiolitic suture. Next, foreland basin sedimentation was interrupted by fold-thrust deformation and final ophiolite emplacement. The observed sequence of events compares favourably with theoretical models of rifted margin to overthrust belt transition and shows that initial phases of continental collision and obduction were completed within 10 to 15 My, with formation of a proto-Himalayan chain by the end of the middle Eocene.

Quantitative integration of geophysical and hydraulic data : from the local towards the regional scale range

Des progrès significatifs ont été réalisés dans le domaine de l'intégration quantitative des données géophysique et hydrologique l'échelle locale. Cependant, l'extension à de plus grandes échelles des approches correspondantes constitue encore un défi majeur. Il est néanmoins extrêmement important de relever ce défi pour développer des modèles fiables de flux des eaux souterraines et de transport de contaminant. Pour résoudre ce problème, j'ai développé une technique d'intégration des données hydrogéophysiques basée sur une procédure bayésienne de simulation séquentielle en deux étapes. Cette procédure vise des problèmes à plus grande échelle. L'objectif est de simuler la distribution d'un paramètre hydraulique cible à partir, d'une part, de mesures d'un paramètre géophysique pertinent qui couvrent l'espace de manière exhaustive, mais avec une faible résolution (spatiale) et, d'autre part, de mesures locales de très haute résolution des mêmes paramètres géophysique et hydraulique. Pour cela, mon algorithme lie dans un premier temps les données géophysiques de faible et de haute résolution à travers une procédure de réduction déchelle. Les données géophysiques régionales réduites sont ensuite reliées au champ du paramètre hydraulique à haute résolution. J'illustre d'abord l'application de cette nouvelle approche dintégration des données à une base de données synthétiques réaliste. Celle-ci est constituée de mesures de conductivité hydraulique et électrique de haute résolution réalisées dans les mêmes forages ainsi que destimations des conductivités électriques obtenues à partir de mesures de tomographic de résistivité électrique (ERT) sur l'ensemble de l'espace. Ces dernières mesures ont une faible résolution spatiale. La viabilité globale de cette méthode est testée en effectuant les simulations de flux et de transport au travers du modèle original du champ de conductivité hydraulique ainsi que du modèle simulé. Les simulations sont alors comparées. Les résultats obtenus indiquent que la procédure dintégration des données proposée permet d'obtenir des estimations de la conductivité en adéquation avec la structure à grande échelle ainsi que des predictions fiables des caractéristiques de transports sur des distances de moyenne à grande échelle. Les résultats correspondant au scénario de terrain indiquent que l'approche d'intégration des données nouvellement mise au point est capable d'appréhender correctement les hétérogénéitées à petite échelle aussi bien que les tendances à gande échelle du champ hydraulique prévalent. Les résultats montrent également une flexibilté remarquable et une robustesse de cette nouvelle approche dintégration des données. De ce fait, elle est susceptible d'être appliquée à un large éventail de données géophysiques et hydrologiques, à toutes les gammes déchelles. Dans la deuxième partie de ma thèse, j'évalue en détail la viabilité du réechantillonnage geostatique séquentiel comme mécanisme de proposition pour les méthodes Markov Chain Monte Carlo (MCMC) appliquées à des probmes inverses géophysiques et hydrologiques de grande dimension . L'objectif est de permettre une quantification plus précise et plus réaliste des incertitudes associées aux modèles obtenus. En considérant une série dexemples de tomographic radar puits à puits, j'étudie deux classes de stratégies de rééchantillonnage spatial en considérant leur habilité à générer efficacement et précisément des réalisations de la distribution postérieure bayésienne. Les résultats obtenus montrent que, malgré sa popularité, le réechantillonnage séquentiel est plutôt inefficace à générer des échantillons postérieurs indépendants pour des études de cas synthétiques réalistes, notamment pour le cas assez communs et importants où il existe de fortes corrélations spatiales entre le modèle et les paramètres. Pour résoudre ce problème, j'ai développé un nouvelle approche de perturbation basée sur une déformation progressive. Cette approche est flexible en ce qui concerne le nombre de paramètres du modèle et lintensité de la perturbation. Par rapport au rééchantillonage séquentiel, cette nouvelle approche s'avère être très efficace pour diminuer le nombre requis d'itérations pour générer des échantillons indépendants à partir de la distribution postérieure bayésienne. - Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending corresponding approaches beyond the local scale still represents a major challenge, yet is critically important for the development of reliable groundwater flow and contaminant transport models. To address this issue, I have developed a hydrogeophysical data integration technique based on a two-step Bayesian sequential simulation procedure that is specifically targeted towards larger-scale problems. The objective is to simulate the distribution of a target hydraulic parameter based on spatially exhaustive, but poorly resolved, measurements of a pertinent geophysical parameter and locally highly resolved, but spatially sparse, measurements of the considered geophysical and hydraulic parameters. To this end, my algorithm links the low- and high-resolution geophysical data via a downscaling procedure before relating the downscaled regional-scale geophysical data to the high-resolution hydraulic parameter field. I first illustrate the application of this novel data integration approach to a realistic synthetic database consisting of collocated high-resolution borehole measurements of the hydraulic and electrical conductivities and spatially exhaustive, low-resolution electrical conductivity estimates obtained from electrical resistivity tomography (ERT). The overall viability of this method is tested and verified by performing and comparing flow and transport simulations through the original and simulated hydraulic conductivity fields. The corresponding results indicate that the proposed data integration procedure does indeed allow for obtaining faithful estimates of the larger-scale hydraulic conductivity structure and reliable predictions of the transport characteristics over medium- to regional-scale distances. The approach is then applied to a corresponding field scenario consisting of collocated high- resolution measurements of the electrical conductivity, as measured using a cone penetrometer testing (CPT) system, and the hydraulic conductivity, as estimated from electromagnetic flowmeter and slug test measurements, in combination with spatially exhaustive low-resolution electrical conductivity estimates obtained from surface-based electrical resistivity tomography (ERT). The corresponding results indicate that the newly developed data integration approach is indeed capable of adequately capturing both the small-scale heterogeneity as well as the larger-scale trend of the prevailing hydraulic conductivity field. The results also indicate that this novel data integration approach is remarkably flexible and robust and hence can be expected to be applicable to a wide range of geophysical and hydrological data at all scale ranges. In the second part of my thesis, I evaluate in detail the viability of sequential geostatistical resampling as a proposal mechanism for Markov Chain Monte Carlo (MCMC) methods applied to high-dimensional geophysical and hydrological inverse problems in order to allow for a more accurate and realistic quantification of the uncertainty associated with the thus inferred models. Focusing on a series of pertinent crosshole georadar tomographic examples, I investigated two classes of geostatistical resampling strategies with regard to their ability to efficiently and accurately generate independent realizations from the Bayesian posterior distribution. The corresponding results indicate that, despite its popularity, sequential resampling is rather inefficient at drawing independent posterior samples for realistic synthetic case studies, notably for the practically common and important scenario of pronounced spatial correlation between model parameters. To address this issue, I have developed a new gradual-deformation-based perturbation approach, which is flexible with regard to the number of model parameters as well as the perturbation strength. Compared to sequential resampling, this newly proposed approach was proven to be highly effective in decreasing the number of iterations required for drawing independent samples from the Bayesian posterior distribution.

Tectonic evolution of the High Himalaya in upper Lahul (NW Himalaya, India)

The Upper Lahul region in the NW Himalaya is located in the transition zone between the High Himalayan Crystalline (HHC) to the SW and the Tethyan Zone sedimentary series to the NE. The tectonic evolution of these domains during the Himalayan Orogeny is the consequence of a succession of five deformation events. An early D1 phase corresponds to synmetamorphic, NE verging folding. This deformation created the Tandi Syncline, which consists of Permian to Jurassic Tethyan metasediments cropping out in the core of a large-scale synformal fold within the HHC paragneiss. This tectonic event is interpreted as related to a NE directed nappe stacking (Shikar Beh Nappe), probably during the late Eocene to the early Oligocene. A subsequent D2a phase caused SW verging folding in the HHC. This deformation is interpreted as contemporaneous with late Oligocene to early Miocene SW directed thrusting along the Main Central Thrust. In the Tethyan Zone, a D2b phase is marked by a decollement thrust, a system of reverse faults, and gentle folds, associated with SW directed tectonic movements. This deformation is related to an imbricate structure, characteristic of a shallow structural level, and developed in the frontal part of a nappe affecting the Tethyan Zone units of SE Zanskar (Nyimaling-Tsarap Nappe). A later D3 phase generated the Chandra Dextral Shear Zone (CDSZ), a large-scale, ductile, dextral strike-slip shear zone, located in the transition zone between the HHC and the Tethyan Himalaya. The CDSZ most likely represents a part of a system of early Miocene extensional and/or dextral, strike-slip shear zones-observed at the HHC-Tethyan Zone contact along the entire Himalaya. A final D4 phase induced large-scale doming and NE:verging back folding.

Thermomechanical analysis of strain localization in a ductile detachment zone

At the latitude of the Thor-Odin dome (British Columbia) the Columbia River Detachment defines the eastern margin of the Shuswap metamorphic core complex and localizes in a 1 km thick muscovite-bearing quartzite mylonite. We present a combined Ar-40/Ar-39, (micro) structural, and oxygen isotope study of the deformation history in the detachment and evaluate the spatial and temporal relationships between microstructure formation and localization of strain. High-precision Ar-40/Ar-39 geochronology from different levels in the mylonite delineates a pattern of increasingly younger (49.0 to 47.9 Ma) deformation ages in deeper levels of the mylonitic footwall. The correlation of Ar-40/Ar-39 ages with decreasing deformation temperatures (similar to 550 degrees-400 degrees C) in the top 200 m of the mylonite indicates that deformation migrated downward from the contact with the hanging wall. Strain localization was diachronous in progressively deeper levels of the footwall and was likely controlled by fluid-assisted strain hardening due to advective heat removal and contemporaneous reaction weakening due to dissolution-reprecipitation of white mica. The observed constant high-stress microstructures across the entire detachment indicate that flow stress was buffered by the interplay of strain rate and temperature, where high strain rates at elevated temperature produced the same microstructure as lower strain rates under decreasing temperature conditions. The combined data suggest that the complex interplay among temporally nonuniform rates of footwall exhumation, heat advection, and embrittlement by meteoric fluids strongly determines the thermomechanical behavior of extensional detachments.

Exhumation history of eastern Ladakh revealed by 40Ar/39Ar and fission-track ages: the Indus River-Tso Morari transect, NW Himalaya

Fission-track and (40)Ar/(39)Ar ages place time constraints on the exhumation of the North Himalayan nappe stack, the Indus Suture Zone and Molasse, and the Transhimalayan Batholith in eastern Ladakh (NW India). Results from this and previous studies on a north-south transect passing near Tso Morari Lake suggest that the SW-directed North Himalayan nappe stack (comprising the Mata, Tetraogal and Tso Morari nappes) was emplaced and metamorphosed by c. 50-45 Ma, and exhumed to moderately shallow depths (c. 10 km) by c. 45-40 Ma. From the mid-Eocene to the present, exhumation continued at a steady and slow rate except for the root zone of the Tso Morari nappe, which cooled faster than the rest of the nappe stack. Rapid cooling occurred at c. 20 Ma and is linked to brittle deformation along the normal Ribil-Zildat Fault concomitant with extrusion of the Crystalline nappe in the south. Data from the Indus Molasse suggest that sediments were still being deposited during the Miocene.