Retrieving the paleoclimatic signal from the deeper part of the EPICA Dome C ice core

An important share of paleoclimatic information is buried within the lowermost layers of deep ice cores. Because improving our records further back in time is one of the main challenges in the near future, it is essential to judge how deep these records remain unaltered, since the proximity of the bedrock is likely to interfere both with the recorded temporal sequence and the ice properties. In this paper, we present a multiparametric study (δD-δ18Oice, δ18Oatm, total air content, CO2, CH4, N2O, dust, high-resolution chemistry, ice texture) of the bottom 60 m of the EPICA (European Project for Ice Coring in Antarctica) Dome C ice core from central Antarctica. These bottom layers were subdivided into two distinct facies: the lower 12 m showing visible solid inclusions (basal dispersed ice facies) and the upper 48 m, which we will refer to as the "basal clean ice facies". Some of the data are consistent with a pristine paleoclimatic signal, others show clear anomalies. It is demonstrated that neither large-scale bottom refreezing of subglacial water, nor mixing (be it internal or with a local basal end term from a previous/initial ice sheet configuration) can explain the observed bottom-ice properties. We focus on the high-resolution chemical profiles and on the available remote sensing data on the subglacial topography of the site to propose a mechanism by which relative stretching of the bottom-ice sheet layers is made possible, due to the progressively confining effect of subglacial valley sides. This stress field change, combined with bottom-ice temperature close to the pressure melting point, induces accelerated migration recrystallization, which results in spatial chemical sorting of the impurities, depending on their state (dissolved vs. solid) and if they are involved or not in salt formation. This chemical sorting effect is responsible for the progressive build-up of the visible solid aggregates that therefore mainly originate "from within", and not from incorporation processes of debris from the ice sheet's substrate. We further discuss how the proposed mechanism is compatible with the other ice properties described. We conclude that the paleoclimatic signal is only marginally affected in terms of global ice properties at the bottom of EPICA Dome C, but that the timescale was considerably distorted by mechanical stretching of MIS20 due to the increasing influence of the subglacial topography, a process that might have started well above the bottom ice. A clear paleoclimatic signal can therefore not be inferred from the deeper part of the EPICA Dome C ice core. Our work suggests that the existence of a flat monotonic ice–bedrock interface, extending for several times the ice thickness, would be a crucial factor in choosing a future "oldest ice" drilling location in Antarctica.

Fluid-related inclusions in Alpine high-pressure peridotite reveal trace element recycling during subduction-zone dehydration of serpentinized mantle (Cima di Gagnone, Swiss Alps).

Serpentinites release at sub-arc depths volatiles and several fluid-mobile trace elements found in arc magmas. Constraining element uptake in these rocks and defining the trace element composition of fluids released upon serpentinite dehydration can improve our understanding of mass transfer across subduction zones and to volcanic arcs. The eclogite-facies garnet metaperidotite and chlorite harzburgite bodies embedded in paragneiss of the subduction melange from Cima di Gagnone derive from serpentinized peridotite protoliths and are unique examples of ultramafic rocks that experienced subduction metasomatism and devolatilization. In these rocks, metamorphic olivine and garnet trap polyphase inclusions representing the fluid released during high-pressure breakdown of antigorite and chlorite. Combining major element mapping and laser-ablation ICP-MS bulk inclusion analysis, we characterize the mineral content of polyphase inclusions and quantify the fluid composition. Silicates, Cl-bearing phases, sulphides, carbonates, and oxides document post-entrapment mineral growth in the inclusions starting immediately after fluid entrapment. Compositional data reveal the presence of two different fluid types. The first (type A) records a fluid prominently enriched in fluid-mobile elements, with Cl, Cs, Pb, As, Sb concentrations up to 10(3) PM (primitive mantle), similar to 10(2) PM Tit Ba, while Rb, B, Sr, Li, U concentrations are of the order of 10(1) PM, and alkalis are similar to 2 PM. The second fluid (type B) has considerably lower fluid-mobile element enrichments, but its enrichment patterns are comparable to type A fluid. Our data reveal multistage fluid uptake in these peridotite bodies, including selective element enrichment during seafloor alteration, followed by fluid-rock interaction along with subduction metamorphism in the plate interface melange. Here, infiltration of sediment-equilibrated fluid produced significant enrichment of the serpentinites in As, Sb, B, Pb, an enriched trace element pattern that was then transferred to the fluid released at greater depth upon serpentine dehydration (type A fluid). The type B fluid hosted by garnet may record the composition of the chlorite breakdown fluid released at even greater depth. The Gagnone study-case demonstrates that serpentinized peridotites acquire water and fluid-mobile elements during ocean floor hydration and through exchange with sediment-equilibrated fluids in the early subduction stages. Subsequent antigorite devolatilization at subarc depths delivers aqueous fluids to the mantle wedge that can be prominently enriched in sediment-derived components, potentially triggering arc magmatism without the need of concomitant dehydration/melting of metasediments or altered oceanic crust.

Apatite as an indicator of fluid salinity: An experimental study of chlorine and fluorine partitioning in subducted sediments

In order to constrain the salinity of subduction zone fluids, piston-cylinder experiments have been conducted to investigate the partitioning behaviour of Cl and F in subducted sediments. These experiments were performed at H2O-undersaturated conditions with a synthetic pelite starting composition containing 800 ppm Cl, over a pressure and temperature range of 2.5–4.5 GPa and 630–900 °C. Repetitive experiments were conducted with 1900 ppm Cl + 1000 ppm F, and 2100 ppm Cl. Apatite represents the most Cl-abundant mineral phase, with Cl concentration varying in the range 0.1–2.82 wt%. Affinity for Cl decreases over the following sequence: aqueous fluid > apatite ⩾ melt > other hydrous minerals (phengite, biotite and amphibole). It was found that addition of F to the Cl-bearing starting composition significantly lowers the Cl partition coefficients between apatite and melt (DClAp–melt) and apatite and aqueous fluid (DClAp–aq). Cl–OH exchange coefficients between apatite and melt (KdCl–OHAp–melt) and apatite and aqueous fluid (KdCl–OHAp–aq) were subsequently calculated. KdCl–OHAp–melt was found to vary from 1 to 58, showing an increase with temperature and a decrease with pressure and displaying a regular decrease with increasing H2O content in melt. Mole fractions of Cl and OH in melt were calculated based on an ideal mixing model for H2O, OH, O, Cl and F. The Cl contents of other hydrous minerals (phengite, biotite and amphibole) fall between 200 and 800 ppm, with resultant Cl partition coefficients from 0.02 to 0.49, appearing independent of the bulk Cl and F content. Preliminary data from this study show that the partitioning behaviour of F is strongly in favour of apatite relative to melt and phengite, with DFAp–melt = 15–51. Apatites from representative eclogite facies metasediments were examined and found to have low Cl contents close to ∼100 ppm. Calculations using our experimentally determined KdCl–OHAp–aq of 0.004 at 2.5 GPa, 630 °C indicate a low salinity character (0.5–2 wt% NaCleq) for the fluid formed during dehydration of subducted oceanic sediment at ∼80 km depth.

Contrasting PT paths within the Barchi-Kol UHP terrain (Kokchetav Complex): Implications for subduction and exhumation of continental crust

The Barchi-Kol terrain is a classic locality of ultrahigh-pressure (UHP) metamorphism within the Kokchetav metamorphic belt. We provide a detailed and systematic characterization of four metasedimentary samples using dominant mineral assemblages, mineral inclusions in zircon and monazite, garnet zonation with respect to major and trace elements, and Zr-in-rutile and Ti-in-zircon temperatures. A typical diamond-bearing gneiss records peak conditions of 49 ± 4 kbar and 950–1000 °C. Near isothermal decompression of this rock resulted in the breakdown of phengite associated with a pervasive recrystallization of the rock. The same terrain also contains mica schists that experienced peak conditions close to those of the diamond-bearing rocks, but they were exhumed along a cooler path where phengite remained stable. In these rocks, major and trace element zoning in garnet has been completely equilibrated. A layered gneiss was metamorphosed at UHP conditions in the coesite field, but did not reach diamond-facies conditions (peak conditions: 30 kbar and 800–900 °C). In this sample, garnet records retrograde zonation in major elements and also retains prograde zoning in trace elements. A garnet-kyanite-micaschist that reached significantly lower pressures (24 ± 2 kbar, 710 ± 20 °C) contains garnet with major and trace element zoning. The diverse garnet zoning in samples that experienced different metamorphic conditions allows to establish that diffusional equilibration of rare earth element in garnet likely occurs at ~900–950 °C. Different metamorphic conditions in the four investigated samples are also documented in zircon trace element zonation and mineral inclusions in zircon and monazite. U-Pb geochronology of metamorphic zircon and monazite domains demonstrates that prograde (528–521 Ma), peak (528–522 Ma), and peak to retrograde metamorphism (503–532 Ma) occurred over a relatively short time interval that is indistinguishable from metamorphism of other UHP rocks within the Kokchetav metamorphic belt. Therefore, the assembly of rocks with contrasting P-T trajectories must have occurred in a single subduction-exhumation cycle, providing a snapshot of the thermal structure of a subducted continental margin prior to collision. The rocks were initially buried along a low geothermal gradient. At 20–25 kbar they underwent near isobaric heating of 200 °C, which was followed by continued burial along a low geothermal gradient. Such a step-wise geotherm is in good agreement with predictions from subduction zone thermal models.

Constraints on the thermal evolution of the Adriatic margin during Jurassic continental break-up: U–Pb dating of rutile from the Ivrea–Verbano Zone, Italy

The Ivrea–Verbano Zone (IVZ), northern Italy, exposes an attenuated section through the Permian lower crust that records high-temperature metamorphism under lower crustal conditions and a protracted history of extension and exhumation associated partly with the Jurassic opening of the Alpine Tethys ocean. This study presents SHRIMP U–Pb geochronology of rutile from seven granulite facies metapelites from the base of the IVZ, collected from locations spanning ~35 km along the strike of Paleozoic fabrics. Rutile crystallised during Permian high-temperature metamorphism and anatexis, yet all samples give Jurassic rutile U–Pb ages that record cooling through 650–550 °C. Rutile age distributions are dominated by a peak at ~160 Ma, with a subordinate peak at ~175 Ma. Both ~160 and ~175 Ma age populations show excellent agreement between samples, indicating that the two distinctive cooling stages they record were synchronous on a regional scale. The ~175 Ma population is interpreted to record cooling in the footwall of rift-related faults and shear zones, for which widespread activity in the Lower Jurassic has been documented along the western margin of the Adriatic plate. The ~160 Ma age population postdates the activity of all known rift-related structures within the Adriatic margin, but coincides with extensive gabbroic magmatism and exhumation of sub-continental mantle to the floor of the Alpine Tethys, west of the Ivrea Zone. We propose that this ~160 Ma early post-rift age population records regional cooling following episodic heating of the distal Adriatic margin, likely related to extreme lithospheric thinning and associated advection of the asthenosphere to shallow levels. The partial preservation of the ~175 Ma age cluster suggests that the post-rift (~160 Ma) heating pulse was of short duration. The regional consistency of the data presented here, which is in contrast to many other thermochronometers in the IVZ, demonstrates the value of the rutile U–Pb technique for probing the thermal evolution of high-grade metamorphic terrains. In the IVZ, a significant decoupling between Zr-in-rutile temperatures and U–Pb ages of rutile is observed, with the two systems recording events ~120 Ma apart.

A systematic evaluation of the Zr-in-rutile thermometer in ultra-high temperature (UHT) rocks

The Zr-in-rutile geothermometer is potentially a widely applicable tool to estimate peak metamorphic temperatures in rocks from diverse geological settings. In order to evaluate its usefulness and reliability to record and preserve high temperatures in granulite facies rocks, rutile from UHT rocks was investigated to assess different mechanisms of Zr (re-)distribution following cooling from high temperature. Granulite facies paragneisses from the lowermost part of the Ivrea Zone, Italy, incorporated as thin sheets into the extensive basaltic body of the Mafic Complex were selected for this study. The results show that Zr-in-rutile thermometry, if properly applied, is well suited to identify and study UHT terranes as it preserves a record of temperatures up to 1190 °C, although the thermometer is susceptible to partial post-peak metamorphic resetting by Zr diffusion. Texturally homogeneous rutile grains preserve Zr concentrations corresponding to temperatures of prograde rutile growth. Diverse rutile textures and relationships between some rutile host grains and included or adjacent Zr-bearing phases bear testimony to varying mechanisms of partial redistribution and resetting of Zr in rutile during cooling and link Zr-in-rutile temperatures to different steps of the metamorphic evolution. Rutile grains that equilibrated their Zr concentrations at temperatures above 1070 °C (i.e. 1.1 wt% Zr) could not retain all Zr in the rutile structure during cooling and exsolved baddeleyite (ZrO2). By subsequent reaction of baddeleyite exsolution lamellae with SiO2, zircon needles formed before the system finally closed at 650–700 °C without significant net loss of Zr from the whole host rutile grain. By reintegration of zircon exsolution needles, peak metamorphic temperatures of up to 1190 °C are derived for the studied rocks, which demonstrates the suitability of this solution thermometer to record UHT conditions and also confirms the extraordinary geological setting of the lowermost part of the Ivrea Zone.