Thermal and redox state of the subcontinental lithospheric mantle of NE Spain from thermobarometric data on mantle xenoliths

Mantle xenoliths in within-plate Cenozoic alkaline mafic lavas from NE Spain are used to assess the local subcontinental lithospheric mantle geotherm and the influence of melting and metasomatism on its oxidation state. The xenoliths are mainly anhydrous spinel lherzolites and harzburgites and gradations between, with minor pyroxenites. Most types show protogranular textures, but transitional protogranular–porphyroclastic and equigranular lherzolites also exist. Different thermometers used in the estimates provide higher subsolidus equilibrium temperatures for harzburgites (1,062 ± 29°C) than for lherzolites (972 ± 89°C), although there is overlap; the lowest temperatures correspond to porphyroclastic lherzolites, whereas pyroxenites give the highest temperatures (up to 1,257°C). Maximum pressures for subsolidus equilibrium of peridotites are at 2.0–1.8 GPa. Later they followed adiabatic decompression and harzburgites registered lower pressures (1.02 ± 0.19 GPa) than lherzolites (1.41 ± 0.27 GPa). One pyroxenite gives values consistent with the spinel lherzolite field (1.08 GPa). The shallowest barometric data are in agreement with the highest local conductive geotherms, which implies that the lithosphere–asthenosphere boundary is at 70–60 km minimum depth. Higher equilibrium temperatures for the harzburgites could be explained by the existence of mafic magma bodies or dykes at the lower crust–mantle boundary. Paleo-fO₂ conditions during partial melting as inferred from the covariation between V and MgO concentrations are mainly between QFM−1 and QFM−2 in log units. However, most thermobarometric fO₂ estimates are between QFM−1 and QFM+1, suggesting oxidation caused by later metasomatism during uplift and cooling.     

Evaluating restoration of man‐made slopes: a threshold approach balancing vegetation and rill erosion

The management of reclaimed slopes derived from industrial and civil activities (e.g. surface mining and road construction) requires the development of practical stability analysis approaches that integrate the processes and mechanisms that rule the dynamics of these ubiquitous emerging ecosystems. This work describes a new modelling approach focused on stability analysis of water‐limited reclaimed slopes, where interactive relationships between rill erosion and vegetation regulate ecosystem stability. Our framework reproduces two main groups of possible trends along the temporal evolution of reclaimed slopes: successful trends, characterized by widespread vegetation development and the effective control of rill erosion processes; and gullying trends, characterized by the progressive loss of vegetation and a sharp logistic increase in erosion rates. Furthermore, this analytical approach allows the determination of threshold values for the state variables (i.e. vegetation cover and rill erosion) that drive the system's stability, facilitating the identification of critical situations that require specific human intervention (e.g. revegetation or, in very problematic cases, revegetation combined with rill network destruction) to ensure the long‐term sustainability of the restored ecosystem. The application of our threshold analysis framework in Mediterranean‐dry reclaimed slopes derived from surface coal mining (the Teruel coalfield in central‐eastern Spain) showed a good field‐based performance. Therefore, we believe that this model is a valuable contribution for the management of water‐limited reclaimed systems, including those associated with rill erosion, as it provides a tool for the evaluation of restoration success and can play an important role in decision‐making during ecosystem restoration in severely disturbed landscapes. Copyright © 2011 John Wiley & Sons, Ltd.     

Two-dimensional geostatistical modeling and prediction of the fracture system in the Albala Granitic Pluton, SW Iberian Massif, Spain

We use a fracture index distribution method of geostatistical modelling and prediction to characterize quantitatively the fracture system in two-dimensions (2D) in the Mina Ratones area, located in the Albalá Granitic Pluton (SW Iberian Massif). The fracture index (FI) is a quantitative estimate of the fracture density in discrete domains. To validate the results of geostatistical modeling a detailed structural map of the area was also made on a scale of 1:1000. The resulting grids, expressed as pixel-maps, describe the continuous value of the FI in 2D for the whole Mina Ratones area. Based on the modelled distribution of the FI and their correlation with mapped faults, we distinguish two structural domains in the studied area: elongated bands of fracture zones with high FI values and romboidal blocks located between them with low FI values. The separation between both domains is gradual. Though a threshold value of the FI that separate both structural domains is not clearly defined, the fracture zones generally present FI>1 and the individualized blocks FI<0.50. As a consequence, the obtained grid of the FI permits the quantitative structural classification of the granitic massif in 2D and understanding fault zone architecture in the Mina Ratones area.     

Histogram and variogram inference in the multigaussian model

Several iterative algorithms are proposed to improve the histogram and variogram inference in the framework of the multigaussian model. The starting point is the variogram obtained after a traditional normal score transform. The subsequent step consists in simulating many sets of gaussian values with this variogram at the data locations, so that the ranking of the original values is honored. The expected gaussian transformation and the expected variogram are computed by an averaging operation over the simulated datasets. The variogram model is then updated and the procedure is repeated until convergence. Such an iterative algorithm can adapt to the case of tied data and despike the histogram. Two additional issues are also examined, referred to the modeling of the empirical transformation function and to the optimal pair weighting when computing the sample variogram.     

Beam Mode Expansion of Corrugated Conical Horns with Phase Correcting Lens: Application to Radioastronomy Receivers

A classical radioastronomy receiver is fed with a corrugated horn and an independent lens, both placed in a cryostat to lower the noise temperature. The beam is focused and directed using a combination of elliptical and plane mirrors. This paper proposes modifying the initial feeding system by placing the lens onto the horn aperture, thereby allowing a size reduction of the horn and lens, and a simplification of their mechanical design. The profiled lens is shaped to correct the phase error on the horn aperture. A quasi-optical model of the horn-plus-lens system has been developed using a Beam Mode Expansion (BME). Results using both a hyperbolic-planar lens and a spherical-elliptical lens, as well as results obtained by using Geometrical Optics (GO) with a Kirchoff–Huygens integration to get the far-field pattern, have been compared with measurements. As a direct application, a full focusing system for the new 40-m radiotelescope at the “Centro Astronómico de Yebes” is presented for the 22, 30 and 45 GHz bands. This paper has developed a QO model for a corrugated conical horn with a phase-correcting lens. This revised version was published online in July 2006 with corrections to the Cover Date.     

Polarized single-crystal FTIR-spectra of natural staurolite

Polarized infrared absorption spectra of thin single-crystal slabs parallel to (010) and (001) of a staurolite from Pizzo Forno, Ticino, with analyzed composition (Fe_2.9Mg_0.9Zn_0.1Mn_0.1)Al_17.5Ti_0.1(Si_7.7Al_0.3)O₄₈H₃ have been measured in the range of 3000–4000 cm^?1. From the pleochroitic behaviour of the OH-vibrations three groups of bands can be distinguished: the bands of group I, a strong band at 3445 cm^?1 plus a weak shoulder at 3358 cm^?1, and the bands of group II, a weak band centered at 3677 cm^?1 plus a shoulder at 3635 cm^?1, are assigned to the H1 and H2 protons, respectively. The bands of group III, a weak band at 3577 cm^?1 plus a shoulder, cannot be interpreted on the basis of the proton positions known so far. We assign them to an additional proton H3, which is bonded to O1 and shows a bifurcated hydrogen bridge to two O5 in a vacant T2 site.     

Assessment of the hydraulic slope flow approach using a mesoscale model

The simplified hydraulic two-layer model for a katabatic flow is analysed using the outputs from a high-resolution mesoscale simulation. A stably stratified night is simulated for the Duero basin, a complex terrain area located in the northern Spanish plateau, with large vertical and horizontal spatial resolution. Well-defined katabatic flows on the basin slopes are generated by the simulation, that are relatively stationary and quasi-bidimensional for some areas in the central part of the night. The bulk quantities used in the two-layer approach as well as the different terms in the equations are computed from the three-dimensional information provided by the mesoscale simulation. This method allows to inspect how well the simplified approach represents the katabatic flow generated by the mesoscale model. The study shows that the hydraulic model allows for a comprehensive analysis of the basic mechanisms of the slope flows but is not able to close the budget equations, since the residuals are large.     

Monitoring of landslide displacements using UAS and control methods based on lines

This study describes a new methodology for monitoring landslide displacements using a set of linear elements included in the terrain. These linear elements are digitized using the products obtained from several photogrammetric flights implemented by means of unmanned aircraft systems in several dates. The control of the displacements between dates is carried out using some new metrics and methods based on lines. The results of the control of the lines are used to analyse the general and particular behaviour of the landslide. Using this type of elements allows us to study some singular behaviour of linear elements included in the landslides (e.g. an affected road can react as a block). The methodology has been tested with success using a real case of an active slope with several landslides which was studied on three separate dates. The results have demonstrated the viability of using this methodology for studying the evolution of the landslides and their consequences regarding certain elements included in the terrain.     

Weak Contrast PP Wave Displacement R/T Coefficients in Weakly Anisotropic Elastic Media

—Approximate PP plane wave displacement coefficients of reflection and transmission for weak contrast interfaces separating weakly but arbitrarily anisotropic elastic media are presented. The PP reflection coefficient for such an interface has been derived recently by Vavry?uk and P?en?ík (1997). The PP transmission coefficient presented in this paper was derived by the same approach. The coefficients are given as a sum of the coefficient for the weak contrast interface separating two nearby isotropic media and a term depending linearly on contrasts of the so-called weak anisotropy (WA) parameters (parameters specifying deviation of properties of the medium from isotropy), across the interface. While the reflection coefficient depends only on 8 of the complete set of the WA parameters describing P-wave phase velocity in weakly anisotropic media, the transmission coefficient depends on their complete set. The PP reflection coefficient depends on "shear-wave splitting parameter" γ. Tests of accuracy of the approximate formulae are presented on several models.