Continental subduction and exhumation of UHP rocks. Structural and geochronological insights from the Dabieshan (East China)

In the Dabieshan, the available models for exhumation of ultrahigh-pressure (UHP) rocks are poorly constrained by structural data. A comprehensive structural and kinematic map and a general cross-section of the Dabieshan including its foreland fold belt and the Northern Dabieshan Domain (Foziling and Luzenguang groups) are presented here. South Dabieshan consists from bottom to top of stacked allochtons: (1) an amphibolite facies gneissic unit, devoid of UHP rocks, interpreted here as the relative autochton; (2) an UHP allochton; (3) a HP rock unit (Susong group) mostly retrogressed into greenschist facies micaschists; (4) a weakly metamorphosed Proterozoic slate and sandstone unit; and (5) an unmetamorphosed Cambrian to Early Triassic sedimentary sequence unconformably covered by Jurassic sandstone. All these units exhibit a polyphase ductile deformation characterized by (i) a NW–SE lineation with a top-to-the-NW shearing, and (ii) a southward refolding of early ductile fabrics.

The Central Dabieshan is a 100-km scale migmatitic dome. Newly discovered eclogite xenoliths in a Cretaceous granitoid dated at 102 Ma by the U–Pb method on titanite demonstrate that migmatization post-dates HP–UHP metamorphism. Ductile faults formed in the subsolidus state coeval to migmatization allow us to characterize the structural pattern of doming. Along the dome margins, migmatite is gneissified under post-solidus conditions and mylonitic–ultramylonitic fabrics commonly develop. The north and west boundaries of the Central Dabieshan metamorphics, i.e. the Xiaotian–Mozitan and Macheng faults, are ductile normal faults formed before Late Jurassic–Early Cretaceous. A Cretaceous reworking is recorded by synkinematic plutons.

North of the Xiaotian–Mozitan fault, the North Dabieshan Domain consists of metasediments and orthogneiss (Foziling and Luzenguang groups) metamorphosed under greenschist to amphibolite facies which never experienced UHP metamorphism. A rare N–S-trending lineation with top-to-the-south shearing is dated at 260 Ma by the ⁴⁰Ar/³⁹Ar method on muscovite. This early structure related to compressional tectonics is reworked by top-to-the-north extensional shear bands.

The main deformation of the Dabieshan consists of a NW–SE-stretching lineation which wraps around the migmatitic dome but exhibits a consistently top-to-the-NW sense of shear. The Central Dabieshan is interpreted as an extensional migmatitic dome bounded by an arched, top-to-the-NW, detachment fault. This structure may account for a part of the UHP rock exhumation. However, the abundance of amphibolite restites in the Central Dabieshan migmatites and the scarcity of eclogites (found only in a few places) argue for an early stage of exhumation and retrogression of UHP rocks before migmatization. This event is coeval to the N–S extensional structures described in the North Dabieshan Domain. Recent radiometric dates suggest that early exhumation and subsequent migmatization occurred in Triassic–Liassic times. The main foliation is deformed by north-verging recumbent folds coeval to the south-verging folds of the South Dabieshan Domain. An intense Cretaceous magmatism accounts for thermal resetting of most of the ⁴⁰Ar/³⁹Ar dates.

A lithosphere-scale exhumation model, involving continental subduction, synconvergence extension with inversion of southward thrusts into NW-ward normal faults and crustal melting is presented.     

Evidence for fullerenes in solid bitumen from pillow lavas of Proterozoic age from Mítov (Bohemian Massif, Czech Republic)

Andesitic pillow lavas containing biogenic, solid bitumen (SB) are a constituent of a Neoproterozoic volcanosedimentary sequence (Teplá-Barrandian unit, Bohemian Massif) in the Mítov area of the Czech Republic. A black shale formation that is crosscut by these andesitic basalts is 565 Ma old. Carbon disulfide extracts of two powdered samples of SB contain 0.2 and 0.3 ppm of C₆₀, respectively, as determined by high-pressure liquid chromatography. The peak assignment based on retention time is fully supported by high-resolution electron ionization mass spectrometry (EI-MS). No C₇₀ was detected, nor was C₆₀ found in two other SB samples from this locality. Other investigated carbonaceous samples from Bohemia (coals and anthracites of Upper Paleozoic age and anthraxolite, graphitoids, and graphite of Upper Proterozoic age) did not contain fullerenes at concentrations above the detection limit of 0.01 ppm. The absence of C₆₀ in these samples was confirmed by EI-MS. The proposed mechanism of fullerene formation involves a primary algal phase, generation of a hydrocarbonaceous mixture in the course of thermal evolution of the sedimentary series, and their high-temperature transformation related to the extrusion of basalt. An important feature for fullerene conservation was the enclosure of fullerenes in SB with a structure similar to glasslike carbon, where the fullerene was protected against oxidation.     

A model for the Am (Km) planetary geomagnetic activity index and application to prediction

We studied the existence of dynamical stochastic relations in the evolution of the am index. A first analysis of the autocorrelation functions showed evidence of several seasonalities. We first used linear (ARMA) models, and it was found that these do not account for the whole internal dynamics of the data series. We then used various non-linear models to provide a better fit to reality. The forecast performances of the non-linear models are not significantly different from those of the linear model. We give a tentative explanation for the failure of the non-linear predictions. Finally, ARCH models were used in order to take into account the fact that the confidence interval for the predicted value depends on past observations.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

CO2 solubility in dacitic melts equilibrated with H2O-CO2 fluids: Implications for modeling the solubility of CO2 in silicic melts

The solubility of CO₂ in dacitic melts equilibrated with H₂O-CO₂ fluids was experimentally investigated at 1250°C and 100 to 500 MPa. CO₂ is dissolved in dacitic glasses as molecular CO₂ and carbonate. The quantification of total CO₂ in the glasses by mid-infrared (MIR) spectroscopy is difficult because the weak carbonate bands at 1430 and 1530 cm⁻¹ can not be reliably separated from background features in the spectra. Furthermore, the ratio of CO_2,mol/carbonate in the quenched glasses strongly decreases with increasing water content. Due to the difficulties in quantifying CO₂ species concentrations from the MIR spectra we have measured total CO₂ contents of dacitic glasses by secondary ion mass spectrometry (SIMS).At all pressures, the dependence of CO₂ solubility in dacitic melts on x^fluid_CO2,total shows a strong positive deviation from linearity with almost constant CO₂ solubility at x_CO2fluid > 0.8 (maximum CO₂ solubility of 795 ± 41, 1376 ± 73 and 2949 ± 166 ppm at 100, 200 and 500 MPa, respectively), indicating that dissolved water strongly enhances the solubility of CO₂. A similar nonlinear variation of CO₂ solubility with x_CO2fluid has been observed for rhyolitic melts in which carbon dioxide is incorporated exclusively as molecular CO₂ (Tamic et al., 2001). We infer that water species in the melt do not only stabilize carbonate groups as has been suggested earlier but also CO₂ molecules.A thermodynamic model describing the dependence of the CO₂ solubility in hydrous rhyolitic and dacitic melts on T, P, f_CO2 and the mol fraction of water in the melt (x_water) has been developed. An exponential variation of the equilibrium constant K₁ with x_water is proposed to account for the nonlinear dependence of x_CO2,totalmelt on x_CO2fluid. The model reproduces the CO₂ solubility data for dacitic melts within ±14% relative and the data for rhyolitic melts within 10% relative in the pressure range 100-500 MPa (except for six outliers at low x_CO2fluid). Data obtained for rhyolitic melts at 75 MPa and 850°C show a stronger deviation from the model, suggesting a change in the solubility behavior of CO₂ at low pressures (a Henrian behavior of the CO₂ solubility is observed at low pressure and low H₂O concentrations in the melt). We recommend to use our model only in the pressure range 100-500 MPa and in the x_CO2fluid range 0.1-0.95. The thermodynamic modeling indicates that the partial molar volume of total CO₂ is much lower in rhyolitic melts (31.7 cm³/mol) than in dacitic melts (46.6 cm³/mol). The dissolution enthalpy for CO₂ in hydrous rhyolitic melts was found to be negligible. This result suggests that temperature is of minor importance for CO₂ solubility in silicic melts.