New LA-ICP-MS cryocell and calibration technique for sub-millimeter analysis of ice cores

Ice cores provide a robust reconstruction of past climate. However, development of timescales by annual-layer counting, essential to detailed climate reconstruction and interpretation, on ice cores collected at low-accumulation sites or in regions of compressed ice, is problematic due to closely spaced layers. Ice-core analysis by laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) provides sub-millimeter-scale sampling resolution (on the order of 100μm in this study) and the low detection limits (ng L–1) necessary to measure the chemical constituents preserved in ice cores. We present a newly developed cryocell that can hold a 1m long section of ice core, and an alternative strategy for calibration. Using ice-core samples from central Greenland, we demonstrate the repeatability of multiple ablation passes, highlight the improved sampling resolution, verify the calibration technique and identify annual layers in the chemical profile in a deep section of an ice core where annual layers have not previously been identified using chemistry. In addition, using sections of cores from the Swiss/Italian Alps we illustrate the relationship between Ca, Na and Fe and particle concentration and conductivity, and validate the LA-ICP-MS Ca profile through a direct comparison with continuous flow analysis results.

Terrestrial Radar Interferometric Measurement of Hillslope Deformation and Atmospheric Disturbances in the Illgraben Debris-Flow Catchment, Switzerland

We describe a method for rapid identification and precise quantification of slope deformation using a portable radar interferometer. A rockslide with creep-like behavior was identified in the rugged and inaccessible headwaters of the Illgraben debris-flow catchment, located in the Central Swiss Alps. The estimated volume of the moving rock mass was approximately 0.5 x 10(6) m(3) with a maximum daily (3-D) displacement rate of 3 mm. Fast scene acquisition in the order of 6 s/scene led to uniquely precise mapping of spatial and temporal variability of atmospheric phase delay. Observations led to a simple qualitative model for prediction of atmospheric disturbances using a simple model for solar radiation, which can be used for advanced campaign planning for short observation periods (hours to days).

Snow shielding factors for cosmogenic nuclide dating inferred from long-term neutron detector monitoring

The depth-dependent attenuation of the secondary cosmic-ray particle flux due to snow cover and its effects on production rates of cosmogenic nuclides constitutes a potential source of uncertainty for studies conducted in regions characterized by frequent seasonal snow burial. Recent experimental and numerical modelling studies have yielded new constraints on the effect of hydrogen-rich media on the production rates of cosmogenic nuclides by low- and high-energy neutrons (<10(-3) MeV and >10(2) MeV, respectively). Here we present long-term neutron-detector monitoring data from a natural setting that we use to quantify the effect of snow cover on the attenuation of fast neutrons (0.1-10 MeV), which are responsible for the production of Ne-21 from Mg and Cl-36 from K. We use data measured between July 2001 and May 2008 at seven stations located throughout the Ecrins-Pelvoux massif (French Western Alps) and its surroundings, at elevations ranging from 200 to 2500 m a.s.l. From the cosmic-ray fluxes recorded during summer, when snow is absent, we infer an apparent attenuation length of 148 g cm(-2) in the atmosphere at a latitude of similar to 45 degrees N and for altitudes ranging from similar to 200 to 2500 m a.s.l. Using snow water-equivalent (SWE) values obtained through snow-coring campaigns that overlap in time the neutron monitoring for five stations, we show that fast neutrons are much more strongly attenuated in snow than predicted by a conventional mass-shielding formulation and the attenuation length estimated in the atmosphere. We suggest that such strong attenuation results from boundary effects at the atmosphere/snow interface induced by the high efficiency of water as a neutron moderator. Finally, we propose an empirical model that allows calculating snow-shielding correction factors as a function of SWE for studies using Ne-21 and Cl-36 analyses in Mg- and K-rich minerals, respectively. This empirical model is of interest for studies with a focus on cosmic-ray exposure dating, particularly if the target rocks are made up of mafic to ultramafic units where seasonal snow-cover is a common phenomenon.

A chrysophyte-based quantitative reconstruction of winter severity from varved lake sediments in NE Poland during the past millennium and its relationship to natural climate variability

Chrysophyte cysts are recognized as powerful proxies of cold-season temperatures. In this paper we use the relationship between chrysophyte assemblages and the number of days below 4 °C (DB4 °C) in the epilimnion of a lake in northern Poland to develop a transfer function and to reconstruct winter severity in Poland for the last millennium. DB4 °C is a climate variable related to the length of the winter. Multivariate ordination techniques were used to study the distribution of chrysophytes from sediment traps of 37 low-land lakes distributed along a variety of environmental and climatic gradients in northern Poland. Of all the environmental variables measured, stepwise variable selection and individual Redundancy analyses (RDA) identified DB4 °C as the most important variable for chrysophytes, explaining a portion of variance independent of variables related to water chemistry (conductivity, chlorides, K, sulfates), which were also important. A quantitative transfer function was created to estimate DB4 °C from sedimentary assemblages using partial least square regression (PLS). The two-component model (PLS-2) had a coefficient of determination of View the MathML sourceRcross2 = 0.58, with root mean squared error of prediction (RMSEP, based on leave-one-out) of 3.41 days. The resulting transfer function was applied to an annually-varved sediment core from Lake Żabińskie, providing a new sub-decadal quantitative reconstruction of DB4 °C with high chronological accuracy for the period AD 1000–2010. During Medieval Times (AD 1180–1440) winters were generally shorter (warmer) except for a decade with very long and severe winters around AD 1260–1270 (following the AD 1258 volcanic eruption). The 16th and 17th centuries and the beginning of the 19th century experienced very long severe winters. Comparison with other European cold-season reconstructions and atmospheric indices for this region indicates that large parts of the winter variability (reconstructed DB4 °C) is due to the interplay between the oscillations of the zonal flow controlled by the North Atlantic Oscillation (NAO) and the influence of continental anticyclonic systems (Siberian High, East Atlantic/Western Russia pattern). Differences with other European records are attributed to geographic climatological differences between Poland and Western Europe (Low Countries, Alps). Striking correspondence between the combined volcanic and solar forcing and the DB4 °C reconstruction prior to the 20th century suggests that winter climate in Poland responds mostly to natural forced variability (volcanic and solar) and the influence of unforced variability is low.

Pollen dispersal properties of Poaceae and Cyperaceae: First estimates of their absolute pollen productivities

Pollen-trap results from the Swiss Alps 1996–2009 were used to assess the pollen dispersal–deposition properties of Poaceae (grasses) and Cyperaceae (sedges). Dispersal parameter values were investigated for a modified version of the Prentice–Sugita pollen dispersal–deposition model. Appropriate values (i.e. realistic in the field and allowing realistic modelling results) for wind speed are suggested to be in the range of 3–7 m s− 1 and for pollen an injection height of 0.03–0.1 m above the ground. The appropriate range of pollen injection height values for grasses and sedges differs from that of trees in the same area, suggesting different pollen dispersal properties between herbs and trees. In addition, logarithmic weighting of the vegetation was tested as an alternative to the modified Prentice–Sugita model. This yielded very similar results, suggesting that the use of such much simpler approximations of the pollen–vegetation relationship is a plausible alternative. Based on the modified Prentice–Sugita model, absolute pollen productivity for Poaceae was estimated to 7300 ± 400 grains cm− 2 year− 1 (1 SE). The data basis for Cyperaceae is smaller than for Poaceae, but the dispersal parameter values determined as appropriate for Poaceae yield good results. Absolute pollen productivity for Cyperaceae was estimated to 6300 ± 1100 grains cm− 2 year− 1 (1 SE).

Linking megathrust earthquakes to brittle deformation in a fossil accretionary complex

Seismological data from recent subduction earthquakes suggest that megathrust earthquakes induce transient stress changes in the upper plate that shift accretionary wedges into an unstable state. These stress changes have, however, never been linked to geological structures preserved in fossil accretionary complexes. The importance of coseismically induced wedge failure has therefore remained largely elusive. Here we show that brittle faulting and vein formation in the palaeo-accretionary complex of the European Alps record stress changes generated by subduction-related earthquakes. Early veins formed at shallow levels by bedding-parallel shear during coseismic compression of the outer wedge. In contrast, subsequent vein formation occurred by normal faulting and extensional fracturing at deeper levels in response to coseismic extension of the inner wedge. Our study demonstrates how mineral veins can be used to reveal the dynamics of outer and inner wedges, which respond in opposite ways to megathrust earthquakes by compressional and extensional faulting, respectively.

Relating orogen width to shortening, erosion, and exhumation during Alpine collision

We investigate along-strike width changes of the thickened, accreted lower plate (TALP) in the Central and in the Eastern Alps. We set the width of the TALP in relation to the inferred amount of collisional shortening and exhumation along six orogen-scale cross sections. Taking the present-day, along-strike gradients in the amount of collisional shortening to represent the temporal evolution of the collisional wedge, it may be concluded that the cross-sectional area of the TALP diminishes during ongoing shortening, indicating that the erosional flux outpaced the accretionary flux. Higher amounts of collisional shortening systematically coincide with smaller widths of the TALP and dramatic increases of the reconstructed eroded rock column. Higher amounts of shortening also coincide with larger amplitudes of orogen-scale, upright folds, with higher exhumation and with higher exhumation rates. Hence, erosion did play a major role in reducing by >30 km the vertical crustal thickness in order to accommodate and allow shortening by folding. Long-term climate differences cannot explain alternating changes of width by a factor of almost 2 along straight segments of the orogen on length scales less than 200 km, as observed from the western Central Alps to the easternmost Eastern Alps. Sedimentary or paleontological evidences supporting such paleo-climatic differences are lacking, suggesting that erosional processes did not directly control the width of the orogen.

Interpretation of fluid inclusions in quartz deformed by weak ductile shearing: reconstruction of differential stress magnitudes and pre-deformation fluid properties

A well developed theoretical framework is available in which paleofluid properties, such as chemical composition and density, can be reconstructed from fluid inclusions in minerals that have undergone no ductile deformation. The present study extends this framework to encompass fluid inclusions hosted by quartz that has undergone weak ductile deformation following fluid entrapment. Recent experiments have shown that such deformation causes inclusions to become dismembered into clusters of irregularly shaped relict inclusions surrounded by planar arrays of tiny, new-formed (neonate) inclusions. Comparison of the experimental samples with a naturally sheared quartz vein from Grimsel Pass, Aar Massif, Central Alps, Switzerland, reveals striking similarities. This strong concordance justifies applying the experimentally derived rules of fluid inclusion behaviour to nature. Thus, planar arrays of dismembered inclusions defining cleavage planes in quartz may be taken as diagnostic of small amounts of intracrystalline strain. Deformed inclusions preserve their pre-deformation concentration ratios of gases to electrolytes, but their H2O contents typically have changed. Morphologically intact inclusions, in contrast, preserve the pre-deformation composition and density of their originally trapped fluid. The orientation of the maximum principal compressive stress (σ1σ1) at the time of shear deformation can be derived from the pole to the cleavage plane within which the dismembered inclusions are aligned. Finally, the density of neonate inclusions is commensurate with the pressure value of σ1σ1 at the temperature and time of deformation. This last rule offers a means to estimate magnitudes of shear stresses from fluid inclusion studies. Application of this new paleopiezometer approach to the Grimsel vein yields a differential stress (σ1–σ3σ1–σ3) of ∼300 MPa∼300 MPa at View the MathML source390±30°C during late Miocene NNW–SSE orogenic shortening and regional uplift of the Aar Massif. This differential stress resulted in strain-hardening of the quartz at very low total strain (<5%<5%) while nearby shear zones were accommodating significant displacements. Further implementation of these experimentally derived rules should provide new insight into processes of fluid–rock interaction in the ductile regime within the Earth's crust.

A thermodynamic model for di-trioctahedral chlorite from experimental and natural data in the system MgO-FeO-Al2O3-SiO2-H2O. Applications to P-T sections and geothermometry

We present a new thermodynamic activity-composition model for di-trioctahedral chlorite in the system FeO–MgO–Al2O3–SiO2–H2O that is based on the Holland–Powell internally consistent thermodynamic data set. The model is formulated in terms of four linearly independent end-members, which are amesite, clinochlore, daphnite and sudoite. These account for the most important crystal-chemical substitutions in chlorite, the Fe–Mg, Tschermak and di-trioctahedral substitution. The ideal part of end-member activities is modeled with a mixing-on-site formalism, and non-ideality is described by a macroscopic symmetric (regular) formalism. The symmetric interaction parameters were calibrated using a set of 271 published chlorite analyses for which robust independent temperature estimates are available. In addition, adjustment of the standard state thermodynamic properties of sudoite was required to accurately reproduce experimental brackets involving sudoite. This new model was tested by calculating representative P–T sections for metasediments at low temperatures (<400 °C), in particular sudoite and chlorite bearing metapelites from Crete. Comparison between the calculated mineral assemblages and field data shows that the new model is able to predict the coexistence of chlorite and sudoite at low metamorphic temperatures. The predicted lower limit of the chloritoid stability field is also in better agreement with petrological observations. For practical applications to metamorphic and hydrothermal environments, two new semi-empirical chlorite geothermometers named Chl(1) and Chl(2) were calibrated based on the chlorite + quartz + water equilibrium (2 clinochlore + 3 sudoite = 4 amesite + 4 H2O + 7 quartz). The Chl(1) thermometer requires knowledge of the (Fe3+/ΣFe) ratio in chlorite and predicts correct temperatures for a range of redox conditions. The Chl(2) geothermometer which assumes that all iron in chlorite is ferrous has been applied to partially recrystallized detrital chlorite from the Zone houillère in the French Western Alps.

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.

The first major incision of the Swiss Deckenschotter landscape

The Swiss Deckenschotter (“cover gravels”) is the oldest Quaternary units in the northern Swiss Alpine Foreland. They are a succession of glaciofluvial gravel layers intercalated with glacial and/or overbank deposits. This lithostratigraphic sequence is called Deckenschotter because it “covers” Molasse or Mesozoic bedrock and forms mesa-type hill-tops. Deckenschotter occurs both within and beyond the extent of the Last Glacial Maximum glaciers. The Swiss Deckenschotter consist of two sub-units: Höhere (Higher) and Tiefere (Lower) Deckenschotter. Although the Höhere Deckenschotter sub-unit (HDS) is topographically higher than the Tiefere Deckenschotter, it is older. The only available age for the Swiss Deckenschotter is 2.5–1.8 Ma based on mammal remains found in HDS at the Irchel site. In this study, we present an exposure age for the topographically lowest HDS, calculated from a cosmogenic 10Be depth-profile. Our results show that the first phase of the Deckenschotter glaciations in the Swiss Alps terminated at least 1,020+80−120 ka ago, which is indicated by a significant fluvial incision. This line of evidence seems to be close to synchronous with the beginning of the Mid-Pleistocene Revolution, when the frequency of the glacial-interglacial cyclicity changed from 41 to 100 ka and the amplitude from low to high, between marine isotope stages 23 and 22.

The triplite–triploidite supergroup: structural modulation in wagnerite, discreditation of magniotriplite, and the new mineral hydroxylwagnerite

Electron-microprobe analysis, single-crystal X-ray diffraction with an area detector, and high-resolution transmission electron microscopy show that minerals related to wagnerite, triplite and triploidite, which are monoclinic Mg, Fe and Mn phosphates with the formula Me2+ 2PO4(F,OH), constitute a modulated series based on the average triplite structure. Modulation occurs along b and may be commensurate with (2b periodicity) or incommensurate but generally close to integer values (∼3b, ∼5b, ∼7b, ∼9b), i.e. close to polytypic behaviour. As a result, the Mg- and F-dominant minerals magniotriplite and wagnerite can no longer be considered polymorphs of Mg2PO4F, i.e., there is no basis for recognizing them as distinct species. Given that wagnerite has priority (1821 vs. 1951), the name magniotriplite should be discarded in favour of wagnerite. Hydroxylwagnerite, end-member Mg2PO4OH, occurs in pyrope megablasts along with talc, clinochlore, kyanite, rutile and secondary apatite in two samples from lenses of pyrope–kyanite–phengite–quartz-schist within metagranite in the coesite-bearing ultrahigh-pressure metamorphic unit of the Dora-Maira Massif, western Alps, Vallone di Gilba, Val Varaita, Piemonte, Italy. Electron microprobe analyses of holotype hydroxylwagnerite and of the crystal with the lowest F content gave in wt%: P2O5 44.14, 43.99; SiO2 0.28, 0.02; SO3 –, 0.01; TiO2 0.20, 0.16; Al2O3 0.06, 0.03; MgO 48.82, 49.12; FeO 0.33, 0.48; MnO 0.01, 0.02; CaO 0.12, 0.10; Na2O 0.01, –; F 5.58, 4.67; H2O (calc) 2.94, 3.36; –O = F 2.35, 1.97; Sum 100.14, 99.98, corresponding to (Mg1.954Fe0.007Ca0.003Ti0.004Al0.002Na0.001)Σ=1.971(P1.003Si0.008)Σ=1.011O4(OH0.526F0.474)Σ=1 and (Mg1.971Fe0.011Ca0.003Ti0.003Al0.001)Σ=1.989(P1.002Si0.001)Σ=1.003O4(OH0.603F0.397)Σ=1, respectively. Due to the paucity of material, H2O could not be measured, so OH was calculated from the deficit in F assuming stoichiometry, i.e., by assuming F + OH = 1 per formula unit. Holotype hydroxylwagnerite is optically biaxial (+), α 1.584(1), β 1.586(1), γ 1.587(1) (589 nm); 2V Z(meas.) = 43(2)°; orientation Y = b. Single-crystal X-ray diffraction gives monoclinic symmetry, space group P21/c, a = 9.646(3) Å, b = 12.7314(16) Å, c = 11.980(4) Å, β = 108.38(4) , V = 1396.2(8) Å3, Z = 16, i.e., hydroxylwagnerite is the OH-dominant analogue of wagnerite [β-Mg2PO4(OH)] and a high-pressure polymorph of althausite, holtedahlite, and α- and ε-Mg2PO4(OH). We suggest that the group of minerals related to wagnerite, triplite and triploidite constitutes a triplite–triploidite super-group that can be divided into F-dominant phosphates (triplite group), OH-dominant phosphates (triploidite group), O-dominant phosphates (staněkite group) and an OH-dominant arsenate (sarkinite). The distinction among the three groups and a potential fourth group is based only on chemical features, i.e., occupancy of anion or cation sites. The structures of these minerals are all based on the average triplite structure, with a modulation controlled by the ratio of Mg, Fe2+, Fe3+ and Mn2+ ionic radii to (O,OH,F) ionic radii.

Valle d’Aosta section of the Sesia Zone: multi-stage HP metamorphism and assembly of a rifted continental margin

The spectrum characteristic of the EMC ranges from eclogites (containing omphacite and/or jadeite, garnet, phengite, glaucophane, zoisite, chloritoid, rutile) to phengite schists, calcschists, and marbles, as well as a variety of orthogneisses. Despite the intense polyphase deformation and HP-metamorphic recrystallization, it is possible in some locations to recognize pre-Alpine characteristics in some of the protoliths. For instance, two types of felsic orthogneiss can be distinguished in the Aosta Valley, one derived from Permian granitoids (with local preservation of intrusive contacts, magmatic inclusions, leucocratic veins and other magmatic structures; Stop 3), the other derived from pre-Variscan leuco-monzogranite, such as the building stone mined at the “Argentera” quarry near Settimo Vittone / Montestrutto (Stop 2; so-called “Verde Argento” contains jadeite, phengite, K-feldspar, quartz). Polycyclic and more rarely monocyclic metasediments contain evidence of a complex Alpine PTDt-evolution, locally including relics of their prograde history from blueschist, one or more stages at eclogite facies. Recent petrochronological studies have dated this HP-evolution of the Sesia Zone in some detail. In the area visited, clear evidence of HP-cycling has been identified in one km-size tectonic slice (Stop 1), but not in adjacent parts of the EMC, indicating “yo-yo tectonics”. Partial retrogression and attendant ductile to brittle deformation of the HP-rocks is evident in one of the outcrops (Stop 4). Apart from the four localities in the Sesia Zone, a final outcrop introduces HP-rocks of the adjacent Piemonte oceanic unit, specifically calc-schists and ophiolite members of the “Zermatt-Saas” zone. The hilltop outcrop (Stop 5) displays foliated antigorite schist with peridotite relics (clinopyroxene, spinel) containing lenses derived from doleritic dykes. These fine-grained metarodingites and the folded veins containing Mg-chlorite and titanoclinohumite within serpentinite once again indicate equilibration under low-temperature eclogite facies conditions. However, these units reached that HP stage more than 20 Ma after the youngest eclogite facies imprint recognized in the Sesia Zone. Despite nearly half a century of intense study in the Sesia Zone, the complex assembly of its HP-terranes and their relation to more external parts of the Western Alps remains incompletely understood. This field guide merely introduces a few of the classic outcrops and discusses some of the critical evidence they contain, but it could not incorporate details on each stage of the evolution recognized so far.

Modelling the system behaviour of wet snow avalanches using an expert system approach for risk management on high alpine traffic roads

The presented approach describes a model for a rule-based expert system calculating the temporal variability of the release of wet snow avalanches, using the assumption of avalanche triggering without the loading of new snow. The knowledge base of the model is created by using investigations on the system behaviour of wet snow avalanches in the Italian Ortles Alps, and is represented by a fuzzy logic rule-base. Input parameters of the expert system are numerical and linguistic variables, measurable meteorological and topographical factors and observable characteristics of the snow cover. Output of the inference method is the quantified release disposition for wet snow avalanches. Combining topographical parameters and the spatial interpolation of the calculated release disposition a hazard index map is dynamically generated. Furthermore, the spatial and temporal variability of damage potential on roads exposed to wet snow avalanches can be quantified, expressed by the number of persons at risk. The application of the rule base to the available data in the study area generated plausible results. The study demonstrates the potential for the application of expert systems and fuzzy logic in the field of natural hazard monitoring and risk management.

Climatology of Vb cyclones, physical mechanisms and their impact on extreme precipitation over Central Europe

Cyclones, which develop over the western Mediterranean and move northeastward are a major source of extreme weather and known to be responsible for heavy precipitation over the northern side of the Alpine range and Central Europe. As the relevant processes triggering these so-called Vb events and their impact on extreme precipitation are not yet fully understood, this study focuses on gaining insight into the dynamics of past events. For this, a cyclone detection and tracking tool is applied to the ERA-Interim reanalysis (1979–2013) to identify prominent Vb situations. Precipitation in the ERA-Interim and the E-OBS data sets is used to evaluate case-to-case precipitation amounts and to assess consistency between the two data sets. Both data sets exhibit high variability in precipitation amounts among different Vb events. While only 23 % of all Vb events are associated with extreme precipitation, around 15 % of all extreme precipitation days (99 percentile) over the northern Alpine region and Central Europe are induced by Vb events, although Vb cyclones are rare events (2.3 per year). To obtain a better understanding of the variability within Vb events, the analysis of the 10 heaviest and lowest precipitation Vb events reveals noticeable differences in the state of the atmosphere. These differences are most pronounced in the geopotential height and potential vorticity field, indicating a much stronger cyclone for heavy precipitation events. The related differences in wind direction are responsible for the moisture transport around the Alps and the orographical lifting along the northern slopes of the Alps. These effects are the main reasons for a disastrous outcome of Vb events, and consequently are absent in the Vb events associated with low precipitation. Hence, our results point out that heavy precipitation related to Vb events is mainly related to large-scale dynamics rather than to thermodynamic processes.