Subduction-zone structure and magmatic processes beneath Costa Rica constrained by local earthquake tomography and petrological modelling

A high-quality data set of 3790 earthquakes were simultaneously inverted for hypocentre locations and 3-D P-wave velocities in Costa Rica. Tests with synthetic data and resolution estimates derived from the resolution matrix indicate that the velocity model is well constrained in central Costa Rica to a depth of 70 km; northwestern and southeastern Costa Rica are less well resolved owing to a lack of seismic stations and seismicity. Maximum H2O content and seismic wave speeds of mid-ocean ridge basalt and harzburgite were calculated for metamorphic phase transformations relevant to subduction. Both the 3-D P-wave velocity structure and petrological modelling indicate the existence of low-velocity hydrous oceanic crust in the subducting Cocos Plate beneath central Costa Rica. Intermediate-depth seismicity correlates well with the predicted locations of hydrous metamorphic rocks, suggesting that dehydration plays a key role in generating intermediate-depth earthquakes beneath Costa Rica. Wadati-Benioff zone seismicity beneath central Costa Rica shows a remarkable decrease in maximum depth toward southeastern Costa Rica. The presence of asthenosphere beneath southeastern Costa Rica, which entered through a proposed slab window, may explain the shallowing of seismicity due to increased temperatures and associated shallowing of dehydration of the slab. Tomographic images further constrain the existence of deeply subducted seamounts beneath central Costa Rica. Large, low P-wave velocity areas within the lower crust are imaged beneath the southeasternmost volcanoes in central Costa Rica. These low velocities may represent anomalously hot material or even melt associated with active volcanism in central Costa Rica. Tomographic images and petrological modelling indicate the existence of a shallow, possibly hydrated mantle wedge beneath central Costa Ric

A lithospheric cross-section through the Swiss Alps—I. Thermokinematic modelling of the Neoalpine orogeny

In this paper we develop a forward 2-D thermokinematic model to investigate the Neoalpine 35-0 Ma phase of orogeny along the European Geotraverse (EGT) through the Swiss Alps on a crustal and lithospheric scale. Using a divergence-free kinematic model (div v=0), we define mass displacements, which subsequently serve as input to a transient thermal model. the thermal model uses critically assessed material prorameters and accounts for the depth dependence of the thermal properties in processes such as crustal thickening and mantle-lithospheric subduction. Based on the presentday density pattern of the deep seismic image and estimated exhumation and shortening rates, we derive, in a first modelling step, a mass-displacement field describing the Neoalpine orogeny as a uniform process in time. In a second—thermal—modelling step, this kinematic scenario is further refined by modelling the non-uniform cooling histories of the southern Lepontine in the Penninic domain. For that purpose we adopt lithospheric shortening rates—and consequently exhumation rates—to agree with total Neoalpine shortening, while keeping the geometry of the kinematic model fixed. the resultant thermokinematic model reflects the main characteristics of Neoalpine tectonics, and shows a good overall agreement with combined geological and geophysical data. the asymmetric feature of the present-day tectonic structure along the profile is strongly reflected in the thermal structure of the lithosphere. This demonstrates the need for a kinematic model to investigate the deep-temperature field in active tectonic provinces. For further refinement of the model, the amounts of shortening have to be more precisely estimated, and a higher spatial density in geochronological and metamorphic data is required. Furthermore, surface heat-flow values are, up to now, too uncertain to constrain the predicted surface heat flow. In summary, our results show that we need, in particular, data constraining the horizontal component of the tectonic and thermal evolution. the results of the Neoalpine orogeny modelling demonstrate that the presented thermokinematic procedure yields a good first-order approximation to investigate crustal-scale and lithospheric processes. We conclude. therefore, that the approach presented provides the potential for application not only to continent-continent collision zones, but also to any active tectonic provinc

Novel microsatellite loci for Sebaea aurea (Gentianaceae) and cross-amplification in related species.

[Premise of the study] Microsatellite loci were developed in Sebaea aurea (Gentianaceae) to investigate the functional role of diplostigmaty (i.e., the presence of additional stigmas along the style).[Methods and Results] One hundred seventy-four and 180 microsatellite loci were isolated through 454 shotgun sequencing of genomic and microsatellite-enriched DNA libraries, respectively. Sixteen polymorphic microsatellite loci were characterized, and 12 of them were selected to genotype individuals from two populations. Microsatellite amplifi cation was conducted in two multiplex groups, each containing six microsatellite loci. Cross-species amplifi cation was tested in seven other species of Sebaea . The 12 novel microsatellite loci amplifi ed only in the two most closely related species to S. aurea (i.e., S. ambigua and S. minutifl ora ) and were also polymorphic in these two species.[Conclusions] These results demonstrate the usefulness of this set of newly developed microsatellite loci to investigate the mating system and population genetic structure in S. aurea and related species.We acknowledge grants to J.K. from the Swiss National Science Foundation (PA00P3_129140) and the Velux Stiftung (project no. 679) and to J.G.S.-M. from a postdoctoral research contract “Ramón y Cajal” from the Ministerio de Ciencia e Innovación (MICINN), Spain.Peer Reviewe

NEW MOHO MAP OF ITALY

In complex tectonics regions, seismological, geophysical, and geodynamic modeling require accurate definition of the Moho geometry. Various active and passive seismic experiments performed in the central Mediterranean region revealed local information on the Moho depth, in some cases used to produce interpolated maps. In this paper, we present a new and original map of the 3-D Moho geometry obtained by integrating selected high-quality controlled source seismic and teleseismic receiver function data. The very small cell size makes the retrieved model suitable for detailed regional studies, crustal corrections in teleseismic tomography, advanced 3-D ray tracing in regional earthquake location, and local earthquake tomography. Our results show the geometry of three different Moho interfaces: the European, Adriatic-Ionian, and Tyrrhenian. The three distinct Moho are fashioned following the Alpine and Apennines subduction, collision, and back-arc spreading and show medium- to high-frequency topographic undulations reflecting the complexity of the geodynamic evolution

Strong paleoclimatic legacies in current plant functional diversity patterns across Europe

Numerous studies indicate that environmental changes during the late Quaternary have elicited long‐term disequilibria between species diversity and environment. Despite its importance for ecosystem functioning, the importance of historical environmental conditions as determinants of FD (functional diversity) remains largely unstudied. We quantified the geographic distributions of plant FD (richness and dispersion) across Europe using distribution and functional trait information for 2702 plant species. We then compared the importance of historical and contemporary factors to determine the relevance of past conditions as predictors of current plant FD in Europe. For this, we compared the strength of the relationships between FD with temperature and precipitation stability since the LGM (Last Glacial Maximum), accessibility to LGM refugia, and contemporary environmental conditions (climate, productivity, soil, topography, and land use). Functional richness and dispersion exhibited geographic patterns with strong associations to the environmental history of the region. The effect size of accessibility to LGM refugia and climate stability since the LGM was comparable to that of the contemporary predictors. Both functional richness and dispersion increased with temperature stability since the LGM and accessibility to LGM refugia. Functional richness' geographic pattern was primarily associated with accessibility to LGM refugia growing degree‐days, land use heterogeneity, diversity of soil types, and absolute minimum winter temperature. Functional dispersion's geographic pattern was primarily associated with accessibility to LGM refugia growing degree‐days and absolute minimum winter temperature. The high explained variance and model support of historical predictors are consistent with the idea that long‐term variability in environmental conditions supplements contemporary factors in shaping FD patterns at continental scales. Given the importance of FD for ecosystem functioning, future climate change may elicit not just short‐term shifts in ecosystem functioning, but also long‐term functional disequilibria

Automatic seismic phase picking and consistent observation error assessment: application to the Italian seismicity

Accuracy of seismic phase observation and consistency of timing error assessment define the quality of seismic waves arrival times. High-quality and large data sets are prerequisites for seismic tomography to enhance the resolution of crustal and upper mantle structures. In this paper we present the application of an automated picking system to some 600 000 seismograms of local earthquakes routinely recorded and archived by the Italian national seismic network. The system defines an observation weighting scheme calibrated with a hand-picked data subset and mimics the picking by an expert seismologist. The strength of this automatic picking is that once it is tuned for observation quality assessment, consistency of arrival times is strongly improved and errors are independent of the amount of data to be picked. The application to the Italian local seismicity documents that it is possible to automatically compile a precise, homogeneous and large data set of local earthquake Pg and Pn arrivals with related polarities. We demonstrate that such a data set is suitable for high-precision earthquake location, focal mechanism determination and high-resolution seismic tomograph

Global late Quaternary megafauna extinctions linked to humans, not climate change

The late Quaternary megafauna extinction was a severe global-scale event. Two factors, climate change and modern humans, have received broad support as the primary drivers, but their absolute and relative importance remains controversial. To date, focus has been on the extinction chronology of individual or small groups of species, specific geographical regions or macroscale studies at very coarse geographical and taxonomic resolution, limiting the possibility of adequately testing the proposed hypotheses. We present, to our knowledge, the first global analysis of this extinction based on comprehensive country-level data on the geographical distribution of all large mammal species (more than or equal to 10 kg) that have gone globally or continentally extinct between the beginning of the Last Interglacial at 132 000 years BP and the late Holocene 1000 years BP, testing the relative roles played by glacial–interglacial climate change and humans. We show that the severity of extinction is strongly tied to hominin palaeobiogeography, with at most a weak, Eurasia-specific link to climate change. This first species-level macroscale analysis at relatively high geographical resolution provides strong support for modern humans as the primary driver of the worldwide megafauna losses during the late Quaternary