On the reliability of the Th-U dating method applied to young volcanic rocks — reply

Hemond and Condomines (1985) claim that our results on zero age volcanic materials of Vesuvius and Etna are incorrect, that no disequilibrium exists in those materials and that the validity of the ²³⁰Th-²³⁸U dating method is thus confirmed. We show that other researchers obtained similar results in zero age rocks. The validity of the ²³⁰Th-²³⁸U dating method is simply ruled out by a comparison of the ages obtained by this method with those obtained on the same samples by different methods.     

Origin and timing of the post-Variscan gabbro–granite complex of Porto (Western Corsica)

The post-Variscan complex of Porto consists of metaluminous to slightly peraluminous A-type biotite granites mingled with gabbro-dioritic rocks, and late dykes with basaltic to trachyandesitic composition. U-Pb zircon dating by LA-ICP-MS on two mafic intrusive samples constrains the time of the gabbro–granite crystallisation at 281 ± 3 Ma and 283 ± 2 Ma. Hornblende ⁴⁰Ar-³⁹Ar ages from a late trachyandesite dyke date the dyking event at 280 ± 2 Ma, which is within error the U-Pb zircon ages of the intrusives. Biotite granites show variable major and trace element compositions and similar initial ε_Nd (−0.3 to +0.9). Whole rock chemistry variations and trace element compositions of plagioclase and allanite indicate that the granites are genetically linked, essentially through fractional crystallisation of feldspars and minor allanite. On the basis of whole-rock chemistry e.g. initial ε_Nd +4.9 to +1.7 and trace element clinopyroxene compositions, we have ascertained that the mafic intrusives and basic dykes formed from isotopically depleted mantle source-derived melts with similar trace element signature. These basic melts experienced slightly different evolutionary histories, controlled by fractional crystallisation and crustal contamination, mainly by the acid magma that gave rise to the associated biotite granites, but also by the enclosing older Variscan granitoids. U-Pb zircon data suggest that the Porto complex was affected by hydrothermal fluid circulation at 259 ± 9 Ma.     

Petrologic evidence of a complex plumbing system feeding the July–August 2001 eruption of Mt. Etna,Sicily, Italy

After the major 1991–1993 eruption, Mt. Etna resumed flank activity in July 2001 through a complex system of eruptive fissures cutting the NE and the S flanks of the volcano and feeding effusive activity, fire fountains, Strombolian and minor phreatomagmatic explosions. Throughout the eruption, magmas with different petrography and composition were erupted. The vents higher than 2,600 m a.s.l. (hereafter Upper vents, UV) erupted porphyritic, plagioclase-rich trachybasalt, typical of present-day summit and flank activity. Differently, the vents located at 2,550 and 2,100 m a.s.l. (hereafter Lower vents, LV) produced slightly more primitive trachybasalt dominated by large clinopyroxene, olivine and uncommon minerals for Etna such as amphibole, apatite and orthopyroxene and containing siliceous and cognate xenoliths. Petrologic investigations carried out on samples collected throughout the eruption provided insights into one of the most intriguing aspects of the 2001 activity, namely the coeval occurrence of distinct magmas. We interpret this evidence as the result of a complex plumbing system. It consists in two separate magma storage systems: a shallow one feeding the activity of the UV and a deeper and more complex storage related to the activity of LV. In this deep storage zone, which is thermally and compositionally zoned, the favourable conditions allow the crystallization of amphibole and the occurrence of cognate xenoliths representing wall cumulates. Throughout 2001 eruption, UV and LV magmas remain clearly distinct and ascended following different paths, ruling out the occurrence of mixing processes between them. Furthermore, integrating the 2001 eruption in the framework of summit activity occurring since 1995, we propose that the 2001 magma feeding the vents lower than 2,600 m a.s.l. is a precursor of a refilling event, which reached its peak during the 2002–2003 Etna flank eruption.     

Lead isotope variations across terrane boundaries of the Tien Shan and Chinese Altay

The Altaid orogen was formed by aggregation of Paleozoic subduction–accretion complexes and Precambrian basement blocks between the Late Proterozoic and the Early Mesozoic. Because the Altaids are the site of abundant granitic plutonism and host some of the largest gold deposits in the world, understanding their formation has important implications on the comprehension of Phanerozoic crustal growth and metallogeny. In this study, we present the first extensive lead isotope data on magmatic and metasedimentary rocks as well as ore deposits of the southern part of the Altaids, including the Tien Shan (Tianshan) and southern Altay (Altai) orogenic belts. Our results show that each terrane investigated within the Tien Shan and southern Altay is characterized by a distinct Pb isotope signature and that there is a SW–NE Pb isotope gradient suggesting a progressive transition from a continental crust environment in the West (the Kyzylkum and Kokshaal segments of the Southern Tien Shan) to an almost 100% juvenile (MORB-type mantle-derived) crust environment in the East (Altay). The Pb isotope signatures of the studied ore deposits follow closely those of magmatic and metasedimentary rocks of the host terranes, thus supporting the validity of lead isotopes to discriminate terranes. Whereas this apparently suggests that no unique reservoir has been responsible for the huge gold concentration in this region, masking of a preferential Pb-poor Au-bearing reservoir by mixing with Pb-rich crustal reservoirs during the mineralizing events cannot be excluded.     

Experimental Sensitivity Analysis of Oxygen Transfer in the Capillary Fringe

Oxygen transfer in the capillary fringe (CF) is of primary importance for a wide variety of biogeochemical processes occurring in shallow groundwater systems. In case of a fluctuating groundwater table two distinct mechanisms of oxygen transfer within the capillary zone can be identified: vertical predominantly diffusive mass flux of oxygen, and mass transfer between entrapped gas and groundwater. In this study, we perform a systematic experimental sensitivity analysis in order to assess the influence of different parameters on oxygen transfer from entrapped air within the CF to underlying anoxic groundwater. We carry out quasi two‐dimensional flow‐through experiments focusing on the transient phase following imbibition to investigate the influence of the horizontal flow velocity, the average grain diameter of the porous medium, as well as the magnitude and the speed of the water table rise. We present a numerical flow and transport model that quantitatively represents the main mechanisms governing oxygen transfer. Assuming local equilibrium between the aqueous and the gaseous phase, the partitioning process from entrapped air can be satisfactorily simulated. The different experiments are monitored by measuring vertical oxygen concentration profiles at high spatial resolution with a noninvasive optode technique as well as by determining oxygen fluxes at the outlet of the flow‐through chamber. The results show that all parameters investigated have a significant effect and determine different amounts of oxygen transferred to the oxygen‐depleted groundwater. Particularly relevant are the magnitude of the water table rise and the grain size of the porous medium.     

Gradual Fault Weakening with Seismic Slip: Inferences from the Seismic Sequences of L’Aquila, 2009, and Northridge, 1994

We estimate seismological fracture energies from two subsets of events selected from the seismic sequences of L’Aquila (2009), and Northridge (1994): 57 and 16 selected events, respectively, including the main shocks. Following Abercrombie and Rice (Geophys J Int 162: 406–424, 2005), we postulate that fracture energy (G) represents the post-failure integral of the dynamic weakening curve, which is described by the evolution of shear traction as a function of slip. Following a direct-wave approach, we compute mainshock-/aftershock-source spectral ratios, and analyze them using the approach proposed by Malagnini et al. (Pure Appl. Geophys., this issue, 2014) to infer corner frequencies and seismic moment. Our estimates of source parameters (including fracture energies) are based on best-fit grid-searches performed over empirical source spectral ratios. We quantify the source scaling of spectra from small and large earthquakes by using the MDAC formulation of Walter and Taylor (A revised Magnitude and Distance Amplitude Correction (MDAC2) procedure for regional seismic discriminants, 2001). The source parameters presented in this paper must be considered as point-source estimates representing averages calculated over specific ruptured portions of the fault area. In order to constrain the scaling of fracture energy with coseismic slip, we investigate two different slip-weakening functions to model the shear traction as a function of slip: (i) a power law, as suggested by Abercrombie and Rice (Geophys J Int 162: 406–424, 2005), and (ii) an exponential decay. Our results show that the exponential decay of stress on the fault allows a good fit between measured and predicted fracture energies, both for the main events and for their aftershocks, regardless of the significant differences in the energy budgets between the large (main) and small earthquakes (aftershocks). Using the power-law slip-weakening function would lead us to a very different situation: in our two investigated sequences, if the aftershock scaling is extrapolated to events with large slips, a power law (a la Abercrombie and Rice) would predict unrealistically large stress drops for large, main earthquakes. We conclude that the exponential stress evolution law has the advantage of avoiding unrealistic stress drops and unbounded fracture energies at large slip values, while still describing the abrupt shear-stress degradation observed in high-velocity laboratory experiments (e.g., Di Toro et al., Fault lubrication during earthquakes, Nature 2011).     

Geochemical heterogeneities and dynamics of magmas within the plumbing system of a persistently active volcano: evidence from Stromboli

We report here the most complete dataset for major and trace elements, as well as Sr isotopic compositions, of magmas erupted by Stromboli since the onset of present-day activity 1,800 years ago. Our data relate to both porphyritic scoria and lava originating in the uppermost parts of the feeding system, plus crystal-poor pumice produced by paroxysmal explosive eruption of deep-seated, fast ascending, magma. The geochemical variations recorded by Stromboli’s products allow us to identify changes in magma dynamics affecting the entire plumbing system. Deep-seated magmas vary in composition between two end-members having different key ratios in strongly incompatible trace elements and Sr isotopes. These features may be ascribed to mantle source processes (fluid/melt enrichment, variable degrees of melting) and occasional contamination by deep, mafic, cumulates. Temporal trends reveal three phases during which magmas with distinct geochemical signatures were erupted. The first phase occurred between the third and fourteenth centuries AD and was characterised by the eruption of evolved magmas sharing geochemical and Sr isotopic compositions similar to those of earlier periods of activity (<12 ka—Neostromboli and San Bartolo). The second phase, which began in the sixteenth century and lasted until the first half of the twentieth century, produced more primitive, less radiogenic, magmas with the lowest Ba/La and Rb/Th ratios of our dataset. The last phase is ongoing and is marked by a magma having the lowest Sr isotopic composition and highest Rb/Th ratio of the dataset. While this new magma can be clearly identified in the pumice erupted during the last two paroxysmal eruptions of 2003 and 2007, shallow degassed magma extruded during this time span records significant geochemical and isotopic heterogeneities. We thus suggest that the shallow reservoir has been only partially homogenised by this new magma influx. We conclude that compositional variations within the shallow magma system of a persistently active volcano provide only a biassed signal of ongoing geochemical changes induced by deep magma refilling. We argue that source changes can only be identified by interpreting the geochemistry of pumice, because it reliably represents magma transferred directly from deep portions of the plumbing system to the surface.