Origin and implications of mafic xenoliths associated with Cenozoic extension-related volcanism in the Válencia Trough, NE Spain

Summary The Válencia Trough of northeastern Spain represents a poorly understood region of Cenozoic extension-related magmatism in western Europe. We present chemical data on mafic xenoliths and their host basalt from the Roca Negra cinder cone in the Catalán Volcanic Zone in the Válencia Trough. The xenoliths consist of augite+hornblende+ oxides±plagioclase±apatite±olivine::Lbiotitelorthopyroxene, and range from clinopyroxenites and hornblendites to gabbroic rocks. The major and trace element compositions of the xenoliths are compatible with formation as cumulates from a range of olivine tholeiitic to strongly alkaline magmas, and show a strong affinity to the Middle Miocene to Recent rift-related magmatism in the Cátalan Volcanic Zone. The xenoliths formed at temperatures of 1110–1180 °C, pressures of 2–7 kb (with 10 kb as the upper limit), and oxygen fugacities corresponding to 1–2 log units above the NNO oxygen buffer. The estimated pressure range implies formation in magma chambers within the crust, possibly also near the crust/mantle boundary (underplating). The presence of crustal cumulates indicates that the total volume of magmas emplaced into the crust in the Cátalan Volcanic Zone may be significantly larger than indicated by surface exposures. Isotopic compositions of the xenoliths are⁸⁷Sr/⁸⁶Sr: 0.703537–0.703647,¹⁴³Nd/¹⁴⁴Nd: 0.512732–0.512766,²⁰⁶Pb/²⁰⁴Pb: 18.097–19.106,²⁰⁷Pb/²⁰⁴Pb: 15.522–15.579, and²⁰⁸Pb/²⁰⁴Pb: 37.850–38.794. This range suggests the involvement of two mantle source components. One of these may be the low-velocity component recognized in Cenozoic magmatism in the eastern Atlantic Ocean to central Europe and western Mediterranean. The other component, which differs from those involved in other domains of Cenozoic volcanism in Europe and adjacent areas, may be located in the lithospheric mantle underlying NE Spain, which was metasomatized during late Variscan to early Alpine deformation.

---

Herkunft und Bedeutung mafischer Xenolithe aus känozoiscben Extensions-bezogenen Vulkaniten aus dem Valencia Trog, Nordost-Spanien

Zusammenfassung Der Válencia Trog in Nordost-Spanien stellt eine noch wenig verstandene Region mit känozoischen Extensions-bezogenern Magmatismus in Westeuropa dar. Wir geben hier chemische Daten über mafische Xenolithe und ihre basaltischen Muttergesteine aus dem Roca Negra Kegei in der vulkanischen Zone von Catalan im Valencia Trog. Die Xenolithe bestehen aus Augit+Hornblende+Oxiden±Plagioklas±Apatit±Olivin±Biotit±Orthopyroxen und reichen in ihrer Zusammensetzung von Klinopyroxeniten und Homblenditen his zu gabbroischen Gesteinen. Die Haupt- und Spurenelementzusammensetzung der Xenolithe entspricht einer Kumulatbildung aus Olivin-tholéitischen bis stark alkalischen Magmen and zeigt eine deutliche Affinität mit dem mittelmiozanen bis rezenten Rift-bezogenen Magmatismus in der vulkanischen Zone von Catalan. Die Xenolithe bildeten sich bei Temperaturen von 1.110 bis 1.180 °C und bei Drucken von 2–7 kb (mit 10 kb als Obergrenze) und bei Sauerstoff-Fugazitäten, die 1 bis 2 log Einheiten über dem NNO Sauerstoff-Puffer liegen. Die abgeschätzten Drucke weisen auf Bildung in Magmakammern innerhalb der Kruste bin, womölich auch in der Nähe der Kruste-Mantelgrenze (Underplating). Das Vorkommen von krustalen Kumulaten zeigt, daß das gesamte Volumen von Magma, das in die Kruste der vulkanischen Zone von Catalan eingedrungen ist, wesentlich größer ist als das, das man aus Oberflächen-Ausbissen entnehmen k6nnte. Die Isotopen-Zusammensetzung der Xenolithe sind⁸⁷Sr/⁸⁶Sr: 0.703537–0.703647,¹⁴³Nd/¹⁴⁴Nd: 0.512732–0.512766,²⁰⁶Pb/204 Pb: 18.097–19.106,²⁰⁷Pb/²⁰⁴Pb: 15.522–15.579, and²⁰⁸Pb/²⁰⁴Pb: 37.850–38.794. Dies läßt erkennen, daß es sich bier um zwei verschiedene Mantel-Quellen handelt. Eine davon könnte die Niedrig-Geschwindigkeitskomponente sein, die man in känozoiscben Magmatiten des östlichen Atlantischen Ozeans, sowie in Zentraleuropa und im westlichen Mittelmeer erkannt hat. Die andere Komponete unterscheidet sich von jenen, die in anderen Bereichen von känozoischem Vulkanismus in Europa und benachbarten Gebieten vorkommen. Sie könnte im lithosphärischen Mantel unterhalb von Nordost-Spanien liegen, der während der variszischen bis zur frühalpidischen Deformation metasomatisch beeinflußt wurde.

---

---

With 10 Figures     

Structural and kinematic analysis of the Corredoiras detachment: evidence for early Variscan synconvergent extension in the Ordenes Complex, NW Spain

 An important detachment is described in the allochthonous Ordenes Complex, in the NW Iberian Massif, and its meaning is related to the kinematics of contemporaneous convergent structures. The Corredoiras Detachment (CD) separates a hangingwall unit, characterised by a medium-pressure metamorphic gradient, from a footwall high-pressure and high-temperature unit and an underlying ophiolitic unit. An associated ductile shear zone, nearly 2000 m thick, developed in the lower part of the hangingwall unit, where the Corredoiras Orthogneiss, a Lower Ordovician metagranite, was progressively transformed into augengneisses, mylonitic and ultramylonitic gneisses. The attitude of the stretching and mineral lineation in the mylonites varies due to late refolding at map scale, but the sense of movement can be estimated, being roughly top to the SE. According to crosscutting relationships, the CD developed subsequent to the thrusting of the high-pressure/high-temperature unit onto the ophiolitic unit, and prior to younger extensional detachments, upright folding and strike-slip tectonics. The geometric relationships of the CD with the previous structures in the footwall unit, the subtractive character of the metamorphic gap between its hangingwall and footwall, and the available isotopic data suggest that the CD is an early Variscan, ductile extensional detachment, the movement of which was roughly simultaneous with the onset of thrusting of the allochthonous complexes over their relative autochthon. Received: 17 November 1998 / Accepted: 4 April 1999     

Physico-chemical evolution of groundwater in tectonically active areas. Application to the Leana hot spring (Murcia Region,SE Spain)

Seismic events can affect the physico-chemical characteristics of groundwater. These anomalies are of a pre-seismic, co-seismic and post-seismic nature and correspond to pulse variations, sudden increases and decreases without return to initial values and upward or downward changes in trend. Continuous and in situ conductivity and temperature monitoring and periodic water sampling at a hot spring associated with neotectonic activity are of great interest for establishing predictive methods. This method is limited to the seismic activity affecting the fracturing system with which the hot spring is associated. The Region of Murcia and surroundings (southeast Spain) was selected as the study area for exploring the nature of these influences on groundwater. A hot spring in the Leana spa (Murcia) was equipped and monitored during the period 2006–2008, allowing for the in situ determination of conductivity and temperature as well as of major and minor constituents at the laboratory. Due to its proximity and related with fault network, we suggest that 86 % of earthquakes located between 0 and 10 km may affect in situ parameters of groundwater, and 75 % may affect laboratory determinations. This percentage drops in more distant zones. Of all earthquakes that seem to influence groundwater, 55 % of the in situ parameter anomalies and 53 % of laboratory were of a pre-seismic nature.     

Smoothing and high risk areas detection in space-time disease mapping: a comparison of P-splines,autoregressive, and moving average models

Recently, several models have been proposed for smoothing risks in disease mapping. These models consider different ways of introducing both spatial and temporal dependence as well as spatio-temporal interactions. In this work, a comparison among some autoregressive, moving average, and P-spline models is performed. Firstly, brain cancer mortality data are used to analyze the degree of smoothness introduced by these models. Secondly, two separate simulation studies (one model-based and the other model-free) are carried out to evaluate the model performance in terms of bias, variability, sensitivity, and specificity. We conclude that P-spline models seem to be a good alternative to autoregressive and moving average models when analyzing highly sparse disease mapping data.     

Uncertainty analysis and probabilistic segmentation of electrical resistivity images: the 2D inverse problem

In this paper, we present the uncertainty analysis of the 2D electrical tomography inverse problem using model reduction and performing the sampling via an explorative member of the Particle Swarm Optimization family, called the Regressive‐Regressive Particle Swarm Optimization. The procedure begins with a local inversion to find a good resistivity model located in the nonlinear equivalence region of the set of plausible solutions. The dimension of this geophysical model is then reduced using spectral decomposition, and the uncertainty space is explored via Particle Swarm Optimization. Using this approach, we show that it is possible to sample the uncertainty space of the electrical tomography inverse problem. We illustrate this methodology with the application to a synthetic and a real dataset coming from a karstic geological set‐up. By computing the uncertainty of the inverse solution, it is possible to perform the segmentation of the resistivity images issued from inversion. This segmentation is based on the set of equivalent models that have been sampled, and makes it possible to answer geophysical questions in a probabilistic way, performing risk analysis.     

A geostatistical protocol to optimize spatial sampling of domestic drinking water supplies in remote environments

This paper deals with the design of optimal spatial sampling of water quality variables in remote regions, where logistics are complicated and the optimization of monitoring networks may be critical to maximize the effectiveness of human and material resources. A methodology that combines the probability of exceeding some particular thresholds with a measurement of the information provided by each pair of experimental points has been introduced. This network optimization concept, where the basic unit of information is not a single spatial location but a pair of spatial locations, is used to emphasize the locations with the greatest information, which are those at the border of the phenomenon (for example contamination or a quality variable exceeding a given threshold), that is, where the variable at one of the locations in the pair is above the threshold value and the other is below the threshold. The methodology is illustrated with a case of optimizing the monitoring network by optimal selection of the subset that best describes the information provided by an exhaustive survey done at a given moment in time but which cannot be repeated systematically due to time or economic constrains.     

Modelling regional effects of artificial groundwater recharge in a multilayer aquifer characterized by perched water tables

This paper presents an approach to estimate the effects of a managed recharge experiment in a multilayer aquifer characterized by the presence of perched water tables in the Medina del Campo groundwater body, Douro basin, central Spain. A numerical model was developed to evaluate the effect of artificial recharge on the shallow sector of a regional-scale aquifer and on formerly active wetlands. The model was developed in the Visual MODFLOW Pro v.2011.1 environment in order to represent and analyse the regional impact of this artificial recharge event. Results suggest that the assumption of a single perched system may prove useful in regional contexts where data is limited. From a study site perspective, managed recharge is observed to increase shallow storage along the riverbanks, which is considered valuable for environmental purposes. However, downstream wetlands are unlikely to experience a significant recovery. Furthermore, only a small percentage of artificial recharge is expected to reach the deep regional aquifer. This method can be exported to settings characterized by the presence of perched aquifers and associated groundwater dependent ecosystems.     

Role of time variant and invariant terms on soil moisture spatio-temporal variability in a subhumid coastal wetland

The knowledge of soil moisture spatio-temporal variability is highly relevant for water resources management. This paper reports an analysis of the spatial–temporal variability of soil moisture data for a small to medium-scale soil-sensors network in a coastal wetland of southwestern Spain. Measurements were taken from five sites located in the Doñana National Park over the time-period of one hydrological year from September 2017 to September 2018. The total area of the soil-sensors network shows an extension about 25 × 3 km. Soil moisture data was separated into time invariant (the temporal mean of the whole period at each site) and time-variant terms (the deviations of soil moisture from the mean, or anomalies). The time-invariant component was generally the main contributor to the total spatial variance of soil moisture and it was mostly controlled by the groundwater levels in the area. Nevertheless, the time variant terms have a huge effect on soil moisture variability in very dry states. Characteristic convex time-dependent patterns for this field site were found between spatially averaged soil moisture and its variability. This information could be used for the up and downscaling of soil moisture from satellite data. Those patterns of relation between spatial mean and variability of soil moisture were only affected by heavy rainfalls giving rise to hysteretic behaviour. This study shows that even though groundwater level is a time-variant variable, it significantly affects soil moisture's time-variant but also time-invariant terms due to the different average groundwater level depths at the different sites.     

Implementation of Solute Transport in the Vadose Zone into the “HYDRUS Package for MODFLOW”

The “HYDRUS package for MODFLOW” is an existing MODFLOW package that allows MODFLOW to simultaneously evaluate transient water flow in both unsaturated and saturated zones. The package is based on incorporating parts of the HYDRUS-1D model (to simulate unsaturated water flow in the vadose zone) into MODFLOW (to simulate saturated groundwater flow). The coupled model is effective in addressing spatially variable saturated-unsaturated hydrological processes at the regional scale. However, one of the major limitations of this coupled model is that it does not have the capability to simulate solute transport along with water flow and therefore, the model cannot be employed for evaluating groundwater contamination. In this work, a modified unsaturated flow and transport package (modified HYDRUS package for MODFLOW and MT3DMS) has been developed and linked to the three-dimensional (3D) groundwater flow model MODFLOW and the 3D groundwater solute transport model MT3DMS. The new package can simulate, in addition to water flow in the vadose zone, also solute transport involving many biogeochemical processes and reactions, including first-order degradation, volatilization, linear or nonlinear sorption, one-site kinetic sorption, two-site sorption, and two-kinetic sites sorption. Due to complex interactions at the groundwater table, certain modifications of the pressure head (compared to the original coupling) and solute concentration profiles were incorporated into the modified HYDRUS package. The performance of the newly developed model is evaluated using HYDRUS (2D/3D), and the results indicate that the new model is effective in simulating the movement of water and contaminants in the saturated-unsaturated flow domains.     

Environmental factors influencing aggregation of manta rays (Manta birostris) off the northeastern coast of the Yucatan Peninsula

Manta rays inhabit tropical, subtropical and temperate waters. Aggregation sites of manta rays have been recognized worldwide, but the reasons for this behavior are still poorly understood. This study describes environmental factors influencing aggregation sites of the giant manta ray (Manta birostris) off the northeastern coast of the Yucatan Peninsula. Observations of manta rays were obtained from scientific surveys conducted during 2006–2011. Environmental data were obtained from satellite imagery. The maximum entropy (Maxent) method for habitat modeling was used to determine the effects of environmental conditions on the species and predict suitable habitat for manta rays in this region. Primary productivity and distance to the coast were the most influential variables, suggesting that aggregation occurs in highly productive coastal waters. The distribution of manta rays predicted by the Maxent model showed that the most suitable habitat within the study area is located off the northeastern coast of the Yucatán Peninsula, more precisely, northeast of Isla Holbox and northwest of Isla Contoy. Seasonal patterns of distribution suggest that the most suitable conditions are present from July through September.