Processes in the Unsaturated Zone by Reliable Soil Water Content Estimation: Indications for Soil Water Management from a Sandy Soil Experimental Field in Central Italy

: Reliable soil moisture data are essential for achieving sustainable water management. In this framework, the performance of devices to estimate the volumetric moisture content by means dielectric properties of soil/water system is of increasing interest. The present work evaluates the performance of the PR2/6 soil moisture proﬁle probe with implications on the understanding of processes involving the unsaturated zone. The calibration at the laboratory scale and the validation in an experimental ﬁeld in Central Italy highlight that although the shape of the moisture proﬁle is the same, there are essential differences between soil moisture values obtained by the calibrated equation and those obtained by the manufacturer one. These differences are up to 10 percentage points for ﬁne-grained soils containing iron oxides. Inaccurate estimates of soil moisture content do not help with understanding the soil water dynamic, especially after rainy periods. The sum of antecedent soil moisture conditions (the Antecedent Soil moisture Index (ASI)) and rainfall related to different stormﬂow can be used to deﬁne the threshold value above which the runoff signiﬁcantly increases. Without an accurate calibration process, the ASI index is overestimated, thereby affecting the threshold evaluation. Further studies on other types of materials and in different climatic conditions are needed to implement an effective monitoring network useful to manage the soil water and to support the validation of remote sensing data and hydrological soil models.


Introduction
In the effort to achieve most of the United Nation (UN) Sustainable Development Goals (SDGs), it is necessary to move towards a sustainable use and management of the soil-water system [1]. For achieving sustainable water management, the monitoring of soil moisture content is mandatory in the understanding of processes involving the unsaturated zone (e.g., infiltration, migration of pollutants, etc.) [2]. As reported by Vanclooster [3], the monitoring is linked to the modeling of soil processes and space-based environmental monitoring systems can give benefits, such as rapid and effective adaptation of the water sector [4]. In this framework, sustainable soil and water management objectives require reliable data. According to Escorihuela and Quintana-Seguí [5], reliable soil moisture data are essential in the Mediterranean region because they affect the characteristics of land processes such as droughts and floods. The ongoing climate change affecting this area decreased rainfall during the recharge period (autumn to spring period, e.g., [6][7][8]) and intensified extreme rainfall events [9]. As recently reported by Mimeau et al. [10], future scenarios indicate a further change of precipitation patterns in the coming years. These changes are associated with an intensification of extreme rainfall events alternating with long periods of drought, increasing the pressure on land with progressive degradation processes affecting soil productivity, irrigation, groundwater recharge, etc. [11]. The moni-toring of the soil-water system is crucial in the efficiency of water use, while also improving the agricultural sustainability [12].
Soil moisture can be obtained by direct methods (thermo-gravimetric laboratory method, [13]) or indirect methods, most of which are based on the measurement of the dielectric properties of the soil-water system. The two methods are used locally, have advantages and disadvantages, and are complementary to each other.
The gravimetric method is the standard technique to measure the gravimetric water content (θ g ). It is also useful to calibrate the estimates of the indirect methods. It requires the soil sampling (destructive method), is time-consuming, and cannot be used for the monitoring in the field [14,15].
Indirect methods allow real-time and continuous monitoring of the volumetric water content (θ) and are more and more used in the latest years. Several instruments are available with different accuracy, costs (including the maintenance), the difficulty of installation, and response time [16]. Compared to the Time Domain Reflectometry (TDR), the Capacitance and Frequency Domain Reflectometry (FDR) allow, by means of the profile probe (PP), the estimation of spatial-temporal evolution of θ [17][18][19][20][21][22][23]. Moreover, PP can contribute to the implementation of an integrated monitoring system, also aiming to map the spatial variation of the soil water content coupled with other electromagnetic methods, such as the Ground Penetrating Radar, GPR [24,25]. Field investigations also help the calibration and validation of remote sensing data, enhancing the understanding of processes at large scales [26,27]. The quality of locally monitored data is the basis for the re-analysis of the soil water content data from satellites [28].
Although the using of capacitance and FDR profile probes is increasing, the calibration problem and performance of the instrumentations cannot be neglected, especially when they are to be applied in large sites [17,23,[29][30][31][32][33][34][35]. In other words, the use of manufacturer equations to convert the electromagnetic properties of soils into θ may lead to incorrect assessments of the soil moisture content, with inevitable effects on soil water management and in the understanding of processes occurring in the unsaturated zone [23,36,37]. In this framework, the present work aims to contribute in evaluating the performance of the PR2/6 soil moisture profile probe with implications on devices operating at a frequency of 100 MHz. An experimental field in Central Italy is selected and the monitored data are presented and discussed, taking into account the equation used to estimate the soil water content, the antecedent rainfalls, and the modeled parameters of the unsaturated soil.

Experimental Field and Soil Characteristics
For the study, an experimental plot 18 m 2 wide was selected. The site is located in Central Italy in the left bank of the Tiber River alluvial plain (43 • 07 33.60" N-12 • 26 04.76" E, Figure 1a), where a fine-grained soil widely outcrops (soil S B ). According to the USDA (United States Department of Agriculture) classification, this soil is a loamy-sand deriving from the erosion of marly-calcareous-siliceous rocks (Flysch) widely outcropping in the Tiber River catchment (the mineralogical composition and organic matter content is shown in Table 1). The lithological model of the site has been reconstructed with the GPR survey [25], with the Dynamic Probe Super Heavy test (DPSH) and with sampling at different depths (0.15, 0.45, 0.60, 0.80, and 0.95 cm) aiming to obtain the main physical parameters of the soil (Figure 1b).
The site is equipped by a weather station (rainfall and temperature) managed by the Servizio Idrografico della Regione Umbria (Figure 1a) and by a PR2/6 profile probe (Delta-T Devices Ltd., Cambridge, UK), which constitutes an integrated monitoring system. The area is characterized by a mean precipitation of about 860 mm/year and mean annual temperature of 14.5 • C (1992-2019 period).
As shown in Figure 1c, the density of sands tends to increase with depth. Within the investigation range of the infiltration process (in the first-meter depth), there is homogeneous loose sand with dry unit weight (γ d ) of about 14.33 kN/m 3 and void index (e) of 0.82 (data obtained on the soil samples collected from S 6 and S 5 ). These are young alluvial sediments related to the two important floods of the Tiber River occurred in November 2005 and in November-December 2010. The geotechnical characteristics of soil S B , together with some other non-plastic soils (sands and loamy sands of different nature), are shown in Table 1. As presented in the introduction section, to check the performance of measurements by indirect methods, the comparison of empirical equations to estimate θ values on similar materials can be useful. The results of this analysis will be presented in Section 3.1.  The PR2/6 capacitance profile probe allows the estimation of the soil volumetric water content (θ) at different depths. The device requires the installation of an access tube of 1.10 m length into the soil. A 9-volt battery generates a signal of 100 MHz (similar to FM radio), which is applied to six pairs of stainless steel rings placed at different depths (0.10, 0.20, 0.30, 0.40, 0.60, and 1.00 m). An electromagnetic field extends all around the device ( Figure 2) and a data logger records the voltage output (V), which is converted into the square root of the dielectric constant ( √ ε ) by a six-order polynomial equation (Equation (1), [40]). The volumetric water content is then obtained by a soil calibration equation (Equation (2)), which contains two parameters, a 0 (soil offset) and a 1 (slope), as proposed by [41]. The dielectric constant (ε) obtained by the PR2/6 probe is the sum of soil real (ε') and imaginary (ε", dielectric loss) permittivity. It is known that for devices working at low frequencies (<1 GHz), the ε" can be affected by the water salinity and therefore influencing the estimates of θ. This problem cannot be neglected for effluent waters with electrical conductivity values higher than 1500-5000 µS/cm [42,43]. According to the PR2/6 manual [40], the probe has a very low intrinsic sensitivity to the salinity of effluent water when it is of an order of magnitude lower than that reported by the above reference studies. In these cases, the accuracy of measurements mainly depends on the proper choice of parameters a 0 and a 1 of Equation (2): in general, the manufacturer suggests to set a 0 = 1.6 and a 1 = 8.4 for mineral soils. The same parameters are suggested for other commercially available electromagnetic devices such as SM200, SM300, and Theta Probe ML2x. All these probes are produced by the Delta-T Devices, work at the same frequency of PR2/6 (100 MHz), and θ values are estimated from voltage outputs of sensors [44].
In this work, the column test has been used for the calibration of parameters a 0 and a 1 as already proposed by other studies [23,37]. Figure 2 summarizes the step-by-step procedure used for calibrating the PR2/6 probe. By using Equation (2), actual θ values are computed and related to √ ε values for obtaining the best a 0 and a 1 parameters for the investigated soil.
where  Step-by-step procedure for the calibration of PR2/6 probe used in this research. The compaction procedure used to compact soil in the column refers to [23].

Double Ring Infiltrometer
The Double Ring Infiltrometer (DRI) is a standardized technique [45] used to estimate the infiltration rate of water in fine-grained soils. The instrument is composed of two steel rings of different diameters: in this study, the inner ring is 30 cm wide while the outer ring is 60 cm. By using a Mariotte bottle, the volume of water needed to keep the water level in the inner rings constant (constant head) is measured over time. The water passing from the inner ring to the soil during a time-step (Δt) represents the infiltration rate (i, Equation (4)). The steady-state infiltration can be assumed as the saturated hydraulic conductivity of the soil (ks), even if the value obtained may be influenced by entrapped air in the voids, biological activity, size of rings, etc. [46]. where: i = infiltration rate (cm/s); V = volume of added water (cm 3 ) in each time-step (Δt); Δt = time step (s); r = radius of the inner ring (cm).

The Soil Water Characteristic Curve (SWCC)
Most of the information for the understanding of infiltration processes within the unsaturated soils is contained in the Soil Water Characteristic Curve (SWCC). The SWCC relates the matric suction (the difference between pore air pressure, ua, and pore water Step-by-step procedure for the calibration of PR2/6 probe used in this research. The compaction procedure used to compact soil in the column refers to [23].

Double Ring Infiltrometer
The Double Ring Infiltrometer (DRI) is a standardized technique [45] used to estimate the infiltration rate of water in fine-grained soils. The instrument is composed of two steel rings of different diameters: in this study, the inner ring is 30 cm wide while the outer ring is 60 cm. By using a Mariotte bottle, the volume of water needed to keep the water level in the inner rings constant (constant head) is measured over time. The water passing from the inner ring to the soil during a time-step (∆t) represents the infiltration rate (i, Equation (4)). The steady-state infiltration can be assumed as the saturated hydraulic conductivity of the soil (k s ), even if the value obtained may be influenced by entrapped air in the voids, biological activity, size of rings, etc. [46]. where: i = infiltration rate (cm/s); V = volume of added water (cm 3 ) in each time-step (∆t); ∆t = time step (s); r = radius of the inner ring (cm).

The Soil Water Characteristic Curve (SWCC)
Most of the information for the understanding of infiltration processes within the unsaturated soils is contained in the Soil Water Characteristic Curve (SWCC). The SWCC relates the matric suction (the difference between pore air pressure, u a , and pore water pressure, u w ) and θ [47]. As reported by Lu and Likos [48], direct measurements of the matric suction are time-consuming, and often only a few data points of the SWCC are obtained [49]. In this work to reconstruct the SWCC of the soil S B outcropping in the experimental field, "The Soil Water Characteristics software (version 6.02.75)" is used (download page: https://hrsl.ba.ars.usda.gov/soilwater/Index.htm). This software, which has been developed by the USDA [50], estimates several hydrological soil properties based on the texture class, gravel fraction, organic matter, salinity, and compaction (dry unit weight). In detail, the program-by using the model developed by [51]-allows the estimation of the matric suction and some other essential parameters (residual volumetric water content, θ r , and saturated volumetric water content, θ s ) that are useful to estimate the suction stress (σ s ). As reported by Lu and Likos [48], σ s can be expressed in terms of normalized volumetric water content (or effective degree of saturation, S e ) and matric suction (Equation (5)): where: σ s = suction stress (kN/m 2 ) θ = volumetric water content (m 3 /m 3 ) θ r = residual volumetric water content (m 3 /m 3 ) θ s = saturated volumetric water content (m 3 /m 3 ) u a = pore air pressure (kN/m 2 ) u w = pore water pressure (kN/m 2 ) (u a − u w ) = matric suction (kN/m 2 ) S e = effective degree of saturation (m 3 /m 3 ).

Calibration of PR2/6 Probe on Soil S B and Comparison with Other Sands
The step-by-step procedure presented in Figure 2 allowed the obtaining of reliable calibration curves for accurate estimates of θ values by the PR2/6 profile probe, an essential step before using the device in the field. Figure 3 shows the calibration curves with parameters a 0 and a 1 of Equation (2) for the different sands listed in Table 1, also referring to previous published data [23,30]. Figure 3 also shows the data of a new sand (soil S C ) collected in a quarry located in the Paglia River alluvial plain, Central Italy (42 • 44 14.92" N-12 • 06 38.13" E, [52]). It is interesting to point out that the equations differ from that suggested by the manufacturer, which tends to overestimate the θ values.  Table 1). θ values are obtained on soil samples by the gravimetric method using Equation (3), following the procedure in Figure 2. (1) data obtained with the PR2/6 probe (taken from [23]); (2) unpublished data obtained with the PR2/6 probe (soil from Paglia River, Central Italy); (3) data obtained with the Theta probe (taken from [30]). PR2/6 and Theta probes work at the same frequency (100 MHz), and θ values are estimated from voltage outputs of sensors.

Soil Water Characteristic Curve of Soil SB and Infiltration Test
The Soil Water Characteristics software has been used to model hydrological soil properties of soil SB outcropping in the experimental field. The input parameters which has been set into the program are: • sand and clay fraction of 82% and 3%, respectively (Table 1); • organic matter content of 1.53% (Table 1); • electrical conductivity of the effluent fluids, 440 μS/cm [23]; • dry unit weight (14.33 kN/m 3 ). Figure 4a shows the SWCC and the hydraulic conductivity curve of soil SB. The modeled residual volumetric water content, θr, and saturated volumetric water content, θs, were 3.0% and 44.9%, respectively. The estimated θs agrees with that calculated by applying the unit-weight relationship (Equation (6)), which provides a value of 45.5%. In the computation, the following data have been considered: γd = 14.33 kN/m 3 ; e = 0.82; specific gravity, Gs = 2.63; γw = 9.81 kN/m 3 .
The modeled saturated hydraulic conductivity (ks) resulted in 3.4 × 10 −3 cm/s ( Figure  4a). Figure 4b shows the result of the DRI test carried out in the field at the end of July 2018. The initial θ value in the first 30 cm depth was about 5%. The infiltration curve is described by the Horton equation [53] and approaches to quasi-steady-state conditions after about 150 min, reaching a value of about 8.0 × 10 −3 cm/s. Although the enclosure of air bubbles in the voids during infiltration prevents maximum saturation, the ks value  Table 1). θ values are obtained on soil samples by the gravimetric method using Equation (3), following the procedure in Figure 2. (1) data obtained with the PR2/6 probe (taken from [23]); (2) unpublished data obtained with the PR2/6 probe (soil from Paglia River, Central Italy); (3) data obtained with the Theta probe (taken from [30]). PR2/6 and Theta probes work at the same frequency (100 MHz), and θ values are estimated from voltage outputs of sensors.

Soil Water Characteristic Curve of Soil S B and Infiltration Test
The Soil Water Characteristics software has been used to model hydrological soil properties of soil S B outcropping in the experimental field. The input parameters which has been set into the program are: • sand and clay fraction of 82% and 3%, respectively (  Figure 4a shows the SWCC and the hydraulic conductivity curve of soil S B . The modeled residual volumetric water content, θ r , and saturated volumetric water content, θ s , were 3.0% and 44.9%, respectively. The estimated θ s agrees with that calculated by applying the unit-weight relationship (Equation (6)), which provides a value of 45.5%. In the computation, the following data have been considered: γ d = 14.33 kN/m 3 ; e = 0.82; specific gravity, G s = 2.63; γ w = 9.81 kN/m 3 .
The modeled saturated hydraulic conductivity (k s ) resulted in 3.4 × 10 −3 cm/s (Figure 4a). Figure 4b shows the result of the DRI test carried out in the field at the end of July 2018. The initial θ value in the first 30 cm depth was about 5%. The infiltration curve is described by the Horton equation [53] and approaches to quasi-steady-state conditions after about 150 min, reaching a value of about 8.0 × 10 −3 cm/s. Although the enclosure of air bubbles in the voids during infiltration prevents maximum saturation, the k s value obtained by the infiltrometer is of the same order of magnitude of the k s predicted by the SWC software.

Precipitation Analyses and Soil Water Profiles
The monitoring of the spatial-temporal variations of soil water content with depth is important in different practical applications (hydrological, hydrogeological, and pedological). To check the response of the PR2/6 probe, the soil water profile dynamic has been analyzed during April-May 2018. Figure 5a shows the daily rainfall and temperature data recorded in the period by the weather station located less than 70 m from the experimental site (Figure 1a). In detail, March 2018 experienced very high precipitations, which were found to be more than three times higher than the historical mean ( Figure  5b). In the subsequent months, rainfall data were in line with the historical mean, with a not negligible increase also found in May. During the observation period (late winter-spring months) the air temperature gradually increased, except for the two rainy periods of March and May 2018. In this framework, the temporal variations of soil water content with depth allow the understanding of processes during wetting and drying periods. Figure 5a also shows the location of four selected comprehensive water content measurements taken before and after important meteorological events. In order to have accurate water content data, each measure has been repeated at least seven times by setting the PR2/6 datalogger with both calibrated parameters and manufacturer ones (mineral soil).

Precipitation Analyses and Soil Water Profiles
The monitoring of the spatial-temporal variations of soil water content with depth is important in different practical applications (hydrological, hydrogeological, and pedological). To check the response of the PR2/6 probe, the soil water profile dynamic has been analyzed during April-May 2018. Figure 5a shows the daily rainfall and temperature data recorded in the period by the weather station located less than 70 m from the experimental site ( Figure 1a). In detail, March 2018 experienced very high precipitations, which were found to be more than three times higher than the historical mean (Figure 5b). In the subsequent months, rainfall data were in line with the historical mean, with a not negligible increase also found in May. During the observation period (late winter-spring months) the air temperature gradually increased, except for the two rainy periods of March and May 2018. In this framework, the temporal variations of soil water content with depth allow the understanding of processes during wetting and drying periods. Figure 5a also shows the location of four selected comprehensive water content measurements taken before and after important meteorological events. In order to have accurate water content data, each measure has been repeated at least seven times by setting the PR2/6 datalogger with both calibrated parameters and manufacturer ones (mineral soil).  Figure 6 shows the selected soil moisture profiles. Following we present the results obtained by setting the HHD datalogger of the PR2/6 probe with the calibrated parameters, moving the discussion about the differences obtained with the manufacturer's parameters for mineral soils to the next section.
The first measurement refers to 26 April 2018 (supplementary), after about two weeks from the end of a very high rainfall period (Figure 6a). Due to the effect of the antecedent rainfalls (about 250 mm in two months), the moisture profile shows the typical shape of a wetting front moving downwards e.g., [54]. The low values in the shallow zone (the first 30 cm depth) are also due to the effect of the evapotranspiration (the air temperature increased of about 10 °C in the week before the measurement). The succession of rainy and dry periods in the following weeks introduced a lot of disturbing effects on the soil water profile. The precipitations occurred at the beginning and middle of May increased the volumetric water content in the first 30 cm depth of about 12 percentage points (Figure 6b,c). The water distribution in the first-meter depth has gradually assumed a constant profile with an effective degree of saturation (Se) of about 36%. As shown in Figure 6d, after few days of no rainfall-coupled with a rise of evapotranspiration phenomena linked to the increase in temperature-the water profile tended to re-assume a shape similar to that recorded on 12 May 2018. Figure 6a,d also display the matric suction profiles during the same observation period, as obtained from the SWC software. In each measurement, the matric suction decreases with depth with maximum changes in the first 30 cm as affected by meteo-climatic conditions. For instance, between 16 April and 12 May, a mean decrease of 150 kPa is observed in the shallow zone.  Figure 6 shows the selected soil moisture profiles. Following we present the results obtained by setting the HHD datalogger of the PR2/6 probe with the calibrated parameters, moving the discussion about the differences obtained with the manufacturer's parameters for mineral soils to the next section.
The first measurement refers to 26 April 2018 (Supplementary Materials), after about two weeks from the end of a very high rainfall period (Figure 6a). Due to the effect of the antecedent rainfalls (about 250 mm in two months), the moisture profile shows the typical shape of a wetting front moving downwards e.g., [54]. The low values in the shallow zone (the first 30 cm depth) are also due to the effect of the evapotranspiration (the air temperature increased of about 10 • C in the week before the measurement). The succession of rainy and dry periods in the following weeks introduced a lot of disturbing effects on the soil water profile. The precipitations occurred at the beginning and middle of May increased the volumetric water content in the first 30 cm depth of about 12 percentage points (Figure 6b,c). The water distribution in the first-meter depth has gradually assumed a constant profile with an effective degree of saturation (S e ) of about 36% (Figure 6g). As shown in Figure 6d, after few days of no rainfall-coupled with a rise of evapotranspiration phenomena linked to the increase in temperature-the water profile tended to re-assume a shape similar to that recorded on 12 May 2018. Figure 6i-l also display the matric suction profiles during the same observation period, as obtained from the SWC software. In each measurement, the matric suction decreases with depth with maximum changes in the first 30 cm as affected by meteo-climatic conditions. For instance, between 26 April and 12 May, a mean decrease of 150 kPa is observed in the shallow zone.

Discussions
The performed analysis confirms that there is no unique relationship between dielectric properties of the medium (√ε) and θ, even when the same instrument is used on similar materials (e.g., non-plastic soils). Most of the sands analyzed here tend to align on the calibration curve having a0 = 1.7 and a1 = 9.5, which differs from that of the manufacturer (mineral soils, a0 = 1.6 and a1 = 8.4). For this sand types (mainly composed by quartz and/or calcium carbonate, Table 1) the use of the manufacture parameters leads to a not excessive overestimation of θ values from 2 to 5 percentage points. In contrast, SB sands-coming from the experimental field-shows the highest parameters (a0 = 1.9 and a1 = 10.6), which in turn give θ values from 5 to 10 percentage points lower than those estimated by using the manufacturer's parameters. The probe calibration is always necessary, but it is particularly important for capacitance probes working at low frequency, especially when used in a medium having high iron content. The iron content affects the imaginary component (ε") of the dielectric constant of the soil/water system [55]. For these soils, the difference between actual and estimated values by the manufacturer parameters tends to increase as the degree of saturation of soil increases. Soil SB has an iron oxides content of about 6.5% [56] with minor contents of pyrite and chlorite. As reported

Discussion
The performed analysis confirms that there is no unique relationship between dielectric properties of the medium ( √ ε ) and θ, even when the same instrument is used on similar materials (e.g., non-plastic soils). Most of the sands analyzed here tend to align on the calibration curve having a 0 = 1.7 and a 1 = 9.5, which differs from that of the manufacturer (mineral soils, a 0 = 1.6 and a 1 = 8.4). For this sand types (mainly composed by quartz and/or calcium carbonate, Table 1) the use of the manufacture parameters leads to a not excessive overestimation of θ values from 2 to 5 percentage points. In contrast, S B sands-coming from the experimental field-shows the highest parameters (a 0 = 1.9 and a 1 = 10.6), which in turn give θ values from 5 to 10 percentage points lower than those estimated by using the manufacturer's parameters. The probe calibration is always necessary, but it is particularly important for capacitance probes working at low frequency, especially when used in a medium having high iron content. The iron content affects the imaginary component (ε") of the dielectric constant of the soil/water system [55]. For these soils, the difference between actual and estimated values by the manufacturer parameters tends to increase as the degree of saturation of soil increases. Soil S B has an iron oxides content of about 6.5% [56] with minor contents of pyrite and chlorite. As reported by Robinson et al. [57], the iron oxide content, coupled with the length of sensor rods used, play an important role in capacitance probe working at 100 MHz, such as the PR2/6 probe. However, recent research has shown that the effects of iron content are more evident for instruments working at even lower frequencies (cf., a WET sensor working at 50 MHz, [55]). The effects of an accurate calibration are also evident when analyzing the temporal variations of soil water content with depth in the experimental field. The use of the PR2/6 probe without an accurate calibration does not help the understanding of what is really happening in the soil. As shown in Figure 6a-d, although the shape of the moisture profiles is the same, there are essential differences between θ values obtained by the calibrated equation and those by the manufacturer one. Referring to Figure 6c, the difference between the two soil moisture profiles is exacerbated after rainy periods. Following to 40 mm of rain in three days (from 14 May to 16 May), the hydrogeological processes that occurred in the shallow soil (first 30 cm depth) increased the θ values from about 18% (calibrated parameters) to 27% (manufacturer parameters). In other words, the fraction of water actually occupied by free water into voids is much higher if the manufacturer equation is used. As presented by Haga et al. [58], the subsurface water dynamics also depends on the antecedent soil moisture conditions, which can be evaluated by using the Antecedent Soil moisture Index (ASI). The sum of ASI and precipitation (P) related to different stormflows allows the identification of threshold values (ASI + P, in mm) above which runoff significantly increases [59,60]. In general, the ASI index is estimated at catchment scale or on plot scale by using remote sensing data or indirect methods, such as TDR or capacitance probes: as illustrated, the calibration process is very important in the latter. An overestimation of the soil moisture affects the ASI index computation and therefore, the definition of reliable runoff thresholds. For example, referring to the middle May rainfall event, the ASI index calculated in the first 30 cm depth is 36 or 51 mm depending on whether the calibrated or factory parameters are used. As a result of rainfalls and processes occurring in the shallow soil (runoff, evapotranspiration, etc.), the water gradually infiltrates, producing the decrease of matric suction and the increase of the effective degree of saturation (S e ). In detail, S e in the first 30 cm depth is overestimated of about 20 percentage points by using manufacturer parameters (Figure 6g). Although the information was obtained on a small plot, the results are important for soil water management. As S e increases, the soil tends to absorb much less water and thus increases the runoff, with practical implications in different disciplines, including soil science, geomorphology, hydrogeology, etc. [61]. The comprehension of landscape processes requires quantifying water and sediment fluxes, i.e., through a connectivitybased approach [61]. As illustrated, the antecedent soil moisture content regulates the partitioning of rainfall into infiltration and runoff, even though other components have to be carefully evaluated. As discussed by Zhao et al. [62], at a local scale, micro-relief acts similar to the vegetation cover, reducing the surface wash velocity and increasing the infiltration rate. These aspects are not considered in this paper since it focused on the effect of soil moisture reliability's estimates. To improve the soil-water management, the concepts discussed in this study should be included. To achieve the SDGs, soil hydraulic properties and soil monitoring need to be assessed, indicating that better measurements and modeling approaches are crucial. According to [63,64], modern technological solutions can provide the needed data and information, but the harmonization of monitoring programs and measurement standards and protocols still needs to be completed. The present work contributes to these topics: good data and sustainable soil management practices contribute in mitigating the land degradation, moving towards a Land Degradation-Neutrality target (LDN), an important goal to increase sustainable land use [65,66].

Conclusions
The achieving of sustainable water management requires reliable soil moisture data, which are increasingly obtained by using electromagnetic methods. The laboratory and field-scale investigations provided useful recommendations on the use of the PR2/6 profile probe, and in general capacitance probes working at 100 MHz, to obtain reliable data for the understanding of some processes in the unsaturated zone made by fine-grained soils. The calibration parameters for the different soils considered are useful in various fields of application where the infiltration play a key role (groundwater recharge, erosion and runoff, irrigation, etc.). Particularly interesting is the results for the flyschoid sand (soil S B ), which widely outcrops along the Central-Northern Apennines, confirming that the iron content deeply affects the soil moisture estimations. This aspect, which is often not carefully taken into account, inevitably leads to erroneous assessments, impacting the soil water management. After the calibration, the probe, coupled with the modeling of non-saturated parameters, can accurately describe the dynamic of water in soils. In this way, further studies on other types of materials and in different pedoclimatic conditions are needed to implement an effective monitoring network that is useful for managing the soil water and for supporting the validation of remote sensing data and hydrological soil models. The work contributes to the achievement of LDN goals, as electromagnetic tools are becoming increasingly popular in modern agricultural practices of sustainable soil management.