Discussion of “Water Content Ratio: An Effective Substitute for Liquidity Index for Prediction of Shear Strength of Clays” by Beshy Kuriakose, Benny Mathews Abraham, A. Sridharan & Babu T. Jose

The discussers read the article by Kuriakose et al. (2017) with great interest, having completed another study on the variation of undrained strength with liquidity index (Vardanega and Haigh 2014). The authors draw attention to previous studies which use the w/w L (e/e L ) ratio to study the determination of the liquid limit (Nagaraj and Jayadeva 1981), and soil compressibility characteristics (e.g., Nagaraj and Murthy 1983, 1986; Griffiths and Joshi 1988). The discussers have some comments related to (a) the new relationships proposed by the authors and (b) whether the new method is sufficient to capture all the features of undrained strength variation, with changes in water content.


List of symbols a
An intercept b A slope c u Undrained shear strength c L Undrained shear strength at liquid limit e Void ratio e L Void ratio at liquid limit I L Liquidity index n Number of data points used to generate a bestfit line p pvalue R 2 Coefficient of determination (where R = the correlation coefficient) RD Relative deviation = 100(1 -R 2 ) 0.5 (see Waters and Vardanega 2009) R MW The computed ratio of strengths at the liquid limit and plastic limit SE Standard error of a regression WCR Water content ratio = w/w L w Water content w L Liquid limit a w L /w p The discussers read the article by Kuriakose et al. (2017) with great interest, having completed another study on the variation of undrained strength with liquidity index (Vardanega and Haigh 2014). The authors draw attention to previous studies which use the w/w L (e/e L ) ratio to study the determination of the liquid limit (Nagaraj and Jayadeva 1981), and soil compressibility characteristics (e.g., Murthy 1983, 1986;Griffiths and Joshi 1988). The discussers have some comments related to (a) the new relationships proposed by the authors and (b) whether the new method is sufficient to capture all the features of undrained strength variation, with changes in water content.

Author's Proposed Relationship
Equation (2) of the original paper appears to contain a small error. The discussers digitised the data from the author's paper (shown in Fig. 1 in the original paper) and fitted the following trend-line to the data: log 10 c u ð Þ ¼ 2:720 À 2:585 Eq. 1 has the following statistical measures: R 2 = 0.935; n = 27; SE = 0.105; RD = 25.5%; p \ 0.001. Perhaps Eq.
(2) in the original paper should have been: This equation would agree more closely with the equations for individual soils given in Fig. 2 of the original paper. Could the authors comment? Both Eqs.
(1) and (2) are plotted along with the digitised data on Fig. 1. We can see a good agreement and match with the data-set presented in the article under discussion. It should be noted that any discrepancy between Eqs. (1) and (2) is almost certainly due to errors that inevitably creep in during the digitisation process. The authors give the following regression equation in Fig. 3 of the original paper: Vardanega and Haigh (2014) (the discussers) assembled a database of 641 fall cone tests on 101 soils from 12 countries and produced the following equation: In Eq. (4) c L = 1.7 kPa (the assumption from Wroth and Wood (1978) of the undrained strength at the liquid limit). Equation 3 can be rearranged to the same form as shown in Wood (1990) to yield Eq. (3a). The assumed ratio of undrained strength between liquid and plastic limit is calculated to be about 30.6 which is more comparable with the factor of 34.3 reported by Vardanega and Haigh (2014) than the commonly assumed factor of 100 (e.g., Schofield andWroth 1968 andWood 1990).
In Eq. (3a) c L is computed to be about 1.36 kPa (this reduction may be explained by the fact that in the article under discussion vane shear data is reported, and that as observed from the original paper's Fig. 3, Eq. (3a) underestimates the measured undrained strengths at liquid limit (I L = 1). As the authors reiterate, Wroth and Wood (1978) reported a range of undrained strength at liquid limit from 0.7 to 2.65 kPa. Figure 2 shows Eqs. (3, 3a) and (4) plotted alongside the relationships of Leroueil et al. (1983) and Wroth and Wood (1978). The good agreement of the Cochin clay data-set with Eq. (4) is encouraging and further confirms that for semi-logarithmic trends linking undrained shear strength and liquidity index the often used factor of 100 between computed undrained strengths at plastic and liquid limits is too high. Now we turn to the suggestion of Kuriakose et al. (2017) of replacing liquidity index with water content ratio as a predictor of undrained soil strength. The discussers re-analysed the large database of fall-cone data that was reported in Vardanega and Haigh (2014) using the authors' method. For the sake of brevity the substantive details about the collection, sources of and analysis of the database will not be repeated here (readers should refer to Vardanega and Haigh 2014 for these details).
The best linear fit to the database (which can be rearranged into the form suggested by Kuriakose et al. 2017) relating values of undrained strength with WCR is given by Eq. (5).  Leroueil et al. (1983) and Wroth and Wood (1978) and Eqs. 3, 3a and 4 This can be re-arranged into the form of Eqs.
(1) and (2) to give: log 10 c u ð Þ ¼ 2:662 À 2:432 w w L ð6Þ Figure 3 shows the fitting of Eq. 5 to the database from Vardanega and Haigh (2014). Figure 4 shows that Eqs. (1) and (6) are functionally equivalent for the range of data shown. Figure 5 shows the accompanying predicted versus measured plot for Eq. (6). For the whole database (reported in Vardanega and Haigh 2014) an acceptable R 2 value of 0.874 is found for the method proposed by Kuriakose et al. (2017) (see Fig. 3). However, this is not as high as the R 2 value of 0.948 achieved from a fit using liquidity index in place of WCR (Vardanega and Haigh 2014). In summary, the water content ratio is somewhat effective in characterising the changes in undrained strength with water content. Lab vane data from Kuriakose et al. (2017) was used to develop a relationship (given as Eq. 1 here-which is similar to Eq. 2) which compares well with that obtained from regression analysis of the fall-cone database of Vardanega and Haigh (2014) (Eq. 6) (see Fig. 4). However, the fit of Eq. (6) to the database is not as good as that given using Eq. (4) as more data-points are observed to migrate outside the ±50% bounds on the prediction versus measured plot (cf. Fig. 5; Vardanega and Haigh 2014).
The increased uncertainty (shown on Fig. 5) originates from the variation of slopes seen for the relationships between WCR and undrained shear strength for different soils, which will now be demonstrated. Given that: and defining, a ¼ w L w P and WCR ¼ w w L we can then write: Therefore: Equation 9 shows clearly how the fitting coefficients in Fig. 7 of the original paper (Kuriakose et al. 2017) can be determined and highlights the point the authors' made in the original paper that WCR is related directly to I L via the a parameter. This said there must be a small error in the calculation of the 'R' values in Fig. 6 of the original paper as the best fit lines for a given soil on the two axes are functionally identical and should thus yield exactly the same value of the correlation coefficient (R), not the small variation shown. The WCR relationship is, marginally poorer in general owing to the greater variability of slope observed. Rearranging Eq. (4) we can say: and after substituting in (9): and then we can write: which has the form: These expressions reveal that the slope of the water content ratio versus the logarithm of undrained strength relationship is implicitly linked to a as shown in Fig. 8  In conclusion, while for a given soil the relationship between WCR and the logarithm of undrained shear strength is indeed linear, unless the relationship is being used solely for the extrapolation of measured undrained strength data for a given soil, there is no statistical basis to use a WCR relationship for undrained strength in preference to one based on liquidity index. One other point must also be made on extrapolation of undrained soil strength data using WCR. The linear relationships between the logarithm of undrained strength and water content break down close to the plastic limit. While this is a problem both for WCR and LI relationships, when using liquidity index it is obvious whether the water content of interest is approaching or drier than the plastic limit. The use of WCR obscures this trend, as the plastic limit is unknown and hence how close soil at a particular water content is to the onset of brittleness (see Haigh et al. 2013) is unclear. However, if extrapolations of fall cone or vane shear strength data close to the liquid limit are being performed then the approach suggested in Kuriakose et al. (2017) will probably suffice. When studying soil behaviour across the plastic range the discussers argue that liquidity index (or logarithmic liquidity index: see: Koumoto and Houlsby 2001) is still the preferred parameter.