Preliminary Validation of a Novel Method for the Assessment of Effective Stress State in Partially Saturated Soils by Cone Penetration Tests

A proper assessment of the soil effective stress state is crucial in many cases to identify a potential geological/geotechnical hazard as shallow landslides or failure of levees that may have a significant impact on human activities and development. This paper is aimed at validating a methodology for the expeditious and economic determination of effective stress state in the vadose zone recently proposed by Lo Presti et al. in 2016. The method is based on the interpretation of cone penetration tests (CPTu). Its validation was carried out by comparing the CPTu predicted values of suction against the measured ones in a well-documented and monitored site. The comparison also includes the prediction of suction that was obtained by using the so-called Modified Kovacs model (MK). Moreover, additional data of water content and saturation degree from another site were used to predict the suction by using the MK model. These values of suction were compared with those inferred by the CPTu.


Introduction
Traditional soil mechanics mainly considers fully saturated or dry soils.The interest in partially saturated soils or soils saturated by capillarity is quite recent and still under development.Indeed, notwithstanding the early formulation by Bishop [1] of the effective stress principle for partially saturated soils, in practice, the effects of partial saturation are mainly neglected or rarely applied.
The difficulty in assessing the in situ saturation degree and/or to obtain reliable measurements of suction, as well as the continuous variation of these parameters is the main reason of the limited use in practice of the huge knowledge that has been gained in recent years on the "soil mechanics of partially saturated soils".
Anyway, a correct estimate of suction or in situ effective stress state is extremely important in many cases.As an example, suction can improve significantly the response of a laterally loaded pile foundation, which is mainly affected by the shallower soil layers.To capture the actual load-deflection curves of some lateral load tests on single pile, Stacul et al. [2] developed an analysis method that considers the increase of both stiffness and resistance in partially saturated soil conditions.Other examples are represented by shallow landslides and failures of levee systems: before the triggering event, soils of interest are usually unsaturated.The quick increase of saturation degree and the development of positive values of pore pressure are the most important causes for shallow landslide and triggering of levee-failure [3,4].
Lo Presti et al. [5] proposed a method for assessing the effective stress state for clayey soils by an appropriate interpretation of CPTu tests.Such a method is briefly described in Section 2.
The method is based on the consideration that the tip resistance strongly depends on the water content in the case of fine-grained soils.Therefore, as known, the tip resistance may be very high in the desiccated crust and even in the vadose zone saturated by capillarity.
The suction effect, in fact, has a significant influence on the cone penetration resistance through suction hardening and its contribution to the effective stress state, as observed recently in the work of Yang and Russell [6].In [6] Lyell silty sand was tested in a calibration chamber, and Yang and Russell highlighted that further studies are necessary for a wide range of unsaturated soil types to generalise CPT interpretation methods.
The typical trend of the tip resistance within a partially saturated clayey soil is shown in Figure 1.This figure shows the tip resistance for two tests carried out at the same location with two different water table depths [5,7].The deposit was homogeneous and on the basis of laboratory testing on undisturbed samples retrieved from the first three meters was mainly classified as CL (low plasticity clay) to CH (high plasticity clay) according to the USCS classification system.The water table was found at a depth of 3.5 m during the wet season (June 2001) and at a depth of 5 m during the dry season (September 2001).Despite the homogeneity of the deposit, it can be observed that the tip resistance (q c ) is influenced to a great extent by the water table depth (due to suction) so that q c increases from 1-2 MPa to 3-4 MPa in the vadose zone above the water table (Figure 1).As can be observed, such an increase is higher during the dry season.
degree and the development of positive values of pore pressure are the most important causes for shallow landslide and triggering of levee-failure [3,4].
Lo Presti et al. [5] proposed a method for assessing the effective stress state for clayey soils by an appropriate interpretation of CPTu tests.Such a method is briefly described in Section 2.
The method is based on the consideration that the tip resistance strongly depends on the water content in the case of fine-grained soils.Therefore, as known, the tip resistance may be very high in the desiccated crust and even in the vadose zone saturated by capillarity.
The suction effect, in fact, has a significant influence on the cone penetration resistance through suction hardening and its contribution to the effective stress state, as observed recently in the work of Yang and Russell [6].In [6] Lyell silty sand was tested in a calibration chamber, and Yang and Russell highlighted that further studies are necessary for a wide range of unsaturated soil types to generalise CPT interpretation methods.
The typical trend of the tip resistance within a partially saturated clayey soil is shown in Figure 1.This figure shows the tip resistance for two tests carried out at the same location with two different water table depths [5,7].The deposit was homogeneous and on the basis of laboratory testing on undisturbed samples retrieved from the first three meters was mainly classified as CL (low plasticity clay) to CH (high plasticity clay) according to the USCS classification system.The water table was found at a depth of 3.5 m during the wet season (June 2001) and at a depth of 5 m during the dry season (September 2001).Despite the homogeneity of the deposit, it can be observed that the tip resistance (qc) is influenced to a great extent by the water table depth (due to suction) so that qc increases from 1-2 MPa to 3-4 MPa in the vadose zone above the water table (Figure 1).As can be observed, such an increase is higher during the dry season.The effect of suction is evident even in terms of soil classification index (I c ) values (Figure 2).I c can be inferred by CPTu tests data as described in Robertson [8], Jefferies and Davies [9] and Robertson and Wride [10].In particular, the I c index of Figure 2, was computed according to the normalized classification system by Robertson and Wride [10].As far as the I c index is concerned, the values decrease from about 3.0 at the water table depth to about 2.0 at a depth of 50 cm.In terms of SBTn (i.e., normalized Soil Behavior Type [10]) classes, silts and sand mixtures become predominant instead of OC stiff clay.
Geosciences 2018, 8, 30 3 of 13 The effect of suction is evident even in terms of soil classification index (Ic) values (Figure 2).Ic can be inferred by CPTu tests data as described in Robertson [8], Jefferies and Davies [9] and Robertson and Wride [10].In particular, the Ic index of Figure 2, was computed according to the normalized classification system by Robertson and Wride [10].As far as the Ic index is concerned, the values decrease from about 3.0 at the water table depth to about 2.0 at a depth of 50 cm.In terms of SBTn (i.e., normalized Soil Behavior Type [10]) classes, silts and sand mixtures become predominant instead of OC stiff clay.This paper is aimed at validating this recently proposed methodology.Its preliminary validation was carried out by comparing the CPTu predicted values of suction against the measured ones in a well-documented and monitored site.

Method
The method applicability requires that the soil is homogeneous, and the water table known.Moreover, the CPTu should be performed down to a depth below the water table.Were these hypotheses respected, the Ic value below the water table is assumed as a target one and the lower Ic values above the water table are corrected by changing the vertical effective stress according to Equations (1)-(4) [8][9][10].This paper is aimed at validating this recently proposed methodology.Its preliminary validation was carried out by comparing the CPTu predicted values of suction against the measured ones in a well-documented and monitored site.

Method
The method applicability requires that the soil is homogeneous, and the water table known.Moreover, the CPTu should be performed down to a depth below the water table.Were these hypotheses respected, the I c value below the water table is assumed as a target one and the lower I c values above the water table are corrected by changing the vertical effective stress according to Equations (1)-(4) [8][9][10].
where Q tn = normalized tip resistance; q t = total tip resistance; σ v0 = total vertical stress of the soil; σ atm = reference pressure, equal to (1) atm (=98 kPa); n = stress exponent expressed as a function of the classification index; σ' v0 = effective vertical stress of the soil; F = normalized friction ratio; fs = sleeve friction.
In case of fine-grained soils, the stress exponent (n) is approximately equal to (1), thus once defined the I c target (I c,t ) value the corresponding effective vertical stress (σ' v0 ) can be easily evaluated according to Equation (5).
Once obtained the effective vertical stress value, the relationship (Equation ( 6)) proposed by Bishop [1] can be used to estimate the matric suction.
where u a = air pressure; u w = water pressure; χ = effective stress parameter that is expressed as a function of the degree of saturation (S r ); and suction ψ is equal to u a − u w .Suction can be evaluated according to Equation ( 6) rewritten in the following way (Equation ( 7)): In order to properly use Equation ( 7) the value of χ should be assessed on the basis of the actual degree of saturation of the soil layers above the water table.The relationship between χ and S r (Figure 3) is available in literature for different soil types (e.g., [11][12][13]).More recently, Konrad and Lebeau [14] related the effective stress parameter (χ) to an effective degree of saturation, which is linked to the capillary component of water retention.Nevertheless, the applicability of the methodology developed in [14] has been validated for non-plastic soils.
In case of fine-grained soils, the stress exponent (n) is approximately equal to (1), thus once defined the Ic target (Ic,t) value the corresponding effective vertical stress (σ'v0) can be easily evaluated according to Equation ( 5).

(
) Once obtained the effective vertical stress value, the relationship (Equation ( 6)) proposed by Bishop [1] can be used to estimate the matric suction.

(
) ( ) where ua = air pressure; uw = water pressure; χ = effective stress parameter that is expressed as a function of the degree of saturation (Sr); and suction ψ is equal to ua − uw.
Suction can be evaluated according to Equation ( 6) rewritten in the following way (Equation ( 7)): In order to properly use Equation ( 7) the value of χ should be assessed on the basis of the actual degree of saturation of the soil layers above the water table.The relationship between χ and Sr (Figure 3) is available in literature for different soil types (e.g., [11][12][13]).More recently, Konrad and Lebeau [14] related the effective stress parameter (χ) to an effective degree of saturation, which is linked to the capillary component of water retention.Nevertheless, the applicability of the methodology developed in [14] has been validated for non-plastic soils.The MK-model (Modified-Kovacs model), as proposed by Aubertin et al. [15,16] (see also Mbonimpa et al. [17,18]), was also used for estimating suction (ψ) on the basis of the measured liquid limit (w L ), void ratio (e) and degree of saturation (S r ) in the shallower soil layers.The MK-model makes use of a parameter defined as the equivalent capillary rise (h c0 ) in the porous medium and is calculated using the expression for the rise of water in a capillary tube with a diameter d.
For the sake of convenience, the expression to estimate the equivalent capillary rise in cohesive/plastic soils is reported below (Equation ( 8), [16]): where ρ s = soil grain density (kg/m 3 ); e = void ratio (−); w L = liquid limit (%).The MK-model uses h c0 as a reference parameter to define the relationship between the degree of saturation (S r ) and the matric suction (ψ).The model considers that water is held by capillary forces, responsible for capillary saturation (S c ), and by adhesive forces, causing saturation by adhesion (S a ).The relationship proposed in the MK-model is written as in Equation ( 9) [16]: The contribution of the capillary and the adhesion components to the total degree of saturation (S r ) is defined as a function of h c0 and ψ using Equations ( 10)-( 12) [16].10) where for cohesive/plastic soils m and a c = model constants that can be assumed equal to 3 × 10 −5 and 7 × 10 −4 , respectively; ψ r = suction at residual water content (Equation ( 12)); ψ 0 = suction value corresponding approximately to complete dryness, in the MK-model is assumed equal to 10 7 cm; ψ n = a normalization parameter and equal to 1 cm when ψ is given in cm, and corresponding to 10 −3 atm.The assessment of suction (ψ) is carried out following these steps: 1. compute the equivalent capillary rise (h c0 ) using ρ s , e and w L measured in situ (Equation ( 8)); 2.
compute S c and S a assuming a trial value for the unknown suction ψ using Equations ( 10) and (11); 3.
repeat steps 2 and 3 assuming different values for the suction ψ in order to obtain a degree of saturation (S r ) approximately equal to that measured in the test site.

Test Sites
Two test sites were considered.Both are located in Oltrepò Pavese, a hilly region of northern Apennines, in northwestern Italy (Figure 4a).The analyzed soils consist of colluvial soils deriving from the weathering of clay bedrock.
The first site (Figure 4b) is named Costa Cavalieri (CC) and it is located on a low-medium topographic gradient slope (7-18 • ).The slope elevation ranges from 460 to 500 m a.s.l.A cultivation of alfalfa grass entirely covers the slope.The bedrock is composed of a clayey melange, with scattered calcareous and marly blocks.Due to the very low permeability of these materials, deep water tables are only seasonally present along the hillslopes inside bedrock levels; perched water table also develops seasonally.
Five main soil horizons were identified in the colluvial cover, thick 2 m (Table 1).These levels were named with subsequent letters, from A to E, according to their position in the soil profile.
The weathered bedrock (WB) is composed of more compacted clays and was identified at 2 m from the ground surface.All the tests for the soil characterization were carried out according to the Recommendation by the Italian Geotechnical Society [19] and the ASTM Standards.
are only seasonally present along the hillslopes inside bedrock levels; perched water table also develops seasonally.
Five main soil horizons were identified in the colluvial cover, thick 2 m (Table 1).These levels were named with subsequent letters, from A to E, according to their position in the soil profile.The weathered bedrock (WB) is composed of more compacted clays and was identified at 2 m from the ground surface.All the tests for the soil characterization were carried out according to the Recommendation by the Italian Geotechnical Society [19] and the ASTM Standards.

Type D (m) Sa (%) S (%) C (%) CaCO3 (%) wL (%) PI (%)
A-cc 0. In the first 1.44 m from the surface, the soil is characterized by a high amount of carbonate in the profile (9.8-13.7%).Moreover, there is the presence of a calcic horizon, rich in carbonate concretions (D-cc level), where the carbonates reach values of 26.7% (Table 1).Instead, E-cc level, below the calcic horizon, does not present carbonate in the soil matrix.An X-ray diffraction analysis on tout-venant revealed that all the soil levels have a high content in mica and in smectite, both generally in the range between 25 and 40%.D-cc level showed an abundance in carbonate minerals (37%), confirming the presence of the calcic horizon.
On the basis of the grain-size distribution (Table 1), soil horizons are very similar.Clay is the most abundant grain size class, ranging in a narrow path between 47.5 and 57.5%.Silt is also in a  In the first 1.44 m from the surface, the soil is characterized by a high amount of carbonate in the profile (9.8-13.7%).Moreover, there is the presence of a calcic horizon, rich in carbonate concretions (D-cc level), where the carbonates reach values of 26.7% (Table 1).Instead, E-cc level, below the calcic horizon, does not present carbonate in the soil matrix.An X-ray diffraction analysis on tout-venant revealed that all the soil levels have a high content in mica and in smectite, both generally in the range between 25 and 40%.D-cc level showed an abundance in carbonate minerals (37%), confirming the presence of the calcic horizon.
On the basis of the grain-size distribution (Table 1), soil horizons are very similar.Clay is the most abundant grain size class, ranging in a narrow path between 47.5 and 57.5%.Silt is also in a narrow range between 39.7 and 46.9%.Sand and gravel contents are very low in all the soil levels, keeping between 0 and 3.2%.According to USCS classification, soil horizons are highly plastic soils (CH), with liquid limit (w L ) ranging from 65.5 to 73.4% and plasticity index (PI) between 45.6 and 53.6% (Table 1).Unit weight (γ) is rather steady along the soil profile, ranging between 18.6 and 19.1 kN/m 3 .
Oedometric tests, performed on undisturbed samples collected in B-cc and C-cc horizons, showed that all the soil levels are over-consolidated.Peak shear strength parameters were obtained for C-cc horizon through a direct-shear test.Peak friction angle (ϕ ) is about 14 • , while effective cohesion (c ) is of 15 kPa.
In this site a hydrological monitoring station (Figure 5) was installed since 27 November 2015 (Bordoni et al. [20]).Data were recorded every 10 min.The monitoring station included water content probes, tensiometers and dielectric sensors.GS3 probes (Decagon Devices Inc., Pullman, WA, USA) were installed at 0.2, 0.4, 0.6, 0.9 and 1.7 m from the ground surface for the measure of the water content.A combination of different devices was used for measuring pore water pressure: two tensiometers (T4e UMS GmbH) and three dielectric sensors (MPS-6 Decagon Devices) were installed at 0.4, 0.6 and 0.9 m from ground level.Tensiometers were installed to measure pore water pressures higher than −101 kPa, while the other pore water pressure devices allow to measure pore water pressures lower than −101 kPa (Bittelli et al. [21]).The monitoring system was equipped with a datalogger (DL-6te, EM-50) powered by batteries.
narrow range between 39.7 and 46.9%.Sand and gravel contents are very low in all the soil levels, keeping between 0 and 3.2%.According to USCS classification, soil horizons are highly plastic soils (CH), with liquid limit (wL) ranging from 65.5 to 73.4% and plasticity index (PI) between 45.6 and 53.6% (Table 1).Unit weight (γ) is rather steady along the soil profile, ranging between 18.6 and 19.1 kN/m 3 .
Oedometric tests, performed on undisturbed samples collected in B-cc and C-cc horizons, showed that all the soil levels are over-consolidated.Peak shear strength parameters were obtained for C-cc horizon through a direct-shear test.Peak friction angle (φ′) is about 14°, while effective cohesion (c′) is of 15 kPa.
In this site a hydrological monitoring station (Figure 5) was installed since 27 November 2015 (Bordoni et al. [20]).Data were recorded every 10 min.The monitoring station included water content probes, tensiometers and dielectric sensors.GS3 probes (Decagon Devices Inc., Pullman, WA, USA) were installed at 0.2, 0.4, 0.6, 0.9 and 1.7 m from the ground surface for the measure of the water content.A combination of different devices was used for measuring pore water pressure: two tensiometers (T4e UMS GmbH) and three dielectric sensors (MPS-6 Decagon Devices) were installed at 0.4, 0.6 and 0.9 m from ground level.Tensiometers were installed to measure pore water pressures higher than −101 kPa, while the other pore water pressure devices allow to measure pore water pressures lower than −101 kPa (Bittelli et al. [21]).The monitoring system was equipped with a datalogger (DL-6te, EM-50) powered by batteries.The second site (Figure 4c) is named Santa Maria della Versa (QUV); it is a western-facing slope also characterized by low-medium topographic gradient (8°).The slope elevation ranges from 260 to 270 m a.s.l.Vineyards entirely cover the slope.
The bedrock consists of a flysch i.e., medium to thickly stratified limestones and marly limestones.The water table is at 3.5 m depth.Soil classification (Table 2) was carried out on samples extracted from two open pits made in the test-site (named QUV1 and QUV2), at depth not greater than about 1.7 m.
Colluvial soils above the bedrock levels were characterized as Costa Cavalieri soils.The soil is rich in clayey minerals, as confirmed by the abundance in clay particles and by a series of shrinkage cracks present during the samplings.
On the basis of the grain-size distribution (Table 2), soil horizons are very similar in depth and between the two sites.Clay is generally the most abundant grain size class, ranging in a path between 46.5 and 66.8%.Silt is also in a range between 25.6 and 50.0%.Sand and gravel contents are very low in all the soil levels, keeping between 0 and 6.0%.The second site (Figure 4c) is named Santa Maria della Versa (QUV); it is a western-facing slope also characterized by low-medium topographic gradient (8 • ).The slope elevation ranges from 260 to 270 m a.s.l.Vineyards entirely cover the slope.
The bedrock consists of a flysch i.e., medium to thickly stratified limestones and marly limestones.The water table is at 3.5 m depth.Soil classification (Table 2) was carried out on samples extracted from two open pits made in the test-site (named QUV1 and QUV2), at depth not greater than about 1.7 m.
Colluvial soils above the bedrock levels were characterized as Costa Cavalieri soils.The soil is rich in clayey minerals, as confirmed by the abundance in clay particles and by a series of shrinkage cracks present during the samplings.
On the basis of the grain-size distribution (Table 2), soil horizons are very similar in depth and between the two sites.Clay is generally the most abundant grain size class, ranging in a path between 46.5 and 66.8%.Silt is also in a range between 25.6 and 50.0%.Sand and gravel contents are very low in all the soil levels, keeping between 0 and 6.0%.
According to USCS classification, soil horizons are highly plastic soils (CH), with liquid limit (w L ) ranging from 71.3 to 84.3% and plasticity index (PI) between 51.6 and 61.2% (Table 2).Unit weight (γ) is rather steady along the soil profile, ranging between 17.9 and 20.1 kN/m 3 .Tables 3 and 4 respectively summarize the saturation degree with depth at CC and QUV, as estimated at the time CPTu were performed.For the test sites "CC" and "QUV" the average degree of saturation (at the time of CPTu tests) in the first 2 m was equal to 55% and 70%, respectively.Considering that the soils were classified as silty clays the best estimate of χ (Figure 3) led to values ranging in between 0.2 and 0.25 (depending on the saturation degree).These values were used in the Bishop equation (Equation ( 7)).

Results
Therefore, it was possible to estimate the values of suction according to the previously described methods and to compare these estimated values to the measured ones (Figure 7).For the test sites "CC" and "QUV" the average degree of saturation (at the time of CPTu tests) in the first 2 m was equal to 55% and 70%, respectively.Considering that the soils were classified as silty clays the best estimate of χ (Figure 3) led to values ranging in between 0.2 and 0.25 (depending on the saturation degree).These values were used in the Bishop equation (Equation ( 7)).
Therefore, it was possible to estimate the values of suction according to the previously described methods and to compare these estimated values to the measured ones (Figure 7).The measured values of suction and those inferred by means of the MK method and from reinterpretation of CPTu are in very good agreement.
Some increase in suction values around depth of 1.6 m in CC1 and CC2 (Figure 7) is predicted using the CPTu-based method.Nevertheless, this is simply due to the locally high tip resistance (qc) values in CPTu-CC1 and CPTu-CC2 at a depth approximately equal to 1.6 m (heterogeneous layer).The same does not happen in CPTu-CC3 case.
Moreover, both the considered methods are capable of perfectly miming the trend of suction with an important reduction in the shallower part because of the raining during few days before CPTu testing.
Similar observations can be done as far as the CPTu tests at QUV site are concerned.In this case, the CPTu-based suction values mimic those predicted by the MK model but exhibit some discrepancies especially at depths in between 0.3 and 1.0 m (Figures 8 and 9).Unfortunately, the experimental measurement of suction at QUV site was not available and therefore it was not possible to establish if such discrepancies depend on limitations in the use of the MK model or of the proposed method.The measured values of suction and those inferred by means of the MK method and from re-interpretation of CPTu are in very good agreement.Some increase in suction values around depth of 1.6 m in CC1 and CC2 (Figure 7) is predicted using the CPTu-based method.Nevertheless, this is simply due to the locally high tip resistance (q c ) values in CPTu-CC1 and CPTu-CC2 at a depth approximately equal to 1.6 m (heterogeneous layer).The same does not happen in CPTu-CC3 case.
Moreover, both the considered methods are capable of perfectly miming the trend of suction with an important reduction in the shallower part because of the raining during few days before CPTu testing.
Similar observations can be done as far as the CPTu tests at QUV site are concerned.In this case, the CPTu-based suction values mimic those predicted by the MK model but exhibit some discrepancies especially at depths in between 0.3 and 1.0 m (Figures 8 and 9).Unfortunately, the experimental measurement of suction at QUV site was not available and therefore it was not possible to establish if such discrepancies depend on limitations in the use of the MK model or of the proposed method.The measured values of suction and those inferred by means of the MK method and from reinterpretation of CPTu are in very good agreement.Some increase in suction values around depth of 1.6 m in CC1 and CC2 (Figure 7) is predicted using the CPTu-based method.Nevertheless, this is simply due to the locally high tip resistance (qc) values in CPTu-CC1 and CPTu-CC2 at a depth approximately equal to 1.6 m (heterogeneous layer).The same does not happen in CPTu-CC3 case.
Moreover, both the considered methods are capable of perfectly miming the trend of suction with an important reduction in the shallower part because of the raining during few days before CPTu testing.
Similar observations can be done as far as the CPTu tests at QUV site are concerned.In this case, the CPTu-based suction values mimic those predicted by the MK model but exhibit some discrepancies especially at depths in between 0.3 and 1.0 m (Figures 8 and 9).Unfortunately, the experimental measurement of suction at QUV site was not available and therefore it was not possible to establish if such discrepancies depend on limitations in the use of the MK model or of the proposed method.It is important to stress that the proposed CPTu-based method directly assesses the effective stress state.As an alternative Khalili and Khabbaz [23] proposed a relationship for the estimate of χ parameter as a function of the ratio between the matric suction and the air entry value.On the other hand, Lu et al. [24] proposed a closed-form equation for the evaluation of the effective stress as a function of the matric suction and the Van Genuchten's model parameters [25].Unfortunately, these methods could not be applied to the present case studies because the necessary parameters were not available.Additionally, soil types validating the methods in [23,24] are very different than the silty-clay mixtures of the present study.

Conclusions
This paper presented a method for the determination of effective stress state in the vadose zone based on the interpretation of piezocone tests (CPTu).The preliminary validation was carried out by comparing the CPTu predicted values of suction against those measured in a monitored site and those predicted by using the so-called Modified Kovacs model (MK).
The values of suction obtained using the CPTu method were assessed based on the soil type and the measured degree of saturation, applying the effective stress equation for partly saturated soils proposed by Bishop.Despite that future research appears necessary to fully validate the proposed methodology.The obtained results are really encouraging and suggest that the proposed method based on the CPTu reinterpretation could be an expeditious and economic tool for the estimate of the in situ effective stress state in fine-grained soil deposits.A proper assessment of the latter is crucial in many cases to identify potential geological/geotechnical hazards as shallow landslides or failure of levees.In particular, the method is proposed for the shallow part of homogeneous clayey deposits.
If the proposed method is applied to evaluate the soil effective stress state the knowledge of the saturation degree is not required.In this case only the water table position and the results from CPTu tests (that can be carried out very easily and quickly) are required.Nevertheless, the knowledge of the saturation degree is necessary to estimate suction by using the "effective stress parameter (χ) vs. degree of saturation" relationships available in literature for specific soil types.
At the same time, the good performance of both the MK model and the relationship between the degree of saturation and the effective stress parameter (χ) of the Bishop equation are confirmed by the data shown.It is important to stress that the proposed CPTu-based method directly assesses the effective stress state.As an alternative Khalili and Khabbaz [23] proposed a relationship for the estimate of χ parameter as a function of the ratio between the matric suction and the air entry value.On the other hand, Lu et al. [24] proposed a closed-form equation for the evaluation of the effective stress as a function of the matric suction and the Van Genuchten's model parameters [25].Unfortunately, these methods could not be applied to the present case studies because the necessary parameters were not available.Additionally, soil types validating the methods in [23,24] are very different than the silty-clay mixtures of the present study.

Conclusions
This paper presented a method for the determination of effective stress state in the vadose zone based on the interpretation of piezocone tests (CPTu).The preliminary validation was carried out by comparing the CPTu predicted values of suction against those measured in a monitored site and those predicted by using the so-called Modified Kovacs model (MK).
The values of suction obtained using the CPTu method were assessed based on the soil type and the measured degree of saturation, applying the effective stress equation for partly saturated soils proposed by Bishop.Despite that future research appears necessary to fully validate the proposed methodology.The obtained results are really encouraging and suggest that the proposed method based on the CPTu re-interpretation could be an expeditious and economic tool for the estimate of the in situ effective stress state in fine-grained soil deposits.A proper assessment of the latter is crucial in many cases to identify potential geological/geotechnical hazards as shallow landslides or failure of levees.In particular, the method is proposed for the shallow part of homogeneous clayey deposits.
If the proposed method is applied to evaluate the soil effective stress state the knowledge of the saturation degree is not required.In this case only the water table position and the results from CPTu tests (that can be carried out very easily and quickly) are required.Nevertheless, the knowledge of the saturation degree is necessary to estimate suction by using the "effective stress parameter (χ) vs. degree of saturation" relationships available in literature for specific soil types.
At the same time, the good performance of both the MK model and the relationship between the degree of saturation and the effective stress parameter (χ) of the Bishop equation are confirmed by the data shown.

Figure 1 .
Figure 1.Measured tip resistance at Broni (PV, Italy).The highlighted layer shows the zone of influence of seasonal changes of the water table.Water table depths at 3.5 m and 5 m during dry and wet seasons respectively.

Figure 1 .
Figure 1.Measured tip resistance at Broni (PV, Italy).The highlighted layer shows the zone of influence of seasonal changes of the water table.Water table depths at 3.5 m and 5 m during dry and wet seasons respectively.

Figure 2 .
Figure 2. I c index from CPTu at Broni (PV, Italy) (Water table depths at 3.5 and 5 m during dry and wet seasons respectively) [5].

Figure 4 .
Figure 4. (a) Location of the two test-sites in northwestern Italy; (b) landslides map and main geomorphological features of Costa Cavalieri site; (c) landslides map and main geomorphological features of Santa Maria della Versa site.

Figure 4 .
Figure 4. (a) Location of the two test-sites in northwestern Italy; (b) landslides map and main geomorphological features of Costa Cavalieri site; (c) landslides map and main geomorphological features of Santa Maria della Versa site.

Figure 5 .
Figure 5. Soil profile and sketch of the Costa Cavalieri (CC) monitoring station.

Figure 5 .
Figure 5. Soil profile and sketch of the Costa Cavalieri (CC) monitoring station.

Figure 6
Figure 6 summarizes the results of CPTu at CC site.The Figure also shows the selected Ic target and Ic values as obtained from CPTu interpretation.

Figure 6 .
Figure 6.CPTu test results at CC site.

Figure 6 .
Figure 6.CPTu test results at CC site.

Figure 7 .
Figure 7. Measured and estimated values of suction at CC site.

Figure
Figure Measured and estimated values of suction at CC site.

Figure 9 .
Figure 9.Estimated values of suction at QUV site.

Figure 9 .
Figure 9.Estimated values of suction at QUV site.

Table 1 .
Selected geotechnical and physical properties of Costa Cavalieri test-site soil and weathered bedrock 1 .

Table 1 .
Selected geotechnical and physical properties of Costa Cavalieri test-site soil and weathered bedrock 1 .