Zircon geochronology and geochemistry to constrain the youngest eruption events and magma evolution of the Mid-Miocene ignimbrite flare-up in the Pannonian Basin,eastern central Europe

A silicic ignimbrite flare-up episode occurred in the Pannonian Basin during the Miocene, coeval with the syn-extensional period in the region. It produced important correlation horizons in the regional stratigraphy; however, they lacked precise and accurate geochronology. Here, we used U–Pb (LA-ICP-MS and ID-TIMS) and (U–Th)/He dating of zircons to determine the eruption ages of the youngest stage of this volcanic activity and constrain the longevity of the magma storage in crustal reservoirs. Reliability of the U–Pb data is supported by (U–Th)/He zircon dating and magnetostratigraphic constraints. We distinguish four eruptive phases from 15.9 ± 0.3 to 14.1 ± 0.3 Ma, each of which possibly includes multiple eruptive events. Among these, at least two large volume eruptions (>10 km³) occurred at 14.8 ± 0.3 Ma (Demjén ignimbrite) and 14.1 ± 0.3 Ma (Harsány ignimbrite). The in situ U–Pb zircon dating shows wide age ranges (up to 700 kyr) in most of the crystal-poor pyroclastic units, containing few to no xenocrysts, which implies efficient recycling of antecrysts. We propose that long-lived silicic magma reservoirs, mostly kept as high-crystallinity mushes, have existed in the Pannonian Basin during the 16–14 Ma period. Small but significant differences in zircon, bulk rock and glass shard composition among units suggest the presence of spatially separated reservoirs, sometimes existing contemporaneously. Our results also better constrain the time frame of the main tectonic events that occurred in the Northern Pannonian Basin: We refined the upper temporal boundary (15 Ma) of the youngest counterclockwise block rotation and the beginning of a new deformation phase, which structurally characterized the onset of the youngest volcanic and sedimentary phase.     

Offshore experiments on styrene spillage in marine waters for risk assessment

Within the context of risk evaluation of chemical spillages into the marine environment, this paper reports on an offshore experiment to study the behaviour of styrene spilled into sea under natural conditions and discusses theoretical approaches. Floating structures were used to enclose the spillage and the gaseous cloud formation, and dissolution processes were in situ monitored. The identification of spill risks for man and marine environment through GESAMP's hazard profile is described for styrene: Styrene is rated as a chemical with a significant health hazard that will float but also evaporate. However, monitoring of the water column in the experiments showed that the concentration of styrene in water during the first hour represents 50% of the product spilled. For the potentially exposed public, the GESAMP hazard rating recommends the closure of beaches and evacuation. The risk assessment developed from experimental data confirms this safety advice.     

The respective roles of surface temperature driven feedbacks and tropospheric adjustment to CO2 in CMIP5 transient climate simulations

An overview of radiative climate feedbacks and ocean heat uptake efficiency diagnosed from idealized transient climate change experiments of 14 CMIP5 models is presented. Feedbacks explain about two times more variance in transient climate response across the models than ocean heat uptake efficiency. Cloud feedbacks can clearly be identified as the main source of inter-model spread. Models with strong longwave feedbacks in the tropics feature substantial increases in cloud ice around the tropopause suggestive of changes in cloud-top heights. The lifting of the tropical tropopause goes together with a general weakening of the tropical circulation. Distinctive inter-model differences in cloud shortwave feedbacks occur in the subtropics including the equatorward flanks of the storm-tracks. Related cloud fraction changes are not confined to low clouds but comprise middle level clouds as well. A reduction in relative humidity through the lower and mid troposphere can be identified as being the main associated large-scale feature. Experiments with prescribed sea surface temperatures are analyzed in order to investigate whether the diagnosed feedbacks from the transient climate simulations contain a tropospheric adjustment component that is not conveyed through the surface temperature response. The strengths of the climate feedbacks computed from atmosphere-only experiments with prescribed increases in sea surface temperatures, but fixed CO₂ concentrations, are close to the ones derived from the transient experiment. Only the cloud shortwave feedback exhibits discernible differences which, however, can not unequivocally be attributed to tropospheric adjustment to CO₂. Although for some models a tropospheric adjustment component is present in the global mean shortwave cloud feedback, an analysis of spatial patterns does not lend support to the view that cloud feedbacks are dominated by their tropospheric adjustment part. Nevertheless, there is positive correlation between the strength of tropospheric adjustment processes and cloud feedbacks across different climate models.     

A thermodynamic model for Fe–Mg dioctahedral K white micas using data from phase-equilibrium experiments and natural pelitic assemblages

The purpose of this study is to derive a solid-solution model for potassic white micas (KWM) encountered in rocks of various bulk compositions, over a wide range of P-T conditions. A compilation of phengite compositions lead us to propose a seven-thermodynamic-component (muscovite, Fe²⁺-Al-celadonite, Mg-Al-celadonite, annite, phlogopite, pyrophyllite and paragonite) ionic solid-solution model which accounts for the Tschermak, Fe-Mg, di/trioctahedral, pyrophyllitic and paragonitic substitutions observed in nature. A four-site mixing model with symmetric Margules parameters to model the Tschermak substitutions, asymmetric Margules parameters to model the other substitutions, and ideal intersite interaction has been adopted. In contrast to previous models, the relevant thermodynamic data and solid-solution properties are calibrated with independent sets of published experiments conducted for the KMASH, KFASH, KFMASH, and KNASH systems, as well as about 200 natural data involving KWM assemblages. The constraints span a wide range of pressure and temperature conditions (150 to 750 °C, 0.5 to 30 kbar), so that our model does not need to be extrapolated outside the calibration range to be used for P-T thermobarometric purposes. The calculated thermodynamic data are interconsistent with the TWQ thermodynamic database and solid-solution models, including that recently published for chlorites.     

Asymptotic Analysis and Layer Decomposition for the Couplex Exercise

This document explores the possibilities of multiscale expansions and domain decomposition to solve part of the Couplex 1 exercise. We concentrate on the hydrostatic pressure and show that the slenderness of the domain and the large variations of the Darcy constants allows an analytical approximation which our test reveals to be true to relative errors smaller than 10^–3. The numerical tests are done in 2D with freefem+ and in 3D with freefem3D. Some considerations are also given for Iodine transport.     

Constraints in 3D gravity inversion

A three-dimensional (3D) inversion program is developed to interpret gravity data using a selection of constraints. This selection includes minimum distance, flatness, smoothness and compactness constraints, which can be combined using a Lagrangian formulation. A multigrid technique is also implemented to resolve separately large and short gravity wavelengths. The subsurface in the survey area is divided into rectangular prismatic blocks and the problem is solved by calculating the model parameters, i.e. the densities of each block. Weights are given to each block depending on depth, a priori information on density and the density range allowed for the region under investigation. The present computer code is tested on modelled data for a dipping dike and multiple bodies. Results combining different constraints and a weight depending on depth are shown for the dipping dike. The advantages and behaviour of each method are compared in the 3D reconstruction. Recovery of geometry (depth, size) and density distribution of the original model is dependent on the set of constraints used. From experimentation, the best combination of constraints for multiple bodies seems to be flatness and a minimum volume for the multiple bodies. The inversion method is tested on real gravity data from the Rouyn-Noranda (Quebec) mining camp. The 3D inversion model for the first 10 km is in agreement with the known major lithological contacts at the surface; it enables the determination of the geometry of plutons and intrusive rocks at depth.     

Coseismic and post-seismic signatures of the Sumatra 2004 December and 2005 March earthquakes in GRACE satellite gravity

The GRACE satellite mission has been measuring the Earth's gravity field and its temporal variations since 2002 April. Although these variations are mainly due to mass transfer within the geofluid envelops, they also result from mass displacements associated with phenomena including glacial isostatic adjustment and earthquakes. However, these last contributions are difficult to isolate because of the presence of noise and of geofluid signals, and because of GRACE's coarse spatial resolution (>400 km half-wavelength). In this paper, we show that a wavelet analysis on the sphere helps to retrieve earthquake signatures from GRACE geoid products. Using a wavelet analysis of GRACE geoids products, we show that the geoid variations caused by the 2004 December ( M _w= 9.2) and 2005 March ( M _w= 8.7) Sumatra earthquakes can be detected. At GRACE resolution, the 2004 December earthquake produced a strong coseismic decrease of the gravity field in the Andaman Sea, followed by relaxation in the area affected by both the Andaman 2004 and the Nias 2005 earthquakes. We find two characteristic timescales for the relaxation, with a fast variation occurring in the vicinity of the Central Andaman ridge. We discuss our coseismic observations in terms of density changes of crustal and upper-mantle rocks, and of the vertical displacements in the Andaman Sea. We interpret the post-seismic signal in terms of the viscoelastic response of the Earth's mantle. The transient component of the relaxation may indicate the presence of hot, viscous material beneath the active Central Andaman Basin.