Modelling changes in stable isotope compositions of minerals during net transfer reactions in a contact aureole: Wollastonite growth at the northern Hunter Mountain Batholith (Death Valley National Park, USA)

One of the world's largest wollastonite deposits was formed at the contact of the northern Hunter Mountain Batholith (California, USA) in Paleozoic sediments. Wollastonite occurs as zones of variable thickness surrounding layers or nodules of quartzite in limestones. A minimum formation temperature of 650 degrees C is estimated from isolated periclase-bearing lenses in that area. Contact metamorphism of siliceous carbonates has produced mineral assemblages that are consistent with heterogeneous, and partly limited infiltration of water-rich fluids, compatible with O-18/O-16 and C-13/C-12 isotopic patterns recorded in carbonates. Oxygen isotope compositions of wollastonites in the study area may also not require infiltration of large quantities of externally-derived fluids that were out of equilibrium with the rocks. 8180 values of wollastonite are high (14.8 parts per thousand to 25.0 parts per thousand; median: 19.7 parts per thousand) and close to those of the host limestone (19.7 parts per thousand to 28 parts per thousand; median: 24.9 parts per thousand) and quartz (18.0 parts per thousand. to 29.1 parts per thousand; median: 22.6 parts per thousand). Isotopic disequilibrium exists at quartz/wollastonite and wollastonite/calcite boundaries. Therefore, classical batch/Rayleigh fractionation models based on reactant and product equilibrium are not applicable to the wollastonite rims. An approach that relies on local instantaneous mass balance for the reactants, based on the wollastonite-forming reaction is suggested as an alternative way to model wollastonite reaction rims. This model reproduces many of the measured delta O-18 values of wollastonite reaction rims of the current study to within +/- 1 parts per thousand, even though the wollastonite compositions vary by almost 10 parts per thousand. (C) 2011 Elsevier B.V. All rights reserved.

Implications of the serpentine phase transition on the behaviour of beryllium and lithium-boron of subducted ultramafic rocks

The Totalp-Platta-Malenco ophiolites in the Eastern Central Alps offer a unique opportunity to study the behaviour of Li, Be and B in ultramafic rocks in response to serpentinization and to progressive Alpine metamorphism. These units represent the remnants of a former ocean-continent transition that was intensely serpentinized during exposure on the Jurassic seafloor of the Ligurian Tethys. From north to the south, three isograd reactions (lizardite double right arrow antigorite + brucite; lizardite + talc double right arrow antigorite; lizardite + tremolite double right arrow antigorite + diopside) have been used to quantify the evolution of the light element content of metamorphic minerals. We determined the Li, Be and B concentrations in major silicate minerals from the ultramafic bodies of Totalp, Platta and Malenco by secondary ion mass spectrometry. Mantle minerals have Be concentrations (e.g. <0.001-0.009 mu g/g in olivine) similar to the metamorphic minerals that replace them (e.g. <0.001-0.016 mu g/g in serpentine). The mantle signature of Be is thus neither erased during seafloor alteration nor by progressive metamorphism from prehnite-pumpellyite to epidote-amphibolite facies. In contrast, the Li and B inventories of metamorphic minerals are related to the lizardite-to-antigorite transition. Both elements display higher concentrations in the low-temperature serpentine polymorph lizardite (max. 156 mu/g Li, max. 318 mu g/g B) than in antigorite (max. 0.11 mu g/g Li, max. 12 mu g/g B). Calculated average B/Li ratios for lizardite (similar to 1395) and antigorite (similar to 115) indicate that Li fractionates from B during the lizardite-to-antigorite transition during prograde metamorphism in ultramafic rocks. In subduction zones, this signature is likely to be recorded in the B-rich nature of forearc fluids. Relative to oceanic mantle the Be content of mantle clinopyroxene is much higher, but similar to Be values from mantle xenoliths and subduction-related peridotite massifs. These data support previous hypothesis that the mantle rocks from the Eastern Central Alps have a subcontinental origin. We conclude that Be behaves conservatively during subduction metamorphism of ultramafic rocks, at least at low-temperature, and thus retains the fingerprint of ancient subduction-related igneous events in mantle peridotites. (C) 2010 Elsevier Ltd. All rights reserved.

Linking temperature estimates and microstructures in deformed polymineralic mantle rocks

To constrain deformation temperatures of mantle shear zones, we studied a strike-slip shear zone (Hilti massif, Semail ophiolite, Oman) and focused on the interaction between microstructural mechanisms and chemical equilibration processes. Quantitative microfabric analysis on harzburgites with different deformation intensity (porphyroclastic tectonite, mylonite, and ultramylonite) was combined with orthopyroxene geothermometry. The average grain size of all phases decreases with decreasing shear zone thickness. Dynamic recrystallization of porphyroclasts in combination with dissolution-precipitation and nucleation result in small-sized, chemically equilibrated pyroxenes. The composition of orthopyroxene was used to calculate deformation temperatures. In the case of the porphyroclastic tectonites, the chemical composition of orthopyroxene has been reset by diffusion yielding temperature estimates of 880-900 degrees C. The mylonites were deformed by dislocation creep of olivine and show a broad range of calculated temperatures, which result from a combination of grain size reduction and inheritance of equilibrium compositions from earlier high-temperature events and diffusion. In mylonites, diffusion profiles combined with geothermometry and grain size analysis indicate a mylonitic deformation temperature of 800-900 degrees C possibly followed by diffusion. In ultramylonites, the smallest grains (<30 mu m) reveal equilibration at temperatures of similar to 700 degrees C during the last stages of ductile deformation, which was dominated by diffusion creep of olivine. Our results provide a crucial link between temperature and evolution of microstructures from dislocation creep to diffusion creep in mantle shear zones.

Seismic attenuation and velocity dispersion in heterogeneous partially saturated porous rocks

Using a numerical approach, we explore wave-induced fluid flow effects in partially saturated porous rocks in which the gas-water saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighbouring regions can exhibit different levels of saturation. To determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. The multiscale nature of the fluid saturation is taken into account by locally computing the physical properties of an effective fluid, which are then used for the larger-scale simulations. We consider two sets of numerical experiments to analyse such effects in heterogeneous partially saturated porous media, where the saturation field is determined by variations in porosity and clay content, respectively. In both cases we also evaluate the seismic responses of corresponding binary, patchy-type saturation patterns. Our results indicate that significant attenuation and modest velocity dispersion effects take place in this kind of media for both binary patchy-type and spatially continuous gas saturation patterns and in particular in the presence of relatively small amounts of gas. The numerical experiments also show that the nature of the gas distribution patterns is a critical parameter controlling the seismic responses of these environments, since attenuation and velocity dispersion effects are much more significant and occur over a broader saturation range for binary patchy-type gas-water distributions. This analysis therefore suggests that the physical mechanisms governing partial saturation should be accounted for when analysing seismic data in a poroelastic framework. In this context, heterogeneities associated with the dry frame properties, which do not play important roles in wave-induced fluid flow processes per se, should be taken into account since they may determine the kind of gas distribution pattern taking place in the porous rock.

Seismic attenuation in heterogeneous partially saturated rocks

We study wave-induced fluid flow effects in porous rocks partially saturated with gas and water, where the saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighboring regions can exhibit different levels of saturation. In order to determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. We consider numerical experiments to analyze such effects in heterogeneous partially saturated porous media, where the saturation field is determined by realistic variations in porosity. Our results indicate that the spatially continuous nature of gas saturation inherent to this study is a critical parameter controlling the seismic response of these environments, which in turn suggests that the physical mechanisms governing partial saturation should be accounted for when analyzing seismic data in a poro-elastic context.

Nucleation and growth of metamorphic minerals in contact aureoles

THESIS ABSTRACT Nucleation and growth of metamorphic minerals are the consequence of changing P-T-X-conditions. The thesis presented here focuses on processes governing nucleation and growth of minerals in contact metamorphic environments using a combination of geochemical analytics (chemical-, isotope-, and trace element composition), statistical treatments of spatial data, and numerical models. It is shown, that a combination of textural modeling and stable isotope analysis allows a distinction between several possible reaction paths for olivine growth in a siliceous dolomite contact aureole. It is suggested that olivine forms directly from dolomite and quartz. The formation of olivine from this metastable reaction implies metamorphic crystallization far from equilibrium. As a major consequence, the spatial distribution of metamorphic mineral assemblages in a contact aureole cannot be interpreted as a proxy for the temporal evolution of a single rock specimen, because each rock undergoes a different reaction path, depending on temperature, heating rate, and fluid-infiltration rate. A detailed calcite-dolomite thermometry study was initiated on multiple scales ranging from aureole scale to the size of individual crystals. Quantitative forward models were developed to evaluate the effect of growth zoning, volume diffusion and the formation of submicroscopic exsolution lamellae (<1 µm) on the measured Mg-distribution in individual calcite crystals and compare the modeling results to field data. This study concludes that Mg-distributions in calcite grains of the Ubehebe Peak contact aureole are the consequence of rapid crystal growth in combination with diffusion and exsolution. The crystallization history of a rock is recorded in the chemical composition, the size and the distribution of its minerals. Near the Cima Uzza summit, located in the southern Adamello massif (Italy), contact metamorphic brucite bearing dolomite marbles are exposed as xenoliths surrounded by mafic intrusive rocks. Brucite is formed retrograde pseudomorphing spherical periclase crystals. Crystal size distributions (CSD's) of brucite pseudomorphs are presented for two profiles and combined with geochemistry data and petrological information. Textural analyses are combined with geochemistry data in a qualitative model that describes the formation periclase. As a major outcome, this expands the potential use of CSD's to systems of mineral formation driven by fluid-infiltration. RESUME DE LA THESE La nucléation et la croissance des minéraux métamorphiques sont la conséquence de changements des conditions de pression, température et composition chimique du système (PT-X). Cette thèse s'intéresse aux processus gouvernant la nucléation et la croissance des minéraux au cours d'un épisode de métamorphisme de contact, en utilisant la géochimie analytique (composition chimique, isotopique et en éléments traces), le traitement statistique des données spatiales et la modélisation numérique. Il est montré que la combinaison d'un modèle textural avec des analyses en isotopes stables permet de distinguer plusieurs chemins de réactions possibles conduisant à la croissance de l'olivine dans une auréole de contact riche en Silice et dolomite. Il est suggéré que l'olivine se forme directement à partir de la dolomie et du quartz. Cette réaction métastable de formation de l'olivine implique une cristallisation métamorphique loin de l'équilibre. La principale conséquence est que la distribution spatiale des assemblages de minéraux métamorphiques dans une auréole de contact ne peut pas être considérée comme un témoin de l'évolution temporelle d'un type de roche donné, puisque chaque type de roche suit différents chemins de réactions, en fonction de la température, la vitesse de réchauffement et le taux d'infiltration du fluide. Une étude thermométrique calcite-dolomite détaillée a été réalisée à diverses échelles, depuis l'échelle de l'auréole de contact jusqu'à l'échelle du cristal. Des modèles numériques quantitatifs ont été développés pour évaluer l'effet des zonations de croissance, de la diffusion volumique et de la formation de lamelles d'exsolution submicroscopiques (<1µm) sur la distribution du magnésium mesuré dans des cristaux de calcite individuels. Les résultats de ce modèle ont été comparés ä des échantillons naturels. Cette étude montre que la distribution du Mg dans les grains de calcite de l'auréole de contact de l'Ubehebe Peak (USA) résulte d'une croissance cristalline rapide, associée aux processus de diffusion et d'exsolution. L'histoire de cristallisation d'une roche est enregistrée dans la composition chimique, la taille et la distribution de ses minéraux. Près du sommet Cima Uzza situé au sud du massif d'Adamello (Italie), des marbres dolomitiques à brucite du métamorphisme de contact forment des xénolithes dans une intrusion mafique. La brucite constitue des pseudomorphes rétrogrades du périclase. Les distributions de taille des cristaux (CSD) des pseudomorphes de brucite sont présentées pour deux profiles et sont combinées aux données géochimiques et pétrologiques. Les analyses textorales sont combinées aux données géochimiques dans un modèle qualitatif qui décrit la formation du périclase. Ceci élargit l'utilisation potentielle de la C5D aux systèmes de formation de minéraux controlés par les infiltrations fluides. THESIS ABSTRACT (GENERAL PUBLIC) Rock textures are essentially the result of a complex interaction of nucleation, growth and deformation as a function of changing physical conditions such as pressure and temperature. Igneous and metamorphic textures are especially attractive to study the different mechanisms of texture formation since most of the parameters like pressure-temperature-paths are quite well known for a variety of geological settings. The fact that textures are supposed to record the crystallization history of a rock traditionally allowed them to be used for geothermobarometry or dating. During the last decades the focus of metamorphic petrology changed from a static point of view, i.e. the representation of a texture as one single point in the petrogenetic grid towards a more dynamic view, where multiple metamorphic processes govern the texture formation, including non-equilibrium processes. This thesis tries to advance our understanding on the processes governing nucleation and growth of minerals in contact metamorphic environments and their dynamic interplay by using a combination of geochemical analyses (chemical-, isotope-, and trace element composition), statistical treatments of spatial data and numerical models. In a first part the thesis describes the formation of metamorphic olivine porphyroblast in the Ubehebe Peak contact aureole (USA). It is shown that not the commonly assumed succession of equilibrium reactions along a T-t-path formed the textures present in the rocks today, but rather the presence of a meta-stable reaction is responsible for forming the olivine porphyroblast. Consequently, the spatial distribution of metamorphic minerals within a contact aureole can no longer be regarded as a proxy for the temporal evolution of a single rock sample. Metamorphic peak temperatures for samples of the Ubehebe Peak contact aureole were determined using calcite-dolomite. This geothermometer is based on the temperature-dependent exchange of Mg between calcite and dolomite. The purpose of the second part of this thesis was to explain the interfering systematic scatter of measured Mg-content on different scales and thus to clarify the interpretation of metamorphic temperatures recorded in carbonates. Numerical quantitative forward models are used to evaluate the effect of several processes on the distribution of magnesium in individual calcite crystals and the modeling results were then compared to measured field. Information about the crystallization history is not only recorded in the chemical composition of grains, like isotope composition or mineral zoning. Crystal size distributions (CSD's) provide essential information about the complex interaction of nucleation and growth of minerals. CSD's of brucite pseudomorphs formed retrograde after periclase of the southern Adamello massif (Italy) are presented. A combination of the textural 3D-information with geochemistry data is then used to evaluate reaction kinetics and to constrain the actual reaction mechanism for the formation of periclase. The reaction is shown to be the consequence of the infiltration of a limited amount of a fluid phase at high temperatures. The composition of this fluid phase is in large disequilibrium with the rest of the rock resulting in very fast reaction rates. RESUME DE LA THESE POUR LE GRAND PUBLIC: La texture d'une roche résulte de l'interaction complexe entre les processus de nucléation, croissance et déformation, en fonction des variations de conditions physiques telles que la pression et la température. Les textures ignées et métamorphiques présentent un intérêt particulier pour l'étude des différents mécanismes à l'origine de ces textures, puisque la plupart des paramètres comme les chemin pression-température sont relativement bien contraints dans la plupart des environnements géologiques. Le fait que les textures soient supposées enregistrer l'histoire de cristallisation des roches permet leur utilisation pour la datation et la géothermobarométrie. Durant les dernières décennies, la recherche en pétrologie métamorphique a évolué depuis une visualisation statique, c'est-à-dire qu'une texture donnée correspondait à un point unique de la grille pétrogénétique, jusqu'à une visualisation plus dynamique, où les multiples processus métamorphiques qui gouvernent 1a formation d'une texture incluent des processus hors équilibre. Cette thèse a pour but d'améliorer les connaissances actuelles sur les processus gouvernant la nucléation et la croissance des minéraux lors d'épisodes de métamorphisme de contact et l'interaction dynamique existant entre nucléation et croissance. Pour cela, les analyses géochimiques (compositions chimiques en éléments majeurs et traces et composition isotopique), le traitement statistique des données spatiales et la modélisation numérique ont été combinés. Dans la première partie, cette thèse décrit la formation de porphyroblastes d'olivine métamorphique dans l'auréole de contact de l'Ubehebe Peak (USA). Il est montré que la succession généralement admise des réactions d'équilibre le long d'un chemin T-t ne peut pas expliquer les textures présentes dans les roches aujourd'hui. Cette thèse montre qu'il s'agirait plutôt d'une réaction métastable qui soit responsable de la formation des porphyroblastes d'olivine. En conséquence, la distribution spatiale des minéraux métamorphiques dans l'auréole de contact ne peut plus être interprétée comme le témoin de l'évolution temporelle d'un échantillon unique de roche. Les pics de température des échantillons de l'auréole de contact de l'Ubehebe Peak ont été déterminés grâce au géothermomètre calcite-dolomite. Celui-ci est basé sur l'échange du magnésium entre la calcite et la dolomite, qui est fonction de la température. Le but de la deuxième partie de cette thèse est d'expliquer la dispersion systématique de la composition en magnésium à différentes échelles, et ainsi d'améliorer l'interprétation des températures du métamorphisme enregistrées dans les carbonates. Des modèles numériques quantitatifs ont permis d'évaluer le rôle de différents processus sur la distribution du magnésium dans des cristaux de calcite individuels. Les résultats des modèles ont été comparés aux échantillons naturels. La composition chimique des grains, comme la composition isotopique ou la zonation minérale, n'est pas le seul témoin de l'histoire de la cristallisation. La distribution de la taille des cristaux (CSD) fournit des informations essentielles sur les interactions entre nucléation et croissance des minéraux. La CSD des pseudomorphes de brucite retrograde formés après le périclase dans le sud du massif Adamello (Italie) est présentée dans la troisième partie. La combinaison entre les données textorales en trois dimensions et les données géochimiques a permis d'évaluer les cinétiques de réaction et de contraindre les mécanismes conduisant à la formation du périclase. Cette réaction est présentée comme étant la conséquence de l'infiltration d'une quantité limitée d'une phase fluide à haute température. La composition de cette phase fluide est en grand déséquilibre avec le reste de la roche, ce qui permet des cinétiques de réactions très rapides.

The magmatic to hydrothermal evolution of the intrusive Mont Saint-Hilaire Complex: Insights into the late-stage evolution of peralkaline rocks

The Cretaceous Mont Saint-Hilaire complex (Quebec, Canada) comprises three major rock units that were emplaced in the following sequence: (I) gabbros; (II) diorites; (III) diverse partly agpaitic foid syenites. The major element compositions of the rock-forming minerals, age-corrected Nd and oxygen isotope data for mineral separates and trace element data of Fe-Mg silicates from the various lithologies imply a common source for all units. The distribution of the rare earth elements in clinopyroxene from the gabbros indicates an ocean island basalt type composition for the parental magma. Gabbros record temperatures of 1200 to 800 degrees C, variable silica activities between 0 center dot 7 and 0 center dot 3, and f(O2) values between -0 center dot 5 and +0 center dot 7 (log delta FMQ, where FMQ is fayalite-magnetite-quartz). The diorites crystallized under uniform a(SiO2) (a(SiO2) = 0 center dot 4-0 center dot 5) and more reduced f(O2) conditions (log delta FMQ similar to-1) between similar to 1100 and similar to 800 degrees C. Phase equilibria in various foid syenites indicate that silica activities decrease from 0 center dot 6-0 center dot 3 at similar to 1000 degrees C to < 0 center dot 3 at similar to 550 degrees C. Release of an aqueous fluid during the transition to the hydrothermal stage caused a(SiO2) to drop to very low values, which results from reduced SiO(2) solubilities in aqueous fluids compared with silicate melts. During the hydrothermal stage, high water activities stabilized zeolite-group minerals. Fluid inclusions record a complex post-magmatic history, which includes trapping of an aqueous fluid that unmixed from the restitic foid syenitic magma. Cogenetic aqueous and carbonic fluid inclusions reflect heterogeneous trapping of coexisting immiscible external fluids in the latest evolutionary stage. The O and C isotope characteristics of fluid-inclusion hosted CO(2) and late-stage carbonates imply that the surrounding limestones were the source of the external fluids. The mineral-rich syenitic rocks at Mont Saint-Hilaire evolved as follows: first, alkalis, high field strength and large ion lithophile elements were pre-enriched in the (late) magmatic and subsequent hydrothermal stages; second, percolation of external fluids in equilibrium with the carbonate host-rocks and mixing processes with internal fluids as well as fluid-rock interaction governed dissolution of pre-existing minerals, element transport and precipitation of mineral assemblages determined by locally variable parameters. It is this hydrothermal interplay between internal and external fluids that is responsible for the mineral wealth found at Mont Saint-Hilaire.

Prograde mica 40Ar/39Ar growth ages recorded in high pressure rocks (Syros, Cyclades, Greece)

Phengites from the eclogite and blueschist-facies sequences of the Cycladic island of Syros (Greece) have been dated by the in situ UV-laser ablation Ar-40/Ar-39 method. A massive, phengite-rich eclogite and an omphacite-rich metagabbro were investigated. The phengites are eubedral and coarse-grained (several 100 mum), strain-free and exhibit no evidence for late brittle deformation or recrystallization. Apparent ages in these samples range from 43 to 50 Ma for the phengite-rich eclogite and 42 to 52 Ma for the ompbacitic metagabbro. This large spread of ages is visible at all scales-within individual grains as well as in domains of several 100 mum and across the entire sample (ca. 2 cm). Such variations have been traditionally attributed to metamorphic cooling or the incorporation of excess argon. However, the textural equilibrium between the phengites and other high pressure phases and the subtle compositional variations within the phengites, especially the preservation of growth textures, alternatively suggest that the observed range in ages may reflect variations of radiogenic argon acquired during phengite formation and subsequent growth, thus dating a discrete event on the prograde path. This implies that the oldest phengite 40Ar/39Ar ages provide the best estimate of a minimum crystallization age, which is in agreement with recently reported U-Pb and Lu-Hf geochronological data. Our results are consistent with available stable isotope data and further suggest that, under fluid-restricted conditions, both stable and radiogenic isotopic systems can survive without significant isotopic exchange during subduction and exhumation from eclogite-facies P-T conditions. (C) 2004 Elsevier B.V. All rights reserved.

Emplacement of ultramafic rocks into the continental crust monitored by light and other trace elements: An example from the Geisspfad body (Swiss-Italian Alps)

In order to evaluate the influence of continental crustal rocks on trace element budgets of serpentinized peridotites incorporated into the continental crust, we have analyzed the chemical composition of whole rock samples and minerals of the Geisspfad ultramafic complex (Swiss-Italian Alps). This complex represents a relict oceanic succession composed of serpentinites, ophicarbonates and metabasic rocks, emplaced into crustal gneisses during Alpine collision. Following peak metamorphic amphibolite facies conditions, fluid flow modified some of the trace element contents of ophicarbonates and deformed serpentinites close to the contact with country rocks. The fluid originated from the surrounding continental crustal rocks as documented by the increase of Pb in the serpentinites, and by the strongly negative all) values (-112 parts per thousand) of some ultramafic rocks close to the contact with surrounding gneisses. Little or no modification of the fluid mobile elements Li, B or U was observed in the serpentinite. In-situ analysis of light elements of serpentinite minerals indicate redistribution of light elements coupled to changes of mineral modes towards the outer 100-150 m of the massif. In the centre of the massif, Li is preferentially concentrated in olivine, while Be and B are hosted by tremolite. In contrast, at the outer rim of the massif, Li and Be are preferentially incorporated into diopside, and B into antigorite. This redistribution of light elements among the different minerals is visible in the serpentinite, at a maximum distance of -100-150 m from the ophicarbonate-metabasite contact. Our results show that interaction of ultramafic rocks and crust-derived fluids can be easily detected by studies of Pb and partial derivative D in whole rocks. We argue that small ultramafic bodies potentially record an emplacement-related trace element signature, and that crustal light element values in ultramafic rocks are not necessarily derived from a subducting slab. (C) 2008 Elsevier B.V. All rights reserved.

Late Oligocene Larger Foraminifera from Nosara (Nicoya Península, Costa Rica) and Windward (Carriacou, Lesser Antilles), calibrated by 87Sr/86Sr isotope stratigraphy.

ABSTRACT: The Late Oligocene first occurrence of Miogypsina gunteri Cole and Miogypsina tani Drooger is verified in the light of larger foraminiferal assemblages from Nosara (Nicoya Peninsula, Costa Rica) and Windward (Carriacou, Grenada, Lesser Antilles). At Windward, they co-occur with planktonic foraminifera and nannofossils studied by earlier workers. 87Sr/86Sr isotope ratios were measured in calcite of larger foraminifera at both sites to determine independent absolute ages of the two outcrops. Late Oligocene shallow water formations unconformably overlie Paleocene-Eocene distal turbidites and siliceous shales in the outcrops in Costa Rica. The shallow-water lithostratigraphy includes near shore volcanic sandstones that alternate with sandy bioclastic limestones. The latter have yielded rich assemblages of larger foraminifera that have been studied in oriented sections, SEM for split material and polished rock thin sections for transmitted light and cathodoluminescence (CL) observation. The association of larger foraminifera includes Heterostegina antillea Cushman, Miogypsina tani Drooger, Miogypsina gunteri Cole, Miogypsina c.f (Miolepidocyclina) panamensis (Cushman), Miogypsina sp., Lepidocyclina (nephrolepidina) vaughani Cushman, Lepidocyclina yurnagurensis Cushman, Lepidocyclina undosa Cushman. At Windward, scattered outcrops of the Belvedere Formation expose mass flow deposits rich in larger foraminifera associated with an early late Oligocene (NP24) planktonic foraminifa assemblage. Planktonic foraminifera and nanno-fossils recovered in rocks slightly downsection indicate an upper Middle Oligocene age. Larger foraminifera species collected just north of Windward include Lepidocyclina undosa Cushman and Miogypsina gunteri Cole. The 87Sr/86Sr ratio was measured in 11 specimens of megalosphaeric Lepidocyclina spp., mechanically extracted from 2 rock samples collected at Punta Peladas (Costa Rica) Ratios range from 0.709088 to 0.708196, which correspond to a model age of 27.4 to 24.94 Ma. Very similar ratios, ranging from 0.708150 to 0.708167 were obtained from Windward Village (Carriacou). All measured 87Sr/86Sr ratios fall within a Chattian (Late Oligocene) age. Nannofossils, planktonic foraminifera and identical 87Sr/86Sr ratios from Punta Peladas and Windward clearly indicate a Chattian (Upper Oligocene) first occurrence of Miogypsina gunteri and Miogypsina tani and confirm this age range for the other larger foraminifera cited above. The first occurrence of Miogypsina gunteri in the lower Athe upper Aquitanian as proposed by European workers cannot be confirmed.

Radiolarian biochronology and Pacific origin of Jurassic-Cretaceous oceanic terranes in Costa Rica and Puerto Rico

AbstractAs demonstrated during several recent geological conferences, there is still a large debate concerning the origins of the Mesozoic oceanic remnants on the Caribbean Plate. The geodynamic models describing the Mesozoic history of the Caribbean realm can be divided into two main categories based on the origin of the Caribbean Plate: 1) An in situ origin between the Americas; 2) A Pacific origin and an eastward transport relative to the Americas. The study of the ribbon-bedded radiolarite is a key in determining the origins of associated Mesozoic oceanic terranes and may help to achieve a general agreement regarding the basic principles on the evolution of the Caribbean Plate. The Early Jurassic to early Late Cretaceous Bermeja Complex of Puerto Rico, witch contains serpentinized peridotite, altered basalt, amphibolite, and chert (Mariquita Chert Formation), and the contemporaneous Santa Rosa Accretionary Complex, which crops out in several half-windows along the south shores of the Santa Elena Peninsula in northwestern Costa Rica, are two of these little-known and crucial ophiolitic mélanges. The Manzanillo and Matambú fore-arc Terranes of the Nicoya Peninsula in the northwestern Costa Rica, which contain Late Cretaceous to Early Paleogene radiolarian-bearing siliceous mudstones and cherts associated with arc-derived mafic to intermediate volcaniclastics, bring important information on the history of the western active margin of the Caribbean Plate. A systematic radiolarian study of these three regions is presented herein in three different articles.The radiolarian biochronology of the Mariquita Chert Formation of the Bermeja Complex presented in this work indicate an early Middle Jurassic to early Late Cretaceous (late Bajocian-early Callovian to middle Albian-middle Cenomanian) age for the Mariquita Chert Formation. The illustrated assemblages contain 150 species, of which 3 are new (Pantanellium karinae, Loopus bermejaense, and L. boricus), and belonging to 59 genera. A review of the previous radiolarian published works on this formation and the results of this study suggest that the Bermeja Complex ranges in age from Middle Jurassic to early Late Cretaceous (late Aalenian to middle Cenomanian) and also reveal a possible feature of the complex, which is the youngling of radiolarian cherts from north to south, evoking a polarity of accretion. On the basis of a currently exhaustive inventory of the ribbonbedded radiolaritic facies on the Caribbean Plate, a re-examination of the distribution of Middle Jurassic sediments associated with oceanic crust from the Caribbean realm, and a paleoceanographical argumentation on the water currents, we come to the conclusion that the radiolarite and associated Mesozoic oceanic terranes of the Caribbean Plate are of Pacific origin. The previous argument for a Pacific origin of the Bermeja Complex presented by Montgomery et al. (1994a), based on their radiolarian age and their estimation of the oldest Proto-Caribbean oceanic crust, is nowadays seriously questionable, owing to the recent progresses in radiolarian biostratigraphy and new discoveries on the age of the first oceanic crust spreading between the Americas. Furthermore, we interpret the radiolarian Parvicingulidae-rich assemblages in the low-latitude Caribbean context as potential indicators of upwelling or land nutrients inputs, instead of indicators of paleolatitudes,as firstly stated by Pessagno and Blome (1986). Eventually, a discussion on the origin of the cherts of the Mariquita Formation illustrated by Middle Jurassic to middle Cretaceous geodynamic models of the Pacific and Caribbean realms bring up the possibility that the rocks of the Bermeja Complex are remnants of two different oceans.The Santa Rosa Accretionary Complex contains various oceanic assemblages of alkaline basalt, radiolarite and polymictic breccias. The radiolarian biochronology (19 illustrated assemblages, 232 species belonging to 63 genera) presented in this work indicate an Early Jurassic to early Late Cretaceous (early Pliensbachian to earliest Turonian) age for the sediments associated with oceanic basalts or recovered from blocks in breccias or megabreccias from the Santa Rosa Accretionary Complex. This study brings to light the Early Jurassic age of a sequence of ribbon-bedded radiolarite, which was previously thought to be of Cretaceous age, intruded by alkaline basalts sills. The presence of Early Jurassic large reworked blocks of radiolarite in a polymictic megabreccia, firstly reported by De Wever et al. (1985) is confirmed. Therefore, the alkaline basalt associated with these radiolarites could be of Jurassic age. In the Carrizal tectonic window, Middle Jurassic radiolarian chert blocks and Early Cretaceous brick-red ribbon-bedded radiolarites overlying pillow basalts are interpreted as fragments of a Middle Jurassic oceanic basement accreted to an Early Cretaceous oceanic plate, in an intra-oceanic subduction context. Whereas, knobby radiolarites and black shale at Playa Carrizal are indicative of a shallower middle Cretaceous paleoenvironment. Other younger oceanic remnants documented the rapid approach of the site of sedimentation to a subduction trench during the late Early Cretaceous (AlbianCenomanian), maybe early Late Cretaceous (Turonian).In total, 60 species belonging to 34 genera were present in relatively well-preserved radiolarian faunas from volcaniclastics and associated pelagic and hemipelagic rocks of the Matambú and Manzanillo terranes, ranging in age from Late Cretaceous to Early Paleogene (middle Turonian-Santonian to late Thanetian-Ypresian). This study shows that radiolarians can provide significant biostratigraphic control in the Nicoya Peninsula where very similar lithologies of different ages are present. Two radiolarian samples directly date the Berrugate Formation for the first time (middle Turonian-Santonian and Coniacian-Santonian). These ages allow to determine a volcanic arc activity on the western edge of the future Caribbean Plate at least since the Santonian that could have lasted through the middle Turonian-early Campanian interval by stratigraphic superposition. Moreover on the basis of these radiolarian ages, the Loma Chumico Formation of Albian age, and the Berrugate Formation of middle Turonian-early Maastrichtian age, can now be clearly differentiated. Two samples from the Sabana Grande Formation give a Coniacian-Santonian age and a Coniacian-Campanian age and indicate that there is a stratigraphic gap of ~10 million years between this formation and the underlying Albian Loma Chumico Formation.RésuméComme cela a pu se vérifier à plusieurs reprises lors de conférences géologiques récentes, le débat sur l'origine des terrains océaniques mésozoïques de la Plaque Caraïbes est toujours d'actualité. Les modèles géodynamiques décrivant l'histoire de la région caraïbes peuvent être classés en deux catégories basées sur l'origine de la Plaque Caraïbes : 1) Une origine in situ entre les Amériques ; 2) Une origine Pacifique et un transport vers l'est, par rapport aux Amériques. L'étude des radiolarites rubanées est capitale pour la détermination de l'origine des terrains océaniques allochtones du Mésozoïque et peut être utile pour parvenir à un compromis général concernant les principes basiques de l'évolution de la Plaque Caraïbes. Le complexe de Bermeja à Porto Rico qui est constitué de péridotites serpentinisées, de basaltes altérés, d'amphibolites et de cherts (Formation des Cherts de Mariquita), et le Complexe d'Accrétion de Santa Rosa qui affleure dans plusieurs demi-fenêtres tectoniques au sud de la Péninsule de Santa Elena au nord-ouest du Costa Rica sont deux de ces mélanges ophiolitiques peu décrits et déterminants. Les terrains de fore-arc de Manzanillo et de Matambu dans la Péninsule de Nicoya au nord-ouest du Costa Rica qui sont composés de calcaires siliceux et de cherts riches en radiolaires associés à du matériel volcanique d'arc mafique à intermédiaire, apportent d'importantes informations sur l'histoire de la marge active occidentale de la Plaque Caraïbe. Une étude systématique des radiolaires de ces trois régions est présentée dans ce travail sous forme de trois articles.La biochronologie des radiolaires de la Formation des Cherts de Mariquita du Complexe d'Accrétion de Santa Rosa présentée dans ce travail indique un âge Jurassique Moyen inférieur à Crétacé Supérieur inférieur (Bajocien supérieur-Callovien inférieur à Albien moyen-Cénomanien moyen) pour la Formation des Cherts de Mariquita. Les assemblages illustrés contiennent 150 espèces, parmis lesquelles 3 sont nouvelles (Pantanellium karinae, Loopus bermejaense et L. boricus), et appartenant à 59 genres différents. Une révision des travaux publiés précédemment sur les radiolaires de cette formation, ainsi que les résultats de cette étude suggèrent que le Complexe de Bermeja a un âge allant du Jurassique moyen au Crétacé Supérieur inférieur (Aalénien supérieur à Cénomanien moyen) et révèle aussi une caractéristique éventuelle du complexe qui est le rajeunissement des radiolarites du nord au sud, évoquant une polarité d'accrétion. Sur la base d'un inventaire actuellement exhaustif du facies radiolaritique rubané sur la Plaque Caraïbes, d'un nouvel examen de la distribution globale des sédiments du Jurassique Moyen associés à de la croûte océanique et d'une argumentation paléocéanographique sur les courants, nous arrivons à la conclusion que les radiolarites et les unités tectoniques océaniques du Mésozoïque associées de la Plaque Caraïbes sont d'origine pacifique. L'argument antérieur pour une origine pacifique du Complexe de Bermeja présenté par Montgomery et al. (1994a), basé sur leur âge à radiolaire et leur estimation de l'âge de la plus vieille croûte océanique des Proto-Caraïbes, est sérieusement remis en question aujourd'hui, en raison des progrès récents de la biostratigraphie des radiolaires et des nouvelles découvertes concernant l'âge du début de l'océanisation entre les Amériques. En outre, dans le contexte de basses latitudes des Caraïbes, nous interprétons les assemblages à radiolaires riches en Parvicingulidae comme étant des indicateurs potentiels d'apports en nutriments des zones d'uppwelling ou des terres, plutôt que des indicateurs de paléolatitudes, comme exposer pour la première fois par Pessagno et Blome (1986). Finalement, une discussion sur l'origine des cherts de la Formation de Mariquita illustrée par des modèles géodynamiques du Jurassique Moyen au Crétacé moyen des régions pacifique et caraïbes, fait poindre la possibilité que les roches du Complexe de Bermeja proviennent de deux océans différents.Le Complexe d'Accrétion de Santa Rosa contient plusieurs assemblages océaniques différents de basaltes alcalins, radiolarites et brèches polymictes. La biochronologie des radiolaires (19 assemblages illustrés, 232 espèces appartenant à 63 genres) présentée dans ce second travail indique un âge Jurassique Inférieur à Crétacé Supérieur inférieur (Pliensbachien inférieur à Turonien initial) pour les sédiments associés aux basaltes océaniques ou provenant de blocs dans des brèches ou des mégabrèches du Complexe d'Accrétion de Santa Rosa. Cette étude met en évidence l'âge Jurassique Inférieur d'une séquence de radiolarites rubanées entrecoupée de sills de basaltes alcalins, dont l'âge estimé était précédemment le Crétacé.La présence de blocs plurimétriques de radiolarites d'âge Jurassique Inférieur remaniés dans une mégabrèche polymicte, dont la présence avait été signalée par De Wever et al. (1985), est confirmée. Par conséquent, les basaltes alcalins associés à ces radiolarites pourraient aussi être d'âge Jurassique. Dans la fenêtre tectonique de Carrizal, des blocs de radiolarites d'âge Jurassique Moyen et des radiolarites du Crétacé Inférieur recouvrant des basaltes en coussins sont interprétés comme des fragments d'une croûte océanique d'âge Jurassique Moyen accrétés à une plaque océanique d'âge Crétacé Inférieur, dans un contexte de subduction intra-océanique. Alors que dans la même zone, les radiolarites « noueuses » et les argiles noires associées sont interprétées comme des indicateurs d'un milieu peu profond au Crétacé. D'autres fragments océaniques plus jeunes documentent une approche rapide du lieu de sédimentation vers une fosse de subduction pendant le Crétacé Inférieur supérieur (Albien-Cénomanien), peut-être Crétacé Supérieur (Turonien).Au total, 60 espèces appartenant à 34 genres ont été déterminées à partir de faunes à radiolaires relativement bien préservées, extraites de roches volcanoclastiques et pélagiques à hémipélagiques associées, provenant des terrains de Matambu et Manzanillo et ayant des âges compris entre le Crétacé Supérieur et le Paléogène Inférieur (Turonien moyen-Santonien à Thanétien supérieur-Yprésien). Cette étude montre que les radiolaires peuvent fournir un contrôle stratigraphique significatif dans la Péninsule de Nicoya, où des lithologies similaires, mais d'âges différents sont présentes. Deux échantillons à radiolaires permettent de dater la Formation de Berrugate pour la première fois (Turonien moyen-Santonien et Coniacien-Santonien). Ces âges permettent d'établir une activité volcanique d'arc le long de la marge occidentale de la futur Plaque Caraïbes au moins depuis le Santonien et qui pourrait avoir durée jusqu'au Turonien moyen-Campanien inférieur. De plus, sur la base de ces âges à radiolaires, la Formation de Loma Chumico d'âge Albien, et la Formation de Berrugate d'âge Turonien moyen-Maastrichtien inférieur, peuvent maintenant être différenciées. Deux échantillons de la Formation de Sabana Grande donnent des âges Coniacien-Santonien et Coniacien-Campanien et indiquent qu'il existe une lacune stratigraphique d'environ 10 millions d'années entre cette formation et la Formation de Loma Chumico sous-jacente d'âge Albien.