Preliminary insights from hydrological field monitoring for the evaluation of landslide triggering conditions over large areas

Rainfall-induced landslides represent a major threat to human activities, and thus an improved understanding of their triggering mechanisms is needed. The paper reports some preliminary inferences on this topic, based on the data recorded over a 2-year period by a multi-parametric monitoring station located on one of the slopes of the Monterosso catchment (Cinque Terre, north-western Italy). This catchment has experienced multiple, concurrent shallow landslides after intense rainfall events. After defining a soil hydraulic model through data interpretation and numerical simulations, slope stability analyses were performed to elucidate several aspects related to shallow landslide occurrence. Both long-term climate conditions and single rainfall events were simulated via physically based approaches. The findings from these simulations enabled us to assume the pattern of infiltration and quantify the impact of soil hydraulic behavior on landslide triggering conditions. In this regard, various analyses were carried out on the same triggering event both at local scale and in the overall catchment, with a view to highlighting the role of initial soil moisture and soil hysteretic behavior in slope stability.


Introduction
Rainfall-induced landslides are the most common form of landslide worldwide (e.g., Dai et al. 2003;Park et al. 2013;Chen et al. 2014;Sajinkumar and Anbazhagan 2015;Marc et al. 2018) and can cause significant economic losses and human casualties (Petley 2012;Haque et al. 2016). In the past few years, both the frequency and the number of landslide and extreme rainfall events have increased (Kirschbaum et al. 2012;Haque et al. 2019); likely due to global climate change (Froude and Petley 2018;Gariano et al. 2020). Hence, the temporal and spatial prediction of rainfall-induced landslides is an issue of growing relevance, and major efforts have been undertaken to develop and/or improve methods and tools designed for this purpose. In particular, a number of real-time early warning systems (EWSs) have recently been developed (e.g., Froese et al. 2005;Badoux et al. 2009;Anderson et al. 2011;Maskrey 2011;Genevois et al. 2018;Abraham et al. 2022) thanks to the availability of improved measuring and monitoring networks (Alfieri et al. 2012). EWSs are generally based on data collected by systems that monitor landslide triggering factors, such as rainfall and/or pore water pressure (Thiebes et al. 2014). The growing interest in EWSs is justified by the fact that, in many cases, they offer a valid solution for risk mitigation with a low environmental and economic impact (Intrieri et al. 2012). Alarm thresholds are generally defined through empirical (e.g., Aleotti 2004;Guzzetti et al. 2008;Brunetti et al. 2010;Mathew et al. 2014;Huang et al. 2015;Dikshit and Satyam 2019;Gariano et al. 2019;Tufano et al. 2019;Valenzuela et al. 2019) or physically based models (e.g., Baum et al. 2008;Liao et al. 2011;Papa et al. 2011;An et al. 2016;Mirus et al. 2018;Hsu and Liu 2019). However, such empirical thresholds have their limitations, mainly related to the need for site-specific landslide inventories and the difficulty of properly considering the complex hydrological processes occurring in slopes (e.g., effective infiltration, evapotranspiration) prior to potential triggering rainfall events (Bogaard and Greco 2018). In this regard, many authors have used antecedent precipitation measurements to determine when landslides are likely to occur, establishing a threshold based on the amount of antecedent rainfall over a specific period (e.g., Govi et al. 1985;Pasuto and Silvano 1998;Cardinali et al. 2006;Abraham et al. 2021). Nevertheless, a key difficulty lies in the definition of the period during which precipitation accumulates. According to Guzzetti et al. (2007), the relevant literature shows a significant scatter in the considered period, that is mostly ascribable to different geological-morphological and climate features, as well as to heterogeneity and incompleteness of rainfall and landslide data used to determine thresholds. Physically based models can provide more reliable results, despite the difficulty of obtaining detailed information on hydrological, morphological, and soil characteristics, especially over large areas (Merritt et al. 2003;Rosso et al. 2006;Berti et al. 2012). Indeed, by simulating the dynamic response of soil wetting, these models can determine the amount of precipitation that is needed to trigger slope failures, directly accounting for antecedent rainfall (Wilson and Wieczorek 1995;Crosta and Frattini 2003). This is why the integration of monitoring data into simplified numerical models can represent a starting point for designing efficient EWSs, as also suggested by Stähli and Bartelt (2007). With this perspective, distributed monitoring can be crucial for a suitable physically based quantification of rainfall-induced landslide scenarios.
In our study, we carried out a series of analyses focused on the assessment of landslide triggering conditions after intense rainfall events. The concept was to make use of the data recorded by an in situ monitoring station in order to assess the triggering process over a wide area. Starting from field monitored data, we first defined several hydraulic parameters 1 3 of the soil, with particular emphasis on those relating to its unsaturated portion, which can play a crucial role in the onset of instabilities within shallow deposits (e.g., Baum et al. 2010;Lu and Godt 2013;Bordoni et al. 2015;Jeong et al. 2017;Rahardjo et al. 2019;Giannini et al. 2022;Martino et al. 2022). Based on our results, we analyzed the slope stability conditions at a local scale via a simplified physically based approach, considering both long-term climate conditions and a specific representative rainfall event. The outcome of these analyses was then used to carry out similar analyses over a larger area using TRIGRS (Baum et al. 2008), a physically based model specifically designed for slope stability analyses at catchment scale. In this regard, the area surrounding our monitoring station (Cinque Terre, north-western Italy) is renowned for numerous landslides that have occurred over the years due to the unique physiographic and meteo-climatic features of its small coastal catchments (e.g., Cevasco et al. 2013;Brandolini et al. 2018;Raso et al. 2019;Di Napoli et al. 2021). Hence, the objective of our work was to evaluate the importance of characterizing the hydraulic behavior of soil covers, and the impact of extrapolating local results to a larger scale. At the same time, we checked whether the information inferred from hydrological field monitoring could contribute to improving basin-scale slope stability analyses.

The study area
The Monterosso catchment (Fig. 1a), which extends over approximately 5.5 km 2 , lies inside Cinque Terre National Park, along the eastern Ligurian coast (north-western Italy). The area features steep (30-40°) slopes, typically carved by deep braided ephemeral streams directly discharging into the sea. These features, which are typical of recently uplifted areas (Faccini et al. 2013), are found in many coastal basins of the easternmost part of Liguria. The combination of this complex orography with the effect of the sea strongly affects the local climate. The latter is typically Mediterranean, with hot and dry summers, and mild winters, fall being the rainiest season. Nevertheless, mean annual precipitation (1033 mm) tends to increase with altitude when moving from the coast to the watershed (Pedemonte 2005). For instance, at the highest elevations of the northern part of the catchment, mean annual precipitation reaches about 1200 mm (Cevasco et al. 2015). Extreme precipitation may also occur, generally between late summer and mid-fall, owing to self-regenerating storm cells or persistent cyclonic Tyrrhenian circulation (Crosta 1998;Cevasco et al. 2009). These rainstorms typically have a short duration (< 24 h), but rainfall intensities may reach or exceed values of 100 mm/h (Cevasco et al. 2015). From a geological point of view, the local bedrock is chiefly composed of a sandstone-claystone flysch and of a pelitic complex, which is overlain by a thin (0.5-3.5 m) regolith cover. The local shallow soils and debris covers have been largely reworked over the centuries for agricultural purposes, resulting into widespread terraced areas (Fig. 1c). These represent the most iconic feature of Cinque Terre National Park (Giordan et al. 2020;Pepe et al. 2020).
However, over the last century, the terraces have been progressively abandoned and covered by Mediterranean scrub, becoming landslide-prone areas. Additionally, also in currently cultivated terraced areas a limited maintenance of the drainage systems can be observed, making the hillslopes even more susceptible to landslide occurrence (Schilirò et al. 2017). In this respect, mention may be made of at least two major rainfall events. The first occurred on the night between October 2 and 3, 1966 when heavy rainfall caused severe injuries to people and damage to buildings and infrastructure, especially in the village of Monterosso. Based on available information, the central roads 1 3 of Monterosso were flooded by water and debris, which spread to basements, shops, and cafes, and swept away boats lying on the beach (Guzzetti et al. 1994;Guzzetti and Tonelli, 2004). During the event, the neighboring Levanto monitoring station recorded 142 mm of rainfall, with maximum intensities of 85 mm/6 h, 127 mm/12 h, and 152 mm/24 h. The second event took place on October 25, 2011, namely between 7:00 and 17:00 UTC. The Monterosso monitoring station recorded a cumulative rainfall of 382 mm, with maximum intensities of 90 mm/h, 195 mm/3 h, and 350 mm/6 h (Regione Liguria 2012). As a consequence of this extreme event, approximately 260 shallow landslides occurred, initiating as debris slides (Fig. 2), and then, in most cases, evolving into debris avalanches or debris flows (Hungr et al. 2014). Huge amounts of earth and debris engulfed the central roads of Monterosso, reaching their maximum thickness (approximately 3 m) in the center of the village (Schilirò et al. 2018).

The Monterosso field monitoring system
After the extreme 2011 weather event, a multi-parametric monitoring system was installed in the Monterosso catchment area, as part of an agreement between Dipartimento per gli Affari Regionali (DARA) of the Italian Government and Dipartimento di Scienze della Terra of Sapienza University of Rome, concerning research on and design of methods to quantify susceptibility to landslides in mountain areas and ability to predict them at an early stage.
The field monitoring station (Fig. 1b) is located on a south-facing vineyard terrace in the central part of the Monterosso catchment, west of the SP38 provincial road, at an elevation of 185-190 m above sea level (asl). The average slope of the terrace is 27°, but it reaches 45° on its flanks (Fig. 2b). This geometry is typical of a particular type of agricultural terrace (i.e., cuiga in local dialect), whose flanks are made up of a mixture of compacted soil and rocks, in place of the more common dry-stone wall. This mixture enhances the formation of a grassy edge which, in turn, reduces runoff and erosion of the terrace flank (Martini et al. 2004). As regards soil cover characteristics, laboratory test results showed a coarse-grained (mostly gravelly) material, with minor components of low-plasticity silt and clay ( Table 1). The unsaturated unit weight of soil was, instead, determined directly in the field via the sand cone method. These features are consistent with an eluvial-colluvial deposit resulting from the weathering of the underlying flysch bedrock and subsequently reworked by farming activities.
The monitoring network consists of ( Fig. 3):  Both the tensiometers and the soil water content sensors were installed at different locations and depths in order to assess spatial changes of hydraulic conditions over time (Fig. 3). In particular, two different depths (50 and 100 cm) were selected to investigate the vertical infiltration process within a depth range consistent with the average soil thickness estimated for terraced areas in the Monterosso catchment (Fig. 4b). In terms of planimetric position, we placed most of the sensors in the outermost part of the terrace, while three sensors (W1, T3, and T4) were installed in the part farthest from the terrace edge, with the aim of observing the hydraulic response of the soil to rainfall events along the slope dip direction.
Data was collected (sampling rate: 5 min) by a CR1000 datalogger (Campbell Scientific Inc.), powered by a photovoltaic panel equipped with a lithium battery. The power unit also feeds a data transmission system based on the GSM (Global System for Mobile Communications)/GPRS (General Packet Radio Service) technology. Specifically, data is automatically transmitted at set intervals to an external server, which can be cloudaccessed from the web. In this way, the application makes it possible not only to display and download data remotely, but also to check the operating status of the monitoring system, detecting potential failures in almost real time.

Physically based modeling
In this work, use was made of different physically based models to reconstruct the hydraulic behavior of the agricultural terrace and its stability conditions.
For hydraulic behavior, we employed HYDRUS-1D version 4.0 (Šimůnek et al. 2008), a software package developed by USDA (United States Department of Agriculture) for simulating one-dimensional, uniform water movement in variably saturated porous media. This model, which is widely used in the literature to analyze water, solutes, and heat propagation in soil (e.g., Kanzari et al. 2018;Shekhar et al. 2019;Davis and Ekwue 2022), uses a modified form of the Richards equation for simulating water flow: where θ is the volumetric water content, t is time, h is the pressure head, S is the sink term, δ is the angle between the flow direction and the vertical axis, and K is the unsaturated hydraulic conductivity given by: where K r and K s are the relative and saturated hydraulic conductivity, respectively. For the analytical hydraulic model to simulate water flow, we chose the van Genuchten-Mualem model (van Genuchten 1980), which may be summarized as follows: where θ is the volumetric water content, θ s is the saturated water content, θ r is the residual water content, α is a fitting parameter depending on pore-air pressure, n is a fitting parameter depending on pore size distribution, u a is the pore-air pressure, and u w is the pore water pressure.
As concerns slope stability, we relied on two different models depending on the scale of our analysis, namely a local-scale model and a catchment-scale model. In the former instance, we resorted to the simplified approach proposed by Lu and Godt (2008). The latter is an extended version of the conventional infinite slope model, which enables calculation of slope stability, in terms of factor of safety (FS), under partially saturated conditions according to: where φ′ is the soil effective friction angle, β is the slope angle, c′ is the soil effective cohesion, γ is the unit weight of soil, H is the potentially unstable soil thickness, and σ S is the suction stress, which may be defined in terms of matric suction (u au w ) for variably saturated materials as (Lu and Likos 2006;Lu et al. 2010): Lu et al. (2010) also proposed a formula for calculating suction stress considering the equivalent degree of saturation (S e ): Thus, taking into account the two formulas for calculating σ S , the factor of safety can be predicted on the basis of both soil water content and pore water pressure values, considering parametrically two different slope angles, 27° and 38°, representing the slope angle of the monitored terrace and the mean slope angle of the 2011 failed slope (Fig. 4c), respectively.
For the catchment-scale analysis, we used TRIGRS (Baum et al. 2008), a model developed by the United States Geological Survey (USGS) that is largely used to predict rainfall-induced landslides over large areas (Weidner et al. 2018;Park et al. 2019;Ávila et al. 2021;Schilirò et al. 2021;Tran et al. 2022). TRIGRS is a model that combines two calculations: (i) a one-dimensional analytical solution for pore pressure response to rainfall infiltration, and (ii) an infinite slope stability calculation. The infiltration model, based on a linearized solution of the Richards equation proposed by Iverson (2000), can account for different types of rainfall inputs on the ground surface. In response to these inputs, TRIGRS can compute the transient pressure head response ψ(Z,t) based on various input parameters, such as slope, soil thickness, depth of the initial steady-state water table, and saturated hydraulic conductivity. At each time step, TRIGRS simulates the rise of the water table if the amount of infiltrating water reaching the surface exceeds the amount that gravity can drain. This analysis is performed by comparing the excess water volume with the pore space directly above the water table, and by applying the resulting water weight at the initial top of the saturated zone at each time step (Baum et al. 2010). Then, the factor of safety is determined at a specific depth Z using the following formula: The model can also account for infiltration into a partially unsaturated surface layer above the water table. In this instance, TRIGRS linearizes the Richards equation by using a number of hydrodynamic parameters, as suggested by Gardner's hydraulic model (1958). In the unsaturated configuration, the pressure head takes into account an effective stress parameter, which is equal to S e (Vanapalli and Fredlund 2000).
To use the TRIGRS model in our study, we relied on a prior analysis (Schilirò et al. 2018), representing the instability scenario arising from the 2011 event in the overall Monterosso catchment. We updated this analysis keeping the input parameters unaltered, except for those concerning the unsaturated phase, whose values were replaced with those obtained in our study. By so doing, we were able to check whether calibrating these parameters with field data could help better define landslide occurrence over large areas. For the spatial resolution of the input data, we selected 4 m × 4 m, which turned out to be the optimal resolution for performing this type of analysis (Schilirò et al. 2018).
Simulation results, expressed in terms of factor of safety (FS) map, were compared with the source areas of the landslides triggered during the event (Fig. 2). To quantify the predictive performance of the model, we carried out an receiver operating characteristic (ROC) curve analysis for different thresholds (i.e., different FS values). This analysis relates the true positive rate (TPR), which is the number of "landslide" cells that are correctly predicted, to the false positive rate (FPR), which is the number of "stable" cells that are erroneously simulated as unstable. The area under the curve (AUC) is the parameter that typically quantifies the performance of the model: the larger the value, the better the prediction.

Analysis of field monitoring data
We analyzed the data recorded by the monitoring station from August 1, 2018 to September 30, 2020. During this period, data was continuously acquired, except for breakdowns of individual sensors and short interruptions due to temporary outage of the power supply system. As regards the rainfall regime (Fig. 5a), several intense, short-duration events were recorded during summer and fall, corresponding to the Mediterranean climate. The most significant event occurred between October 27 and 29, 2018, when approximately 225 mm of rainfall were recorded in just three days (159 mm on October 29 itself). This event also showed extremely high hourly rainfall peaks, with a reached maximum value of 66.8 mm.
Rainfall directly induced sudden drops in the air and soil temperatures (Fig. 5b); however, a clear seasonal trend was observed. It is interesting to note that the air temperature had more significant daily changes than soil temperature, which reached approximately 30 °C in summer and 10 °C in winter.
As regards hydraulic monitoring, the data recorded by tensiometers (Fig. 5c) showed very high matric suction values (i.e., negative soil water pressure) under hot dry summer conditions (drying periods). These values, reaching up to − 900 hPa at the deepest sensor (T1), became approximately equal to zero immediately after the first fall rains, and were mostly unchanged until the next summer (wetting periods). Moreover, soil moisture records showed a clear relation with rainfall peaks (Fig. 5d). However, the deepest sensor (i.e., W1) recorded lower values (0.05 to 0.15) than the others (0.10 to 0.3).
If the spatial distribution of sensors is considered (Fig. 3), we can notice that those located downslope (i.e., W5 and W3) generally recorded higher water content values than the sensor located upslope at the same depth (i.e., W2-W4 and W1).
During one of the most important rainfall events (July 27-28, 2019), T2 and T4 were the first sensors to detect a significant pore water pressure change, approximately one hour after the first rainfall peak, given their position closer to the surface (Fig. 6). In contrast, T1 showed a sharp decrease in matric suction only at 3:00 am on July 28, when both T1 and T2 recorded their maximum values. T3, the deepest sensor upslope of the terrace, started to record pore water pressure changes only at 6:00 am on July 28, reaching maximum values approximately six hours later.

Reconstruction of the soil hydraulic behavior via numerical modeling
The data recorded over a two-year timespan made it possible to reconstruct the hydraulic behavior of the investigated soil in a comprehensive way, considering extremely variable atmospheric conditions. As mentioned above, this behavior was strongly related not only to the saturated phase but also to the unsaturated one. The analysis of unsaturated hydraulic properties is usually performed by reconstructing soil water characteristic curves (SWCCs), which describe the relationship between volumetric water content and matric suction. In this study, different SWCCs were determined in the upper (0-50 cm) and lower (50-100 cm) parts of the soil profile by coupling the records of T2-W2 and T1-W3, respectively. In particular, use was made of the nonlinear least-square method (Kemmer and Keller 2010) for fitting Eq. 3 to soil moisture-pore water pressure values, in order to determine the unknown parameters α and n. However, SWCCs can exhibit a hysteretic behavior depending on soil wetting and drying (Fredlund et al. 2011;Likos et al. 2013;Yang et al. 2019). In wetting periods, the volumetric water content induced by rainfall infiltration corresponded to a specific water pressure value, which was lower than that resulting from the same soil moisture level in drying periods. This type of process can notably affect the unsaturated shear strength of soils (Kristo et al. 2019). Thus, for each pair of sensors, we reconstructed a drying and wetting curve in relation to the period of the year in which the data had been recorded (Fig. 7). Some data scatter was observed; as also pointed out by Bordoni et al. (2015), this scatter can be correlated with very rapid hysteresis processes occurring after more intense rainfall events. Fitting parameters (Table 2), especially those concerning the wetting process of the lower soil portion, showed particularly high values, which can be explained by the presence of gravel in the investigated soil (Yang et al. 2019). To assess the effectiveness of the fitting procedure, we calculated the root mean square  (Table 2), thus substantiating the reliability of the fitted curves.
After defining the SWCCs of the upper and lower soil portions, we used the resulting fitting parameters as input data to the HYDRUS-1D model. More particularly, we first carried out a numerical simulation over a six-month period (i.e., August 1, 2018-January 31, 2019) to calibrate the α values, together with the saturated hydraulic conductivity under drying (K sd ) and wetting (K sw ) conditions. A subsequent simulation considered the remaining period (i.e., February 1, 2019-September 30, 2020), in order to validate all of the hydraulic parameters, with a view to using them in the following slope stability analyses (Sect. 4.3). In the latter simulation, we used W2 and W3 as control points of volumetric water content, while in the first simulation we employed both water pressure and water content data for calibrating the above-mentioned parameters via inverse optimization. A 100-cm soil profile inclined by 27° was chosen as geometric configuration of the investigated soil profile. A free drainage condition was considered along the bottom boundary, while hourly data was used as rainfall input. Evapotranspiration was accounted for by using daily maximum and minimum temperature records within the Hargreaves equation (Jensen et al. 1997).
The results of parameter calibration revealed that the lower part of the soil had a higher K sd value than the upper one (1.16 × 10 -5 vs. 9.49 × 10 -6 m/s), while K sw was equal

Slope stability analyses
After validating the hydraulic parameters of the soil on a preliminary basis, we reconstructed slope stability conditions through the simplified approach described in Sect. 3.2. As mechanical parameters of the soil, we used a friction angle of 31° (Schilirò et al. 2018), while effective cohesion was assumed to be equal to zero, considering the high degree of soil remolding associated with farming practices. As a first step, we calculated FS over the entire monitoring period in order to check the potential onset of instability conditions at depths of 50 and 100 cm. To define suction stress (σ S ), we generally used Eq. 7 with data recorded by W2 and W3. For near-saturation conditions only, we directly used positive water pressure values recorded by T1 and T2. It is worth pointing out that, in line with the above, we used the "wetting" parameters of van Genuchten-Mualem model during wetting periods, and vice versa. The resulting trends ( Fig. 9) showed that FS always lays above 1, which is consistent with field evidence during the two years of observation. The trend was more stable at a depth of 100 cm, whereas the upper part of the soil had larger FS changes depending on its rapid water content fluctuations. In other words, individual rainfall events were supposed to cause sudden and significant FS drops even during spring-summer periods, when the highest FS values were reached. In winter, the FS trend at a depth of 50 cm was variable, but within a smaller range.
The same approach was taken to reconstruct slope stability conditions at local scale during the rainfall event that took place on October 25, 2011 (see Sect. 2). In this instance, hydraulic conditions were reconstructed with HYDRUS simulations, employing the calibrated model described in the preceding section. We first evaluated the initial soil moisture conditions (i.e., prior to the intense rainfall event) through a numerical simulation using daily rainfall data in the previous month. Afterward, FS was calculated during the event (7:00-17:00 UTC) using hourly rainfall data provided by ARPAL (Agenzia Regionale di Protezione dell'Ambiente Ligure). Results (Fig. 10a) indicate that FS at a depth of 50 cm rapidly diminished after 9:00 UTC (start of the most critical stage of the rainfall event) and reached a value of less than 1 at 11:00 UTC. Conversely, the reduction of FS at a depth of 100 cm was less marked and started about one hour later (10:00 UTC). Also in this instance, the most critical value was attained at 11:00 UTC, but it remained above 1 (i.e., 1.18). At the end of rainfall, FS tended to increase in both cases, with a greater gradient at a depth of 50 cm. Nevertheless, the local conditions of the investigated slope did not represent the overall Monterosso catchment. This is why different results can be achieved when considering a slope angle of 38° (i.e., the average slope angle measured in the 2011 landslide source areas; see Fig. 4c).
By replacing this value with the preceding one (27°), representing the monitored terraced slope, we observed that critical conditions also arose at a depth of 100 cm as early as at 10:00 UTC (Fig. 10b). Taking into account that soil thickness in source areas was equal to or greater than 1 m (Fig. 4d) on average, we concluded that, also in this instance, simulation results were consistent with field evidence.
After local-scale analysis, we decided to simulate the 2011 event in the overall Monterosso catchment using the TRIGRS model. Results highlighted a good agreement between TRIGRS cells at FS < 1 and the location of landslides (Fig. 11). This performance was corroborated by the area under curve (AUC) metrics obtained from the ROC (receiver operator characteristic) curve; the latter returned a value of 79.3% (Table 3), given by a reduced number of false positives and the absence of landslides above FS > 2.

Discussion
The data recorded by tensiometers and water content sensors (Fig. 5c, d) yielded some preliminary indications about the hydrological behavior of the soil. In the first place, the lag time between rainfall events and hydraulic changes was very low, which is consistent with the hydraulic conductivity typical of coarse-grained materials (Table 1). We identified an upper soil portion (0-50 cm) with a K sd value of 9.49 × 10 -6 m/s, whereas the lower part (50-100 cm) was slightly more permeable (1.16 × 10 -5 m/s). As regards the parameters of the unsaturated phase ( Fig. 7 and Table 2), the difference between the two parts proved to depend on the different distribution of pore sizes as a function of depth, which might in turn be correlated with grain-size changes. The high α values obtained in the bottom part suggested the presence of material with larger pores. The latter may be ascribed to the method used to build farming terraces in Cinque Terre, which involves the upslope excavation of the filling material. The geometry of the terrace could also explain the differences in water content trends. The results shown in Figs. 5 and 6 best represent a wetting front infiltration scheme (Schilirò et al. 2019), in which a reduction of matric suction first occurs in the upper part of the soil and, if rainfall persists, positive pore water pressures can also develop (Fig. 12a) due to the accumulation of infiltrating water on the edge of the terrace (Fig. 12b). In this respect, we are aware that a detailed description of this process should include an evaluation of the role of the terrace flank, which probably acts as a low permeability boundary. However, this would require a 2D-3D analysis going beyond the scope of this work. In fact, the main objective of our study was to better define the hydraulic parameters of the soil based on recorded data and then use them in slope stability analyses at catchment scale. Nevertheless, in the light of data collected during specific rainfall events, we deemed it reasonable to assume the above-described infiltration scheme.
With regard to slope stability, the analysis performed over the two-year monitoring period showed strong FS changes upon significant rainfall peaks and/or prolonged rainfall periods (Fig. 9). These changes were more marked in the shallowest part of the soil, where the contribution of partial saturation to soil shear strength was higher. During the investigated period, no landslide phenomena were recorded in the area, which agrees with the absence of FS values below 1 in the modeling results. By contrast, during the 2011 event, the slope accommodating the monitoring station experienced very shallow failures. The latter evidence was consistent with the outcome of numerical simulations, which confirmed that FS was below 1 at a depth of 50 cm, upon the 11:00 UTC rainfall peak, whereas the FS reduction at a depth of 100 cm was not sufficient to induce instability (Fig. 10a). After these analyses, we decided to simulate the 2011 event once again, with the slope conditions of the investigated terrace, but assuming a change in initial soil moisture, with a view to gaining further insight into its role in the onset of failure under the same triggering event. Considering that HYDRUS simulations returned an initial S e equal to 0.45 at the beginning of the October 25, 2011 event, we selected drier (S e : 25%) and wetter (S e : 75%) initial conditions in slope stability calculation, referring to a mean slope angle of 27°. At the same time, we used both "drying" and "wetting" retention curve parameters (Table 2) to also quantify the effect of the hysteretic behavior of the soil. The resulting trends showed, once again, that the highest FS changes occurred in the upper part of the soil, where the initial FS value dropped by about 1, passing from S e 25% to 75% (Fig. 13a). It is also worth noting that the use of "wetting" parameters significantly reduced stability, especially with an initial S e of 25%. In the latter instance, the reduction was sufficient to induce failure at 11:00 UTC, as in the other instances featured by a higher S e . This finding stresses the importance of considering not only existing soil moisture conditions, but also the proper soil hysteresis path in relation to the season in which the rainfall event occurs.
As regards FS at a depth of 100 cm (Fig. 13b), the differences with variable initial soil moisture and "drying"/ "wetting" parameters were less pronounced, and failure Fig. 12 Wetting front infiltration scheme considering the presence of: a soil water content sensors, and b tensiometers. The blue arrows represent rainfall input conditions were never reached. This point can be explained by the infiltration scheme inferred from recorded data and replicated in HYDRUS simulations i.e., the advance of the wetting front. As stated above, the presence of a wetting front moving downward during intense rainfall events develops positive pore water pressures mostly in the shallowest part of the soil. These pressures become progressively lower with depth, where the instability process mainly depends on the loss of matric suction. This means that the lower part of the soil is less susceptible to failure, unless other types of hydraulic processes take place (i.e., temporary water table rising). The data recorded by the monitoring station did not allow us to hypothesize this type of process. However, it is worthwhile emphasizing that only rainfall data was available for the most significant rainfall event that occurred in the investigated period (October 27-29, 2018). Indeed, the hydraulic data of the soil was not collected during the event due to a temporary failure of the sensors. Hence, the formation of a perched water table cannot be excluded a priori.
As regards the reconstruction of the 2011 landslide events in the overall Monterosso catchment via TRIGRS, the good performance of our simulation was confirmed by the AUC value, which was higher than that obtained in a similar analysis carried out by Schilirò et al. 2018 (Table 3). The improvement was related to the FPR value, which was much lower in the new analysis. Specifically, our analysis predicted a reasonable percentage of source areas (42.4%) with a very low number of false positives. Conversely, although the previous analysis correctly simulated a slightly higher percentage of landslides (49.5%), the extent of areas erroneously predicted as "unstable" was much larger (Fig. 14a). This means that the use of updated hydraulic parameters improves the performance of the model. Furthermore, a strong reduction in false alarms is certainly welcomed in view of the potential application of the model in EWSs.
Finally, it is important to stress that the above-described insights rely on interpretation of field monitoring data and simplified modeling as part of a wider and still ongoing research activity. The next step will be the confirmation of the hypothesized infiltration scheme (Fig. 12) via 2D-3D numerical analyses, in order to account for spatial changes in the hydraulic properties of the soil across the terrace, with particular regard to its flanks. These analyses will be supported by an in-depth physical and hydraulic characterization of Fig. 13 FS trends at depths of 50 cm (a) and 100 cm (b), calculated during the October 25, 2011 event, assuming different initial soil moisture conditions (in terms of Se) and "drying"/ "wetting" parameters the terrace in different sectors, alongside with a series of laboratory flume tests. The latter will be conducted on a soil slope model reproducing field conditions, following the procedure described in Schilirò et al. 2019. The results of the planned activities will provide useful information to enhance the reliability of analyses over large areas.

Conclusions
In this paper, we have presented and analyzed the data recorded by an in situ monitoring station located in the Monterosso catchment, making part of Cinque Terre National Park, in northern Italy. The aim of our study was to obtain some preliminary insights into the hydraulic behavior of soil covers and shallow landslide triggering mechanisms, to be used as a basis for stability analyses over large areas. Starting from atmospheric data and soil hydraulic data collected over a two-year period, we first built a hydraulic model of the soil, paying particular attention to characterizing its unsaturated phase. The resulting model, consisting of an upper and lower part (both with hydraulic parameters typical of coarse-grained materials), was obtained from numerical simulations with HYDRUS. The simulations also allowed us to infer the advance of a wetting front as a potential infiltration scheme. As regards landslide triggering conditions, local-scale stability analyses showed that the strongest FS changes occurred in the upper part of the soil, where the loss of matric suction was faster and the development of positive pore water pressures more pronounced. Simulation of the October 25, 2011 event emphasized the major role of initial moisture and hysteretic behavior of the soil, as well as local slope angle, in slope stability conditions. The same event was modeled at catchment scale. In the latter modeling exercise, the use of parameters calibrated with field data improved the predictive performance of the TRIGRS model vs. prior analyses carried out in the same area. This evidence stressed the importance of hydraulically characterizing weathered soil covers that are landslide-prone, with a view to designing more reliable future scenarios. In view of the above, future efforts will be devoted to an in-depth physical and hydraulic characterization of the investigated slope, taking into account potential spatial changes in soil parameter values in relation to the geometry of terraces. These efforts will hopefully enable us to refine and improve the hydraulic model of the soil discussed in this paper and to better define potential conditions for rainfall-triggered shallow landslides initially at local scale and then in the overall catchment. Therefore, this work should be considered as a first attempt to couple field monitoring data with simplified physically based models, with a view to developing EWSs specifically calibrated for the morpho-climatic features of a wide area.