AgeSpectraAnalyst: A MATLAB based package to model zircon age distributions in silicic magmatic systems.

In the last decade, improvements in the analytical precision achievable by zircon U-Pb geochronological techniques have allowed to resolve complexities of zircon crystallization histories in magmatic rocks to an unprecedented level. A number of studies have strived to link resolvable dispersion in zircon age spectra of samples from fossil magmatic systems to the physical parameters of their parent magma bodies. However, the methodologies developed have so far been limited to reproduce the effect of simple thermal histories on the final distribution of zircon ages. In this work we take a more nuanced approach, fine-tuning a thermodynamics-based zircon saturation model to predict the relative distribution of zircon ages in samples from silicic magma reservoirs experiencing open-system processes (e.g. heat/mass addition, mechanical mixing). Employing the MATLAB package (AgeSpectraAnalyst) presented in this contribution:•Users can forward model the effect that diverse thermal histories and mechanical mixing processes characteristic of silicic magma bodies have on zircon age distributions as measured by high-precision, chemical abrasion thermal ionization mass spectrometry (CA-ID-TIMS) U-Pb geochronology.•Zircon CA-ID-TIMS datasets from silicic magmatic systems can be easily compared with model output to gain semi-quantitative information on thermo-mechanical history of the system of interest.•We demonstrated (Tavazzani et al., in press) that distribution of high-precision zircon ages in crystallized remnants of shallow (∼ 250 MPa), silicic magma reservoirs can discriminate between systems that experienced catastrophic, caldera-forming eruptions and systems that underwent monotonic cooling histories.

Recent speciation associated with range expansion and a shift to self-fertilization in North American Arabidopsis.

The main processes classically evoked for promoting reproductive isolation and speciation are geographic separation reducing gene flow among populations, divergent selection, and chance genomic change. In a case study, we present evidence that the additional factors of climate change, range expansion and a shift in mating towards inbreeding can initiate the processes leading to parapatric speciation. At the end of the last Pleistocene glaciation cycle, the North American plant Arabidopsis lyrata expanded its range and concomitantly lost its reproductive mode of outcrossing multiple times. We show that in one of the newly colonized areas, the self-fertilizing recolonization lineage of A. lyrata gave rise to selfing A. arenicola, which expanded its range to subarctic and arctic Canada and Greenland, while the parental species remained restricted to temperate North America. Despite the vast range expansion by the new species, mutational load did not increase, probably because of selfing and quasi-clonal selection. We conclude that such peripheral parapatric speciation combined with range expansion and inbreeding may be an important but so far overlooked mode of speciation.

Garnet petrochronology reveals the lifetime and dynamics of phonolitic magma chambers at Somma-Vesuvius.

Somma-Vesuvius is one of the most iconic active volcanoes with historic and archeological records of numerous hazardous eruptions. Petrologic studies of eruptive products provide insights into the evolution of the magma reservoir before eruption. Here, we quantify the duration of shallow crustal storage and document the evolution of phonolitic magmas before major eruptions of Somma-Vesuvius. Garnet uranium-thorium petrochronology suggests progressively shorter pre-eruption residence times throughout the lifetime of the volcano. Residence times mirror the repose intervals between eruptions, implying that distinct phonolite magma batches were present throughout most of the volcano's evolution, thereby controlling the eruption dynamics by preventing the ascent of mafic magmas from longer-lived and deeper reservoirs. Frequent lower-energy eruptions during the recent history sample this deeper reservoir and suggest that future Plinian eruptions are unlikely without centuries of volcanic quiescence. Crystal residence times from other volcanoes reveal that long-lived deep-seated reservoirs and transient upper crustal magma chambers are common features of subvolcanic plumbing systems.

Monitoring air pollution close to a cement plant and in a multi-source industrial area through tree-ring analysis.

Thirty-two trace elements were examined in the tree rings of downy oak to evaluate the pollution levels close to a cement plant isolated in a rural context and an industrial area where multiple sources of air pollution are or were present. Tree cores were collected from trees growing 1 km from both the cement plant and the industrial area that are located 8 km from each other. The analysis of the trace elements was performed on annual tree rings from 1990 to 2016 using laser ablation inductively coupled plasma mass spectrometry. Trace elements Cs, Mg, Mn, S and Zn reflected the emission history of the cement plant. Their values have increased since early 2000s, when the cement plant started its activity. However, the lack of significant trends of pollutants in the tree rings from the industrial area and the possible effect of translocation and volatility of some elements left open questions. The very weak changes of the other trace elements in the period 1990-2016 suggest those elements do not mark any additional effect of the industrial activity on the background pollution. The results confirm that downy oak trees growing close to isolated industrial plants must be considered a pollution forest archive accessible through dendrochemistry.

Demographic Processes Linked to Genetic Diversity and Positive Selection across a Species' Range.

Demography determines the strength of genetic drift, which generally reduces genetic variation and the efficacy of selection. Here, we disentangled the importance of demographic processes at a local scale (census size and mating system) and at a species-range scale (old split between population clusters, recolonization after the last glaciation cycle, and admixture) in determining within-population genomic diversity and genomic signatures of positive selection. Analyses were based on re-sequence data from 52 populations of North American Arabidopsis lyrata collected across its entire distribution. The mating system and range dynamics since the last glaciation cycle explained around 60% of the variation in genomic diversity among populations and 52% of the variation in the signature of positive selection. Diversity was lowest in selfing compared with outcrossing populations and in areas further away from glacial refugia. In parallel, reduced positive selection was found in selfing populations and in populations with a longer route of postglacial range expansion. The signature of positive selection was also reduced in populations without admixture. We conclude that recent range expansion can have a profound influence on diversity in coding and non-coding DNA, similar in magnitude to the shift toward selfing. Distribution limits may in fact be caused by reduced effective population size and compromised positive selection in recently colonized parts of the range.

Deep versus shallow sources of CO2 and Rn from a multi-parametric approach: the case of the Nisyros caldera (Aegean Arc, Greece).

Estimating the quantity of CO2 diffusively emitted from the Earth's surface has important implications for volcanic surveillance and global atmospheric CO2 budgets. However, the identification and quantification of non-hydrothermal contributions to CO2 release can be ambiguous. Here, we describe a multi-parametric approach employed at the Nisyros caldera, Aegean Arc, Greece, to assess the relative influence of deep and shallow gases released from the soil. In April 2019, we measured diffuse soil surface CO2 fluxes, together with their carbon isotope compositions, and at a depth of 80 cm, the CO2 concentration, soil temperature, and the activities of radon and thoron. The contributions of deep CO2 and biogenic CO2 fluxes were distinguished on the basis of their carbon isotope compositions. A Principal Component Analysis (PCA), performed on the measured parameters, effectively discriminates between a deep- and a shallow degassing component. The total CO2 output estimated from a relatively small testing area was two times higher with respect to that observed in a previous survey (October 2018). The difference is ascribed to variation in the soil biogenic CO2 production, that was high in April 2019 (a wet period) and low or absent in October 2018 (a dry period). Accounting for seasonal biogenic activity is therefore critical in monitoring and quantifying CO2 emissions in volcanic areas, because they can partially- or completely overwhelm the volcanic-hydrothermal signal.

How do volatiles escape their shallow magmatic hearth?

Only a small fraction (approx. 1-20%) of magmas generated in the mantle erupt at the surface. While volcanic eruptions are typically considered as the main exhaust pipes for volatile elements to escape into the atmosphere, the contribution of magma reservoirs crystallizing in the crust is likely to dominate the volatile transfer from depth to the surface. Here, we use multiscale physical modelling to identify and quantify the main mechanisms of gas escape from crystallizing magma bodies. We show that most of the outgassing occurs at intermediate to high crystal fraction, when the system has reached a mature mush state. It is particularly true for shallow volatile-rich systems that tend to exsolve volatiles through second boiling, leading to efficient construction of gas channels as soon as the crystallinity reaches approximately 40-50 vol.%. We, therefore, argue that estimates of volatile budgets based on volcanic activity may be misleading because they tend to significantly underestimate the magmatic volatile flux and can provide biased volatile compositions. Recognition of the compositional signature and volumetric dominance of intrusive outgassing is, therefore, necessary to build robust models of volatile recycling between the mantle and the surface. This article is part of the Theo Murphy meeting issue 'Magma reservoir architecture and dynamics'.

The influence of water in silicate melt on aluminium excess in plagioclase as a potential hygrometer.

Measuring water contents of magmas is fundamental to resolving a number of geological questions, such as the mechanisms of silicic magma evolution, the triggering of volcanic eruptions, and the formation of porphyry copper deposits. This study focuses on the correlation between apparent deviations from stoichiometry of plagioclase crystals and high water concentration in the magmatic melt from which they grew. We considered this relationship as a potential geo-hygrometer (water activity indicator). To test and potentially calibrate this new technique, a range of natural and experimental plagioclase crystals were analysed, with particular care taken to identify and avoid analytical bias and artefacts. In contrast to recently published material, we found no systematic aluminium excess in plagioclase, irrespective of the water concentration of the silicate melt it crystallised from. This suggests that aluminium excess in plagioclase cannot serve as a geo-hygrometer. The high likelihood of misinterpreting analytical artefacts (due to alkali migration and imprecise standardisation) as small deviations from stoichiometry, also requires its application as a mineral exploration tool to be treated with caution.

Controls on explosive-effusive volcanic eruption styles.

One of the biggest challenges in volcanic hazard assessment is to understand how and why eruptive style changes within the same eruptive period or even from one eruption to the next at a given volcano. This review evaluates the competing processes that lead to explosive and effusive eruptions of silicic magmas. Eruptive style depends on a set of feedback involving interrelated magmatic properties and processes. Foremost of these are magma viscosity, gas loss and external properties such as conduit geometry. Ultimately, these parameters control the speed at which magmas ascend, decompress and outgas en route to the surface, and thus determine eruptive style and evolution.

LA-ICP-MS and SIMS U-Pb and U-Th zircon geochronological data of Late Pleistocene lava domes of the Ciomadul Volcanic Dome Complex (Eastern Carpathians).

This article provides laser-ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) and secondary ionization mass spectrometry (SIMS) U-Pb and U-Th zircon dates for crystals separated from Late Pleistocene dacitic lava dome rocks of the Ciomadul Volcanic Dome Complex (Eastern Carpathians, Romania). The analyses were performed on unpolished zircon prism faces (termed rim analyses) and on crystal interiors exposed through mechanical grinding an polishing (interior analyses). 206Pb/238U ages are corrected for Th-disequilibrium based on published and calculated distribution coefficients for U and Th using average whole-rock and individually analyzed zircon compositions. The data presented in this article were used for the Th-disequilibrium correction of (U-Th)/He zircon geochronology data in the research article entitled "The onset of the volcanism in the Ciomadul Volcanic Dome Complex (Eastern Carpathians): eruption chronology and magma type variation" (Molnár et al., 2018) [1].

LA-ICP-MS U-Pb zircon geochronology data of the Early to Mid-Miocene syn-extensional massive silicic volcanism in the Pannonian Basin (East-Central Europe).

This article provides LA-ICP-MS in-situ U-Pb zircon dates performed on single crystals from dacitic to rhyolitic ignimbrites of the Bükkalja Volcanic Field (Hungary, East-Central Europe) temporally covering the main period of the Neogene silicic volcanic activity in the Pannonian Basin. The data include drift-corrected, alpha dose-corrected, Th-disequilibrium-corrected, and filtered data for geochronological use. The data presented in this article are interpreted and discussed in the research article entitled "Early to Mid-Miocene syn-extensional massive silicic volcanism in the Pannonian Basin (East-Central Europe): eruption chronology, correlation potential and geodynamic implications" by Lukács et al. (2018) [1].

Long-term magmatic evolution reveals the beginning of a new caldera cycle at Campi Flegrei.

Understanding the mechanisms that control the accumulation of large silicic magma bodies in the upper crust is key to determining the potential of volcanoes to form caldera-forming eruptions. Located in one of the most populated regions on Earth, Camp Flegrei is an active and restless volcano that has produced two cataclysmic caldera-forming eruptions and numerous smaller eruptive events over the past 60,000 years. Here, we combine the results of an extensive petrological survey with a thermomechanical model to investigate how the magmatic system shifts from frequent, small eruptions to large caldera-forming events. Our data reveal that the most recent eruption of Monte Nuovo is characterized by highly differentiated magmas akin to those that fed the pre-caldera activity and the initial phases of the caldera-forming eruptions. We suggest that this eruption is an expression of a state shift in magma storage conditions, whereby substantial amounts of volatiles start to exsolve in the shallow reservoir. The presence of an exsolved gas phase has fundamental consequences for the physical properties of the reservoir and may indicate that a large magma body is currently accumulating underneath Campi Flegrei.

Insights into the Vitis complex in the Danube floodplain (Austria).

European grapevine populations quickly disappeared from most of their range, massively killed by the spread of North American grapevine pests and diseases. Nowadays taxonomic pollution represents a new threat. A large Vitis complex involves escaped cultivars, rootstocks, and wild grapevines. The study aimed to provide insight into the Vitis complex in the Danube region through field and genetic analyses. Among the five other major rivers in Europe which still host wild grapevine populations, the Danube floodplain is the only one benefiting from an extensive protected forest area (93 km²) and an relatively active dynamic flood pulse. The Donau-Auen National Park also regroups the largest wild grapevine population in Europe. Ninety-two percent of the individuals collected in the park were true wild grapevines, and 8% were hybrids and introgressed individuals of rootstocks, wild grapevines, and cultivars. These three groups are interfertile acting either as pollen donor or receiver. Hybrids were established within and outside the dykes, mostly in anthropized forest edges. The best-developed individuals imply rootstock genes. They establish in the most erosive parts of the floodplain. 42% of the true wild grapevines lived at the edges of forest/meadow, 33.3% at the edges forest/channels, and 23.9% in forest gaps. DBH (Diameter Breast Height) varied significantly with the occurrence of flooding. Clones were found in both true wild and hybrids/introgressed grapevines. The process of cloning seemed to be prevented in places where flooding dynamics is reduced. The current global distribution of true wild grapevines shows a strong tendency toward clustering, in sites where forestry practices were the most extensive. However, the reduced flooding activity is a danger for long-term sustainability of the natural wild grapevine population.

Wild European apple (Malus sylvestris (L.) Mill.) population dynamics: insight from genetics and ecology in the Rhine Valley. Priorities for a future conservation programme.

The increasing fragmentation of forest habitats and the omnipresence of cultivars potentially threaten the genetic integrity of the European wild apple (Malus sylvestris (L.) Mill). However, the conservation status of this species remains unclear in Europe, other than in Belgium and the Czech Republic, where it has been declared an endangered species. The population density of M. sylvestris is higher in the forests of the upper Rhine Valley (France) than in most European forests, with an unbalanced age-structure, an overrepresentation of adults and a tendency to clump. We characterize here the ecology, age-structure and genetic diversity of wild apple populations in the Rhine Valley. We use these data to highlight links to the history of this species and to propose guidelines for future conservation strategies. In total, 255 individual wild apple trees from six forest stands (five floodplain forests and one forest growing in drier conditions) were analysed in the field, collected and genotyped on the basis of data for 15 microsatellite markers. Genetic analyses showed no escaped cultivars and few hybrids with the cultivated apple. Excluding the hybrids, the genetically "pure" populations displayed high levels of genetic diversity and a weak population structure. Age-structure and ecology studies of wild apple populations identified four categories that were not randomly distributed across the forests, reflecting the history of the Rhine forest over the last century. The Rhine wild apple populations, with their ecological strategies, high genetic diversity, and weak traces of crop-to-wild gene flow associated with the history of these floodplain forests, constitute candidate populations for inclusion in future conservation programmes for European wild apple.

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 amplification was conducted in two multiplex groups, each containing six microsatellite loci. Cross-species amplification was tested in seven other species of Sebaea. The 12 novel microsatellite loci amplified only in the two most closely related species to S. aurea (i.e., S. ambigua and S. minutiflora) 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.