Streamflow response to regional climate model output in the mountainous watershed: a case study from the Swiss Alps

Regional climate model (RCM) outputs are often used in hydrological modeling, in particular for streamflow forecasting. The heterogeneity of the meteorological variables such as precipitation, temperature, wind speed and solar radiation often limits the ability of the hydrological model performance. This paper assessed the sensitivity of RCM outputs from the PRUDENCE project and their performance in reproducing the streamflow. The soil and water assessment tool was used to simulate the streamflow of the Rhone River watershed located in the southwestern part of Switzerland, with the climate variables obtained from four RCMs. We analyzed the difference in magnitude of precipitation, maximum and minimum air temperature, and wind speed with respect to the observed values from the meteorological stations. In addition, we also focused on the impact of the grid resolution on model performance, by analyzing grids with resolutions of 50 × 50 and 25 × 25 km². The variability of the meteorological inputs from various RCMs is quite severe in the studied watershed. Among the four different RCMs, the Danish Meteorological Institute provided the best performance when simulating runoff. We found that temperature lapse rate is significantly important in the mountainous snow and glacier dominated watershed as compared to other variables like precipitation, and wind speed for hydrological performance. Therefore, emphasis should be given to minimum and maximum temperature in the bias correction studies for downscaling climatic data for impact modeling in the mountainous snow and glacier dominated complex watersheds.     

Behaviour of radiogenic Pb in zircon during ultrahigh-temperature metamorphism: an ion imaging and ion tomography case study from the Kerala Khondalite Belt,southern India

Zircon crystals from a locally charnockitized Paleoproterozoic high-K metagranite from the Kerala Khondalite Belt (KKB) of southern India have been investigated by high-spatial resolution secondary ion mass spectrometry analysis of U–Th–Pb and rare earth elements (REE), together with scanning ion imaging and scanning ion tomography (depth-profiled ion imaging). The spot analyses constrain the magmatic crystallization age of the metagranite to ca. 1,850 Ma, with ultrahigh-temperature (UHT) metamorphism occurring at ca. 570 Ma and superimposed charnockite formation at ca. 520–510 Ma, while the ion imaging reveals a patchy distribution of radiogenic Pb throughout the zircon cores. Middle- to heavy-REE depletion in ca. 570 Ma zircon rims suggests that these grew in equilibrium with garnet and therefore date the UHT metamorphism in the KKB. The maximum apparent ²⁰⁷Pb/²⁰⁶Pb age obtained from the unsupported radiogenic Pb concentrations is also consistent with formation of the Pb patches during this event. The superimposed charnockitization event appears to have caused additional Pb-loss in the cores and recrystallization of the rims. The results of depth-profiling of the scanning ion tomography image stack show that the Pb-rich domains range in size from <5 nm to several 10 nm (diameter if assumed to be spherical). The occurrence of such patchy Pb has previously been documented only from UHT metamorphic zircon, where it likely results from annealing of radiation-damaged zircon. The formation of a discrete, heterogeneously distributed and subsequently immobile Pb phase effectively arrests the normal Pb-loss process seen at lower grades of metamorphism.     

Petrogenesis of a voluminous Quaternary adakitic volcano: the case of Baru volcano

The origin of adakite magmas remains controversial because initially the term adakite had petrogenetic significance implying an origin by direct melting of the eclogitized subducting oceanic crust. Many models have been produced for their origin, and until now there has not been a straightforward method to discriminate between these models in a given adakite suite. Here, we use detailed chronological and geochemical studies of selected adakitic edifices that allows for the determination of the magmatic output rate parameter (Qe), which has been correlated with the rates of magma generation deep within subduction zones. By providing temporal and eruption rate estimates, we provide constraints on the possible petrogenetic processes involved in the generation of adakite-like signatures. Adakite magmas derived from the melting of the subducting slab should be volumetrically insignificant when compared to the adakite-like magmas produced by typical arc magma generation processes. In this study, we use this observation and the extraordinary stratigraphic exposure from Miocene to present in an adakitic volcano in Panama and to study the temporal and chemical variation in erupted magmas to estimate rates of magma generation. Detailed chemical and geochronological analyses of Baru volcano indicate that the volcanic edifice was constructed in its entirety during the Quaternary and magmas display adakite-like features such as steep rare earth elements patterns, pronounced depletions in the heavy rare earth elements, low Y, high Sr, and high Sr/Y. The magmatic output rates (Qe) that we have calculated show that compared to other typical adakitic volcanoes, most of the volcanic edifice of Baru volcano was constructed extremely rapidly (<~213 k.a.) and in time frames that are similar to typical arc volcanoes. The observed chemical and mineralogical variation, coupled with the high magma production rates, indicate that Baru volcano is more representative of a typical arc volcano than a small-volume melt of the subducting oceanic crust. The technique we outline may have broader application in determining the petrogenetic conditions of other adakite suites.     

Origin and significance of poikilitic and mosaic peridotite xenoliths in the western Pannonian Basin: geochemical and petrological evidences

Peridotite xenoliths erupted by alkali basaltic volcanoes in the western Pannonian Basin can be divided into two fundamentally contrasting groups. Geochemical characteristics of the abundant protogranular, porphyroclastic and equigranular nodules suggest that these samples originate from an old consolidated and moderately depleted lithospheric mantle domain. In contrast, the geochemical features of the worldwide rare, but in the Pannonian Basin relatively abundant, poikilitic xenoliths attest to a more complex evolution. It has been argued that the origin of the peculiar texture and chemistry may be intimately linked to melt/rock reactions at successively decreasing liquid volumes in a porous melt flow system. The most likely site where such reactions can take place is the asthenosphere–lithosphere boundary. In this context, poikilitic xenoliths may provide petrological and geochemical evidence for reactions between magmatic liquids issued from the uprising asthenosphere and the solid mantle rocks of the lithosphere. These reactions are important agents of the thermal erosion of the lithosphere; thus, they could have considerably contributed to the thinning of the lithosphere in the Pannonian region. We suggest that in the Pannonian Basin, there could be a strong relation between the unusual abundance of poikilitic mantle xenoliths and the strongly eroded lithosphere.     

Overstepping the garnet isograd: a comparison of QuiG barometry and thermodynamic modeling

The consequences of overstepping the garnet isograd reaction have been investigated by comparing the composition of garnet formed at overstepped P–T conditions (the overstep or “OS” model) with the P–T conditions that would be inferred by assuming garnet nucleated in equilibrium with the matrix assemblage at the isograd (the equilibrium or “EQ” model). The garnet nucleus composition formed at overstepped conditions is calculated as the composition that produces the maximum decrease in Gibbs free energy from the equilibrated, garnet-absent, matrix assemblage for the bulk composition under study. Isopleths were then calculated for this garnet nucleus composition assuming equilibrium with the matrix assemblage (the EQ model). Comparison of the actual P–T conditions of nucleation (the OS model) with those inferred from the EQ model reveals considerable discrepancy between the two. In general, the inferred garnet nucleation P–T conditions (the EQ model) are at a lower temperature and higher or lower pressure (depending on the coexisting calcic phase(s)) than the actual (OS model) nucleation conditions. Moreover, the degree of discrepancy increases with the degree of overstepping. Independent estimates of the pressure of nucleation of garnet were made using the Raman shift of quartz inclusions in garnet (quartz-in-garnet or QuiG barometry). To test the validity of this method, an experimental synthesis of garnet containing quartz inclusions was made at 800 °C, 20 kbar, and the measured Raman shift reproduced the synthesis conditions to within 120 bars. Raman band shifts from three natural samples were then used to calculate an isochore along which garnet was presumed to have nucleated. Model calculations were made at several temperatures along this isochore (the OS model), and these P–T conditions were compared to those computed assuming equilibrium nucleation (the EQ model) to estimate the degree of overstepping displayed by these samples. A sample from the garnet isograd in eastern Vermont is consistent with overstepping of around 10 degrees and 0.6 kbar (affinities of around 2 kJ/mole garnet). A sample from the staurolite–kyanite zone in the same terrane requires overstepping of around 50 °C and 2–5 kbar (affinities of around 10–18 kJ/mole garnet). A similar amount of overstepping was inferred for a blueschist sample from Sifnos, Greece. These results indicate that overstepping of garnet nucleation reactions may be common and pronounced in regionally metamorphosed terranes, and that the P–T conditions and paths inferred from garnet zoning studies may be egregiously in error.     

Thermal and rheological controls on the formation of mafic enclaves or banded pumice

Magma mixing can occur in a fluid manner to produce banded pumice or in a brittle manner to form enclaves. We propose that the critical control on mixing style is a competition between developing networks of crystals in the intruding magma that impart a strength to the magma and melting and disrupting those networks in the host. X-ray computed tomography analysis demonstrates that banded pumice from the 1915 Mt. Lassen eruption lacks crystal networks. In contrast, rhyodacite hosts with mafic enclaves from Chaos Crags contain well-developed networks of large crystals. We present a one-dimensional conductive cooling model that predicts mixing style, either ductile or brittle, as a function of magma compositions, temperatures, and the size of the intruding dike. Our model relies on three assumptions: (1) Mixing is initiated by the injection of a hot dike into a cooler magma body with a yield strength; (2) when magma crystallinity exceeds a critical value, 13 vol% plagioclase, the magma develops a yield strength; and (3) when total crystallinity exceeds 40 vol%, the magma has a penetrative crystal network and is effectively solid. Importantly, because the two magmas are of different compositions, their crystallinities and viscosities do not have the same variations with temperature. As the intruding magma cools, it crystallizes from the outside in, while simultaneously, host magma temperature near the intruder rises. Mixing of the two magmas begins when the host magma is heated sufficiently to (1) disrupt the crystal network and (2) initiate convection. If the shear stress exerted by the convecting host magma on the dike is greater than the yield strength of the dike margin (and dike crystallinity does not exceed 40 %), then fluid mixing occurs, otherwise enclaves form by brittle deformation of the dike. Application of the model to magma compositions representative of Lassen and Chaos Crags shows that emplacement of dikes <1 m thick should produce enclaves, whereas thicker dikes should generate fluid mixing and form banded pumice within days to weeks of emplacement. Similar relationships apply to other modeled magmatic systems, including Pinatubo, Unzen, and Ksudach/Shtuybel’ volcanoes. For all studied systems, the absolute size of the intruding dike, not just its proportion relative to the host, influences mixing style.     

Primary melt from Sannome-gata volcano,NE Japan arc: constraints on generation conditions of rear-arc magmas

The conditions under which rear-arc magmas are generated were estimated using primary basalts from the Sannome-gata volcano, located in the rear of the NE Japan arc. Scoriae from the volcano occur with abundant crustal and mantle xenoliths, suggesting that the magma ascended rapidly from the upper mantle. The scoriae show significant variations in their whole-rock compositions (7.9–11.1 wt% MgO). High-MgO scoriae (MgO > ~9.5 wt%) have mostly homogeneous ⁸⁷Sr/⁸⁶Sr ratios (0.70318–0.70320), whereas low-MgO scoriae (MgO < ~9 wt%) have higher ⁸⁷Sr/⁸⁶Sr ratios (>0.70327); ratios tend to increase with decreasing MgO content. The high-MgO scoriae are aphyric, containing ~5 vol% olivine microphenocrysts with Mg# [100 × Mg/(Mg + Fe²⁺)] of up to 90. In contrast, the low-MgO scoriae have crustal xenocrysts of plagioclase, alkali feldspar, and quartz, and the mineralogic modes correlate negatively with whole-rock MgO content. On the basis of these observations, it is inferred that the high-MgO scoriae represent primary or near-primary melts, while the low-MgO scoriae underwent considerable interaction with the crust. Using thermodynamic analysis of the observed petrological features of the high-MgO scoriae, the eruption temperature of the magmas was constrained to 1,160–1,220 °C. Given that the source mantle was depleted MORB-source mantle, the primary magma was plausibly generated by ~7 % melting of a garnet-bearing spinel peridotite; taking this into consideration, and considering the constraints of multi-component thermodynamics, we estimated that the primary Sannome-gata magma was generated in the source mantle with 0.5–0.6 wt% H₂O at 1,220–1,230 °C and at ~1.8 GPa, and that the H₂O content of the primary magma was 6–7 wt%. The rear-arc Sannome-gata magma was generated by a lower degree of melting of the mantle at greater depths and lower temperatures than the frontal-arc magma from the Iwate volcano, which was also estimated to be generated by ~15 % melting of the source mantle with 0.6–0.7 wt% H₂O at ~1,250 °C and at ~1.3 GPa.     

Settling and compaction of chromite cumulates employing a centrifuging piston cylinder and application to layered mafic intrusions

The time scales and mechanics of gravitationally driven crystal settling and compaction is investigated through high temperature (1,280–1,500 °C) centrifuge-assisted experiments on a chromite-basalt melt system at 100–1,500g (0.5 GPa). Subsequently, the feasibility of this process for the formation of dense chromite cumulate layers in large layered mafic intrusions (LMIs) is assessed. Centrifugation leads to a single cumulate layer formed at the gravitational bottom of the capsule. The experimentally observed mechanical settling velocity of a suspension of ~24 vol% chromite is calculated to be about half (~0.53) of the Stokes settling velocity, with a sedimentation exponent n of 2.35 (3). Gravitational settling leads to an orthocumulate layer with a porosity of 0.52 (all porosities as fraction). Formation times for such a layer from a magma with initial chromite contents of 0.1–1 vol% are 140–3.5 days, equal to a growth rate of 0.007–0.3 m/day for grain sizes of 1–2 mm. More compacted chromite layers form with increasing centrifugation time and acceleration through chemical compaction: An increase of grain contact areas and grain sizes together with a decrease in porosity is best explained by pressure dissolution at grain contacts, reprecipitation and grain growth into the intergranular space and a concomitant expulsion of intergranular melt. The relation between the porosity in the cumulate pile and effective pressure integrated over time (Δρ · h · a · t) is best fit with a logarithmic function, in fact confirming that a (pressure) dissolution–reprecipitation process is the dominant mechanism of compaction. The experimentally derived equation allows calculating compaction times: 70–80 % chromite at the bottom of a 1-m-thick chromite layer are reached after 9–250 years, whereas equivalent compaction times are 0.2–0.9 years for olivine (both for 2 mm grain size). The experiments allow to determine the bulk viscosities of chromite and olivine cumulates to be of magnitude 10⁹ Pa s, much lower than previously reported. As long as melt escape from the compacting cumulate remains homogeneous, fluidization does not play any role; however, channelized melt flow may lead to suspension and upward movement of cumulate crystals. In LMIs, chromitite layers are typically part of a sequence with layers of mafic minerals, compaction occurs under the additional weight of the overlying layers and can be achieved in a few years to decades.     

Zircon evidence for a ~200 k.y. supereruption-related thermal flare-up in the Miocene southern Black Mountains,western Arizona,USA

The Silver Creek caldera (southern Black Mountains, western Arizona) is the source of the 18.8 Ma, >700 km³ Peach Spring Tuff (PST) supereruption, the largest eruption generated in the Colorado River Extensional Corridor (CREC) of the southwestern United States. Within and immediately surrounding the caldera is a sequence of volcanics and intrusions ranging in age from ~19 to 17 Ma. These units offer a record of magmatic processes prior to, during, and immediately following the PST eruption. To investigate the thermal evolution of the magmatic center that produced the PST, we applied a combination of Ti-in-zircon thermometry, zircon saturation thermometry, and high-precision U–Pb CA–TIMS zircon dating to representative pre- and post-supereruption volcanic and intrusive units from the caldera and its environs. Similar to intracaldera PST zircons, zircons from a pre-PST trachytic lava (19 Ma) and a post-PST caldera intrusion (18.8 Ma) yield exceptionally high-Ti concentrations (most >20 ppm, some up to nearly 60 ppm), corresponding to calculated temperatures that exceed 900 °C. In these units, Ti-in-zircon temperatures typically surpass zircon saturation temperatures (ZSTs), suggesting the entrainment of zircon that had grown in hotter environments within the magmatic system. Titanium concentrations in younger volcanic and intrusive units (~18.7–17.5 Ma) decline through time, corresponding to an average cooling rate of 10^?3.5 °C/year. The ~200 k.y. thermal peak evident at Silver Creek caldera is spatially limited: elsewhere in the Miocene record of the northern CREC, Ti-in-zircon concentrations and ZSTs are much lower, suggesting that felsic magmas were generally substantially cooler.     

Deciphering petrogenic processes using Pb isotope ratios from time-series samples at Bezymianny and Klyuchevskoy volcanoes,Central Kamchatka Depression

The Klyuchevskoy group of volcanoes in the Kamchatka arc erupts compositionally diverse magmas (high-Mg basalts to dacites) over small spatial scales. New high-precision Pb isotope data from modern juvenile (1956–present) erupted products and hosted enclaves and xenoliths from Bezymianny volcano reveal that Bezymianny and Klyuchevskoy volcanoes, separated by only 9 km, undergo varying degrees of crustal processing through independent crustal columns. Lead isotope compositions of Klyuchevskoy basalts–basaltic andesites are more radiogenic than Bezymianny andesites (²⁰⁸Pb/²⁰⁴Pb = 37.850–37.903, ²⁰⁷Pb/²⁰⁴Pb = 15.468–15.480, and ²⁰⁶Pb/²⁰⁴Pb = 18.249–18.278 at Bezymianny; ²⁰⁸Pb/²⁰⁴Pb = 37.907–37.949, ²⁰⁷Pb/²⁰⁴Pb = 15.478–15.487, and ²⁰⁶Pb/²⁰⁴Pb = 18.289–18.305 at Klyuchevskoy). A mid-crustal xenolith with a crystallization pressure of 5.2 ± 0.6 kbars inferred from two-pyroxene geobarometry and basaltic andesite enclaves from Bezymianny record less radiogenic Pb isotope compositions than their host magmas. Hence, assimilation of such lithologies in the middle or lower crust can explain the Pb isotope data in Bezymianny andesites, although a component of magma mixing with less radiogenic mafic recharge magmas and possible mantle heterogeneity cannot be excluded. Lead isotope compositions for the Klyuchevskoy Group are less radiogenic than other arc segments (Karymsky—Eastern Volcanic Zone; Shiveluch—Northern Central Kamchatka Depression), which indicate increased lower-crustal assimilation beneath the Klyuchevskoy Group. Decadal timescale Pb isotope variations at Klyuchevskoy demonstrate rapid changes in the magnitude of assimilation at a volcanic center. Lead isotope data coupled with trace element data reflect the influence of crustal processes on magma compositions even in thin mafic volcanic arcs.