Erupted cumulate fragments in rhyolites from Lipari (Aeolian Islands)

Over the last ~267 ky, the island of Lipari has erupted magmas ranging in compositions from basaltic andesites to rhyolites, with a notable compositional gap in the dacite field. Bulk geochemical and isotopic compositions of the volcanic succession, in conjunction with major and trace elemental compositions of minerals, indicate that the rhyolites were dominantly generated via crystal fractionation processes, with subordinate assimilation. Radiogenic (Sr, Nd, and Pb) and stable (O) isotopes independently suggest ≤30 % of crustal contamination with the majority of it occurring in mafic compositions, likely relatively deep in the system. Within the rhyolites, crystal-rich, K₂O-rich enclaves are common. In contrast to previous interpretations, we suggest that these enclaves represent partial melting, remobilization and eruption of cumulate fragments left-over from rhyolite melt extraction. Cumulate melting and remobilization is supported by the presence of (1) resorbed, low-temperature minerals (biotite and sanidine), providing the potassic signature to these clasts, (2) reacted Fo-rich olivine, marking the presence of mafic recharge, (3) An_38–21 plagioclase, filling the gap in feldspar composition between the andesites and the rhyolites and (4) strong enrichment in Sr and Ba in plagioclase and sanidine, suggesting crystallization from a locally enriched melt. Based on Sr-melt partitioning, the high-Sr plagioclase would require ~2300 ppm Sr in the melt, a value far in excess of Sr contents in Lipari and Vulcano magmas (50–1532 ppm) but consistent with melting of a feldspar-rich cumulate. Due to the presence of similar crystal-rich enclaves within the rhyolites from Vulcano, we propose that the eruption of remobilized cumulates associated with high-SiO₂ rhyolites may be a common process at the Aeolian volcanoes, as already attested for a variety of volcanic systems around the world.     

An individual-based population dynamic model for estimating biomass yield and nutrient fluxes through an off-shore mussel (Mytilus galloprovincialis) farm

The fluxes of carbon, nitrogen and phosphorus through an off-shore long-line Mytilus galloprovincialis farm during a typical rearing cycle were estimated by combining a simple population dynamic model, based on a new individual model, and a set of field data, concerning the composition of the seston, as well as that of mussel meat and faeces. The individual model, based on an energy budget, was validated against a set of original field data, which were purposely collected from July 2006 to May 2007 in the North-Western Adriatic Sea (Italy) and was further tested using historical data. The model was upscaled to the population level by means of a set of Monte Carlo simulations, which were used for estimating the size structure of the population. The daily fluxes of C, N and P associated with mussel filtration, excretion and faeces and pseudo-faeces production were integrated over the 10-month-long rearing cycle and compared with the total amount of C, N and P removed by harvesting. The results indicate that the individual model compares well with an existing literature model and provides reliable estimations of the growth of mussel specimen over a range of trophic conditions which are typical of the Northern Adriatic Sea coastal area. The results of the budget calculation indicate that, even though the harvest represents a net removal of phosphorus and nitrogen from the ecosystem, the mussel farm increases the retention time of both nutrients in the coastal area, via the deposition of faeces and pseudo-faeces on the sea-bed. In fact, the amount of nitrogen associated with deposition is approximately twice the harvested one and the amount of phosphorus is approximately five times higher. These findings are in qualitative agreement with the results of literature budget and model calculations carried out in a temperate coastal embayment. This agreement suggests that the proper assessment of the overall effect of long-line mussel farming on both the benthic and pelagic ecosystem asks for an integrated modelling approach, which should include the dynamic of early diagenesis processes, as well as of that of nutrients released from the surface sediment.     

Assessment of rockfall susceptibility by integrating statistical and physically-based approaches

In Val di Fassa (Dolomites, Eastern Italian Alps) rockfalls constitute the most significant gravity-induced natural disaster that threatens both the inhabitants of the valley, who are few, and the thousands of tourists who populate the area in summer and winter.To assess rockfall susceptibility, we developed an integrated statistical and physically-based approach that aimed to predict both the susceptibility to onset and the probability that rockfalls will attain specific reaches. Through field checks and multi-temporal aerial photo-interpretation, we prepared a detailed inventory of both rockfall source areas and associated scree-slope deposits. Using an innovative technique based on GIS tools and a 3D rockfall simulation code, grid cells pertaining to the rockfall source-area polygons were classified as active or inactive, based on the state of activity of the associated scree-slope deposits. The simulation code allows one to link each source grid cell with scree deposit polygons by calculating the trajectory of each simulated launch of blocks. By means of discriminant analysis, we then identified the mix of environmental variables that best identifies grid cells with low or high susceptibility to rockfalls. Among these variables, structural setting, land use, and morphology were the most important factors that led to the initiation of rockfalls.We developed 3D simulation models of the runout distance, intensity and frequency of rockfalls, whose source grid cells corresponded either to the geomorphologically-defined source polygons (geomorphological scenario) or to study area grid cells with slope angle greater than an empirically-defined value of 37° (empirical scenario). For each scenario, we assigned to the source grid cells an either fixed or variable onset susceptibility; the latter was derived from the discriminant model group (active/inactive) membership probabilities.Comparison of these four models indicates that the geomorphological scenario with variable onset susceptibility appears to be the most realistic model. Nevertheless, political and legal issues seem to guide local administrators, who tend to select the more conservative empirically-based scenario as a land-planning tool.     

NanoSr – A New Carbonate Microanalytical Reference Material for In Situ Strontium Isotope Analysis

The in situ measurement of Sr isotopes in carbonates by MC‐ICP‐MS is limited by the availability of suitable microanalytical reference materials (RMs), which match the samples of interest. Whereas several well‐characterised carbonate reference materials for Sr mass fractions > 1000 µg g^?1 are available, there is a lack of well‐characterised carbonate microanalytical RMs with lower Sr mass fractions. Here, we present a new synthetic carbonate nanopowder RM with a Sr mass fraction of ca. 500 µg g^?1 suitable for microanalytical Sr isotope research (‘NanoSr’). NanoSr was analysed by both solution‐based and in situ techniques. Element mass fractions were determined using EPMA (Ca mass fraction), as well as laser ablation and solution ICP‐MS in different laboratories. The ⁸⁷Sr/⁸⁶Sr ratio was determined by well‐established bulk methods for Sr isotope measurements and is 0.70756 ± 0.00003 (2s). The Sr isotope microhomogeneity of the material was determined by LA‐MC‐ICP‐MS, which resulted in ⁸⁷Sr/⁸⁶Sr ratios of 0.70753 ± 0.00007 (2s) and 0.70757 ± 0.00006 (2s), respectively, in agreement with the solution data within uncertainties. Thus, this new reference material is well suited to monitor and correct microanalytical Sr isotope measurements of low‐Sr, low‐REE carbonate samples. NanoSr is available from the corresponding author.     

Temperature and density of the Tyrrhenian lithosphere and slab and new interpretation of gravity field in the Tyrrhenian Basin

The gravity anomaly field of the Tyrrhenian basin and surrounding regions reflects the complex series of geodynamic events active in this area since the Oligocene–Miocene. They can resume in lithospheric thinning and asthenospheric rising beneath the Tyrrhenian Basin, coexisting with the roll-back subduction of the African plate margin westward sinking beneath the Calabrian Arc. The geographic closeness between these processes implies an intense perturbation of the mantle thermal regime and an interference at regional scale between the related gravity effects.A model of the litho-asthenospheric structure of this region is suggested, showing a reasonable agreement with both the evidences in terms of regional gravity anomaly pattern and the results concerning thermal state and petro-physical features of the mantle. The first phase of this study consisted of the computation of the isotherms in the crust–mantle system beneath the Tyrrhenian Basin and, afterwards, of the density distribution within the partially melted upwelling asthenosphere. The second phase consisted of a temperature/density modelling of the slab subducting beneath the Calabrian Arc. Finally, a 2^1 / 2 interpretation of gravity data was carried out by including as constraints the results previously obtained. Thus, the final result depicts a model matching both gravity, thermal and petrographic data. They provide (a) a better definition of the thermal regime of the passive mantle rise beneath the Tyrrhenian basin by means of the estimation of the moderate asthenospheric heating and (b) a model of lithospheric slab subducting with rates that could be smaller than generally suggested in previous works.     

Patterns of fluid flow in the contact aureole of the Late Miocene Monte Capanne pluton (Elba Island,Italy): the role of structures and rheology

Fluid–rock interaction was investigated in the inner aureole of the Late Miocene Monte Capanne pluton on Elba Island (Tuscany, central Italy) by integrating structural, petrological, fluid inclusion, and stable isotope analyses. In the north-western sector of the aureole (Procchio–Spartaia area), calc–silicates alternate with nearly pure carbonate layers at the metre scale. Close to the pluton, the prograde metamorphic sequence includes calc–silicates that transition within a few metres to overlying nearly pure calcite marbles. The calc–silicates are extensively metasomatised to form massive wollastonite-grossular-bearing exoskarn. The mineralogical assemblage found in the marbles and the unshifted carbon and oxygen isotopic ratios in calcite attest that the fluid phase was internally buffered. On the other hand, the calc–silicates constituted channels for infiltration of disequilibrium fluids of magmatic origin. Fluid infiltration was enhanced by hydrofracturing and structurally-controlled by existing planar anisotropies in calc–silicates (layering and lithological boundaries). At the metamorphic peak (∼600°C and 1.5–2 kbar), the marble–calc–silicate interface acted as a barrier to fluids exsolved from the crystallising intrusions, separating two different flow patterns in the inner aureole: a high fluid–flux region on its higher grade side (Wol-zone) and a low fluid–flux region on the lower-grade side (Cpx zone). Results of this study: (1) documented that fluid pathways in the aureole rocks at the top of the pluton were largely horizontal, controlled by the lithological layering and the pluton–host rock contact; and (2) elucidated the primary control exerted by the structural and rheological properties of the host rocks on the geometry of fluid flow during pluton emplacement.     

Mineralogical and geochemical study of granular xenoliths from the Alban Hills volcano,Central Italy: bearing on evolutionary processes in potassic magma chambers

Granular xenoliths (ejecta) from pyroclastic deposits emplaced during the latest stages of activity of the Alban Hills volcano range from ultramafic to salic. Ultramafic types consist of various proportions of olivine, spinel, clinopyroxene and phlogopite. They show low SiO₂, alkalies and incompatible element abundances and very high MgO. However, Cr, Co and Sc are anomalously low, at a few ppm level. Olivine is highly magnesian (up to Fo%=96) and has rather high CaO (1% Ca) and very low Ni (around a few tens ppm) contents. These characteristics indicate a genesis of ultramafic ejecta by thermal metamorphism of a siliceous dolomitic limestone, probably with input of chemical components from potassic magma. The other xenoliths have textures and compositional characteristics which indicate that they represent either intrusive equivalents of lavas or cumulates crystallized from variably evolved ultrapotassic magmas. One sample of the former group has major element composition resembling ultrapotassic rocks with kamafugitic affinity. Some cumulitic rocks have exceedingly high abundances of Th (81–84 ppm) and light rare-earth elements (LREE) (La+Ce=421–498 ppm) and extreme REE fractionation (La/Yb=288–1393), not justified by their modal mineralogy which is dominated by sanidine, leucite and nepheline. Finegrained phases are dispersed through the fractures and within the interstices of the main minerals. Semiquantitative EDS analyses show that Th and LREE occur at concentration levels of several tens of percent in these phases, indicating that their presence is responsible for the high concentration of incompatible trace elements in the whole rocks. The interstitial position of these phases and their association with fluorite support a secondary origin by deposition from fluorine-rich fluids separated from a highly evolved potassic liquid. The Nd isotopic ratios of the cjecta range from 0.51182 to 0.51217. ⁸⁷Sr/⁸⁶Sr ratios range from 0.70900 to 0.71036. With the exception of one sample, these values are lower than those of the outcropping lavas, which cluster around 0.7105±3. This indicates either the occurrence of several isotopically distinct potassic magmas or a variable interaction between magmas and wall rocks. However, this latter hypothesis requires selective assimilation of host rocks in order to explain isotopic and geochemical characteristics of lavas and xenoliths. The new data indicate that the evolutionary processes in the potassic magmas of the Alban Hills were much more complex than envisaged by previous studies. Interaction of magmas with wall rocks may be an important process during magmatic evolution. Element migration by gaseous transfer, often invoked but rarely constrained by sound data, is shown to have occurred during the latest stages of magmatic evolution. Such a process was able to produce selective enrichment of Th, U, LREE and, to a minor degree, Ta and Hf in the wall rocks of potassic magma chamber. Finally, the occurrence of xenoliths with kamafugitic composition points to the existence of this type of ultrapotassic magma at the Alban Hills.     

Comparing predictive capability of statistical and deterministic methods for landslide susceptibility mapping: a case study in the northern Apennines (Reggio Emilia Province, Italy)

Statistical and deterministic methods are widely used in geographic information system based landslide susceptibility mapping. This paper compares the predictive capability of three different models, namely the Weight of Evidence, the Fuzzy Logic and SHALSTAB, for producing shallow earth slide susceptibility maps, to be included as informative layers in land use planning at a local level. The test site is an area of about 450 km² in the northern Apennines of Italy where, in April 2004, rainfall combined with snowmelt triggered hundreds of shallow earth slides that damaged roads and other infrastructure. An inventory of the landslides triggered by the event was obtained from interpretation of aerial photos dating back to May 2004. The pre-existence of mapped landslides was then checked using earlier aerial photo coverage. All the predictive models were run on the same set of geo-environmental causal factors: soil type, soil thickness, land cover, possibility of deep drainage through the bedrock, slope angle, and upslope contributing area. Model performance was assessed using a threshold-independent approach (the ROC plot). Results show that global accuracy is as high as 0.77 for both statistical models, while it is only 0.56 for SHALSTAB. Besides the limited quality of input data over large areas, the relatively poorer performance of the deterministic model maybe also due to the simplified assumptions behind the hydrological component (steady-state slope parallel flow), which can be considered unsuitable for describing the hydrologic behavior of clay slopes, that are widespread in the study area.     

Harkerite from the Alban Hills,Italy

Harkerite, found in metamorphic ejecta of the Alban Hills associated with cuspidine, grossular, phlogopite, vesuvian, biotite, and minor amounts of diopside, aegirinaugite, leucite, magnetite and calcito, shows cubic Laue symmetry m3m, possible space groups Fm3m, F432, F43m, a₀ = 14.82 Å.On the basis of isomorphous replacements suggested by the crystal structure analysis, chemical data may be represented by the formula: Ca₄₈Mg₁₆(AlSi₄O₁₆)₄(BO₃)₁₂(CO₃)₂₀ · 4H₂O. Refractive index n_D = 1.6490.The relations between harkerite from Albano, harkerite from Skye and other known harkerites and sakhaites are discussed.