Triggering processes of deep-seated gravitational slope deformation (DSGSD) in an un-glaciated area of the Cavargna Valley (Central Southern Alps) during the Middle Holocene

Triggering mechanisms and causative processes of deep-seated gravitational slope deformations (DSGSD) in Alpine and high mountain areas include, among others, post-glacial de-buttressing, earthquake-induced ground shaking, or co-seismic surface faulting. Distinguishing between climatic or tectonically driven factors is challenging since faults and fracture systems can play both an active and/or passive role in the process initiation. We applied an integrated approach, including morpho-structural analysis, geologic field survey, a paleoseismological approach applied to trenching, radiocarbon dating and detailed sedimentological analysis, to a DSGSD located in the Cavargna Valley (N Italy), an area that was never occupied by extensive ice tongues in the Upper Pleistocene and Holocene. We were able to document at least two phases of movement from the Middle Holocene onwards and related triggering factors. Thanks to the comparison with dated lacustrine turbidites in the Alpine Lakes, we argue a seismic triggering for the landslide onset, with a possible source located in southern Switzerland, close to the Insubric Line. The later evolution of the landslide was instead driven by climatic predisposing conditions (humid and cold), instead, corresponding to a phase of enhanced slope instability.


Introduction
Triggering mechanisms and/or causative processes of deep-seated gravitational slope deformations (DSGSD) in alpine and high mountain areas include, among others, post-glacial de-buttressing, earthquake-induced ground shaking, or co-seismic surface faulting, but they are generally hard to be identified. Distinguishing between climatic or tectonically driven factors can be particularly challenging since faults and fracture systems can play both an active and/ or passive role in the process initiation, whereas the chronological correlation between DSGSD activation and climatic-driven increase instability requires precise radiometric dating. Nevertheless, this issue surely has significant implications for recent landscape evolution models as much as for paleoseismological studies and risk assessment.
A possible approach to face this issue is the dating of the onset and phases of activity of DSGSDs and a careful comparison with other external evidence of slope instabilities in the same region, looking for spatio-temporal clusters that would point to a triggering mechanism localized both in time and space (e.g., a seismic triggering). Conversely, a temporal cluster of evidence, more widespread over the entire region, possibly indicates a climatic driving factor, instead.
In this line, the surface expression of DSGSDs includes, especially in the uphill sector of the slope movement, morphostructural features such as double-crested ridges, counter slope curvilinear or straight scarps, open fissures, and morphological trenches (Onida 2001;Agliardi et al. 2001;Gutiérrez et al. 2009;Mariani and Zerboni 2020). These morpho-structures, generally characterized by ridge-top depressions bounded by fault scarps with opposite orientations, allow eluvial and/ or colluvial sediments to be deposited, offering the opportunity to constrain the timing of onset and the evolution of gravity-driven slope deformations through a paleoseismological approach (Forcella et al. 2001;Tibaldi et al. 2004;Gutiérrez-Santolalla et al. 2005;Agliardi et al. 2009;Gutiérrez et al. 2009;McCalpin and Corominas 2019).
We applied an integrated approach to a DSGSD in a mountain area of northern Italy that was never occupied by extensive ice tongues in the Upper Pleistocene and Holocene and thus where slope unloading, due to the removal of glacial confinement, cannot be invoked. The case study is in the Cavargna Valley, N Italy ( Fig. 1a), an area located outside the extent of the Alpine ice cap during Last Glacial Maximum (LGM) and only partially interested by small cirque glaciers. Thanks to an integrated approach, including morpho-structural analysis, geologic field survey, a paleoseismological approach applied to trenching, radiocarbon dating, and detailed sedimentological analysis, we were able to identify two phases of slope deformation and discuss between different causative mechanisms of both.

Geological, structural, and geomorphological setting
In the following paragraphs, we will outline the geological and geomorphological framework of the area, particularly focusing on its late Quaternary landscape evolution, to consider all the possible predisposing factors and causative processes of the DSGSD.

Regional setting
The study site is located in the European Southern Alps, close to the Tonale Line, a segment of the Periadriatic Lineament (Schmid et al. 1989) that accommodated through a dextral strike-slip motion the indentation of the Southern Alps against the Central Alps during the Oligo-Miocene.
Bedrock of the studied area belongs to a portion of the upper crustal tectono-metamorphic units of the Southern Alps (i.e., the Dervio -Olgiasca Zone - Fig. 1b; Spalla et al. 2002), located between two regional tectonic boundaries: the Musso Fault, to the N, and the Lugano -Valgrande Fault, to the south (Fig. 1). The Musso Fault is a cataclastic Alpine structure that was probably active at depth, during pre-Alpine tectonics (Siletto et al. 1990;Bertotti et al. 1993). The Lugano -Valgrande Fault is a Lower Jurassic normal fault that was later upturned during Alpine tectonics (Bertotti, 1991) and that is marked by the presence of a thick greenschist-facies mylonitic belt, where the deepest fault zone crops out (Fig. 1b).
Brittle structures crosscutting through the bedrock, as emerging from a lineament analysis (Fig. 1c), show a predominant N110 direction, closely related to the local strike of the Lugano-Valgrande and Musso faults, with secondary modes at N30 and N150, possibly two Riedel systems of the former, formed during the dextral motion of the Tonale Line.

Study site
The study site is located at the head of Piazza valley (Fig. 2), a small ca. NE-SW trending tributary of the main Cavargna Valley. The Piazza valley is almost entirely interested by a DSGSD feature, ~ 1250 m wide and 1300 m long, which deformed the valley slope.
The valley presents a moderate relief energy, with an elevation drop of ~ 650 m (from 1750 to 1100 m asl) and is characterized by a steep down valley sector (~ 35-40°) and a gentler uphill slope (~ < 20°). An overall convex profile is marked by a clear break of slope between 1300 and 1400 m asl.
A considerable slope adjustment must be inferred for this region during the late Pleistocene, deeply interested by extensive glaciations and by periglacial and paraglacial processes. As common in the Alps, paraglacial environment, in particular, was influenced by post-glacial de-buttressing, fall of local base-level and gravitational adjustment of poorly vegetated slopes, following the retreat of glacial bodies (Knight and Harrison 2009). A detailed reconstruction of the recent geological and geomorphological evolution of the area is therefore a fundamental step for the understanding of the evolution and the causative processes of this DSGSD.
The Cavargna Valley DSGSD area and its surrounding were never interested by extensive glaciations in the Upper Pleistocene. According to Schlüchter et al. (2009), the main LGM ice cap reached an elevation of ~ 950 m asl; additionally, small cirque glaciers locally occupied the most elevated sectors of the Alps, typically above 1700 m asl, as highlighted by glacial and periglacial erosional landforms.
The area experienced a considerable regional deformation due to the isostatic answer to glacial retreat (post glacial rebound). Present-day uplift of the Alps is measured in ~ 1 mm/year through geodesy and remote sensing (e.g., Kahle et al. 1997;Gudmundsson 1994) and predicted that post-LGM rebound presumably contributes significantly to this value.
The temporal evolution of glacial withdrawal, bracketed through exposure dating of boulders and bedrock surfaces (Ivy-Ochs et al. 2006) highlights a very fast deglaciation phase along the southern Alps, starting from ~ 21 ka BP and leading to the almost complete dismantling of glaciers before 15 ka BP, thus implying a quick passage from full glacial conditions to a disequilibrium stage.
Finite-element modeling for the post-LGM isostatic rebound (Norton and Hampel 2010) highlights, for the site area, a displacement of ~ 90 m, cumulated between 21 and 13 ka BP and only ~ 10 m from 13 ka BP to present. The temporal evolution is thus characterized by a quick uplift in the first few thousands of years after glacial removal, followed by a long-lasting period of slow movement.
All the processes closely related to post-glacial rebound that could have caused and/or triggered slope mass movements (e.g., glacial seismotectonics and changes in local base-level) had therefore to be active in a time-window close to glacial withdrawal and, consequently, the onset and timing of DSGSD features, as well, have to be consistent and in time-dependency with these causative processes.
From this point of view, this case study is representative of the landscape morphology of high-altitude sectors in the whole Central Southern Alps, and can provide a significant contribution in analyzing slope instabilities non-related to paraglacial dynamics.
The area is covered by an extensive mapping of the active, inactive, and quiescent landslides ( Fig. 1c; IFFI database -"Inventario dei Fenomeni Franosi in Italia" http:// www. proge ttoif fi. ispra mbien te. it/ carta netif fi/ proge tto. asp? lang= EN; Trigila and Iadanza 2008;Trigila et al. 2010) that highlight the high density of instabilities in this high mountain catchment. More than 750 features were mapped in the IFFI database in Cavargna Valley: ~ 700 features are classified as rotational-translational or complex, 20 are rockfall-toppling and only 10 features are of flow type. This class distribution of the mapped slope instabilities characterizes valley flanks that were subject to recent Fig. 1 a Regional geological sketch map of the Alps: the white circles indicate the locations of dated landslides nearby the study area (A, Tibaldi et al. (2004); B, Forcella et al. (2001); C, Poschinger et al. (2006); D, Agliardi et al. (2009); b Regional geological map with emphasis on the tectonic-metamorphic zones of the South Alpine basement in the Lake Como area: DCZ, Domaso-Cortafò Zone; DOZ, Dervio -Olgiasca Zone; MMZ, Mt. Muggio Zone (after, Spalla et al. 2002); c Shaded Digital Elevation Model (DTM 5 m; Regione Lombardia) with digitized lineaments (rose diagram on the right) and mapped landslides from the Italian landslide database (source: IFFI database, http:// www. proge ttoif fi. ispra mbien te. it/ carta netif fi/ proge tto. asp? lang= EN); the yellow star indicates the study site in all the panels. ◂ re-adjustment, with a limited thickness of shallow loose deposits and typically composed by basement lithologies.
Detailed geological field survey, conducted at 1:10.000 scale, and structural data collection allowed to recognize that the local structural setting played a role in predisposing the slope to deep seated deformations (Fig. 3).
The bedrock is widely outcropping along the slopes of the area, whereas landslide blocky deposits, later incised by the local drainage, are outcropping in the lower sectors of the valley. Rock foliation, here, considered also as a primary discontinuity driving possible slope movements, is moderately to highly dipping with an almost constant ca. N-S strike with most of the foliation planes dipping downslope, to the west. Meso-scale folding records a first deformation event (D1 in Fig. 3), ca. E-W directed, later involved into a NW-SE compressional event (D2 in Fig. 3), resulting in moderately expressed superposed folding structures. The area interested by the DSGSD is chaotically deformed, with pervasive jointing disrupting the rock mass.

Materials and methods
We performed a morpho-structural analysis by means of the photointerpretation of aerial stereo-photo pairs (acquisition year 1980, nominal scale 1:20.000) interpreted by means of an analogic SOK-KIA stereo-viewer equipped with a 10 × magnifier.
Interpreted features have been later digitized in a GIS environment and compared with a 5 m spaced digital terrain model (Regione Lombardia, https:// www. carto grafia. servi zirl. it/ arcgis/ servi ces/ wms/ DTM5_ RL_ wms/ MapSe rver/ WMSSe rver). Geological and geomorphological field surveys were performed at 1:10.000 scale in a single survey season.
We excavated an exploratory trench in the niche area of the DSGSD to temporally constrain the onset and evolution of the DSGSD in the release zone by means of a traditional paleoseismological approach (e.g., McCalpin 2009). The trench is ~ 5 m long and 2.5 m deep.
The northern wall of the trench was cleaned and gridded with string, after which we drew a field trench log on graph paper at a scale of 1:20. Stratigraphic units were defined based on color, texture, and sedimentary structures and according to soil-forming processes. Soil horizons were identified according to the guidelines of the Food and Agriculture Organization (FAO 2006) and sampled for physical and chemical analyses, charcoal identification, and dating. The diagnostic soil horizons of buried palaeosols in the sequence were defined according to the international classification systems (Zerboni et al. 2011(Zerboni et al. , 2015WRB 2014;Soil Survey Staff 2014).
Particle size distribution was determined using laser diffraction (Malvern Mastersizer MS-2000) after H 2 O 2 and HCl treatments, according to the procedure described in (Crouvi et al. 2008), for removal of organic matter and precipitated carbonates, respectively.
The age of the polycyclic soil sequence was obtained by dating with radiocarbon (AMS-14 C) three samples of charcoal. AMS-14 C dating results were calibrated (2 s range) using the INTCAL20 curve (Reimer et al. 2020). The results of analyses and dating on soil horizons are reported in Zerboni et al. (2019); the reader is referred to that paper for the discussion of their significance in relation with the formation of the stratigraphic sequence.
Samples for charcoal analysis were wet washed through 4, 2, and 1 mm mesh sieves. The recovered fragments were analyzed using a Nikon Optiphot reflected light microscope at 50-400 × magnification (analyst: Gianalberto Losapio), on transverse, radial and tangential sections. Reference specimens of Laboratorio di Archeobiologia of Musei Civici of Como (Italy) as well as specialized literature (Schweingruber 1978(Schweingruber , 1990 were used.

Photointerpretation and morphostructural analysis
The SE flank of the Piazza Valley is occupied by an almost circular sector showing the typical characteristics of a deep-seated slope deformation (Figs. 4a and 5a). Following a detailed photointerpretation of the area, we re-mapped the extent of the DSGSD, slightly differing from the previously published IFFI one, being larger than previously thought (Fig. 4a).
An upper sector shows a gently dipping and regular profile, resulting from the erosive action of a hanging glacier (Fig. 2) and is crosscut by several ca. N-S trending lineaments and uphill-facing scarps, delimiting trench-like elongated depressions (Figs. 4a and 5b). This upper sector corresponds to the releasing niche area of the DSGSD, indicating an almost E-W directed extension. A major slope break separates this sector from the lower one (Fig. 4a), where secondary and shallower rotational slope deformations are expressed by typical arcuate niche scarps and slightly bulging at their feet. Bodies of terraced blocky landslide deposits, later incised by the drainage network, are present at the toe of the DSGSD (Fig. 5d). We can infer a relatively shallow lower shear plane for the DSGSD (~ 50-80 m below the ground level at maximum), as deduced from the envelope of the deformed sectors (Fig. 4b).
The extent of the deformed area is also confirmed by a structural analysis on the fracture distribution (Fig. 4a) with a drastic increase in the joint frequency observed at the survey stations located inside the interpreted DSGSD. Joint orientation is subparallel to the scarp strike in the upper sector or to the outer boundary of the landslide, whereas stations located outside the deformed area (e.g., VCA8, VCA11, and VCA13) record rare almost north-dipping joint systems.

Trench stratigraphy and geochronology
We excavated an exploratory trench in the upper sector, across an elongated morphological depression (Fig. 5b). We exposed a multiple sequence of slope/colluvial deposits weathered into soil

Original Paper
Landslides horizons that cumulated in a depression bounded, at the sides, by two gravitationally driven normal faults (Fig. 6); subsequent colluvial events show a fining upward trend.
Sediments, deposited through colluvial slope processes, were subsequently weathered, and reorganized into a pedosequence of three stacked soils/paleosols (Fig. 6). For a detailed description of the (paleo) pedological setting of the trench, see a recent companion paper (Zerboni et al. 2019).
The uppermost unit 1 corresponds to the extant soil developed on a silty loam and differentiated into an upper organic A horizon and a lower Bw horizon. Below, unit 2 is a truncated and buried paleosol divided into three horizons: an eluvial 2E horizon occupies the upper position above a rubified 2Bs horizon. Intercalated in unit 2, there is a residual lens of macroscopic charcoal fragments, several centimeters thick, identified as the remains of an open-air hearth. Dating from two charcoal samples taken from the charcoalrich lens gave a result of 2730 ± 43 and 2683 ± 42 years uncal BP (2932-2756 years cal BP and 2865-2743 years cal BP, respectively). The two samples are mainly composed of Abies charcoals, with presence of Laburnum cf. alpinum, Picea/Larix, and Fagus (Table 1). All the species are typical of high-altitude mountain settings, except for beech that is presently occupying areas at altitudes lower than 1400 m asl. Silver fir is the most common species, present at elevations up to 1800 m asl at least until the XVI cent. A.D. but presently almost lacking anywhere (Castelletti and Motella De Carlo 2012).
Below, at the bottom of unit 2, a 2BC horizon follows, with common reddish mottle. Unit 3 is an eluvial 3Et horizon on a silty clay, in which reddish mottles, comparable to those of the units above, are still present. Unit 4 is a weathered rubified 3Bs horizon that overlie a 3C (unit 5) horizon, the latter marking the boundary to the bedrock (unit 6) at about 150 cm below the present-day topographic surface. A charcoal fragment from unit 5 was dated to 6850 ± 20 years uncal BP (7733-7618 years cal BP).
The whole sequence is thickening up inside the fault-bounded depression (graben) and is crosscut by two conjugate sets of normal faults. The kinematic indicators are provided by the orientation of the clast's long axis, along a fault zone with characteristic shear fabric and indicating an almost pure dip slip kinematic, with less than 10% of horizontal component (Fig. 6). The bedding of the slope deposits is slightly dipping to the east, involved in a gentle synclinal folding with the beds getting steeper close to the fault (drag folding). Brittle secondary structures cut through units 4 and 3, displace the base of unit 2 as well, and die out upward inside unit 2 (i.e., fault concealment sensu Bonilla and Lienkaemper 1991). The movement of these faults is therefore postdating both the deposition of the parent material of unit 3 and its successive pedogenesis and the deposition of the lowermost unit 2.

Grain size distribution
The thickening of the units inside the graben reflects an increased sedimentation space that formed as the depression deepened and Landslides widened through time. In order to unravel possible phases of increased movements and recognize events of surface rupture and reworking of slope deposits across the scarp, we proceeded with a high-resolution grain size analysis of the units.
We determined the particle grain size distribution of the uppermost four units. These were sampled (Fig. 6) both outside (H sample series) and inside the fault-bounded depression (G sample series).
Results are shown in Fig. 7: curves are presented coupled per sampled unit and after successive treatments for organic matter and CaCO 3 removal.
Overall, we note a good overlapping of the curves for the same units both in the H and G sampling locations, indicating a direct rework of the sediments from the hanging wall with no further weathering and a very limited transportation. Unit 1 is characterized by a predominant organic fraction, even in coarse aggregates (i.e., up to 50 mm). It shows an increase in the coarse inorganic fraction in the sample collected inside the trench (hundreds of mm), a fraction that is totally lacking in the H sample (located downslope) and that is thus indicating the occurrence of continuous colluvial events whose products were trapped in the depression.
Unit 2 shows a basal part (2 Gb) with a grain size curve perfectly overlapping with the H sample, indicating the same parent material, whereas the upper level (2 Ga) displays a marked shift toward the finer fraction, in all curves, possibly indicating a re-colluvium of sediments from the H sector into the depression, as also indicated by the erosive unconformity at the top of unit 2. Units 3 and 4, conversely, are almost perfectly matching in the H and G sectors.

Interpretation and discussion: characteristics and onset of the Cavargna Valley DSGSD
From the pedostratigraphy described above and considering the interpretation by Zerboni et al. (2019), the following evolution of the Cavargna Valley DSGSD can be inferred, basing on the restoration of the trench log (Fig. 8). Considering that the eastern end of the trench is not well exposed, we restored the trench basing on the cutoff, crosscut relationships, and folding against the western end of the depression. At ~ 7.7 ka BP (step 0 in Fig. 8), unit 5 deposited as a coarsegrained slope/colluvial deposit. This happened in the early Middle Holocene (Davis et al. 2003), a period of climatic stability along the Southern Alps (Magny et al. 2009(Magny et al. , 2012Badino et al. 2018) that recorded, in Southern Europe, slightly colder average temperatures We have no evidence of slope movements at that time or before. Conversely, the inception of the slope deformations happens at some time after, during the deposition of units 4 and 3 but before the development of soil 3 (step 1 in Fig. 8). We can consistently assume that the development of Soil 3 can be ascribed to a late Atlantic period, that ends up close to 5.0 ka BP (Mayewski et al. 2004;Zerboni et al. 2019). This offer a major constrain for the onset of slope failure and thus we can suggest that the landslide inception can be dated back between 7.7 and 5.0 ka BP; we record a displacement of 22 cm during this stage.
Then, deformation slowed down with only 24 cm of displacement cumulated up to step 2 (i.e., ~ 2.7 ka cal. BP; Fig. 8) and with a partial infilling of the fault-bounded depression due to the colluvial event in charge of the accumulation of unit 2.
Later deformation is due to a second pulse of slope movements (step 3) with several decimeters of offset displacing the buried soil 3 and the base of unit 2. Then, soil 2 developed, including an upper organic horizon that was later truncated by erosion.
Finally, a last phase of slow and continuous deformation, coupled with natural or human-driven colluvial activity and pedogenesis can be attested and probably dated back since the Roman warm period up to the present-day (Zerboni et al. 2019).
The evolution of the Cavargna Valley DSGSD records an almost continuous sequence of slope movements from the Middle Holocene onwards, far after the disintegration of the LGM ice cap. This fact and the fact that the area has never been interested by any major ice tongue let use to exclude that a major role in predisposing or triggering the mass movement has been played by glacial and post-glacial dynamics, including glacial debuttressing. We thus consider other processes including a climatic and/

Original Paper
Landslides or tectonic control on the landslide triggering; both of these two, possibly concurrent, triggering mechanisms have been invoked for known large landslides in the Alps, that seem to cluster in the same period (see e.g., Borgatti and Soldati 2010;Carbonel et al. 2013;Strasser et al. 2013;Zerathe et al. 2014;Oswald et al. 2021). Seismic triggering is one of the major causes for large deep-seated landslides (e.g., Keefer 1984Keefer , 2000Guzzetti et al. 2009;Delgado et al. 2011;Wang et al. 2019;Chunga et al. 2019). One of the possible supporting evidence for claiming a seismic triggering is the spatio-temporal clustering of large landslides, as observed also in the Alps (Sanchez et al. 2010;Zerathe et al. 2014) where a clustering of large events is ascribed to the subboreal chronozone (i.e., ~ 4.4-3.0 ka BP; Prager et al. 2008;Viganò et al. 2021;Oswald et al. 2021). Seismic shaking can also act as a weakening and predisposing factor for later landslides triggered by relatively moderate events, as also observed in offshore lacustrine records (Oswald et al. 2021). The onset of the Cavargna Valley landslide is ascribed to the Middle Holocene, a period when we record an apparently increased seismic activity in the whole Alpine chain (Schnellmann et al. 2002;Sanchez et al. 2010;Strasser et al. 2013;Kremer et al. 2017).
If we consider a seismic triggering for the pristine activation of this landslide, we must investigate the possibility that this could have been triggered by a seismogenic source included in the current seismotectonic setting of the chain. We thus consider a circular area, centered on the Cavargna Valley site and with radius of 250 km (Fig. 9) to identify all the possible sources triggering the earthquakes. This area corresponds to the maximum extent of the area affected by seismically-triggered landslides (Michetti et al. 2007), given an earthquake of maximum epicentral intensity X MCS (i.e., the upper bound for earthquakes in the Alpine seismotectonic setting; Dal Zilio et al. 2018). We proceeded with a simple buffer analysis, following the attenuation regressions for the maximum distance of seismically triggered landslides, given the earthquake epicentral intensity (Papadopoulos and Plessa 2000;Prestininzi and Romeo 2000). We conservatively considered an integrated upper bound relation, by adopting the maximum value from both the regressions for any given epicentral intensity. For any recorded earthquake from historical or instrumental catalogues (Fäh et al. 2011;Rovida et al. 2022), we buffered a region, around the epicenter, according to the regressions given above: none of the recorded strong earthquakes would have been potentially capable of triggering a landslide at the site (Fig. 9).
It results that either a presently unknown seismogenic source, closer to the site, must be supposed or a greater seismogenic potential must be invoked for some of the considered seismogenic sources. In the former hypothesis, there are several attested capable faults (ITHACA database; https:// www. ispra mbien te. gov. it/ en/ proje cts/ soil-and-terri tory/ italy-hazar ds-from-capab le-fault ing-1? set_ langu

Original Paper
Landslides age= en) or potentially active faults (Hetényi et al. 2018) in the area not covered by the buffered region (Fig. 9). These include some structures located in the core of the Alps, such as the Alpe di Cauritt Fault, the Bormio Fault, and the Cima Soliva Fault (see locations in Fig. 9) and two major capable faults at the southern margin of the range: the Gonfolite Backthrust (Sileo et al. 2007) and the Albese con Cassano Anticline (Michetti et al. 2012).
In any case, external evidence such as a worldwide relationship between magnitude and plate convergence rates (Dal Zilio et al. 2018) and macroseismic data, as well , both indicate an upper magnitude bound of M L 6.5 (i.e., I 0 IX-X MCS) for the present day seismotectonic setting of the Central Alps.
Conversely, a possible hypothesis for claiming an increased seismotectonic potential in the Alps during past periods is related to a fast and transient crustal deformation induced by postglacial rebound (Sanchez et al. 2010). Similar observations are widespread in the Fennoscandia shield area (Mörner 2004;van Loon et al. 2016) and in continental Europe (Brandes et al. 2015), where Late Glacial large earthquakes occurred with impressive frequency. A recent estimate of the post-glacial rebound in the Alps, based on an analytical solution for a half-space lithosphere of the Alps (Norton and Hampel 2010), predict ~ 100 m for the area close to the Cavargna Valley: 90 m regained during an initial fast phase of rebound (between 21 and 13 Ka BP) followed by a long-tailored slow deforming period when the remaining amount has been recovered. We expect that, if an increase of crustal deformation and, in turn, seismicity happened, this had to be synchronous with the crustal rebound or with a peak of frequency delayed only a few thousands of years, as observed elsewhere (Ballantyne et al. 2014a, b).
A climatic hypothesis must be considered as well; in fact, several authors suggest a Middle/Late Holocene climatic-triggered (re) activation of landslides and more in general slope instability in the southern Alps (Borgatti and Soldati 2010;Pelfini et al. 2014;Zerathe et al. 2014) and along the northern Apennine of Italy (Bertolini and Tellini 2001;Bertolini 2007). Even though the structural, seismic, and lithologic parameters oversee activation and dynamic of slope processes at different scales, it has been acknowledged that a higher frequency of rainfall and increased cooling led to (re)activation of landslides (Trauth et al. 2000(Trauth et al. , 2003McGuire and Maslin 2012). In fact, increased rainfall and discharge and high fluctuations in precipitation, runoff and erosion cause significantly increased destabilization of slopes; moreover, the same factors influence the pore-water pressures lowering the critical thresholds in slopes susceptible to failure (Trauth et al. 2000;Guzzetti et al. 2007;Peruccacci et al. 2012;Segoni et al. 2018).
Records of flooding events are usually helpful in the reconstruction of rainfall fluctuations. Floods are related to the increased Fig. 9 Buffer analysis for a possible seismic triggering of the Cavargna Valley landslide, based on the analysis of instrumental and historical seismicity (Fäh et al. 2011;Rovida et al. 2022); capable (ITHACA database; https:// www. ispra mbien te. gov. it/ en/ proje cts/ soil-andterri tory/ italy-hazar ds-fromcapab le-fault ing-1? set_ langu age= en) and potentially active faults (Hetényi et al. 2018) are drawn in the closeness of the study site; see text for details. Codes: AA, Albese con Cassano Anticline; AdC, Alpe di Cauritt Fault; BF, Bormio Fault; CS, Cima Soliva Fault; GB, Gonfolite Backthrust Fig. 10 Comparison of evidence for a climatic or seismic triggering for the Cavargna Valley DSGSD: a Location of offshore dated lacustrine turbidites (Kremer et al. 2020) and dated landslides (labeled as A, B, and C)within 100 km from the site: locations highlighted in red show a spatio-temporal clustering consistent with the Cavargna Valley site and potentially indicating a seismic triggering mechanism; b Stacked flood records for the Southern Alps spanning the past 10 kyr (Wirth et al. 2013) correlated with the timing of the Cavargna Valley landslide movements (blue-shaded areas indicate increased frequency of floods, including the latest period during the Little Ice Age-LIA); c Offshore evidence of paleoseismic events from dated lake turbidites in the Central Alps (Kremer et al. 2020;Rapuc et al. 2022 for Lake Iseo); see text for details Landslides runoff caused by more precipitations. The same process also drives the start of rhexistasy phases and the development of slope instability. We used one of the longest offshore lacustrine records of flooding for the Southern Alps (Wirth et al. 2013); floods are related to increased runoff under rhexistasy phases (Erhart 1967), and the same setting generally promotes the activation of slope instability. According to this geological record, major flooding events are clustered at the 6.4-5.0 ka BP, 4.2-2.4 ka BP (with a peak at 2.6-2.4 ka BP) intervals, and during the Little Ice Age (Fig. 10). Other geomorphological and hydroclimatic records suggest, for the same intervals, glacial advances, cooler climatic conditions, fluctuating lake levels, and shifts in the timberline (see a review in Magny et al. 2009Magny et al. , 2012Badino et al. 2018).
The seismic and climatic triggers are thus both reliable factors controlling the onset and evolution of the Val Cavargna DSGSD. In order to disentangle between the two hypotheses, we checked a possible spatio-temporal correlation between the onset of the pulses of landslide movements and hydroclimatic variability (Wirth et al. 2013) together with a dataset of offshore lacustrine paleoseismic indicators. Along the Alps, lacustrine turbidites have been identified in many cores and attributed to paleoseismicity (i.e., lacustrine turbidites; Kremer et al. 2020). Figure 10 illustrates the locations of offshore dated lacustrine turbidites and landslides inside a radius of 100 km around the study site.
The inception of the Cavargna Valley landslide, is poorly constrained in a quite long time window (i.e., 7.7-5.0 ka BP), but we noticed that it happened contemporary to a period which records two clusters of paleoseismic evidence (i.e., 7.3-7.0 ka BP and 6.5-5.5 ka BP), closely related in time and space (yellow intervals in Fig. 10b). Notably, other lakes, spreading farther to the NE, do not record a contemporary increase in the turbiditic events. The spatial distribution of the indicators, suggest that the likely location for the potential seismogenic source lies in the middle of the Ticino District (southern Switzerland) close to the Insubric Line; a seismic event at such location may have started the activity of the Cavargna Valley DSGSD. Yet, from the hydroclimatic point of view, the same time window lies across the switch from a warm period to one characterized by an increase in precipitations. We can suppose that the new climatic conditions could have contributed to promoting the movement of the Cavargna Valley DSGSD during the second part of this time window, but the spatial clustering of the evidence strongly suggest the triggering from a localized seismogenic source, instead.
The inception of the second pulse of movement (2.8-2.0 Ka BP) is closely temporally related to the peak in precipitations recorded by the lacustrine archives. Additionally, a widespread evidence of an increasing rate of turbiditic episodes is recorded in several lakes across the entire considered area (Fig. 10b), thus lacking a specific spatial clustering around the Cavargna Valley site. For the same period, a noteworthy contribution of human activities (wood clearance and herding) has possibly amplified the effect of the climate (Zerboni et al., 2019).
The role played by the climate conditions as a predisposing factor driving the evolution of such large landslides has to be carefully taken into account. Remarkably, the two phases of movement documented at the study site are coeval with cold and humid periods reconstructed in the Alpine region, in particular cold event CE3 between 7.5-6.9 ka BP and CE8 at 2.6-2.3 ka BP (Haas et al. 1998). The intensity of the first phase of cooling was relatively low and likely did not play a major role in the activation of the Val Cavargna DSGSD. On the contrary, several factors may have contributed to the reactivation of the DSGSD during the second phase, including lower temperatures, greater precipitations and all the environmental consequences of such processes: decrease in the local plant cover (negative shift of the treeline), increased periglacial conditions and formation of permafrost or even formation of small local glaciers. Such events have a clear effect on the reactivation of shallow slope processes (i.e., formation of colluvium), but their influence on deep movements is less evident.
In this perspective, we can thus propose that the onset of the Cavargna Valley landslide was triggered by seismic event(s), whereas its subsequent evolution is markedly influenced by climatic factors, thus highlighting a major shift in the factors overseeing the DSGSD dynamics.
We argue that, for settings similar to our case study, a combined action of enhanced seismicity and hydroclimatic conditions could be considered as both a triggering and predisposing factor: the progressive increase in the rock fatigue process and the reach of a critical state of the slope could in fact finally result in a major reactivation triggered by a moderate event (Gischig et al. 2016;Oswald et al. 2021).
Nonetheless, based on evidence, we suggest that the pristine activation of the landslide (initial Middle Holocene) occurred as a response to seismic activity, whereas subsequent pulses (late Holocene) of the structure seem to be much more correlated to phases of climatic instability. Likely, the structural weakness created by the initial event offered the proper conditions for subsequent climate-triggered reactivation of the landslide. Oscillations in climatic conditions and in the amount of precipitation are recorded throughout the Holocene at both millennial and centennial time scales (Magny et al. 2013).

Conclusions
In this study, we investigated a DSGSD in the Cavargna Valley (N Italy) through a morphostructural, pedostratigraphical, geomorphological, and paleoseismological approach. Dating and retrodeformation of the sedimentary sequence unearthed from an exploratory trench allow to document at least two phases of movement from the Middle Holocene onwards; the first pulse is bracketed at 7.7-5.0 ka BP, while the inception of the second pulse is dated at 2.8-2.0 ka BP.
We attempted to identify the triggering mechanism for the Cavargna Valley DSGSD by considering the historical earthquake catalog, a dataset of offshore lacustrine paleoseismological indicators and regional flood chronology. We exclude post-glacial debuttressing because the study region was never occupied by extensive glacial covers. Enhanced seismic activity related to postglacial rebound does not match the timing of the movements, so it can be excluded as well.
We conclude that a seismic triggering is likely for the onset of the DSGSD during the initial Middle Holocene, based on the spatiotemporal clustering with a set of offshore evidence, with a possible source located in an area lacking known historical seismicity or active faults, pointing to a possible knowledge gap in the seismotectonics of the Alps. Later evolution and successive pulses in the activity of the Cavargna Valley DSGSD (Late Holocene), instead, seem to be more correlated to regional proxies of climatic changes, showing a possible correspondence with periods of increased surface instability.
Understanding the timing and triggering mechanisms for DSGSD in Cavargna Valley and surrounding regions is paramount for reconstructing the local landscape evolution and eventually refine the implications for seismic hazard and for the cumulative effect of seismicity, structural weakness, and climatic instability. thanked the fruitful discussion and for the scientific coordination of the Interreg 2007-2013 "SITINET Project": during the field campaign under the umbrella of this project we firstly started the investigations here described.

Funding
Open access funding provided by Università degli Studi dell'Insubria within the CRUI-CARE Agreement.

Competing Interests
The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.