Redistribution of HFSE elements during rutile replacement by titanite

Titanite growth at the expense of rutile during retrograde hydration of eclogite into amphibolite is a common phenomenon. We investigated an amphibolite sample from the Tromsø eclogite facies terrain in Northern Norway to gain insight into the trace element distribution between rutile and titanite during incomplete resorption of the rutile by titanite. Patchy compositional zoning of Al, Ti, and F in titanite relates to the presence of a fluid with variable Ti/Al and/or F during its growth. Laser ablation ICP–MS and electron microprobe data for high field strength elements (HFSE: Nb, Zr, Ta, and Hf) of rutile resorbed by titanite indicate a pronounced enrichment of these elements in the rim of a large single rutile crystal (~8 mm) and a systematic decrease towards uniform HFSE contents in the large core. HFSE contents of smaller rutile grains (~0.5 mm) and rutile inclusions (<100 μm) in the titanite overgrowth are similar or higher than in the rims of large rutile crystals. Element profiles from the rim inward demonstrate that HFSE enrichment in rutile is controlled by diffusion. HFSE ratios in diffusion-altered rutile show systematic variations compared with the uniform core composition of the large rutile. Modelling of Zr and Nb diffusion in rutile indicates that diffusion coefficients in rutile in fluid-dominated natural systems must be considerably higher than those determined experimentally at 1 bar in dry systems. Variations of HFSE contents in the newly formed titanite show no systematic spatial distribution. HFSE ratios in titanite and the rims of rutile are different, indicating different solid/fluid distribution coefficients in these minerals. Element fractionation by diffusion into the relict rutile and during fluid-mediated growth of new titanite could substantially change the HFSE budget of these minerals and could affect their use for geochemical tracing and other applications, such as Zr-based geothermobarometry.     

Pre-Cenozoic intra-plate magmatism along the Central Andes (17–34°S): Composition of the mantle at an active margin

Sr–Nd–Pb isotope ratios of alkaline mafic intra-plate magmatism constrain the isotopic compositions of the lithospheric mantle along what is now the eastern foreland or back arc of the Cenozoic Central Andes (17–34°S). Most small-volume basanite volcanic rocks and alkaline intrusive rocks of Cretaceous (and rare Miocene) age were derived from a depleted lithospheric mantle source with rather uniform initial ¹⁴³Nd/¹⁴⁴Nd ( 0.5127–0.5128) and ⁸⁷Sr/⁸⁶Sr ( 0.7032–0.7040). The initial ²⁰⁶Pb/²⁰⁴Pb ratios are variable (18.5–19.7) at uniform ²⁰⁷Pb/²⁰⁴Pb ratios (15.60 ± 0.05). A variety of the Cretaceous depleted mantle source of the magmatic rocks shows elevated Sr isotope ratios up to 0.707 at constant high Nd isotope ratios. The variable Sr and Pb isotope ratios are probably due to radiogenic growth in a metasomatized lithospheric mantle, which represents the former sub-arc mantle beneath the early Palaeozoic active continental margin. Sr–Nd–Pb isotope signatures of a second mantle type reflected in the composition of Cretaceous (one late Palaeozoic age) intra-plate magmatic rocks (¹⁴³Nd/¹⁴⁴Nd  0.5123, ⁸⁷Sr/⁸⁶Sr  0.704, ²⁰⁶Pb/²⁰⁴Pb  17.5–18.5, and ²⁰⁷Pb/²⁰⁴Pb  15.45–15.50) are similar to the isotopic composition of old sub-continental lithospheric mantle of the Brazilian Shield.

Published Nd and Sr isotopic compositions of Mesozoic to Cenozoic arc-related magmatic rocks (18–40°S) represent the composition of the convective sub-arc mantle in the Central Andes and are similar to those of the Cretaceous (and rare Miocene) intra-plate magmatic rocks. The dominant convective and lithospheric mantle type beneath this old continental margin is depleted mantle, which is compositionally different from average MORB-type depleted mantle. The old sub-continental lithospheric mantle did not contribute to Mesozoic to Cenozoic arc magmatism.     

The LMDZ4 general circulation model: climate performance and sensitivity to parametrized physics with emphasis on tropical convection

The LMDZ4 general circulation model is the atmospheric component of the IPSL–CM4 coupled model which has been used to perform climate change simulations for the 4th IPCC assessment report. The main aspects of the model climatology (forced by observed sea surface temperature) are documented here, as well as the major improvements with respect to the previous versions, which mainly come form the parametrization of tropical convection. A methodology is proposed to help analyse the sensitivity of the tropical Hadley–Walker circulation to the parametrization of cumulus convection and clouds. The tropical circulation is characterized using scalar potentials associated with the horizontal wind and horizontal transport of geopotential (the Laplacian of which is proportional to the total vertical momentum in the atmospheric column). The effect of parametrized physics is analysed in a regime sorted framework using the vertical velocity at 500 hPa as a proxy for large scale vertical motion. Compared to Tiedtke’s convection scheme, used in previous versions, the Emanuel’s scheme improves the representation of the Hadley–Walker circulation, with a relatively stronger and deeper large scale vertical ascent over tropical continents, and suppresses the marked patterns of concentrated rainfall over oceans. Thanks to the regime sorted analyses, these differences are attributed to intrinsic differences in the vertical distribution of convective heating, and to the lack of self-inhibition by precipitating downdraughts in Tiedtke’s parametrization. Both the convection and cloud schemes are shown to control the relative importance of large scale convection over land and ocean, an important point for the behaviour of the coupled model.     

Synchronous pattern of moisture availability on the southern Tibetan Plateau since 17.5 cal. ka BP – the Tangra Yumco lake sediment record

A possible asynchronicity of the spatial and temporal moisture availability on the Tibetan Plateau has been a controversial subject of discussion in recent years. Here we present the first attempt to systematically investigate possible spatial and temporal variations in moisture availability by examining two lakes, Tangra Yumco and Nam Co, on an east–west transect on the southern Tibetan Plateau using identical proxies for palaeoenvironmental reconstruction. In this study, an independent record from Tangra Yumco was analysed applying a multi‐proxy approach to reconstruct variations in moisture availability since the Lateglacial. Results were subsequently compared with previously published records from Nam Co and additional records from Tso Moriri (northwestern Himalaya) and Naleng Co (southeastern Tibetan Plateau). Our results show that Tangra Yumco was at least partially ice covered prior to 17.1 cal. ka BP. A temperature rise after 17.1 cal. ka BP probably resulted in thawing of the permafrost. At 16.0 cal. ka BP moisture availability increased, representing an initial monsoonal intensification. Warmer conditions between 13.0 and 12.4 cal. ka BP and cooler conditions between 12.4 cal. ka BP and the onset of the Holocene reflect the Bølling‐Allerød and Younger Dryas. At the onset of the Holocene moisture availability rapidly increased, with moisture highest prior to 8.5 cal. ka BP when temperatures were also highest. After 8.5 cal. ka BP the moisture availability gradually decreased and showed only minor amplitude variations. These findings are consistent with the records from large lakes like Nam Co, Tso Moriri, and Naleng Co, revealing a synchronous pattern of moisture availability on the southern Tibetan Plateau.     

Shallow sediment deformation styles in north-western Campeche Knolls,Gulf of Mexico and their controls on the occurrence of hydrocarbon seepage

During German R/V Meteor M67/2 expedition to Campeche Knolls, southern Gulf of Mexico, a set of 2D high resolution seismic data was acquired to study the near-surface sediment structure and its relationship with hydrocarbon seepages in this salt province. The comprehensive survey covered 20 individual bathymetric highs or ridges and identified three principle structural types: Passive Type, Chaopopte Type, and Asymmetric Flap Type. The first type is the result of passive diapirism, whereas the latter two were initialized by a regional compressional event in the Miocene, but are later differently modified by salt tectonism. Chapopote Type structures appear as symmetrical domes, with uplifted coarse-grained Miocene sediments in their cores and rather thin syn-kinematic sediments covering the crests. Asymmetric Flap Type structures are also first folded as domes or ridges, but one flap later subsided together with salt evacuation, resulting in single uplifted monoclines. With the coarse-grained pre-kinematic sediments as reservoir units, both structural types can focus and accumulate hydrocarbons. The geometries of the structures suggest that hydrocarbons are accumulated in the center of the Chapopote Type structures and in the subsided flaps of the Asymmetric Flap Type structures. Hydrocarbon leakage from these thinly sealed reservoirs is regarded as the principle mechanism for the seepage in the study area, and accordingly the most seepage-prone positions are above these reservoirs. The seep locations suggested by analysis of sea-surface oil slick images of SAR satellite data are also examined in this study. These independently derived seep locations confirm the seepage-prone positions to be above the shallow buried reservoirs. This study suggest that the shallow sediment structures control the distribution of the hydrocarbon seeps of the north-western Campeche Knolls, although the hydrocarbons are sourced from the greater depth.     

Effects of variable wind shear on the mesoscale circulation forced by slab‐symmetric diabatic heating

Abstract

We examine the response of stably stratified airflow to a slab‐symmetric diabatic forcing associated with condensation in long‐lasting precipitation bands. The steady‐state linearized Boussinesq equations are used to model the diagnostic relationship between the vertical motion field, the heating source and the ambient flow. The basic‐state flow is assumed to be horizontally uniform and non‐rotating, but the static stability and wind vary in the vertical. Linear theory shows that the speed of the along‐band wind component is unimportant for slab‐symmetric heating since it cannot contribute towards the advection of buoyancy or vertical motion.

For typical atmospheric stratification and a moving heating source associated with a cloud band, the Taylor‐Goldstein equation is solved numerically. The numerical results show that the cross‐band wind shear tilts the updraft core and broadens it. While the magnitude of the shear is increased, the circulation becomes stronger. The details of the wind profile are also important in determining the intensity and structure of the circulation. When the wind profile indicates a convex bulge (i.e. the low‐level shear is weaker than the upper‐level shear), the circulation becomes slightly weaker in comparison with the linear wind profile. Conversely, the circulation becomes stronger when the wind profile has a concave shape. Increasing the concave bulge tends to enhance the circulation but not in a monotonic fashion. This non‐monotonic relation between the vertical motion and the parabolic wind profile is interpreted in terms of kinetic energy changes of parcels that interchange their altitudes.     

In situ LA-ICPMS Isotopic and Geochronological Studies on Carbonatites and Phoscorites from the Guli Massif,Maymecha-Kotuy,Polar Siberia

In this study we present a fresh isotopic data, as well as U–Pb ages from different REE-minerals in carbonatites and phoscorites of Guli massif using in situ LA-ICPMS technique. The analyses were conducted on apatites and perovskites from calcio-carbonatite and phoscorite units, as well as on pyrochlores and baddeleyites from the carbonatites. The ⁸⁷Sr/⁸⁶Sr ratios obtained from apatites and perovskites from the phoscorites are 0.70308–0.70314 and 0.70306–0.70313, respectively; and 0.70310–0.70325 and 0.70314–0.70327, for the pyrochlores and apatites from the carbonatites, respectively.Furthermore, the in situ laser ablation analyses of apatites and perovskites from the phoscorite yield εNd from 3.6 (±1) to 5.1 (±0.5) and from 3.8 (±0.5) to 4.9 (±0.5), respectively; εNd of apatites, perovskites and pyrochlores from carbonatite ranges from 3.2 (±0.7) to 4.9 (±0.9), 3.9 (±0.6) to 4.5 (±0.8) and 3.2 (±0.4) to 4.4 (±0.8), respectively. Laser ablation analyses of baddeleyites yielded an eHf(t)d of +8.5 (± 0.18); prior to this study Hf isotopic characteristic of Guli massif was not known. Our new in situ εNd, ⁸⁷Sr/⁸⁶Sr and eHf data on minerals in the Guli carbonatites imply a depleted source with a long time integrated high Lu/Hf, Sm/Nd, Sr/Rb ratios.In situ U–Pb age determination was performed on perovskites from the carbonatites and phoscorites and also on pyrochlores and baddeleyites from carbonatites. The co-existing pyrochlores, perovskites and baddeleyites in carbonatites yielded ages of 252.3 ± 1.9, 252.5 ± 1.5 and 250.8 ± 1.4 Ma, respectively. The perovskites from the phoscorites yielded an age of 253.8 ± 1.9 Ma. The obtained age for Guli carbonatites and phoscorites lies within the range of ages previously reported for the Siberian Flood Basalts and suggest essentially synchronous emplacement with the Permian-Triassic boundary.     

Zoning and exsolution in alkali feldspars from Laacher See volcano (Western Germany): constraints on temperature history prior to eruption

Compositional zoning and exsolution patterns of alkali feldspars in carbonatite-bearing cognate syenites from the 6.3 km³ (D.R.E) phonolitic Laacher See Tephra (LST) deposit in western Germany (12.9 ka) are reported. These rocks represent the cooler outer portion and crystal-rich products of a cooling magma reservoir at upper crustal levels. Major and trace-element difference between cores and rims in sanidine crystals represent two generations of crystal growth separated by unmixing of a carbonate melt. Trace-element differences measured by LA–ICP–MS are in accordance with silicate–carbonate unmixing. Across the core–rim boundary, we extracted gray-scale profiles from multiple accumulations of back-scattered electron images. Gray scales directly represent K/Na ratios owing to low concentrations of Ba and Sr (<?30 ppm). Diffusion gradients are modeled to solve for temperature using known pre-eruptive U–Th zircon ages (0–20 ky) of each sample (Schmitt et al., J Petrol 51:1053–1085, 2010). Estimated temperatures range from 630 °C to 670 °C. For the exsolution boundaries, a diffusive homogenization model is constrained by the solvus temperature of ~ 712^_725 °C and gives short time scales of only 15–50 days. Based on our results, we present a model for the temperature–time history of these rocks. The model also constrains the thermal variation across the cooling crystal-rich carapace of the magma reservoir over 20 ka and suggests a thermal reactivation of cumulates, the cooling carapace, and probably the entire system only a few years prior to the explosive eruption of the remaining molten core of the phonolitic magma reservoir.