Monitoring the Hydrological Balance of a Landslide-Prone Slope Covered by Pyroclastic Deposits over Limestone Fractured Bedrock

: Many mountainous areas in Campania, Southern Italy, are characterized by steep slopes covered by loose unsaturated pyroclastic deposits laying upon fractured limestone bedrock. The soil covers are mainly constituted by layers of ashes and pumices. Large and intense rainfall events trigger shallow landslides, often turning into debris ﬂows that cause huge damage and casualties. The slope of Cervinara, around 40 km Northeast of Naples, was involved in a catastrophic ﬂowslide on 16 December 1999, triggered by a rainstorm of 325 mm in 48 h. To capture the main e ﬀ ects of precipitation on the slope stability, hydro-meteorological monitoring activities have been carried out at the slope to assess the water balance for three years (2017–2020). The ﬁeld monitoring data allowed the identiﬁcation of the complex hydrological processes involving the unsaturated pyroclastic soil and the shallow groundwater system developing in the limestone bedrock, which control the conditions that potentially predispose the slope to landslide triggering. Speciﬁcally, late autumn has been identiﬁed as the potentially most critical period, when slope drainage processes are not yet e ﬀ ective, and soil covers already receive large amounts of precipitation.


Introduction
The numerous hydrological processes involving slopes in response to meteorological forces are extremely important to understand the predisposing conditions of rainfall-induced landslide initiation [1][2][3][4][5][6][7], neglected or only indirectly considered even in the most recent landslide early warning systems (e.g., [8]). The hydrological predisposing causes are frequently connected with an area exceeding the failed slope, i.e., they are related to large-scale processes (in time and space). The need to assess the water balance, including large-scale hydrological processes, is true for both shallow and deep-seated landslides. On one hand, deep-seated landslides, commonly referred to as deep-seated gravitational slope deformations, are very slow mass movements of a large volume (usually >0.5 km 3 ) of soil and rock [9][10][11][12][13]. The deformation processes in slopes present various forms [14], and are considered to be precursors of eventual damages [15][16][17]. Several authors have studied the transition process from gravitational deformation to slide, especially from the mechanical viewpoint [18]. Causal factors include: lithology and geologic structure [19], tectonic and topographic stresses and their concentration [20,21], seismicity [22], changes in groundwater regime [23], and human activity [24]. The hydrogeological response of the groundwater regime for the initiation of deep-seated landslides is strongly correlated with antecedent precipitation history [25,26]. Moreover, the groundwater flow in a The slope is covered with shallow deposits of pyroclastic materials generated from several eruptions of two main volcanic complexes of Campania, the Somma-Vesuvius and the Phlegrean Fields, which occurred during the last 40,000 years [63]. The soil cover, usually in an unsaturated condition, is mainly constituted by layers of volcanic ash of high porosity (up to 75%), and pumice lapilli, laying upon fractured calcareous bedrock. The total thickness of the cover is between 1.5 m and 2.5 m in the steepest part of the slope, while at its foot it reaches several meters [54]. The inclination of the slope ranges from 35° to 45° at an altitude between 500 m and 750 m. The climate of the area is Mediterranean, with the precipitation regime characterized by a dry season (from May/June to September), and a wet season (from October to the end of April/May). The mean annual precipitation in the area is around 1600 mm, and the total potential evapotranspiration , determined with the Thornthwaite formula, is between 700 mm (at the altitude of 750 m) and 800 mm (at 360 m) [48]. In the area, extreme daily rainfall heights of about 130 mm and 170 mm correspond to return periods of 10 and 50 years, respectively. The vegetation cover consists of cultivated chestnut The slope is covered with shallow deposits of pyroclastic materials generated from several eruptions of two main volcanic complexes of Campania, the Somma-Vesuvius and the Phlegrean Fields, which occurred during the last 40,000 years [63]. The soil cover, usually in an unsaturated condition, is mainly constituted by layers of volcanic ash of high porosity (up to 75%), and pumice lapilli, laying upon fractured calcareous bedrock. The total thickness of the cover is between 1.5 m and 2.5 m in the steepest part of the slope, while at its foot it reaches several meters [54]. The inclination of the slope ranges from 35 • to 45 • at an altitude between 500 m and 750 m. The climate of the area is Mediterranean, with the precipitation regime characterized by a dry season (from May/June to September), and a wet season (from October to the end of April/May). The mean annual precipitation in the area is around 1600 mm, and the total potential evapotranspiration ET 0 , determined with the Thornthwaite formula, is between 700 mm (at the altitude of 750 m) and 800 mm (at 360 m) [48]. In the area, extreme daily rainfall heights of about 130 mm and 170 mm correspond to return periods of 10 and 50 years, respectively. The vegetation cover consists of cultivated chestnut trees (Castanea sativa). A dense Water 2020, 12, 3309 4 of 23 underbrush develops from late April to September, mainly formed by ferns (Pteridium aquilinum) and other seasonal shrubs.

Hydro-Meteorological Monitoring Station
Hydrological monitoring activities have been carried out at the slope of Cervinara since 2002, by measuring precipitation depth, soil water content and capillary tension [47,54]. Since December 2017, an automatic hydro-meteorological station has been installed at the altitude of about 585 m, immediately near the scarp of the 1999 landslide ( Figure 1). The acquisition, with a time resolution of 1 h, and the data storage are ensured by a data logger, powered by a 12 V battery connected to two solar panels through a charge controller. Hydrological monitoring equipment includes five "Jet-fill" tensiometers equipped with tension transducers, a self-refilling tensiometer, and six probes for time domain reflectometry (TDR), connected through coaxial cables and a multiplexer to a reflectometer. The TDR probes are constituted by three metallic rods, with length of 0.10 m and diameter of 3 mm, spaced 1.5 cm apart. TDR is a widely used indirect technique for field measurements of soil water content [64], providing soil bulk dielectric permittivity, ε r , which can be related to the soil volumetric water content, θ, through specific calibration expressions [65][66][67]. The meteorological station consists of a rain gauge with a resolution of 0.197 mm, thermo-hygrometer, soil temperature sensor, anemometer, pyranometer and barometric sensor. Table 1 summarizes the main characteristics of the installed devices. The solar panels and the meteorological sensors are fixed to a metal pole between 1.50 m and 2.50 m above the ground. The rain gauge is located on top of another metal pole, 2.00 m above the ground surface. Differently, the thermistor is buried under the soil surface at depth of 0.10 m (Figure 2).  For the installation of the soil hydrological monitoring devices, a trench was excavated, and the TDR probes were horizontally pushed into its walls at various depths between 0.30 m and 2.00 m below the ground surface. Nearby, the tensiometers were installed vertically into the soil at the bottom of boreholes, previously excavated by a manual drill, between 0.20 m and 3.0 m below the ground surface.
From the top to the bottom of the soil cover, the following layers were found during the excavation: (i) layer A, organic top soil derived from volcanic ashes, with thickness of about 0.90 m; (ii) layer B, coarse pumices, about 0.40 m thick; (iii) layer C, volcanic ashes, about 0.20 m thick; (iv) layer D, fine pumices, with thickness of about 0.10-0.15 m; (v) layer E, volcanic ashes, more than 0.50 m thick. Unluckily, the total excavation depth (i.e., 2.15 m) did not allow to reach the underlying fractured calcareous bedrock. Figure 3 reports a schematic cross-section of the local stratigraphy with the location along the depth of the tensiometers and the TDR probes.   For the installation of the soil hydrological monitoring devices, a trench was excavated, and the TDR probes were horizontally pushed into its walls at various depths between 0.30 m and 2.00 m below the ground surface. Nearby, the tensiometers were installed vertically into the soil at the bottom of boreholes, previously excavated by a manual drill, between 0.20 m and 3.0 m below the ground surface.
From the top to the bottom of the soil cover, the following layers were found during the excavation: (i) layer A, organic top soil derived from volcanic ashes, with thickness of about 0.90 m; (ii) layer B, coarse pumices, about 0.40 m thick; (iii) layer C, volcanic ashes, about 0.20 m thick; (iv) layer D, fine pumices, with thickness of about 0.10-0.15 m; (v) layer E, volcanic ashes, more than 0.50 m thick. Unluckily, the total excavation depth (i.e., 2.15 m) did not allow to reach the underlying fractured calcareous bedrock. Figure 3 reports a schematic cross-section of the local stratigraphy with the location along the depth of the tensiometers and the TDR probes.
Disturbed samples were collected to retrieve the grain size distributions of the various layers, as shown in Figure 4.
The coarsest soil is the pumice layer B, which can be classified as gravel with sand, while the finest is the ash layer E, which is silty sand.
The ceramic tips of the tensiometers were located at depths of 0.20 m, 0.50 m, and 0.80 m in soil layer A, at 1.60 m (layer D) and at 1.80 m (layer E) of depth below the surface. The deepest one was located at 3.00 m of depth. The TDR probes were buried at depths of 0.30 m, 0.75 m, and 0.95 m in the shallowest layer (A), at 1.10 m in the coarse pumice layer (B), at 1.40 m in the volcanic ash layer (C), and, the last one at 2.00 m depth below the ground (layer E). Water 2020, 12, x FOR PEER REVIEW 6 of 24 Disturbed samples were collected to retrieve the grain size distributions of the various layers, as shown in Figure 4. The coarsest soil is the pumice layer B, which can be classified as gravel with sand, while the finest is the ash layer E, which is silty sand.
The ceramic tips of the tensiometers were located at depths of 0.    Disturbed samples were collected to retrieve the grain size distributions of the various layers, as shown in Figure 4. The coarsest soil is the pumice layer B, which can be classified as gravel with sand, while the finest is the ash layer E, which is silty sand.
The ceramic tips of the tensiometers were located at depths of 0.

Surface and Groundwater Circulation
In the surrounding area, a network of streams is present that intercept runoff and water from several ephemeral and perennial springs located in the surroundings ( Figure 5, Table 2). The largest perennial springs, located at low elevations, are supplied by the deep groundwater system typical of the limestone massifs of the area [39]. Aiming at the assessment of the hydrological balance of the slope, the water levels of two streams located at the foot of the slope have been manually monitored every three or four weeks since December 2017 ( Figure 5). Afterwards, in one of the two streams, located at 550 m a.s.l. (Section 1), a multiparameter probe, measuring water pressure, temperature and electrical conductivity, was installed in March 2019 to monitor the hydrologic regime of the stream with a finer time resolution (one hour), and to identify the origin of the observed discharge, i.e., overland or subsurface runoff, or groundwater flow, which typically present electrical conductivities that differ more than one order of magnitude from each other [68][69][70].
several ephemeral and perennial springs located in the surroundings ( Figure 5, Table 2). The largest perennial springs, located at low elevations, are supplied by the deep groundwater system typical of the limestone massifs of the area [39]. Aiming at the assessment of the hydrological balance of the slope, the water levels of two streams located at the foot of the slope have been manually monitored every three or four weeks since December 2017 ( Figure 5). Afterwards, in one of the two streams, located at 550 m a.s.l. (Section 1), a multiparameter probe, measuring water pressure, temperature and electrical conductivity, was installed in March 2019 to monitor the hydrologic regime of the stream with a finer time resolution (one hour), and to identify the origin of the observed discharge, i.e., overland or subsurface runoff, or groundwater flow, which typically present electrical conductivities that differ more than one order of magnitude from each other [68][69][70].

Monitoring Data
In this section, data from the field monitoring are described. Monitoring started on 7 December 2017 and is still ongoing. During this time, the monitoring has been quite continuous, except for three periods where the data logger was out of order. Moreover, the instrumentation was integrated with

Monitoring Data
In this section, data from the field monitoring are described. Monitoring started on 7 December 2017 and is still ongoing. During this time, the monitoring has been quite continuous, except for three periods where the data logger was out of order. Moreover, the instrumentation was integrated with a thermo-hygrometer in May 2018. The collected data have been used to calculate the various terms of the water balance of the slope, consisting of the precipitation input, the stream discharge, the water exchanges with the atmosphere and with the deep groundwater system, and the water stored in the soil cover and in the perched aquifer temporarily developing in the upper part of the fractured limestone (epikarst). Figure 6 shows the hourly rainfall recorded by the rain gauge during the monitored period. exchanges with the atmosphere and with the deep groundwater system, and the water stored in the soil cover and in the perched aquifer temporarily developing in the upper part of the fractured limestone (epikarst). Figure 6 shows the hourly rainfall recorded by the rain gauge during the monitored period. In the first case, 130 mm of rain fell, with a maximum hourly intensity of 18 mm/h. In the second case, 100 mm and 175 mm of rain were recorded in the two days, fallen after one and a half months of frequent precipitation (around 500 mm from the beginning of November until 19 December). After a maximum hourly intensity value of 41 mm/h occurred at 02:00 PM of 21 December, a shallow landslide was triggered along a slope about one kilometre far from the slope of Cervinara. It turned into a flowslide, and, after reaching the foot, it was channelized through a stream, reaching the nearby town of San Martino Valle Caudina. Moreover, the collected rainfall data are close in agreement with the data from the rain gauge of the regional Civil Protection Agency of Cervinara, not far from the monitoring area.

Evapotranspiration
To assess the hydrologic water balance accurately, actual evapotranspiration has to be estimated, based on the available meteorological data. The potential evapotranspiration has been estimated with the formulas of Hargreaves [71] (1), Blaney-Criddle [72,73] (2), with the Penman-Monteith equation [74] (3), and with the empirical formula of Thornthwaite [75] (4): In the previous equations, = 0.0820 MJ/m min represents the solar constant; (m) is the inverse relative Earth-Sun distance; (rad) is the sunset hour angle; (rad) is the latitude; (rad) is the solar radiation. In the first case, 130 mm of rain fell, with a maximum hourly intensity of 18 mm/h. In the second case, 100 mm and 175 mm of rain were recorded in the two days, fallen after one and a half months of frequent precipitation (around 500 mm from the beginning of November until 19 December). After a maximum hourly intensity value of 41 mm/h occurred at 02:00 PM of 21 December, a shallow landslide was triggered along a slope about one kilometre far from the slope of Cervinara. It turned into a flowslide, and, after reaching the foot, it was channelized through a stream, reaching the nearby town of San Martino Valle Caudina. Moreover, the collected rainfall data are close in agreement with the data from the rain gauge of the regional Civil Protection Agency of Cervinara, not far from the monitoring area.

Evapotranspiration
To assess the hydrologic water balance accurately, actual evapotranspiration has to be estimated, based on the available meteorological data. The potential evapotranspiration has been estimated with the formulas of Hargreaves [71] (1), Blaney-Criddle [72,73] (2), with the Penman-Monteith equation [74] (3), and with the empirical formula of Thornthwaite [75] (4): In the previous equations, G cs = 0.0820 MJ/m min represents the solar constant; d r (m) is the inverse relative Earth-Sun distance; w s (rad) is the sunset hour angle; ϕ (rad) is the latitude; δ (rad) is the solar radiation.
T ai 5 1.514 In the previous equations, I is the heat index, which depends on the mean monthly temperatures T ai ( • C). In this case, as no continuous temperature data were collected, the monthly temperatures Water 2020, 12, 3309 9 of 23 recorded between 1989 and 1999 from Montesarchio station (560 m a.l.s), about 10 km from Cervinara, were used. In Table 3, the parameters required by Equations (1)-(4) are summarized. Table 3. Summary of the parameters (determined using the meteorological data collected by the Cervinara station or from literature or evidence on site).

Blaney-Criddle (mm/day)
RH min = daily minimum relative humidity (%) * U z = daily wind speed at 2 m above the ground (m/s) * r i = ratio between monthly average daily bright sunshine duration n(h) and monthly average maximum daily sunshine duration N(h) The leaf area index (LAI), required for the calculation of ET 0 with Penman-Monteith Equation (3), has been assumed equal to 4 m 2 /m 2 between May and August when foliage and dense underbrush are present (as indicated by [76,77] for Italian Castanea Sativa), while it has been considered negligible from November to February, when the trees are leafless. Smaller values have been assumed when foliage and underbrush are not fully developed, namely 1.5 m 2 /m 2 in March and October, and 3 m 2 /m 2 in April and September. The tree height has been assigned equal to 18 m through visual inspection, in agreement with [78]. Moreover, as the magnitude of daily soil heat flux G over a period of 10-30 days is very small compared to net radiation flux, R n , it has been neglected [79]. The short-wave radiation reflection coefficient α (albedo) has been considered equal to 0.23 according to vegetation cover type [79,80]. Figure 7 reports the average monthly potential evapotranspiration rate estimated by Equations (1)-(4).
underbrush are present (as indicated by [76,77] for Italian Castanea Sativa), while it has been considered negligible from November to February, when the trees are leafless. Smaller values have been assumed when foliage and underbrush are not fully developed, namely 1.5 / in March and October, and 3 / in April and September. The tree height has been assigned equal to 18 m through visual inspection, in agreement with [78]. Moreover, as the magnitude of daily soil heat flux G over a period of 10-30 days is very small compared to net radiation flux, , it has been neglected [79]. The short-wave radiation reflection coefficient α (albedo) has been considered equal to 0.23 according to vegetation cover type [79,80]. Figure 7 reports the average monthly potential evapotranspiration rate estimated by Equations (1)-(4). The main pro of using the simple formulas of Hargreaves (1) and Thornthwaite (4) is the minimum data requirement, i.e., only air temperature data are needed. However, the application of the Blaney-Criddle formula (2) and Penman-Monteith equation (3) provides more reliable estimates, as both take into account the effects of actual insolation time (ratio n/N between actual and maximum possible bright sunshine hours) and all the measured meteorological variables, while the latter also allows combining them with the estimates of stomatal resistance to evapotranspiration by vegetation. The potential evapotranspiration values estimated by Penman-Monteith and Blaney-Criddle equations are in many months larger than the ones provided by expressions (1) and (4), especially from April to September, with the highest evapotranspiration rate (10 mm/day) calculated in May 2020. Indeed, the three monitored years were much warmer than the average. The lowest estimated evapotranspiration rates, namely 0.02 mm/day and 0.1 mm/day, respectively, in December 2019 and January 2020, confirm the interplay of the absence of foliage and underbrush with low solar radiation The main pro of using the simple formulas of Hargreaves (1) and Thornthwaite (4) is the minimum data requirement, i.e., only air temperature data are needed. However, the application of the Blaney-Criddle formula (2) and Penman-Monteith Equation (3) provides more reliable estimates, as both take into account the effects of actual insolation time (ratio n/N between actual and maximum possible bright sunshine hours) and all the measured meteorological variables, while the latter also allows combining them with the estimates of stomatal resistance to evapotranspiration by vegetation. The potential evapotranspiration values estimated by Penman-Monteith and Blaney-Criddle equations are in many months larger than the ones provided by expressions (1) and (4), especially from April to September, with the highest evapotranspiration rate (10 mm/day) calculated in May 2020. Indeed, the three monitored years were much warmer than the average. The lowest estimated evapotranspiration rates, namely 0.02 mm/day and 0.1 mm/day, respectively, in December 2019 and January 2020, confirm the interplay of the absence of foliage and underbrush with low solar radiation during the winter season [73,81]. To consider the specific climatic features of the monitoring years, the monthly values of ET 0 from the Blaney-Criddle formula have been adopted for the assessment of the water balance, except for the periods where the meteorological station was out of order-when the mean estimates of Thornthwaite formula have been used. Indeed, the Penman-Monteith equation suffered from too long periods where there was a lack of data due to the discontinuous functioning of the pyranometer (Figure 7). To obtain actual evapotranspiration ET, the estimated ET 0 must be multiplied by a crop coefficient, depending on the characteristics of the vegetation and on the climatic conditions. The suction measurements (Section 3.4) indicate that the soil only rarely became so dry it induced stress to the vegetation (i.e., this could have occurred only sporadically from late July to September). In this period, ET can therefore be assumed close to potential evapotranspiration ET 0 . Conversely, the dormant vegetation implies that in autumn and winter ET is likely quite smaller than ET 0 . Hence, the crop coefficients have been assumed between 0.3 (from December to February) and 0.9 (from May to July). The estimated annual ET 0 and ET are reported in Table 4.  Figure 8 shows the manually measured water levels in the two monitored stream sections-1 and 2 ( Figure 5).  Figure 8 shows the manually measured water levels in the two monitored stream sections-1 and 2 ( Figure 5). The water level in stream Section 1, located at lower altitude than Section 2, starts growing in autumn (i.e., between late September and early December, depending on the year), and reaches its maximum between late winter and early spring. In late spring and summer, the water level progressively decreases, and the stream becomes dry in late summer. In Section 2, instead, water flow has been observed only sporadically, after extremely rainy periods. Given the small extension of the catchments of the two sections, the delay and persistence of the water level variations, compared to the precipitation, seem to indicate that the discharge in the streams cannot be considered as the result of surface or subsurface runoff. To estimate the stream discharge [L/T] from the water level readings H [L] reported in Figure 8, a monomial stage-discharge relationship has been defined for each hydrological year:

Stream Discharges
where and represent the maximum water level, and the corresponding estimated flow discharge, and α [-] is an exponent assumed to be equal to 2.5. The maximum flow discharge value The water level in stream Section 1, located at lower altitude than Section 2, starts growing in autumn (i.e., between late September and early December, depending on the year), and reaches its maximum between late winter and early spring. In late spring and summer, the water level progressively decreases, and the stream becomes dry in late summer. In Section 2, instead, water flow has been observed only sporadically, after extremely rainy periods. Given the small extension of the catchments of the two sections, the delay and persistence of the water level variations, compared to the precipitation, seem to indicate that the discharge in the streams cannot be considered as the result of surface or subsurface runoff. To estimate the stream discharge Q s [L/T] from the water level readings H [L] reported in Figure 8, a monomial stage-discharge relationship has been defined for each hydrological year: where H max and Q max represent the maximum water level, and the corresponding estimated flow discharge, and α [-] is an exponent assumed to be equal to 2.5. The maximum flow discharge value has been estimated by measuring the cross-sectional area and evaluating the mean flow velocity. Figure 9 reports the estimated discharges of the two monitored streams.
Water 2020, 12, x FOR PEER REVIEW 12 of 24 has been estimated by measuring the cross-sectional area and evaluating the mean flow velocity. Figure 9 reports the estimated discharges of the two monitored streams.     . Estimated flow discharges of the two monitored stream sections (the grey shaded areas indicate periods with a lack of data). Figure 10 shows the water level and electrical conductivity in stream Section 1 ( Figure 5), measured by the multiparametric probe between April 2019 and August 2020. The values of electrical conductivity are mostly between 0.2 and 0.4 mS/cm, confirming that the stream discharge is likely originated by baseflow and not by surface runoff, which would be associated with electrical conductivity two orders of magnitude smaller [69]. During rain events, when the water level peaks, the electrical conductivity quickly decreases, indicating that rainwater mixes with the water flowing through the stream. Some large rainfall events have been chosen to infer the importance of overland and subsurface runoff in the observed discharge generation ( Figure  11). Specifically, for each event, runoff coefficients C [-] have been estimated as the ratio between total specific flood runoff volume (calculated on the basis of the area between the dashed lines, indicating The values of electrical conductivity are mostly between 0.2 and 0.4 mS/cm, confirming that the stream discharge is likely originated by baseflow and not by surface runoff, which would be associated with electrical conductivity two orders of magnitude smaller [69]. During rain events, when the water level peaks, the electrical conductivity quickly decreases, indicating that rainwater mixes with the water flowing through the stream. Some large rainfall events have been chosen to infer the importance of overland and subsurface runoff in the observed discharge generation ( Figure 11). Specifically, for each event, runoff coefficients C [-] have been estimated as the ratio between total specific flood runoff volume (calculated on the basis of the area between the dashed lines, indicating the estimated baseflow trend, and the total discharge line in the graphs of Figure 11) and the precipitation depth of each event (Table 5).  The values reported in Table 5 indicate that the runoff coefficients are in all cases very small, although the selected events are among the largest occurred during the monitoring period. In fact, just a small discharge growth during rainfall is observed, while a stable increase is established after The values reported in Table 5 indicate that the runoff coefficients are in all cases very small, although the selected events are among the largest occurred during the monitoring period. In fact, just a small discharge growth during rainfall is observed, while a stable increase is established after the end of each event, which can be interpreted as a relatively fast baseflow. During rainfall events smaller than the selected ones, no significant increase was observed.
Consequently, the response of stream discharge to precipitations seems mostly supplied by baseflow, which, given the relatively fast response, can be related to groundwater circulation occurring in the shallower part of the limestone bedrock. Hence, both overland and subsurface runoff contributions can be neglected in the assessment of long-term water balance. Figure 12 shows, from top to bottom panels, the hourly rainfall, soil water and suction measurements during the monitored period. the end of each event, which can be interpreted as a relatively fast baseflow. During rainfall events smaller than the selected ones, no significant increase was observed.

Soil Water Content and Suction
Consequently, the response of stream discharge to precipitations seems mostly supplied by baseflow, which, given the relatively fast response, can be related to groundwater circulation occurring in the shallower part of the limestone bedrock. Hence, both overland and subsurface runoff contributions can be neglected in the assessment of long-term water balance. Figure 12 shows, from top to bottom panels, the hourly rainfall, soil water and suction measurements during the monitored period. Despite the lack of data during the autumns in 2018 and 2019, it is possible to recognize seasonal suction fluctuations characterized by trends similar to those observed in other pyroclastic covers in Campania [46,82,83]. Specifically, suction values in the deepest layer always remained lower than that recorded in the uppermost part of the cover, with the exception of short periods during the rainy seasons, when the rainwater infiltration did not yet affect the deepest part of the cover. Furthermore, given the very high porosity of the deposit (i.e., up to 0.75 for the ashes) and the mean water content recorded during the wet seasons, which never exceeded 0.50, the soil was always far from saturation. Figure 13 shows the soil response in terms of water content and suction observed during the wet period (25 March-30 April 2019, characterized by a total rainfall depth of 123 mm). Despite the lack of data during the autumns in 2018 and 2019, it is possible to recognize seasonal suction fluctuations characterized by trends similar to those observed in other pyroclastic covers in Campania [46,82,83]. Specifically, suction values in the deepest layer always remained lower than that recorded in the uppermost part of the cover, with the exception of short periods during the rainy seasons, when the rainwater infiltration did not yet affect the deepest part of the cover. Furthermore, given the very high porosity of the deposit (i.e., up to 0.75 for the ashes) and the mean water content recorded during the wet seasons, which never exceeded 0.50, the soil was always far from saturation. Figure 13 shows the soil response in terms of water content and suction observed during the wet period (25 March-30 April 2019, characterized by a total rainfall depth of 123 mm).

Soil Water Content and Suction
The effects of infiltration in response to rainfall generally depend on depth, as fluctuations of soil water content and suction are delayed and smoothened (less marked) in the deepest layers, while in the shallowest ones, they tend to be rapid and higher (values of suction decrease quickly and immediately after the observed rainfall). Figure 14 shows the soil response observed during the dry period (1 June-31 July 2019).
In this case, it is easy to note that the suction starts to increase earlier and faster in the shallowest layers, and then in the deepest ones. This response can be related to the upward evapotranspiration fluxes due to the water uptake from the roots, which are denser in the upper part of the soil cover due to the flourishing underbrush. Although the soil moisture dynamics are strictly connected to the specific precipitation trends of each hydrologic year, it is possible to define different periods with a similar profile in terms of soil water content and suction ( Table 6). The effects of infiltration in response to rainfall generally depend on depth, as fluctuations of soil water content and suction are delayed and smoothened (less marked) in the deepest layers, while in the shallowest ones, they tend to be rapid and higher (values of suction decrease quickly and immediately after the observed rainfall). Figure 14 shows the soil response observed during the dry period (1 June-31 July 2019). In this case, it is easy to note that the suction starts to increase earlier and faster in the shallowest layers, and then in the deepest ones. This response can be related to the upward evapotranspiration fluxes due to the water uptake from the roots, which are denser in the upper part of the soil cover due to the flourishing underbrush. Although the soil moisture dynamics are strictly connected to the specific precipitation trends of each hydrologic year, it is possible to define different periods with a similar profile in terms of soil water content and suction (Table 6). The effects of infiltration in response to rainfall generally depend on depth, as fluctuations of soil water content and suction are delayed and smoothened (less marked) in the deepest layers, while in the shallowest ones, they tend to be rapid and higher (values of suction decrease quickly and immediately after the observed rainfall). Figure 14 shows the soil response observed during the dry period (1 June-31 July 2019). In this case, it is easy to note that the suction starts to increase earlier and faster in the shallowest layers, and then in the deepest ones. This response can be related to the upward evapotranspiration fluxes due to the water uptake from the roots, which are denser in the upper part of the soil cover due to the flourishing underbrush. Although the soil moisture dynamics are strictly connected to the specific precipitation trends of each hydrologic year, it is possible to define different periods with a similar profile in terms of soil water content and suction ( Table 6).  * As observed by rainfall data in Figure 6, the dry seasons start at the end of May and finish in September. ** Typical profiles observed in this study (see Figure 15): 1 moisture strongly increasing upward, 2 moisture evenly high, 3 moisture decreasing upward, 4 moisture evenly low.
September. ** Typical profiles observed in this study (see Figure 15): 1 moisture strongly increasing upward, 2 moisture evenly high, 3 moisture decreasing upward, 4 moisture evenly low. Figure 15 shows soil water content and suction profiles throughout the soil cover, estimated from the local values at the depths where TDR probes and tensiometers were located, in late autumn, middle winter, late spring and summer from 2017 to 2020.   Figure 15 shows soil water content and suction profiles throughout the soil cover, estimated from the local values at the depths where TDR probes and tensiometers were located, in late autumn, middle winter, late spring and summer from 2017 to 2020.
In late autumn (Figure 15a), the slopes are being rewetted in the upper part (layer A) after the typical dry and hot Mediterranean summer, as the wetting front in response to rainfall is not still able to propagate to further depths (only the shallowest layer has a quick response to precipitation). In these cases, the water content and suction in superficial levels of soil (layer A) reach values between 0.42-0.50 and 5-8 kPa, respectively. On the other hand, during the middle of winter (Figure 15b), the prolonged rainy periods allow the rainwater infiltration to reach the deepest layer (layer E), leading to a relatively uniform wetting of the entire thickness of the pyroclastic cover, with water content ranging between 0.38 and 0.55 and the correspondent suction ranging between 2-3 kPa in the uppermost part and about 15 kPa at greater depth. The base of the soil profile becomes wet, indicating that in this period gravitational leakage from the soil profile may occur. In the late spring (Figure 15c), the soil profile starts being progressively dried, mostly in the uppermost meter, due to the water uptake from the growing vegetation. In fact, the water content ranges between 0.35 and 0.45 (with an exceptional case of 0.5 during 2018) and suction between 8 kPa and 22 kPa. At the base of the profile (layer E), the soil is still wet given the slower downward-directed gravitational drainage towards the underlying fractured bedrock. At the end of summer (Figure 15d), after a prolonged dry period (with only a few isolated summer rainstorms), the soil profile becomes evenly dry. The values of suction reach 60 kPa, while the water content drops below 0.3 (except during 2008).

Discussion
The measured precipitation and soil water content, and the estimated evapotranspiration and stream discharges, have been used to assess the water balance of the hydrological system along three hydrologic years, between 1 September 2017 and 31 August 2020 ( Figure 16): At the right hand side of Equation (6), written in terms of specific storage [L], V s and V b represent the water volume stored in a column of unit area of soil and bedrock (epikarst and groundwater system), respectively. At the left hand side, R is rainfall intensity [L/T], ET is the actual evapotranspiration [L/T], and q s is the specific stream discharge [L/T]. The various terms of Equation (6), cumulated over each monitored hydrological year, are plotted in Figure 16. The measured precipitations were used to assess the rainfall intensity R. During the periods where the rain gauge was out of order (Section 3.1), the data were supplemented with those collected by the nearby rain gauge of Cervinara, managed by the regional Civil Protection Agency, located at the foot of the slope. The empirical formula of Blaney-Criddle was used to estimate the evapotranspiration ET (Section 3.2). The total discharge (sum'of estimated discharges of the two monitored streams, Section 3.3) was divided by the underground watershed extension, estimated around 2.5 km 2 from the geologic map of the area, to obtain the specific discharge q s . The measurements of soil water content, obtained by the TDR probes installed at depths of 0.30 m, 0.95 m, 1.40 m and 2.00 m (Section 3.4), allow the estimation of soil specific storage (V s = θ i ∆z i ) in the whole soil cover, for which the total thickness was assumed to be equal to 2.60 m. Finally, the term V b was obtained from the other measured/estimated terms with the balance Equation (6). To help understand the annual dynamics of the water balance, in Figure 16 both the specific storage terms V s and V b have been set to zero at the beginning of each hydrologic year. From the yearly trend of the various terms of the water balance, it is clearly shown how the monitored streams are almost dry from the end of summer until nearly the end of the year. Afterwards, the shallow groundwater circulation supplying the streams becomes active, with maximum discharge in early April (17)(18), in February (18)(19), in December (19)(20), and returning to dry in September (2018) and in August (2019 and 2020). The slow response of the stream to precipitation indicates that part of the water is temporarily stored within the epikarst, before coming out into the streams. From the plots of Figure 16, it clearly seems that stream discharge appears when water starts to be stored in the fractured bedrock. As the soil cover remains far from saturation during the entire rainy season (from October to April/May), a condition not favouring fast lateral subsurface flow through the soil cover, the stream water likely originates from the upper part of the fractured bedrock, where a temporary perched aquifer likely develops. It is also worth noting the different response of the soil cover to precipitation in different periods of the rainy season. In fact, soil storage strongly increases only in late autumn and early winter, being nearly equal to the total precipitation (e.g., December 2017 and November 2019) and indicating that little water reaches the fracture bedrock in this period of the year, probably limited by the still dry lower part of the soil cover. Afterwards, although significant precipitation normally occurs during the whole of winter and early spring, soil storage exhibits small fluctuations around values still far from saturation (e.g., March 2018, February-May 2019 and April 2020). As evapotranspiration is limited in this period of the year, it can be argued that most infiltrating water crosses the soil cover and reaches the underlying fractured bedrock. Hence, the transfer of water from the soil cover to the fractured bedrock appears an important process to prevent the achievement of high moisture conditions within the soil cover. This result can help understand why many landslides have been recorded in late autumn and early winter in Campania [48,49].

Conclusions
To investigate the hydrologic behaviour of an unsaturated pyroclastic slope cover laying upon a fractured limestone bedrock in response to meteorological forces, an automatic meteo-hydrological monitoring station was installed along the slope of Cervinara, Northeast of Naples, where catastrophic rainfall-induced landslides occurred in December 1999. The acquired set of data, consisting of rainfall input, meteorological variables, stream discharge, and water contents of the soil cover, enabled the assessment of the water balance of the slope for three hydrological years (2017-2018, 2018-2019, 2019-2020) and, in particular, the terms of water storage in the soil cover and in the fractured bedrock, and the exchanges of water with the atmosphere, through actual evapotranspiration, and with the epikarst and the deeper groundwater system through leakage at the bottom of the soil cover. The comparison between the times of generation of the stream discharge and of the bedrock storage suggests that a perched aquifer temporarily develops in the epikarst once the soil cover has become wet. In fact, the leakage from the soil cover to the perched aquifer seems to be related to the moisture conditions at the base of the soil profile, which exhibit recurrent seasonal behaviour. Specifically, varying soil response in terms of water content and suction during wet and dry seasons is observed, and four different periods with recurrent similar soil profiles are defined, i.e., late autumn, middle winter, late spring and summer end.
The observed recurrent soil moisture conditions in early winter can be of use to understand why many landslides occur in this period in Campania. In fact, the still dry conditions of the lowermost part of the soil cover, as well as the still poor accumulation of water in the epikarst, hamper the drainage through the soil-bedrock interface, and thus favouring the fast accumulation of water within the cover during rainfall events, which can eventually lead to slope failure.
The obtained results, while confirming that, in the studied geomorphological context, soil moisture monitoring can be of great help for the improvement of landslide prediction, also indicate that other variables, such as stream water level, can be a proxy of the hydrological conditions, which can predispose the slopes to landslides. This can be of help for the development of effective landslide early warning systems, simple enough to be operated by public authorities in charge of civil protection actions.
Funding: This research was funded by Università degli Studi della Campania "L. Vanvitelli" as part of the Ph.D. project "Modelling hydrological processes affecting rainfall-induced landslides for the development of early warning systems" within the Doctoral Course "A.D.I.". Experimental activities were also funded by Università degli Studi della Campania "L. Vanvitelli" through the programme "VALERE: VAnviteLli pEr la RicErca".