Achieving a more realistic assessment of rockfall hazards by coupling three-dimensional process models and field-based tree-ring data

Sound knowledge of the spatial and temporal patterns of rockfalls is fundamental for the management of this very common hazard in mountain environments. Process-based, three-dimensional simulation models are nowadays capable of reproducing the spatial distribution of rockfall occurrences with reasonable accuracy through the simulation of numerous individual trajectories on highly-resolved digital terrain models. At the same time, however, simulation models typically fail to quantify the ‘real’ frequency of rockfalls (in terms of return intervals). The analysis of impact scars on trees, in contrast, yields real rockfall frequencies, but trees may not be present at the location of interest and rare trajectories may not necessarily be captured due to the limited age of forest stands. In this article, we demonstrate that the coupling of modeling with tree-ring techniques may overcome the limitations inherent to both approaches. Based on the analysis of 64 cells (40 m × 40 m) of a rockfall slope located above a 1631-m long road section in the Swiss Alps, we illustrate results from 488 rockfalls detected in 1260 trees. We illustrate that tree impact data cannot only be used (i) to reconstruct the real frequency of rockfalls for individual cells, but that they also serve (ii) the calibration of the rockfall model Rockyfor3D, as well as (iii) the transformation of simulated trajectories into real frequencies. Calibrated simulation results are in good agreement with real rockfall frequencies and exhibit significant differences in rockfall activity between the cells (zones) along the road section. Real frequencies, expressed as rock passages per meter road section, also enable quantification and direct comparison of the hazard potential between the zones. The contribution provides an approach for hazard zoning procedures that complements traditional methods with a quantification of rockfall frequencies in terms of return intervals through a systematic inclusion of impact records in trees.

Variations of clinopyroxene/melt element partitioning during assimilation of olivine/peridotite by low-Mg diorite magma

The effects of crystal chemistry and melt composition on the control of clinopyroxene/melt element partitioning (D) during the assimilation of olivine/peridotite by felsic magma have been investigated in Mesozoic high-Mg diorites from North China. The assimilation resulted in significant increase of Mg, Cr and Ni and only slight (< 30%) decrease of incompatible elements of the magma, and the compositional variations have been mirrored by the normally and reversely zoned clinopyroxene microphenocrysts formed at the early stage of the magma evolution. The Mg# [100 × Mg / (Mg + Fe)] values of the reversely zoned clinopyroxenes increase from 65 to 75 in the core to 85–90 in the high-Mg midsection, and reduce back to 73–79 at the rim. Trace element profiles across all these clinopyroxene domains have been measured by LA-ICP-MS. The melt trace element composition has been constrained from bulk rock analyses of the fine-grained low- and high-Mg diorites. Clinopyroxene/melt partition coefficients for rare earth elements (REE) and Y in the high-Mg group zonings (Mg# > 73–79, DDy < 1.2) are positively correlated with tetrahedral IVAl and increase by a factor of 3–4 as tetrahedral IVAl increases from 0.01 to 0.1 per formula unit (pfu). These systematic variations are interpreted to be controlled by the clinopyroxene composition. In contrast, partition coefficients for low-Mg group zonings (Mg# < 75–79, DDy > 1.2) are elevated by up to an order of magnitude (for REE and Y) or more (for Zr and Hf) at similar IVAl, indicating dominant control of melt composition/structure. DZr and DHf show a larger sensitivity to the compositional change of crystal and melt than DREE. DTi values for the low- and high-Mg zonings show a uniform dependence on IVAl. DSr and DLi are insensitive to the compositional change of clinopyroxene and melt, resulting in Sr depletions in the clinopyroxene zonings with elevated REE without crystallization of plagioclase. Our observations show that crystal chemistry and melt composition/structure may alternatively control clinopyroxene/melt partitioning during the assimilation of peridotite by felsic magma, and may be useful for deciphering clinopyroxene compositions and related crust–mantle processes.

Contrasting PT paths within the Barchi-Kol UHP terrain (Kokchetav Complex): Implications for subduction and exhumation of continental crust

The Barchi-Kol terrain is a classic locality of ultrahigh-pressure (UHP) metamorphism within the Kokchetav metamorphic belt. We provide a detailed and systematic characterization of four metasedimentary samples using dominant mineral assemblages, mineral inclusions in zircon and monazite, garnet zonation with respect to major and trace elements, and Zr-in-rutile and Ti-in-zircon temperatures. A typical diamond-bearing gneiss records peak conditions of 49 ± 4 kbar and 950–1000 °C. Near isothermal decompression of this rock resulted in the breakdown of phengite associated with a pervasive recrystallization of the rock. The same terrain also contains mica schists that experienced peak conditions close to those of the diamond-bearing rocks, but they were exhumed along a cooler path where phengite remained stable. In these rocks, major and trace element zoning in garnet has been completely equilibrated. A layered gneiss was metamorphosed at UHP conditions in the coesite field, but did not reach diamond-facies conditions (peak conditions: 30 kbar and 800–900 °C). In this sample, garnet records retrograde zonation in major elements and also retains prograde zoning in trace elements. A garnet-kyanite-micaschist that reached significantly lower pressures (24 ± 2 kbar, 710 ± 20 °C) contains garnet with major and trace element zoning. The diverse garnet zoning in samples that experienced different metamorphic conditions allows to establish that diffusional equilibration of rare earth element in garnet likely occurs at ~900–950 °C. Different metamorphic conditions in the four investigated samples are also documented in zircon trace element zonation and mineral inclusions in zircon and monazite. U-Pb geochronology of metamorphic zircon and monazite domains demonstrates that prograde (528–521 Ma), peak (528–522 Ma), and peak to retrograde metamorphism (503–532 Ma) occurred over a relatively short time interval that is indistinguishable from metamorphism of other UHP rocks within the Kokchetav metamorphic belt. Therefore, the assembly of rocks with contrasting P-T trajectories must have occurred in a single subduction-exhumation cycle, providing a snapshot of the thermal structure of a subducted continental margin prior to collision. The rocks were initially buried along a low geothermal gradient. At 20–25 kbar they underwent near isobaric heating of 200 °C, which was followed by continued burial along a low geothermal gradient. Such a step-wise geotherm is in good agreement with predictions from subduction zone thermal models.