Flow-R, a model for susceptibility mapping of debris flows and other gravitational hazards at a regional scale

The development of susceptibility maps for debris flows is of primary importance due to population pressure in hazardous zones. However, hazard assessment by process-based modelling at a regional scale is difficult due to the complex nature of the phenomenon, the variability of local controlling factors, and the uncertainty in modelling parameters. A regional assessment must consider a simplified approach that is not highly parameter dependant and that can provide zonation with minimum data requirements. A distributed empirical model has thus been developed for regional susceptibility assessments using essentially a digital elevation model (DEM). The model is called Flow-R for Flow path assessment of gravitational hazards at a Regional scale (available free of charge under http://www.flow-r.org) and has been successfully applied to different case studies in various countries with variable data quality. It provides a substantial basis for a preliminary susceptibility assessment at a regional scale. The model was also found relevant to assess other natural hazards such as rockfall, snow avalanches and floods. The model allows for automatic source area delineation, given user criteria, and for the assessment of the propagation extent based on various spreading algorithms and simple frictional laws. We developed a new spreading algorithm, an improved version of Holmgren's direction algorithm, that is less sensitive to small variations of the DEM and that is avoiding over-channelization, and so produces more realistic extents. The choices of the datasets and the algorithms are open to the user, which makes it compliant for various applications and dataset availability. Amongst the possible datasets, the DEM is the only one that is really needed for both the source area delineation and the propagation assessment; its quality is of major importance for the results accuracy. We consider a 10 m DEM resolution as a good compromise between processing time and quality of results. However, valuable results have still been obtained on the basis of lower quality DEMs with 25 m resolution.

Geochemical variability in distal and proximal glass from the Youngest Toba Tuff eruption

The Youngest Toba Tuff (YTT, erupted ca. 74 ka ago) is a distinctive and widespread tephra marker across south and southeast Asia. The climatic, human and environmental consequences of the YTT eruption are widely debated. Although a considerable body of geochemical data is available for this unit, there has not been a systematic study of the variability of the ash geochemistry. Intrinsic (magmatic) and extrinsic (post-depositional) chemical variations bring fundamental information regarding the petrogenesis of the magma, the distribution of the tephra and the interaction between the ash and the receiving environment. Considering the importance of the geochemistry of the YTT for stratigraphic correlations and eruptive models, it is central to the YTT debate to quantify and interpret such variations. Here we collate all published geochemical data on the YTT glass, including analyses from 68 sites described in the literature and three new samples. Two principal sources of chemical variation are investigated: (i) compositional zonation of the magma reservoir, and (ii) post-depositional alteration. Post-depositional leaching is responsible for up to ca. 11% differences in Na2O/K2O and ca. 1% differences in SiO2/Al2O3 ratios in YTT glass from marine sites. Continental tephra are 2% higher in Na2O/K2O and 3% higher in SiO2/Al2O3 respect to the marine tephra. We interpret such post-depositional glass alteration as related to seawater induced alkali migration in marine environments, or to site-specific water pH. Crystal fractionation and consequential magmatic differentiation, which produced order-of-magnitude variations in trace element concentrations reported in the literature, also produced major element differences in the YTT glass. FeO/Al2O3 ratios vary by about 50 %, which is analytically significant. These variations represent magmatic fractionation involving Fe-bearing phases. We also compared major element concentrations in YTT and Oldest Toba Tuff (OTT) ash samples, to identify potential compositional differences that could constrain the stratigraphic identity of the Morgaon ash (Western India); no differences between the OTT and YTT samples were observed.

Structural, metamorphic and geochronological relations between the Zanskar Shear Zone and the Miyar Shear Zone (NW Indian Himalaya): Evidence for two distinct tectonic structures and implications for the evolution of the High Himalayan Crystalline of Zanskar

The geologic structures and metamorphic zonation of the northwestern Indian Himalaya contrast significantly with those in the central and eastern parts of the range, where the high-grade metamorphic rocks of the High Himalayan Crystalline (HHC) thrust southward over the weakly metamorphosed sediments of the Lesser Himalaya along the Main Central Thrust (MCT). Indeed, the hanging wall of the MCT in the NW Himalaya mainly consists of the greenschist facies metasediments of the Chamba zone, whereas HHC high-grade rocks are exposed more internally in the range as a large-scale dome called the Gianbul dome. This Gianbul dome is bounded by two oppositely directed shear zones, the NE-dipping Zanskar Shear Zone (ZSZ) on the northern flank and the SW-dipping Miyar Shear Zone (MSZ) on the southern limb. Current models for the emplacement of the HHC in NW India as a dome structure differ mainly in terms of the roles played by both the ZSZ and the MSZ during the tectonothermal evolution of the HHC. In both the channel flow model and wedge extrusion model, the ZSZ acts as a backstop normal fault along which the high-grade metamorphic rocks of the HHC of Zanskar are exhumed. In contrast, the recently proposed tectonic wedging model argues that the ZSZ and the MSZ correspond to one single detachment system that operates as a subhorizontal backthrust off of the MCT. Thus, the kinematic evolution of the two shear zones, the ZSZ and the MSZ, and their structural, metamorphic and chronological relations appear to be diagnostic features for discriminating the different models. In this paper, structural, metamorphic and geochronological data demonstrate that the MSZ and the ZSZ experienced two distinct kinematic evolutions. As such, the data presented in this paper rule out the hypothesis that the MSZ and the ZSZ constitute one single detachment system, as postulated by the tectonic wedging model. Structural, metamorphic and geochronological data are used to present an alternative tectonic model for the large-scale doming in the NW Indian Himalaya involving early NE-directed tectonics, weakness in the upper crust, reduced erosion at the orogenic front and rapid exhumation along both the ZSZ and the MSZ.

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.

Tectono-metamorphic evolution of a nappe stack: A case study of the Swiss Alps

Fold-and-thrust belts are prominent structures that occur at the front of compressional orogens. To unravel the tectonic and metamorphic evolution of such complexes, kinematic investigations, quantitative microstructural analysis and geothermometry (calcite–graphite, calcite–dolomite) were performed on carbonate mylonites from thrust faults of the Helvetic nappe stack in Central Switzerland. Paleo-isotherms of peak temperature conditions and cooling stages (fission track) of the nappe pile were reconstructed in a vertical section and linked with the microstructural and kinematic evolution. Mylonitic microstructures suggest that under metamorphic conditions close to peak temperature, strain was highly localized within thrust faults where deformation temperatures spatially continuously increased in both directions, from N to S within each nappe and from top–down in the nappe stack, covering a temperature range of 180–380 °C. Due to the higher metamorphic conditions, thrusting of the lowermost nappe, the Doldenhorn nappe, was accompanied by a much more pronounced nappe internal ductile deformation of carbonaceous rock types than was the case for the overlying Wildhorn- and Gellihorn nappes. Ongoing thrusting brought the Doldenhorn nappe closer to the surface. The associated cooling resulted in a freezing in of the paleo-isotherms of peak metamorphic conditions. Contemporaneous shearing localized in the basal thrust, initially still in the ductile deformation regime and finally as brittle faulting and cataclasis inducing ultimately an inverse metamorphic zonation. With ongoing exhumation and the formation of the Helvetic antiformal nappe stack, a bending of large-scale tectonic structures (thrusts, folds), peak temperature isotherms and cooling isotherms occurred. While this local bending can directly be attributed to active deformation underneath the section investigated up to times of 2–3 ma, a more homogeneous uplift of the entire region is suggested for the very late and still active exhumation stage.