Petrological evidence for stepwise accretion of metamorphic soles during subduction infancy (Semail ophiolite,Oman and UAE)

Metamorphic soles are tectonic slices welded beneath most large‐scale ophiolites. These slivers of oceanic crust metamorphosed up to granulite facies conditions are interpreted as forming during the first million years of intraoceanic subduction following heat transfer from the incipient mantle wedge towards the top of the subducting plate. This study reappraises the formation of metamorphic soles through detailed field and petrological work on three key sections from the Semail ophiolite (Oman and United Arab Emirates). Based on thermobarometry and thermodynamic modelling, it is shown that metamorphic soles do not record a continuous temperature gradient, as expected from simple heating by the upper plate or by shear heating as proposed in previous studies. The upper, high‐T metamorphic sole is subdivided in at least two units, testifying to the stepwise formation, detachment and accretion of successive slices from the down‐going slab to the mylonitic base of the ophiolite. Estimated peak pressure–temperature conditions through the metamorphic sole, from top to bottom, are 850°C and 1 GPa, 725°C and 0.8 GPa and 530°C and 0.5 GPa. These estimates appear constant within each unit but differing between units by 100–200°C and ~0.2 GPa. Despite being separated by hundreds of kilometres below the Semail ophiolite and having contrasting locations with respect to the ridge axis position, metamorphic soles show no evidence for significant petrological variations along strike. These constraints allow us to refine the tectonic–petrological model for the genesis of metamorphic soles, formed via the stepwise stacking of several homogeneous slivers of oceanic crust and its sedimentary cover. Metamorphic soles result not so much from downward heat transfer (ironing effect) as from progressive metamorphism during strain localization and cooling of the plate interface. The successive thrusts originate from rheological contrasts between the sole, initially the top of the subducting slab, and the peridotite above as the plate interface progressively cools. These findings have implications for the thickness, the scale and the coupling state at the plate interface during the early history of subduction/obduction systems.     

Atomic oxygen on the Venus nightside: Global distribution deduced from airglow mapping

The Visible and Infra-Red Thermal Imaging Spectrometer (VIRTIS) instrument on board the Venus Express spacecraft has measured the O₂(a¹Δ) nightglow distribution at 1.27 μm in the Venus mesosphere for more than two years. Nadir observations have been used to create a statistical map of the emission on Venus nightside. It appears that the statistical 1.6 MR maximum of the emission is located around the antisolar point. Limb observations provide information on the altitude and on the shape of the emission layer. We combine nadir observations essentially covering the southern hemisphere, corrected for the thermal emission of the lower atmosphere, with limb profiles of the northern hemisphere to generate a global map of the Venus nightside emission at 1.27 μm. Given all the O₂(a¹Δ) intensity profiles, O₂(a¹Δ) and O density profiles have been calculated and three-dimensional maps of metastable molecular and atomic oxygen densities have been generated. This global O density nightside distribution improves that available from the VTS3 model, which was based on measurements made above 145 km. The O₂(a¹Δ) hemispheric average density is 2.1 × 10⁹ cm^?3, with a maximum value of 6.5 × 10⁹ cm^?3 at 99.2 km. The O density profiles have been derived from the nightglow data using CO₂ profiles from the empirical VTS3 model or from SPICAV stellar occultations. The O hemispheric average density is 1.9 × 10¹¹ cm^?3 in both cases, with a mean altitude of the peak located at 106.1 km and 103.4 km, respectively. These results tend to confirm the modeled values of 2.8 × 10¹¹ cm^?3 at 104 km and 2.0 × 10¹¹ cm^?3 at 110 km obtained by Brecht et al. [Brecht, A., Bougher, S.W., Gérard, J.-C., Parkinson, C.D., Rafkin, S., Foster, B., 2011a. J. Geophys. Res., in press] and Krasnopolsky [Krasnopolsky, V.A., 2010. Icarus 207, 17–27], respectively. Comparing the oxygen density map derived from the O₂(a¹Δ) nightglow observations, it appears that the morphology is very different and that the densities obtained in this study are about three times higher than those predicted by the VTS3 model.     

Yearly evolution of Euro-Atlantic weather regimes and of their sub-seasonal predictability

It is often assumed that weather regimes adequately characterize atmospheric circulation variability. However, regime classifications spanning many months and with a low number of regimes may not satisfy this assumption. The first aim of this study is to test such hypothesis for the Euro-Atlantic region. The second one is to extend the assessment of sub-seasonal forecast skill in predicting the frequencies of occurrence of the regimes beyond the winter season. Two regime classifications of four regimes each were obtained from sea level pressure anomalies clustered from October to March and from April to September respectively. Their spatial patterns were compared with those representing the annual cycle. Results highlight that the two regime classifications are able to reproduce most part of the patterns of the annual cycle, except during the transition weeks between the two periods, when patterns of the annual cycle resembling Atlantic Low regime are not also observed in any of the two classifications. Forecast skill of Atlantic Low was found to be similar to that of NAO+, the regime replacing Atlantic Low in the two classifications. Thus, although clustering yearly circulation data in two periods of 6 months each introduces a few deviations from the annual cycle of the regime patterns, it does not negatively affect sub-seasonal forecast skill. Beyond the winter season and the first ten forecast days, sub-seasonal forecasts of ECMWF are still able to achieve weekly frequency correlations of r = 0.5 for some regimes and start dates, including summer ones. ECMWF forecasts beat climatological forecasts in case of long-lasting regime events, and when measured by the fair continuous ranked probability skill score, but not when measured by the Brier skill score. Thus, more efforts have to be done yet in order to achieve minimum skill necessary to develop forecast products based on weather regimes outside winter season.

     

Mesozoic to Cenozoic tectono‐metamorphic history of the South Pamir–Hindu Kush (Chitral,NW Pakistan): Insights from phase equilibria modelling,and garnet–monazite petrochronology

The Karakoram–Hindu Kush–Pamir and adjacent Tibetan plateau belt comprise a series of Gondwana‐derived crustal fragments that successively accreted to the Eurasian margin in the Mesozoic as the result of the progressive Tethys ocean closure. These domains provide unique insights into the thermal and structural history of the Mesozoic to Cenozoic Eurasian plate margin, which are critical to inform the initial boundary conditions (e.g. crustal thickness, structure and thermo‐mechanical properties) for the subsequent development of the large and hot Tibetan–Himalaya orogen, and the associated crustal deformation processes. Using a combination of microstructural analyses, thermobarometry modelling and U–Th–Pb monazite and Lu–Hf garnet geochronology, the study reappraises the metamorphic history of exposed mid‐crustal metapelites in the Chitral region of the South Pamir–Hindu Kush (NW Pakistan). This study also demonstrates that trace elements in monazite (especially Y and Dy), combined with thermodynamical modelling and Lu–Hf garnet dating, provides a powerful integrated toolbox for constraining long‐lived and polyphased tectono‐metamorphic histories in all their spatial and temporal complexity. Rocks from the Chitral region were progressively deformed and metamorphosed at sub‐ and supra‐solidus conditions through at least four distinct episodes from the Mesozoic to the Cenozoic. Rocks were first metamorphosed at ~400–500°C and ~0.3 GPa in the Late Triassic–Early Jurassic (210–185 Ma), likely in response to the accretion of the Karakoram during the Cimmerian orogeny. Pressure and temperature subsequently increased by ~0.3 GPa and 100°C in the Early‐ to Mid Cretaceous (140–80 Ma), coinciding with the intrusion of calcalkaline granitic plutons across the Karakoram and Pamir regions. This event is interpreted as the record of crustal thickening and the development of a proto‐plateau within the Eurasian margin due to a long‐lived episode of slab flattening in an Andean‐type margin. Peak metamorphism was reached in the Late Eocene–Early Oligocene (40–30 Ma) at conditions of 580–600°C and ~0.6 GPa and 700–750°C and 0.7–0.8 GPa for the investigated staurolite schists and sillimanite migmatites respectively. This crustal heating up to moderate anatexis likely resulted in the underthrusting of the Indian plate after a NeoTethyan slab‐break off or to the Tethyan Himalaya–Lhasa microcontinent collision and subsequent oceanic slab flattening. Near‐isothermal decompression/exhumation followed in the Late Oligocene (28–23 Ma) as marked by a pressure decrease in excess of ~0.1 GPa. This event was coeval with the intrusion of the 24 Ma Garam Chasma leucogranite. This rapid exhumation is interpreted to be related to the reactivation of the South Pamir–Karakoram suture zone during the ongoing collision with India. The findings of this study confirm that significant crustal shortening and thickening of the south Eurasian margin occurred during the Mesozoic in an accretionary‐type tectonic setting through successive episodes of terrane accretions and probably slab flattening, transiently increasing the coupling at the plate interface. Moreover, they indicate that the south Eurasian margin was already hot and thickened prior to Cenozoic collision with India, which has important implications for orogen‐scale strain‐accommodation mechanisms.     

Timing of metamorphism and deformation in the Swat valley,northern Pakistan:Insight into garnet-monazite HREE partitioning

New metamorphic petrology and geochronology from the Loe Sar dome in the Swat region of northern Pakistan place refined constraints on the pressure, temperature and timing of metamorphism and deformation in that part of the Himalayan orogen. Thermodynamic modelling and monazite petrochronology indicate that metamorphism in the area followed a prograde evolution from ～525 ± 25 ℃and 6 士 0.5 kbar to ～610 ± 25 ℃ and 9 士 0.5 kbar, between ca. 39 Ma and 28 Ma. Partitioning of heavy rare earth elements between garnet rims and 30-28 Ma monazite are interpreted to indicate coeval crystallization at peak conditions. Microtextural relationships indicate that garnet rim growth post-dated the development of the main foliation in the area. The regional foliation is folded about large-scale N-S trending fold axes and overprinting E-W trending folds to form km-scale domal culminations. The textural relationships observed indicate that final dome development must be younger than the 30-28 Ma monazite that grew with garnet rims post-regional foliation development, but pre-doming-related deformation. This new timing constraint helps resolve discrepancy between previous interpretations,which have alternately suggested that N-S trending regional folds must be either pre-or post-early Oligocene. Finally, when combined with existing hornblende and white mica cooling ages, these new data indicate that the study area was exhumed rapidly following peak metamorphism.     

Spatial correlation of OH Meinel and O2 infrared atmospheric nightglow emissions observed with VIRTIS-M on board Venus Express

We present the two-dimensional distribution of the O₂ a¹Δ–X³Σ (0–0) band at 1.27 μm and the OH Δv = 1 Meinel airglow measured simultaneously with the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) on board Venus Express. We show that the two emissions present very similar spatial structures. A cross-correlation analysis indicates that the highest level of correlation is reached with only very small relative shifts of the pairs of images. In spite of the strong spatial correlation between the morphology of the bright spots in the two emissions, we also show that their relative intensity is not constant, in agreement with earlier statistical studies of their limb profiles. We conclude that the two emissions have a common precursor that controls the production of both excited species. We argue that atomic oxygen, which produces O₂ (¹Δ) molecules by three-body recombination and is the precursor of ozone formation, also governs to a large extent the OH airglow morphology through the H + O₃ → OH^* + O₂ reaction.     

Metamorphic sole formation,emplacement and blueschist facies overprint: early subduction dynamics witnessed by western Turkey ophiolites

The largest ophiolite on Earth, in western Turkey, is a key place to study obduction and early subduction dynamics. Ophiolite remnants derived from the same Neotethyan branch (preserved as a result of long‐lived Late Cretaceous continental subduction and later obduction) are underlain by hundred‐metre‐thick extensive metamorphic soles. These soles formed synchronously, at c. 93 Ma, and were welded to the base of the ophiolite, thereby dating the start of intra‐oceanic subduction. This contribution focuses on the structure, petrology and pressure–temperature evolution of the soles and other subduction‐derived units. Peak pressure–temperature conditions were estimated at 10.5 ± 2 kbar and 800 ± 50 °C for the sole by means of pseudosection calculations using Theriak/Domino and at 12 kbar and 425 °C for the unique, enigmatic blueschist facies overprint of the sole. This study provides clues to the mechanisms of sole underplating during early subduction, later cooling, and the nature of the western Turkey ophiolite.     

Atomic oxygen distributions in the Venus thermosphere: Comparisons between Venus Express observations and global model simulations

Nightglow emissions provide insight into the global thermospheric circulation, specifically in the transition region (～70–120 km). The O₂ IR nightglow statistical map created from Venus Express (VEx) Visible and InfraRed Thermal Imaging Spectrometer (VIRTIS) observations has been used to deduce a three-dimensional atomic oxygen density map. In this study, the National Center of Atmospheric Research (NCAR) Venus Thermospheric General Circulation Model (VTGCM) is utilized to provide a self-consistent global view of the atomic oxygen density distribution. More specifically, the VTGCM reproduces a 2D nightside atomic oxygen density map and vertical profiles across the nightside, which are compared to the VEx atomic oxygen density map. Both the simulated map and vertical profiles are in close agreement with VEx observations within a ～30° contour of the anti-solar point. The quality of agreement decreases past ～30°. This discrepancy implies the employment of Rayleigh friction within the VTGCM may be an over-simplification for representing wave drag effects on the local time variation of global winds. Nevertheless, the simulated atomic oxygen vertical profiles are comparable with the VEx profiles above 90 km, which is consistent with similar O₂ (¹Δ) IR nightglow intensities. The VTGCM simulations demonstrate the importance of low altitude trace species as a loss for atomic oxygen below 95 km. The agreement between simulations and observations provides confidence in the validity of the simulated mean global thermospheric circulation pattern in the lower thermosphere.