Comparative phylogeography and population genetic structure of three widespread mollusc species in the Mediterranean and near Atlantic

Although several studies have evaluated the genetic structure and phylogeographic patterns in many species of marine invertebrates, a general model that applies to all of them remains elusive. For example, some species present an admixture of populations with high gene flow, whereas others exhibit more complex patterns characterized by small‐scale unstructured genetic heterogeneity, even at a local scale. These differences are thought to be due to clear biological aspects such as direct versus indirect development, or the presence of lecithotrophic versus planktotrophic larvae, but few studies compare animals with similar distributions and life modes. Here, we explore the phylogeographic and genetic structure patterns in two chiton (Chiton olivaceus and Lepidopleurus cajetanus) and one abalone (Haliotis tuberculata) species co‐occurring in the same habitat. Samples were obtained from shallow rocky bottoms along the Iberian Peninsula (Atlantic and Mediterranean coasts), Italy, Croatia and Greece, and the mitochondrial markers COI and 16S rRNA gene were sequenced. Our data show evidence of admixture and population expansion in C. olivaceus and H. tuberculata, whereas L. cajetanus exhibited a ‘chaotic patchiness’ pattern defined by a high genetic variability with locality‐exclusive haplotypes, high genetic divergence, and a lack of geographic structure. Shared haplotypes were sampled in both coasts of Iberia (for H. tuberculata) and in the Western and Eastern Mediterranean (for C. olivaceus), potentially indicating high dispersal ability and a recent expansion. The processes underlying the fine‐scale structuring in L. cajetanus remain a mystery. These results are especially interesting because the reproductive mode of the two chitons is similar but differs from that of the abalone, with a veliger larva, while instead the genetic structure of C. olivaceus and H. tuberculata are similar, thus contrasting with predictions based on the life history of the three molluscs and showing that the genetic patterns of marine species may be shaped by many factors, including historical ones.     

MVL spatiotemporal analysis for model intercomparison in EPS: application to the DEMETER multi-model ensemble

In a recent paper, Gutiérrez et al. (Nonlinear Process Geophys 15(1):109–114, 2008) introduced a new characterization of spatiotemporal error growth—the so called mean–variance logarithmic (MVL) diagram—and applied it to study ensemble prediction systems (EPS); in particular, they analyzed single-model ensembles obtained by perturbing the initial conditions. In the present work, the MVL diagram is applied to multi-model ensembles analyzing also the effect of model formulation differences. To this aim, the MVL diagram is systematically applied to the multi-model ensemble produced in the EU-funded DEMETER project. It is shown that the shared building blocks (atmospheric and ocean components) impose similar dynamics among different models and, thus, contribute to poorly sampling the model formulation uncertainty. This dynamical similarity should be taken into account, at least as a pre-screening process, before applying any objective weighting method.     

Determination of the optimal size of area to be sampled by use of the moss biomonitoring technique

The aim of the present study was to determine the optimal size of the sampling area that should be used in the moss biomonitoring technique, in relation to the spatial variation in the concentrations of contaminants determined in samples collected within short distances (<30 m). For this, the spatial structure in the concentrations of various metals and metalloids was explored by use of semivariograms calculated by use of a robust estimator of the variance. In each of the 3 sampling sites considered in the study (one sampled on two occasions), 50 samples of the moss Pseudoscleropodium purum were collected and the tissue concentrations of As, Cd, Cu, Hg, K, Ni, Pb, Se and Zn were determined. The results obtained show that in most cases (25 out of 32), there was no spatial structure in the concentrations of the elements and that the semivariograms were subject to the nugget effect. For those elements for which spatial independence or a linear spatial structure was observed, the size of the sampling area did not affect the results, and the dimensions were irrelevant within the range studied. The existence of spatial structure, which can be fitted to a spherical model, enabled a range of autocorrelation to be defined, which corresponded to an area of diameter >16 m. The proposed method must be applied to other species and in other regions, in order to standardize the size of the sampling sites in the moss biomonitoring technique.     

The role of small‐scale fold and fault development in seismogenic zones: example of the western Huércal‐Overa Basin (eastern Betic Cordillera,Spain)

The NW–SE shortening between the African and the Eurasian plates is accommodated in the eastern Betic Cordillera along a broad area that includes large N‐vergent folds and kilometric NE–SW sinistral faults with related seismicity. We have selected the best exposed small‐scale tectonic structures located in the western Huércal‐Overa Basin (Betic Cordillera) to discuss the seismotectonic implications of such structures usually developed in seismogenic zones. Subvertical ESE–WNW pure dextral faults and E–W to ENE–ESW dextral‐reverse faults and folds deform the Quaternary sediments. The La Molata structure is the most impressive example, including dextral ESE–WNW Neogene faults, active southward‐dipping reverse faults and associated ENE–WSW folds. A molar M1 assigned to Mimomys savini allows for precise dating of the folded sediments (0.95–0.83 Ma). Strain rates calculated across this structure give ～0.006 mm a^?1 horizontal shortening from the Middle Pleistocene up until now. The widespread active deformations on small‐scale structures contribute to elastic energy dissipation around the large seismogenic zones of the eastern Betics, decreasing the seismic hazard of major fault zones. Copyright © 2009 John Wiley & Sons, Ltd.     

A modified multiple tension upward infiltration method to estimate the soil hydraulic properties

Determination of saturated hydraulic conductivity, K_s, and the shape parameters α and n of the water retention curve, θ(h), is of paramount importance to characterize the water flow in the vadose zone. This work presents a modified upward infiltration method to estimate K_s, α and n from numerical inverse analysis of the measured cumulative upward infiltration (CUI) at multiple constant tension lower boundary conditions. Using the HYDRUS‐2D software, a theoretical analysis on a synthetic loam soil under different soil tensions (0, 0–10, 0–50 and 0–100 cm), with and without an overpressure step of 10 cm high from the top boundary condition at the end of the upward infiltration process, was performed to check the uniqueness and the accuracy of the solutions. Using a tension sorptivimeter device, the method was validated in a laboratory experiment on five different soils: a coarse and a fine sand, and a 1‐mm sieved loam, clay loam and silt‐gypseous soils. The estimated α and n parameters were compared to the corresponding values measured with the TDR‐pressure cell method. The theoretical analysis demonstrates that K_s and θ(h) can be simultaneously estimated from measured upward cumulative infiltration when high (>50 cm) soil tensions are initially applied at the lower boundary. Alternatively, satisfactory results can be also obtained when medium tensions (<50 cm) and the K_s calculated from the overpressure step at the end of the experiment are considered. A consistent relationship was found between the α (R² = 0.86, p < 0.02) and n (R² = 0.97, p < 0.001) values measured with the TDR‐pressure cell and the corresponding values estimated with the tension sorptivimeter. The error between the α (in logarithm scale) and n values estimated with the inverse analysis and the corresponding values measured with pressure chamber were 3.1 and 6.1%, respectively. Copyright © 2016 John Wiley & Sons, Ltd.     

Evaluation of groundwater availability and sustainable extraction rate for the Upper Santiago Valley Aquifer,Chile

The use of groundwater as water supply has increased dramatically in Santiago, Chile, during the last decades, and there is a need for accurately estimating the actual groundwater supply capacity in the upper Santiago Valley aquifer. The behavior of this aquifer was studied in order to determine the availability of water and the long-range sustainable extraction rate.Water-table depths were simulated using a numerical model with information on recharge from the last 48 years under different extraction policies. With this series of groundwater level data, groundwater level probability distribution functions were determined and extraction statistics were estimated by fitting time series models and by using the crossing theory. With this information, it has been possible to calculate the risk of being unable to supply groundwater demand because the results verify that only 67% of the water rights granted are able to be extracted on a sustained basis with a 90% exceedance probability. Furthermore, the results obtained demonstrate that the method is adequate for determining exceedance probabilities of groundwater flow.

---

Resumen El uso de las aguas subterráneas como fuente de abastecimiento se ha incrementado enormemente en Santiago de Chile durante las últimas décadas, por lo que se necesita estimar de forma adecuada la capacidad actual de suministro del acuífero superior del Valle de Santiago. Se ha estudiado el comportamiento de este acuífero para determinar la disponibilidad de recursos y la tasa de extracción sustentable a largo plazo.Se ha simulado los niveles freáticos del acuífero utilizando un modelo numérico con información de la recarga en los últimos 48 años bajo diferentes regímenes de extracción. Con esta serie de datos de los niveles en el acuífero, se ha determinado la distribución de probabilidad de los niveles de las aguas subterráneas y se ha estimado las estadísticas de las extracciones mediante el ajuste de modelos de series temporales y la teoría de cruce. Con esta información, se ha podido calcular el riego de no satisfacer la demanda de abastecimiento, puesto que los resultados verifican que sólo se podría extraer el 67% de las concesiones de agua existentes de manera sustentable con un 90% de probabilidad de excedencia. Aún más, los resultados obtenidos demuestran que el método es adecuado para determinar las probabilidades de excedencia del flujo de las aguas subterráneas.

---

Resumé Lutilisation de leau souterraine pour leau potable sest considérablement accrue à Santiago (Chili) au cours des derniers lustres; cest pourquoi il est maintenant nécessaire destimer avec précision la capacité réelle dalimentation de la nappe de laquifère supérieur de la vallée de Santiago. Le fonctionnement de cet aquifère a été étudié afin de déterminer la disponibilité en eau et le taux dextraction durable à long terme.Les niveaux de la nappe ont été simulés avec un modèle numérique à partir dinformation sur la recharge au cours des 48 dernières années pour différentes politiques dextraction. A partir de cette série de données piézométriques, des fonctions de distribution de probabilité des niveaux de la nappe ont été déterminées et des statistiques dextraction ont été estimées par ajustement de modèles de séries chronologiques en utilisant la théorie du croisement. Avec ces informations, il a été possible de calculer le risque de ne pas pouvoir satisfaire la demande deau souterraine, puisque les résultats vérifient que seulement 67% des droits deau attribués peuvent être extraits sur une base durable avec une probabilité de dépassement de 90%. En outre, les résultats obtenus démontrent que la méthode est adaptée à la détermination des probabilités de dépassement de lécoulement souterrain.

     

Origin and depositional environment of Ordovician stratiform iron mineralization from Zamora (NW Iberian Peninsula)

The Ordovician stratiform iron deposits at Zamora (NW Iberia) are arranged in several levels ranging between 0.2 and 1.5 m in thickness, which are interstratified in the upper member of the “Pielgo” Quartzites Formation (Arenig). The sandy nature, and trace-fossils corresponding to the ichnogenus Cruziana and Daedalus, together with major and trace element contents of this formation suggest an inter and subtidal, shallow marine depositional environment, which on a global scale formed part of a broad shelf situated in the northern margin of the Gondwana continent. The iron mineralization displays a foliated and banded structure due to the alternation of quartzitic, phosphatic (apatite), chamositic, chamositic-biotitic and ferriferous (magnetite and hematite) beds. They have high TiO₂, Ta, Sc, V, Nb, Co, Zn and Y contents. The magnetite contains unusually large amounts of TiO₂, V, Cr and Ni; there is also a clear depletion in Eu and the (Eu/Sm)_CN ratio is <1 and the (Sm/Yb)_CN ratio is >1. The chamosite contains high concentrations of Cr and V. These results suggest that iron was supplied from the weathering of a continental source, in combination with volcanic activity, such as within-basin basic volcanism or the presence of basic volcanic rocks in the exposed land. The physicochemical conditions of iron mineral crystallization calculated from chamosite compositions are the following: log f O₂: −38.8 to −30.7, log f S₂: −13.2 to −9.5 and T: 200 to 330 °C. These results together with the δ¹⁸O value (∼2‰) of the magnetite suggest that chamosite and magnetite were crystallized during later diagenesis and early low-grade metamorphism under redox conditions below the magnetite-hematite buffer. Received: 25 September 1996 / Accepted: 2 February 1998