Bimodal magmatism as a consequence of the post-collisional readjustment of the thickened Variscan continental lithosphere (Aiguilles Rouges - Mont Blanc Massifs, Western Alps)

High Precision U-Pb zircon and monazite dating in the Aiguilles Rouges-Mont Blanc area allowed discrimination of three short-lived bimodal magmatic pulses: the early 332 Ma Mg-K Pormenaz monzonite and associated 331 Ma peraluminous Montees Pelissier monzogranite; the 307 Ma cordierite-bearing peraluminous Vallorcine and Fully intrusions; and the 303 Fe-K Mont Blanc syenogranite. All intruded syntectonically along major-scale transcurrent faults at a time when the substratum was experiencing tectonic exhumation, active erosion recorded in detrital basins and isothermal decompression melting dated at 327-320 Ma. Mantle activity and magma mixing are evidenced in all plutons by coeval mafic enclaves, stocks and synplutonic dykes. Both crustal and mantle sources evolve through time, pointing to an increasingly warm continental crust and juvenile asthenospheric mantle sources. This overall tectono-magmatic evolution is interpreted in a scenario of post-collisional restoration to normal size of a thickened continental lithosphere. The latter re-equilibrates through delamination and/or erosion of its mantle root and tectonic exhumation/erosion in an overall extensional regime. Extension is related to either gravitational collapse or back-are extension of a distant subduction zone.

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.

Timing of Palaeozoic magmatism in the Maggia and Sambuco nappes and paleogeographic implications (Central Lepontine Alps)

Magmatic rocks from the pre-Mesozoic basements of the Sambuco and Maggia nappes have been dated by U-Pb zircon ages with the LA-ICPMS technique. Several magmatic events have been identified in the Sambuco nappe. The mafic banded calc-alkaline suite of Scheggia is dated at 540 Ma, an age comparable to that of mafic rocks in the Austroalpine Silvretta nappe. The Sasso Nero peraluminous augengneiss has an age of 480-470 Ma, like many other ``older orthogneisses'' in Alpine basement units. It hosts a large proportion of inherited zircons, which were dated around 630 Ma, a Panafrican age indicating the Gondwanan affiliation of the Sambuco basement. The calc-alkaline Matorello pluton yielded ages around 300 Ma, similar to numerous Late Carboniferous intrusions in other basement units of the Lower Penninic (Monte Leone, Antigorio, Verampio) and Helvetic domains (Gotthard and other External Crystalline Massifs). Associated lamprophyric dykes are slightly younger (300-290 Ma), like similar dykes sampled in gneiss blocks included in the sedimentary cover of the underlying Antigorio nappe (290-285 Ma). The Cocco granodiorite and Ruscada leucogranite, both intruding the basement of the neighbouring Maggia nappe, yielded ages of ca. 300-310 Ma, identical within errors to the age of the Matorello pluton. They are significantly older than former age determinations. This age coincidence, coupled with remarkable petrologic similarities between the Cocco and Matorello granodiorites, strongly suggests paleogeographic proximity of the Sambuco and Maggia nappes in Late Carboniferous times. In recent publications these two nappes have been interpreted as belonging to distinct Mesozoic paleogeographic domains: ``European'' for Sambuco and ``Brian double dagger onnais'' for Maggia, separated by the ``Valais'' oceanic basin. In this case, the similarity of the Matorello and Cocco intrusions would demonstrate the absence of any significant transcurrent movement between these two continental domains. Alternatively, according to a more traditional view, Sambuco and Maggia might belong to a single large Alpine tectonic unit.

Synthèse de l'évolution de la plateforme Urgonienne (Barrémien tardif à Aptien précoce) du Sud-Est de la France : facies, micropaléontologie, géochimie, géométries, paléotectonique et géomodélisation

Il y a environ 125 millions d'années, au Crétacé inférieur, la position des continents et le climat terrestre étaient bien différents de ce que l'on connait aujourd'hui. Le Sud-Est de la France, secteur de cette étude, était alors recouvert d'eau, sous un climat chaud et humide. Sur la bordure de cette étendue d'eau (appelée bassin Vocontien), qui correspond aujourd'hui aux régions de la Provence, du Vaucluse, du Gard, de l'Ardèche et du Vercors, des plateformes carbonatées, (telles que les Bahamas), se développaient. Le calcaire, formé à partir des sédiments accumulés sur ces plateformes, est appelé Urgonien. L'objectif de cette étude est de définir les facteurs qui ont influencé le développement de cette plateforme carbonatée dite « urgonienne » et dans quelle mesure. Plusieurs missions de terrain ont permis de récolter de nombreux échantillons de roche en 52 lieux répartis sur l'ensemble du Sud-Est de la France. Les observations réalisées sur le terrain ainsi que les données acquises en laboratoire (microfaune, microfacies et données géo-chimiques) ont permis, de subdiviser chacune des 52 séries urgoniennes en séquences stratigraphiques et cortèges sédimentaires. La comparaison des épaisseurs et des faciès de chaque cortège sédimentaire permet de concevoir la géométrie et l'évolution paléogéographique de la plateforme urgonienne. Les résultats de cette étude démontrent que son organisation est principalement dirigée par des failles qui ont jouées pendant le dépôt des sédiments. Sur la bordure nord du bassin Vocontien, trois failles subméridiennes contrôlent la géométrie et la répartition des environnements de dépôt. Sur sa bordure sud, ces failles synsédimentaires d'orientation N30° et N110° délimitent des blocs basculés. En tête de bloc, des séries d'épaisseurs réduites à faciès de lagon interne se sont déposées alors que les pieds de blocs sont caractérisés par des épaisseurs importantes et la présence de faciès plus externes. Ces concepts ont ensuite été testés en construisant un modèle numérique en trois dimensions de l'Urgonien du Sud-Est de la France. Sa cohérence avec les données acquises tout au long de cette étude d'une part, et sa cohérence géométrique d'autre part, valide les théories avancées. Des formations équivalentes à l'Urgonien sont réparties dans le monde entier et notamment au Moyen-Orient où elles constituent les réservoirs pétroliers les plus importants. Etre capable de caractériser les facteurs ayant influencé son architecture permet par la suite une meilleure exploitation de ses ressources énergétiques. -- Au Crétacé inférieur, l'intense activité magmatique due à la dislocation du super-continent Pangée influence fortement les conditions environnementales globales. Au Barrémien terminal et Aptien basal, période géologique dont fait l'objet cette étude, le bassin Vocontien, puis Bédoulien, recouvre le Sud-Est de la France, sous un climat chaud et humide. Sur les bordures de ces bassins, des plateformes carbonatées se mettent en place. Les sédiments qui se déposent sur ces plateformes sont à l'origine de la formation urgonienne. Afin d'étudier cette formation, une charte biostratigraphique, principalement basée sur les Orbitolinidés, et un modèle de faciès ont été développés. Les assemblages faunistiques, la succession des faciès, les observations de terrain ainsi que l'étude de signaux géochimiques ont permis le découpage séquentiel de la série urgonienne le long de 54 coupes et puis, répartis sur l'ensemble du Sud-Est de la France. Les corrélations induites par cette étude stratigraphique ont mis en évidence d'importantes variations d'épaisseur et d'environnements de dépôt au sein même de la plateforme urgonienne. Ces variations sont expliquées par le jeu de failles syn-sédimentaires qui ont compartimentées la plateforme urgonienne en blocs. Sur la bordure sud du bassin Vocontien, ces failles d'orientation N30° et N110° délimitent six blocs basculés. Au sommet du Barrémien terminal, la subsidence des blocs situés le plus au sud s'amplifie jusqu'à provoquer l'ouverture du bassin de la Bédoule au sud du secteur d'étude. Cette théorie d'évolution a ensuite été testée par l'élaboration d'un modèle numérique en trois dimensions de l'Urgonien du Sud-Est de la France. Sa cohérence avec les données acquises tout au long de cette étude d'une part, et sa cohérence géométrique d'autre part, valide les théories avancées. Des analogues de l'Urgonien sont répartis dans le monde entier et notamment au Moyen-Orient où ils représentent d'importants réservoirs pétroliers. Être capable de caractériser les facteurs ayant influencé l'architecture de l'Urgonien du Sud-Est de la France permet par la suite une meilleure exploitation de ses ressources énergétiques. -- During the Early Cretaceous epoch, intensive magmatic activity due to the dislocation of the super-continent Pangaea, highly influenced global environmental conditions, which were characterized by a warm and generally humic climate. In this context, carbonate platforms were important in tropical and subtropical shallow-water regions, and especially during the late Barremian and early Aptian, platform carbonates of so-called Urgonian affinity are widespread. In southeastern France, the Urgonian platform was part of the northern Tethyan margin and bordered the Vocontian and the Bedoulian basins. The goal of this thesis was the systematic study of the Urgonian Formation in this region, and in order to achieve this goal, a biostratigraphic chart and a facies model were developed. The faunistic assemblages, the facies succession, the field observations and the study of geochemical signals lead to a sequential subdivision of the Urgonian series along 54 sections and wells allocated in five different regions in southeastern France (Gard, Ardèche, Vercors, Vaucluse and Provence). Correlations from this stratigraphic study highlight important variations in thickness and depositional environments of the Urgonian series. These variations are explained by relative movements induced by syn-sedimentary faults, which divided the Urgonian platforms into blocks. On the southern border of the Vocontian basin, these faults, oriented N30° and N110°, delineate six tilted blocks. At the top of the upper Barremian carbonates, subsidence of the two southern blocks accelerated leading to the opening of the Bedoulian basin. The reconstruction of the sequence-stratigraphic and paleoenvironmental evolution of the Urgonian platforms was then tested by the construction of a 3D numerical model of the Urgonian formation of southeastern France. Firstly, its consistency with the data collected during this study, and secondly, its geometrical coherence validate the proposed theory. Urgonian analogs exist all over the world and particularly in Middle East where they constitute important oil reservoirs. The exact reconstruction of the major factors, which influenced the architecture of these formations, will allow for a better exploitation of these energy resources.

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.

The role of basalt replenishment in the generation of basaltic andesites of the ongoing activity at Arenal volcano, Costa Rica: evidence from clinopyroxene and spinel

The bulk composition of magma erupted from Volcan Arenal has remained nearly constant (SiO2 = 53.6-54.9 wt%; MgO = 5.0-4.5 wt%) during almost 30 years of continuous activity (1969-1996). None the less, clinopyroxene (cpx) phenocrysts and their spinel inclusions record a much more complex open-system evolution in which steady-state production of the erupted basaltic andesitic magma is linked to episodic injections of basalt into Arenal's magma conduit/reservoir system. High-resolution major element zoning profiles (electron microprobe) on a large number of phenocrysts (>14,000 analyses), tied to back-scattered electron (BSE) images, have been used to assess the compositional characteristics of the magmatic end members as well as the timing and dynamics of magma replenishment events. No two cpx phenocrysts have exactly the same zoning profile. The vast majority of our analyses record the crystallization of cpx (Cr2O3 < 0.12 wt%; Mg# = 65-79; Al/Ti = 2-7) from a liquid comparable to or more evolved than erupted magma compositions. However, half of all cpx grains are cored by high-Cr cpx (Cr2O3 = 0.2-0.72 wt%) or contain similar basaltic compositions as abrupt growth bands in phenocrysts with and without high-Cr cores; phenocrysts with high-Cr cpx occur throughout the ongoing activity. In a few cases, high-Cr cpx occurs very near the outer margin of the grain without an apparent growth hiatus, particularly in 1968/69 and 1992/93. The main conclusions are: (1) all basaltic andesitic lavas erupted at Arenal during the ongoing activity that began in July, 1968, are the products of magma mixing, (2) clinopyroxenes record multiple replenishment events of basaltic magma in contrast to the near constancy of erupted bulk compositions, (3) some phenocrysts preserve records of multiple interactions with basaltic magmas requiring magmatic processes to operate on time-scales shorter than residence times of some phenocrysts, (4) multiple occurrences of clinopyroxene with high-Cr rims suggest that basalt replenishment events have occurred with sub-decadal frequency and may predate eruption by months or less. From this we infer that Arenal volcano is underlain by a continuously active, small-volume magmatic reservoir maintained in quasi-steady state by basalt recharge over several decades. The monotony of erupting Arenal magmas implies that fractionation, recharge, ascent, and eruption are well balanced in order for magmas to be essentially uniform while containing phenocrysts with vastly different growth histories at the time of eruption.

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.

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.

Low-pressure, water-assisted anatexis of basic dykes in a contact metamorphic aureole, Fuerteventura (Canary Islands): oxygen isotope evidence for a meteoric fluid origin

Migmatites produced by low-pressure anatexis of basic dykes are found in a contact metamorphic aureole around a pyroxenite-gabbro intrusion (PX2), on Fuerteventura. Dykes outside and inside the aureole record interaction with meteoric water, with low or negative delta O-18 whole-rock values (+0.2 to -3.4 parts per thousand), decreasing towards the contact. Recrystallised plagioclase, diopside, biotite and oxides, from within the aureole, show a similar evolution with lowest delta O-18 values (-2.8, -4.2, - 4.4 and -7.6 parts per thousand, respectively) in the migmatite zone, close to the intrusion. Relict clinopyroxene phenocrysts preserved in all dykes, retain typically magmatic delta O-18 values up to the anatectic zone, where the values are lower and more heterogeneous. Low delta O-18 values, decreasing towards the intrusion, can be ascribed to the advection of meteoric water during magma emplacement, with increasing fluid/rock ratios (higher dyke intensities towards the intrusion acting as fluid-pathways) and higher temperatures promoting increasing exchange during recrystallisation.

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.

Late Jurassic oceanic crust and Upper Cretaceous Caribbean plateau picritic basalts exposed in the Duarte igneous complex, Hispaniola

Four distinct rock units have been recognized near El Aguacate, in the Janico-Juncalito-La Vega area of the Duarte complex (Dominican Republic): (1) serpentinites crosscut by numerous diabasic dikes, (2) basalts interbedded with Late Jurassic ribbon cherts, (3) picrites and ankaramites relatively enriched in incompatible trace elements, and (4) amphibolites and gneissic amphibolites chemically similar to Oceanic Plateau Basalts. Similar Ar-Ar ages of late magmatic amphibole from a picrite, and hornblende from an amphibolite (86.1 +/- 1.3 Ma and 86.7 +/- 1.6 Ma, respectively), suggest that the Duarte picrites are contemporaneous with the Deep Sea Drilling Program Leg 15 and Ocean Drilling Program Leg 126 basalts drilled from the Caribbean oceanic plateau. These basalts are associated with sediments containing Late Cretaceous faunas. Sr, Nd, and Pb data show that enriched picrites and amphibolites are isotopically similar to mafic lavas from previously described Caribbean plateau and Galapagos hotspot basalts. Major element, trace element, and lead isotopic features of Late Jurassic basalts and diabases are consistent with those of normal oceanic crust basalt. However, these basalts differ from typical N-MORB because they have lower epsilon Nd ratios that plot within the range of Ocean Island Basalts. These rocks appear to represent remnants of the Caribbean Jurassic oceanic crust formed from an oceanic ridge possibly close to a hotspot. Later, they were tectonically juxtaposed with Late Cretaceous slices of the Caribbean-Colombian plateau.

Enriched basaltic andesites from mid-crustal fractional crystallization, recharge, and assimilation (Pilavo Volcano, Western Cordillera of Ecuador)

The origin of andesite is an important issue in petrology because andesite is the main eruptive product at convergent margins, corresponds to the average crustal composition and is often associated with major Cu-Au mineralization. In this study we present petrographic, mineralogical, geochemical and isotopic data for basaltic andesites of the latest Pleistocene Pilavo volcano, one of the most frontal volcanoes of the Ecuadorian Quaternary arc, situated upon thick (30-50 km) mafic crust composed of accreted Cretaceous oceanic plateau rocks and overlying mafic to intermediate Late Cretaceous-Late Tertiary magmatic arcs. The Pilavo rocks are basaltic andesites (54-57 center dot 5 wt % SiO(2)) with a tholeiitic affinity as opposed to the typical calc-alkaline high-silica andesites and dacites (SiO(2) 59-66 wt %) of other frontal arc volcanoes of Ecuador (e.g. Pichincha, Pululahua). They have much higher incompatible element contents (e.g. Sr 650-1350 ppm, Ba 650-1800 ppm, Zr 100-225 ppm, Th 5-25 ppm, La 15-65 ppm) and Th/La ratios (0 center dot 28-0 center dot 36) than Pichincha and Pululahua, and more primitive Sr ((87)Sr/(86)Sr similar to 0 center dot 7038-0 center dot 7039) and Nd (epsilon(Nd) similar to +5 center dot 5 to +6 center dot 1) isotopic signatures. Pilavo andesites have geochemical affinities with modern and recent high-MgO andesites (e.g. low-silica adakites, Setouchi sanukites) and, especially, with Archean sanukitoids, for both of which incompatible element enrichments are believed to result from interactions of slab melts with peridotitic mantle. Petrographic, mineral chemistry, bulk-rock geochemical and isotopic data indicate that the Pilavo magmatic rocks have evolved through three main stages: (1) generation of a basaltic magma in the mantle wedge region by flux melting induced by slab-derived fluids (aqueous, supercritical or melts); (2) high-pressure differentiation of the basaltic melt (at the mantle-crust boundary or at lower crustal levels) through sustained fractionation of olivine and clinopyroxene, leading to hydrous, high-alumina basaltic andesite melts with a tholeiitic affinity, enriched in incompatible elements and strongly impoverished in Ni and Cr; (3) establishment of one or more mid-crustal magma storage reservoirs in which the magmas evolved through dominant amphibole and clinopyroxene (but no plagioclase) fractionation accompanied by assimilation of the modified plutonic roots of the arc and recharge by incoming batches of more primitive magma from depth. The latter process has resulted in strongly increasing incompatible element concentrations in the Pilavo basaltic andesites, coupled with slightly increasing crustal isotopic signatures and a shift towards a more calc-alkaline affinity. Our data show that, although ultimately originating from the slab, incompatible element abundances in arc andesites with primitive isotopic signatures can be significantly enhanced by intra-crustal processes within a thick juvenile mafic crust, thus providing an additional process for the generation of enriched andesites.

Arrested kinetic Li isotope fractionation at the margin of the llimaussaq complex, South Greenland: Evidence for open-system processes during final cooling of peralkaline igneous rocks

Li contents [Li] and isotopic composition (delta Li-7) of mafic minerals (mainly amphibole and clinopyroxene) from the alkaline to peralkaline Ilimaussaq plutonic complex, South Greenland, track the behavior of Li and its isotopes during magmatic differentiation and final cooling of an alkaline igneous system. [Li] in amphibole increase from < 10 ppm in Caamphiboles of the least differentiated unit to >3000 ppm in Na-amphiboles of the highly evolved units. In contrast, [Li] in clinopyroxene are comparatively low (<85 ppm) and do not vary systematically with differentiation. The distribution of Li between amphibole and pyroxene is controlled by the major element composition of the minerals (Ca-rich and Na-rich, respectively) and changes in oxygen fugacity (due to Li incorporation via coupled substitution with ferric iron) during magmatic differentiation. delta(7) Li values of all minerals span a wide range from + 17 to - 8 parts per thousand, with the different intrusive units of the complex having distinct Li isotopic systematics. Amphiboles, which dominate the Li budget of whole-rocks from the inner part of the complex, have constant delta Li-7 of + 1.8 +/- 2.2 parts per thousand (2 sigma, n = 15). This value reflects a homogeneous melt reservoir and is consistent with their mantle derivation, in agreement with published O and Nd isotopic data. Clinopyroxenes of these samples are consistently lighter, with Delta Li-7(amph-cpx). as large as 8 parts per thousand and are thus not in Li isotope equilibrium. These low values probably reflect late-stage diffusion of Li into clinopyroxene during final cooling of the rocks, thus enriching the clinopyroxene in 6 Li. At the margin of the complex delta(7) Li in the syenites increases systematically, from +2 to high values of + 14 parts per thousand. This, coupled with the observed Li isotope systematics of the granitic country rocks, reflects post-magmatic open-system processes occurring during final cooling of the intrusion. Although the shape and magnitude of the Li isotope and elemental profiles through syenite and country rock are suggestive of diffusion-driven isotope fractionation, they cannot be modeled by one-dimensional diffusive transport and point to circulation of a fluid having a high 67 Li value (possibly seawater) along the chilled contact. In all, this study demonstrates that Li isotopes can be used to identify complex fluid- and diffusion-governed processes taking place during the final cooling of such rocks. (c) 2007 Elsevier B.V All rights reserved.