Deep to shallow lacustrine evaporites in the Libros Gypsum (southern Teruel Basin, Miocene, NE Spain): an occurrence of pelletal gypsum rhythmites

ABSTRACT A number of non‐marine evaporite units composed of primary gypsum were deposited in saline lakes that developed in the southern Teruel Basin (NE Spain) during the Miocene. In the basin depocentre, a continuum of lacustrine evaporite lithofacies influenced by the activity of organisms is displayed. The Libros Gypsum was deposited in a deep lake, in which water stratification became unstable with progressive shoaling. Rhythmites, composed of laminae of pelletal gypsum and laminae of very fine lenticular gypsum crystals mixed with siliceous microorganisms, formed in addition to gypsum turbidites, intraformational gypsum breccias and slump structures. The pelletal laminae originated from the faecal activity of animals (crustaceans?) ingesting gypsum crystallites in the lake water during episodes of maximum evaporation, whereas the laminae of very fine lenticular gypsum mixed with microorganisms accumulated during episodes of relative dilution. In the wide marginal zones of the basin, the Libros Gypsum unit consists of massive to thin‐bedded bioturbated gypsum and thin‐bedded clotted gypsum, which formed in intermediate to very shallow (palustrine) water depths. The bioturbated gypsum lithofacies were produced by the action of diverse organisms, presumably worms and coleopterans, and chironomid larvae to a lesser extent; the massive lithofacies precipitated in very shallow water; and the thin‐bedded lithofacies formed in shallow to deeper settings. The thin‐bedded clotted gypsum is a relatively deep facies that may have diverse origins (e.g. bioturbation, compaction, disruption of soft sediments and early diagenesis). There is a well‐developed metre‐scale cyclicity in the marginal lake sequences, which is not observed in the inner lake deposits. This suggests a depth control in the various lacustrine subenvironments to record cyclic evaporitic processes. The isotopic composition of the gypsum indicates early sulphate‐reducing bacterial activity in the bottom of the lake and suggests that the sulphate was derived from the chemical recycling of Triassic evaporites of the country rocks.     

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.     

Conceptually-based shot noise modeling of streamflows at short time interval

A conceptual-stochastic approach to short time runoff data modelling is proposed, according to the aim of reproducing the hydrological aspects of the streamflow process and of preserving as much as possible the dynamics of the process itself. This latter task implies preservation of streamflow characteristics at higher scales of aggregation and, within a conceptual framework, involves compatibility with models proposed for the runoff process at those scales. At a daily time scale the watershed response to the effective rainfall is considered as deriving from the response of three linear reservoirs, respectively representing contributions to streamflows of large deep aquifers, with over-year response lag, of aquifers which run dry by the end of the dry season and of subsurface runoff. The surface runoff component is regarded as an uncorrelated point process. Considering the occurrences of effective rainfall events as generated by an independent Poisson process, the output of the linear system represents a conceptually-based multiple shot noise process. Model identification and parameter estimation are supported by information related to the aggregated runoff process, in agreement to the conceptual framework proposed, and this allows parameter parsimony, efficient estimation and effectiveness of the streamflow reproduction. Good performances emerged from the model application and testing made with reference to some daily runoff series from Italian basins.     

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.     

Ar / Ar dating of high-pressure rocks from the Ligurian Alps: Evidence for a continuous subduction–exhumation cycle

We present ³⁹Ar–⁴⁰Ar dating of phengite, muscovite and paragonite from a set of mafic and metasedimentary rocks sampled from the high-pressure (HP) metaophiolites of the Voltri Group (Western Alps) and from clasts in the basal layer conglomerates from the Tertiary molasse which overlie the high-pressure basement. The white mica-bearing rocks display peak eclogitic and blueschist-facies parageneses, locally showing complex greenschist-facies replacement textures. The internal discordance of age spectra is proportional to the chemical complexity of the micas. High-Si phengites from eclogite clasts record a ³⁹Ar–⁴⁰Ar age of ca. 49 Ma for the eclogite stage and ca. 43 Ma for the blueschist retrogression; phengites from a blueschist basement sample yield an age of ca. 40 Ma; low-Si muscovite from a metasediment dates the formation of the greenschist paragenesis at ca. 33 Ma. Our data indicate that the analyzed samples reached high-pressure conditions at different times over a time-span of c.a. 10 Ma. Subduction was continuing during exhumation and blueschist retrograde re-equilibration of higher-pressure, eclogite-facies rocks. This process kept the isotherms depressed, allowing the older HP-rocks to escape thermal re-equilibration. Our results, added to literature data, fit a tectonic model of a subduction–exhumation cycle, with different tectonic slices subducted at different times from Early Eocene until the Eocene–Oligocene boundary.