Regional retrograde alteration of sub-greenschist facies chlorite to smectite

Dioctahedral smectite is present as a retrograde alteration product of chlorite in Permian-Triassic red slates of the Malaguide Complex in Sierra de Espuña (Betic Cordillera). Mineral assemblages and textures, illite crystallinity indices, and fluid inclusion data indicate sub-greenschist facies conditions that reached at least 180°C in the higher-grade tectonic unit of the Malaguide Complex, preceding formation of smectite. Smectite, having K as the dominant interlayer cation, occurs ubiquitously intercalated with trioctahedral chlorite as thin packets of layers and as individual layers that commonly change to chlorite along layers. Although some chlorite is typically homogeneous and trioctahedral, much chlorite shows signs of alteration and has compositions corresponding to different degrees of smectite contaimination. The incompatibility of metamorphic grade with the occurrence of smectite, the general association of chlorite and smectite, and the textural relations collectively show that dioctahedral smectite is derived through replacement of trioctahedral chlorite. Such replacement occurs on a regional basis and demonstrates that caution must be used in interpreting the occurrence of smectite in pelites as being due to prograde processes. Alteration of trioctahedral chlorite under oxidizing conditions due to introduction of phreatic water after uplift of the Betic Cordillera is proposed as the cause of formation of smectite.     

Subfilter-scale Fluxes over a Surface Roughness Transition. Part I: Measured Fluxes and Energy Transfer Rates

A wind-tunnel experiment was designed and carried out to study the effect of a surface roughness transition on subfilter-scale (SFS) physics in a turbulent boundary layer. Specifically, subfilter-scale stresses are evaluated that require parameterizations and are key to improving the accuracy of large-eddy simulations of the atmospheric boundary layer. The surface transition considered in this study consists of a sharp change from a rough, wire-mesh covered surface to a smooth surface. The resulting magnitude jump in aerodynamic roughnesses, M = ln(z ₀₁/z ₀₂), where z ₀₁ and z ₀₂ are the upwind and downwind aerodynamic surface roughnesses respectively, is similar to that of past experimental studies in the atmospheric boundary layer. The two-dimensional velocity fields used in this study are measured using particle image velocimetry and are acquired at several positions downwind of the roughness transition as well as over a homogeneous smooth surface. Results show that the SFS stress, resolved strain rate and SFS transfer rate of resolved kinetic energy are dependent on the position within the boundary layer relative to the surface roughness transition. A mismatch is found in the downwind trend of the SFS stress and resolved strain rate with distance from the transition. This difference of behaviour may not be captured by some eddy-viscosity type models that parameterize the SFS stress tensor as proportional to the resolved strain rate tensor. These results can be used as a benchmark to test the ability of existing and new SFS models to capture the spatial variability SFS physics associated with surface roughness heterogeneities.     

Zircon M127 – A Homogeneous Reference Material for SIMS U–Pb Geochronology Combined with Hafnium,Oxygen and,Potentially, Lithium Isotope Analysis

In this article, we document a detailed analytical characterisation of zircon M127, a homogeneous 12.7 carat gemstone from Ratnapura, Sri Lanka. Zircon M127 has TIMS‐determined mean U–Pb radiogenic isotopic ratios of 0.084743 ± 0.000027 for ²⁰⁶Pb/²³⁸U and 0.67676 ± 0.00023 for ²⁰⁷Pb/²³⁵U (weighted means, 2s uncertainties). Its ²⁰⁶Pb/²³⁸U age of 524.36 ± 0.16 Ma (95% confidence uncertainty) is concordant within the uncertainties of decay constants. The δ¹⁸O value (determined by laser fluorination) is 8.26 ± 0.06‰ VSMOW (2s), and the mean ¹⁷⁶Hf/¹⁷⁷Hf ratio (determined by solution ICP‐MS) is 0.282396 ± 0.000004 (2s). The SIMS‐determined δ⁷Li value is ?0.6 ± 0.9‰ (2s), with a mean mass fraction of 1.0 ± 0.1 μg g^?1 Li (2s). Zircon M127 contains ~ 923 μg g^?1 U. The moderate degree of radiation damage corresponds well with the time‐integrated self‐irradiation dose of 1.82 × 10¹⁸ alpha events per gram. This observation, and the (U–Th)/He age of 426 ± 7 Ma (2s), which is typical of unheated Sri Lankan zircon, enable us to exclude any thermal treatment. Zircon M127 is proposed as a reference material for the determination of zircon U–Pb ages by means of SIMS in combination with hafnium and stable isotope (oxygen and potentially also lithium) determination.     

Future regional projections of extreme temperatures in Europe: a nonstationary seasonal approach

This paper analyzes changes of maximum temperatures in Europe, which are evaluated using two state-of-the-art regional climate models from the EU ENSEMBLES project. Extremes are expressed in terms of return values using a time-dependent generalized extreme value (GEV) model fitted to monthly maxima. Unlike the standard GEV method, this approach allows analyzing return periods at different time scales (monthly, seasonal, annual, etc). The study focuses on the end of the 20th century (1961?C2000), used as a calibration/validation period, and assesses the changes projected for the period 2061?C2100 considering the A1B emission scenario. The performance of the regional models is evaluated for each season of the calibration period against the high-resolution gridded E-OBS dataset, showing a similar South-North gradient with larger values over the Mediterranean basin. The inter-RCM changes in the bias pattern with respect to the E-OBS are larger than the bias resulting from a change in the boundary conditions from ERA-40 to ECHAM5 20c3m. The maximum temperature response to increased green house gases, as projected by the A1B scenario, is consistent for both RCMs. Under that scenario, results indicate that the increments for extremes (e.g. 40-year return values) will be two or three times higher than those for the mean seasonal temperatures, particularly during Spring and Summer in Southern Europe.     

An enhanced ensemble Kalman filter scheme incorporating model error in sequential coupling between flow and geomechanics

In this work, we construct a new methodology for enhancing the predictive accuracy of sequential methods for coupling flow and geomechanics while preserving low computational cost. The new computational approach is developed within the framework of the fixed-stress split algorithm procedure in conjunction with data assimilation based on the ensemble Kalman filter (EnKF). In this context, we identify the high-fidelity model with the two-way formulation where additional source term appears in the flow equation containing the time derivative of total mean stress. The iterative scheme is then interlaced with data assimilation steps, which also incorporate the modeling error inherent to the EnKF framework. Such a procedure gives rise to an “enhanced one-way formulation,” exhibiting substantial improvement in accuracy compared with the classical one-way method. The governing equations are discretized by mixed finite elements, and numerical simulation of a 2D slab problem between injection and production wells illustrate the tremendous achievement of the method proposed herein.     

Red micas from basal ignimbrites of Mt. Vulture (Italy): interlayer content appraisal by a multi-methodic approach

The crystal chemistry of red phlogopites from Mt. Vulture (Italy) ignimbrites has been studied by electron microprobe, secondary ion mass spectrometry (SIMS), single crystal structural investigation and Fourier transform infrared (FTIR) spectroscopy. The analysed phlogopite has Fe/(Fe + Mg) ∼ 0.35, TiO₂ (wt%): 2.8–5.0 and H₂O (wt%): 1.24–3.37. Infrared spectra revealed the presence of bands due to the NH₄⁺ and H₂O stretching and bending vibrations. The samples belong to the 1M polytype. The bimodal behaviour of several structural parameters allows red micas to be clustered into two distinct groups: K⁺ ↔ NH₄⁺, H₂O and M³⁺-vacancy substitutions dominate in the first group; M^3+,4+-oxy, in the second group. It has to be pointed out that quantitative analysis of hydrogen (via SIMS) together with the characterization of the local environment of the anionic site (via FTIR) are fundamental in assessing the correct structural formula and the substitution mechanisms in micas. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Time relationships between volcanism–plutonism–alteration–mineralization in Cu-stratabound ore deposits from the Michilla mining district,northern Chile: a <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar geochronological approach

The Michilla mining district comprises one of the most important stratabound and breccia-style copper deposits of the Coastal Cordillera of northern Chile, hosted by the Middle Jurassic volcanic rocks of the La Negra Formation. ⁴⁰Ar/³⁹Ar analyses carried out on igneous and alteration minerals from volcanic and plutonic rocks in the district allow a chronological sequence of several magmatic and alteration events of the district to be established. The first event was the extrusion of a thick lava series of the La Negra Formation, dated at 159.9 ± 1.0 Ma (2σ) from the upper part of the series. A contemporaneous intrusion is dated at 159.6 ± 1.1 Ma, and later intrusive events are dated at 145.5 ± 2.8 and 137.4 ± 1.1 Ma, respectively. Analyzed alteration minerals such as adularia, sericite, and actinolite apparently give valid ⁴⁰Ar/³⁹Ar plateau and miniplateau ages. They indicate the occurrence of several alteration events at ca. 160–163, 154–157, 143–148, and 135–137 Ma. The first alteration event, being partly contemporaneous with volcanic and plutonic rocks, was probably produced in a high thermal gradient environment. The later events may be related either to a regional low-grade hydrothermal alteration/metamorphism process or to plutonic intrusions. The Cu mineralization of the Michilla district is robustly bracketed between 163.6 ± 1.9 and 137.4 ± 1.1 Ma, corresponding to dating of actinolite coexisting with early-stage chalcocite and a postmineralization barren dyke, respectively. More precisely, the association of small intrusives (a dated stock from the Michilla district) with Cu mineralization in the region strongly suggests that the main Michilla ore deposit is related to a magmatic/hydrothermal event that occurred between 157.4 ± 3.6 and 163.5 ± 1.9 Ma, contemporaneous or shortly after the extrusion of the volcanic sequence. This age is in agreement with the Re–Os age of 159 ± 16 Ma obtained from the mineralization itself (Tristá-Aguilera et al., Miner Depos, 41:99–105,2006).     

On the explicit determination of stability constants for linearized geodetic boundary value problems

The theory of GBVPs provide the basis to the approximate methods used to compute global gravity models. A standard approximation procedure is least squares, which implicitly assumes that data, e.g. gravity disturbance and gravity anomaly, are given functions in L ²(S). We know that solutions in these cases exist, but uniqueness (and coerciveness which implies stability of the numerical solutions) is the real problem. Conditions of uniqueness for the linearized fixed boundary and Molodensky problems are studied in detail. They depend on the geometry of the boundary; however, the case of linearized fixed boundary BVP puts practically no constraint on the surface S, while the linearized Molodensky BVP requires the previous knowledge of very low harmonics, for instance up to degree 25, if we want the telluroid to be free to have inclinations up to 60^°.     

A Modulated-Gradient Parametrization for the Large-Eddy Simulation of the Atmospheric Boundary Layer Using the Weather Research and Forecasting Model

The performance of the modulated-gradient subgrid-scale (SGS) model is investigated using large-eddy simulation (LES) of the neutral atmospheric boundary layer within the weather research and forecasting model. Since the model includes a finite-difference scheme for spatial derivatives, the discretization errors may affect the simulation results. We focus here on understanding the effects of finite-difference schemes on the momentum balance and the mean velocity distribution, and the requirement (or not) of the ad hoc canopy model. We find that, unlike the Smagorinsky and turbulent kinetic energy (TKE) models, the calculated mean velocity and vertical shear using the modulated-gradient model, are in good agreement with Monin–Obukhov similarity theory, without the need for an extra near-wall canopy model. The structure of the near-wall turbulent eddies is better resolved using the modulated-gradient model in comparison with the classical Smagorinsky and TKE models, which are too dissipative and yield unrealistic smoothing of the smallest resolved scales. Moreover, the SGS fluxes obtained from the modulated-gradient model are much smaller near the wall in comparison with those obtained from the regular Smagorinsky and TKE models. The apparent inability of the LES model in reproducing the mean streamwise component of the momentum balance using the total (resolved plus SGS) stress near the surface is probably due to the effect of the discretization errors, which can be calculated a posteriori using the Taylor-series expansion of the resolved velocity field. Overall, we demonstrate that the modulated-gradient model is less dissipative and yields more accurate results in comparison with the classical Smagorinsky model, with similar computational costs.