Minding the geotechnical data gap: appraisal of the variability of key soil parameters for slope stability modelling in Saint Lucia

Identification of failure thresholds and critical uncertainties associated with slope stability often requires the specification of geotechnical parameter values for input into a physically-based model. The variation of these parameters (including mechanical soil properties such as effective friction angle and cohesion) can have a significant impact on the computed factor of safety. These uncertainties arise from natural variations in soils, measurement techniques, and lack of reliable information. Researchers may use statistical analysis coupled with numerical simulation to determine possible ranges of slope factors of safety and the relative influence of geotechnical and other parameters, such as topsoil depth and rainfall. This study investigates the variation of geotechnical parameters observed on the island of Saint Lucia in the Eastern Caribbean. A database of particle size distributions, in-situ moisture contents, Atterberg and direct shear box test results is compiled from 91 samples of tropical soils in Saint Lucia. A study of various probability distributions shows that the Weibull distribution may be favoured for the effective friction angle of the Saint Lucian soils considered based on the Akaike information criterion, employed as an estimator of the relative quality of statistical models dealing with the trade-off between goodness-of-fit and simplicity of the model.


Notation c'
cohesion intercept, kPa CF clay fraction k number of parameters in a model L likelihood n number of data points PI plasticity index, % p p-value; probability of rejecting the null hypothesis and thus concluding that no correlation exists between two parameters r correlation coefficient R 2 coefficient of determination SF silt-clay fraction, % (taken here as percent passing the 0.075mm sieve) w nat natural water content, % w L liquid limit, % w P plastic limit, % x a variable β Weibull shape factor δ Weibull scale factor μ statistical mean σ' effective normal stress, kPa τ shear stress, kPa ϕ' effective friction angle, degrees ϕ' crit effective critical state friction angle, degrees ϕ' peak effective peak friction angle, degrees Introduction Landslides and geotechnical data Landslides triggered by rainfall and seismic events affect people, the built environment and economies worldwide. The susceptibility of slopes to landslides is due to a combination of preparatory factors including slope angle, geology, soil properties, groundwater conditions and vegetation cover. Human actions such as the construction of buildings and roads or changes to natural vegetation can affect these preparatory factors and decrease the stability of a slope. Rapid rates of urbanisation and land use change (both planned and informal) in landslide-prone developing countries mean that they are disproportionately affected by landslide hazards and their consequences (Petley, 2009). In such locations, landslides frequently occur due to a combination of localised factors such as slope-cutting, surcharge loading, altered surface water drainage patterns, and leakage from water services, water tanks and latrinesas reported in the city of São José dos Campos, Brazil, by Mendes et al. (2018), for example.
Implementation of landslide hazard assessments at the required resolution to inform landslide exposure, vulnerability and risk assessments, and risk management are often hindered by the limited availability of data in such regions. Acquisition of soil geotechnical information poses a particular challenge due to the expense of soil sampling in relevant locations and laboratory testing, the inherent heterogeneity of soils, the disturbance of in-situ soil structure and moisture content that occurs with sample extraction, and the resulting uncertainties associated with parameter values and sparse data. Furthermore, any such data from previous studies and projects tend to be inaccessible or undocumented. Thus, engineers and slope stability modellers encounter a geotechnical data gap when attempting to use physics-based models to investigate highly localised urban landslide processes (and potential mitigation measures) and when parameterising spatially-distributed landslide hazard models in Geographical Information Systems (GIS) at city or regional scales (e.g., Mutekanga et al. 2010;Ringrose-Voase et al. 2017;Du et al. 2018).

Effective friction angle
Slope stability calculations are sensitive to soil strength parameters such as the effective friction angle (ϕ', degrees), which is a measure of the resistance of the soil sample to shear stress, as defined by the Mohr-Coulomb failure criterion, Eq. 1 (e.g., Parry, 2004): where τ is the shear stress, c' the cohesion intercept, and σ' the effective normal stress. Simple factor of safety calculations (and the associated design charts) used by engineers for rapid slope stability assessments require input values for the soil effective friction angle, effective cohesion intercept, unit weight, slope geometry and groundwater conditions (e.g., Taylor, 1937;Anderson and Lloyd, 1991;Anderson et al. 1997;Michalowski, 2002 andLi et al. 2008). The c' and the ϕ' can also be back-calculated from field observations (e.g., van Asch, 1984). During shearing, a dense soil sample reaches a peak effective angle of friction (ϕ' peak ) at a relatively low strain; as strain increases the angle drops to a critical effective angle of internal friction (ϕ' crit ) (c.f. Craig, 2004, p.103) which should be used if designing against slope failure (Take and Bolton, 2011). However, ϕ' crit is still rarely reported in geotechnical characterisation studies. More sophisticated models representing dynamic hydrological processes of rainfall infiltration, sub-surface flow, pore water pressure and slope stability response have shown that in certain cases the choice of friction angle design values can affect computed factors of safety to the same extent as variations in rainfall for deeply weathered residual soil slopes in the humid tropics (e.g., Holcombe et al. 2016, Beesley et al. 2017and Shepheard et al., 2018a. Physically-based models can also be applied using reliability or stochastic approaches, in which input parameters are described using probability distributions rather than discrete values. Sensitivity analysis of model inputs and outputs may then be performed to identify dominant slope stability mechanisms for different classes of slope, system behavioural thresholds and ultimately data acquisition priorities.

Geotechnical correlations and classification systems
Where extensive laboratory testing is not available or affordable, it is common practice to use comparisons with other soil index properties, such as the plasticity index (PI) as a proxy for friction angle (e.g., Kenney, 1959;Brooker and Ireland, 1965, Ladd et al. 1977and Sorensen and Okkels, 2013. Such relationships allow geotechnical engineers to draw additional information from the routine soil tests carried out as part of civil engineering design and construction of foundations and transport infrastructure. Another source of proxies for soil mechanical properties can be agricultural soil maps and databases developed at national and worldwide scales. Agricultural soil data often include soil indices such as particle size distribution and soil bulk density and aspects of soil chemistry relating to clay particles, though not always the PI. Current engineering and agricultural databases of soil properties are most representative of North America, Europe and other regions with resources for soil testing. The most landslide-prone areas, typically encompassing developing countries in sub-tropical and tropical climatic zones, are less well represented. On the other hand, more recently, studies based in developing countries are becoming available (e.g., Roopnarine et al. 2012, Havaee et al. 2015. Roopnarine et al. (2012) reported the results of a study of Trinidadian soils and used soil physical properties to predict friction angle, reporting that the sand and clay fractions were the key predictors of residual and peak friction angle. Havaee et al. (2015) obtained significant positive correlation between friction angle and gravel content for Central Iran, and they also observed improvement in prediction models of soil shear strength derived using basic soil properties and normalized difference vegetation index.

Study aims
This paper investigates the extent to which the gap between available and required geotechnical data, specifically friction angle, can be reduced using soil data typically available in a landslide-prone humid tropical country. The first step was the collation of non-georeferenced hard copies of soil test results routinely undertaken by the Government of Saint Lucia, including soil classification and indexing, and direct shear test data. The resulting database is analysed with the following research objectives: (i) Based on preliminary work by Shepheard et al. (2018b) and Vardanega et al. (2018), to investigate whether significant soil property correlations can be derived to link soil friction angle to basic soil parameters for the Saint Lucia database (hereafter referred to as the 'SL database'); (ii) To distinguish different classes of Saint Lucia's soils using both a soil type classification framework and the well-known Casagrande approach (Casagrande 1947); (iii) To fit probability distribution curves to the soil properties, the functions that best fit the available data.
The overall aim is to provide local engineers and landslide researchers with information to make a priori estimates of key geotechnical parameters based on available local data instead of (or in addition to) more generic probability distributions derived from regional or global datasets.

Materials and methods
The Saint Lucia soil data were made available by the Government of Saint Lucia materials testing laboratory in 2016 as part of collaboration on the landslide hazard assessment component of a World Bank funded public infrastructure asset risk management project ('Vision 2030'). The database contains the results of a variety of historical laboratory tests carried out on soil samples from across the island. These data represent a subset of the laboratory's historic soil test records, which are currently archived in largely analogue format. Tests typically included the determination of field moisture content, particle size distribution and Atterberg limits as well as direct shear tests on partially-disturbed samples. The direct shear apparatus (DSA) in Saint Lucia has previously been  RRL, 1966, survey). Approximate soil sample locations marked by dots satisfactorily benchmarked against a modified DSA at the University of Bristol by repeating the tests reported in Lings and Dietz (2004), using loosely and densely compacted samples of Leighton Buzzard sand. The DSA testing procedure for Saint Lucia was also reviewed as part of the Vision 2030 project and found to have been consistently applied by the laboratory technicians. The Saint Lucia procedure is to extract samples from the field and place them into the direct shear The first stage of data processing was the creation of a digital database of the soil test records and the estimation of the location from which each sample was obtained. Of the 91 database entries there were 70 samples, taken from 44 sites, for which approximate map co-ordinates could be determined (assigned). Figure 1 indicates the 44 approximate sample locations associated with parent material geologies and the soil groups identified in the 1966 Saint Lucia soil and land-use survey (Regional Research Laboratory, 1966). The next stage was to identify classes of soils that might be expected to show distinct geotechnical properties and that could be distinguished from a combination of field observations, basic (often incomplete) soil descriptions in the original soil test records, and locations of the samples with respect to a particular parent material or soil type. The rationale for this classification is described in the next section along with the interpretation of the Saint Lucia soil data, the 1966 soil survey and local expert knowledge.

Geotechnical correlations
In soil engineering and geotechnics it is common to estimate engineering parameters that are complex or hard to measure by using a simpler property, or combination of properties (e.g., Kulhaway and Mayne, 1990;Phoon, 2014a, 2014b;Havaee et al. 2015;Zolfaghari, et al. 2016;Ahmed (2018), Bayat et al. (2018); Jie et al. (2018); Pham et al. (2018) and Schjønning and Lamandé, 2018). Preliminary statistical analysis of this database has been reported in Shepheard et al. (2018b) and Vardanega et al. (2018). Figure 2 shows the soil friction angles (and apparent cohesion intercepts) in the database plotted against the liquid limit (w L ), plastic limit (w p ), PI, silt-clay fraction (SF) and natural water content (w nat ). As clay fraction (CF) is not routinely available from the Saint Lucia geotechnical laboratory, SF is used in place of this more common parameter, which is defined as the percentage of material passing the 0.075 mm sieve. The effective friction angles are assumed to correspond to a normal stress range of 37.2 kPa to 112.4 kPa based on examination of local working practice in Saint Lucia. This is analogous to reports by Rouse et al. (1986) and Rouse (1990) of a series of effective friction angles of Dominican soils obtained from reversible shear box testing using a normal stress range of 41.8 kPa to 229.8 kPa.
In Fig. 2 the number of data-points available for each subcorrelation (n); the correlation coefficient (r) and the p-value (p) are shown underneath each sub-plot. The p-value represents the probability value estimated to test the null hypothesis that the coefficient of correlation between the variables is zero (no effect). A low p-value leads to the rejection of the null hypothesis, suggesting that the tested variable can be a meaningful addition for the prediction of friction angle. Examination of Fig. 2 shows that the two best correlations with effective friction angle are the natural water content (r = −0.60, n = 52) and the plasticity index (r = −0.39, n = 55), where the water content (w) is characterised by p < 0.001 and the plasticity index (PI) by p < 0.004, making these two parameters good candidates for a regression model for predicting the friction angle. Figure 3 shows the regression of ϕ' peak with PI used to construct Eq. (2). Figure 4 shows the regression of ϕ' peak with w nat to derive Eq. 3.
Equation 2 was found to have an R 2 = 0.15 for n = 55 with p = 0.004.
Equation 3 was found to have an R 2 = 0.36 for n = 52 with p < 0.001.  The reduction of effective friction angle with increasing plasticity index has been identified previously by many researchers (e.g., Kenney, 1959;Brooker and Ireland, 1965, Ladd et al. 1977and Sorensen and Okkels, 2013. The data are highly scattered for all the correlations of basic soil parameters with ϕ' peak and c'. Figure 2 shows only basic soil parameters that do not capture the effects of sample disturbance. Figure 5 illustrates how the regression of ϕ' peak with PI for the SL database compares with some previously quoted by Kenney (1959), Ladd et al. (1977), Brooker and Ireland (1965) and Sorensen and Okkels (2013). The Saint Lucia trend-line correlates with a lower ϕ' peak than other sources. This may be due to sample disturbance effects. The line from Brooker and Ireland (1965) is closest to that for the SL data. Shepheard et al. (2018b) performed multiple-linearregression analysis and found the following relationship to predict ϕ' peak (to within about ±50%) from basic soil properties: Equation 4 was found to have an R 2 = 0.56 for n = 47 with p < 0.001. Vardanega et al. (2018) found that when ϕ' peak was regressed against liquidity index for the SL database a R 2 of 0.43 was obtained.

Classification of soils based on formation and weathering
The natural structure of soil, and thus its geotechnical behaviour, is dependent on its parent material and the climate, topography, biological factors (micro-organisms, plants and animals) and time. The formation process for many finegrained soils can be categorised simply as either 'residual' or 'sedimentary'. Sedimentary soils are formed from minerals and organic materials that have been eroded (or produced by a volcanic eruption), transported and deposited by air, water or ice, and become consolidated (Wesley, 1990). Over time, sedimentary soil horizons may develop, and bonds can form between particles so that the material eventually behaves more like intact rock. For engineers, the stress history of sedimentary soils is recognised to be an important factor in determining their geotechnical behaviour. Conversely, taking the definition of The Geological Society Professional Handbook on Tropical Residual Soils (Fookes, 1997a, p10), residual soil profiles are formed from the in situ physical and chemical weathering of rock, leaching and accumulation of insoluble minerals and movement of fine particles, animal activity, plant root growth, and incorporation of organic materials, typically leading to cohesive soils and in some cases the cementation of the soil (e.g. laterites); essentially the reverse of the sedimentary soil formation process (Wesley, 1990). Large areas of the earth are mantled by residual soils, and deep profiles can form in humid tropical regions, such as the Caribbean, where readily available moisture and high temperatures lead to aggressive weathering. Both the mineralogy and structure of the parent material and the degree of weathering decomposition strongly influence the geotechnical behaviour of residual soils (cf. Wesley, 2009).
The 1966 soil survey of Saint Lucia (Regional Research Laboratory, 1966) primarily focuses on soils regarding their agricultural usage. It names 53 soil types based on parent material, geomorphological or topographical context, estimated particle size distribution, mineralogy, chemistry, nutrient availability, drainage, erodibility and depth. These soils are grouped into six classes according to the prevailing soil science taxonomy of the time (see legend of Fig. 1b). Although soil science terms such as 'latosols' are not often used by engineers (Wesley, 2009), they do encapsulate information on the soil formation process, parent material and structure. Thus, the soil series and classes defined in global and national soil maps, such as those compiled by FAO-UNESCO and reported by Hartemink et al. (2013), for example, can provide a starting point for investigating potential soil strength properties (e.g., Bonilla and Johnson, 2012). When the 1966 Saint Lucia soil survey information was combined with field observations, soil sample descriptions and depths, and local geotechnical engineering knowledge of the soils, it was possible to identify three distinct soil classes for this study: A) residual soils B) agglomerate soils and C) ash-derived soils, the second two classes having been formed through the weathering of sedimentary materials that were originally deposited by volcanic activity. As well as identifying which of the three classes each soil sample in the SL database belonged to, a further distinction was made based on the degree of weathering of the residual soils (A) and the soil matrix from the agglomerate material (B). For residual soils the weathered state of rock is often defined in terms of gradations from rock to soil, numbered from the fresh parent rock, Grade I, to completely weathered residual soil, Grade VI (GCO Geotechnical Control Office, 1982;Fookes, 1997aFookes, , 1997bToll, 2012). In this paper, the term 'residual soil' is used to describe soil samples in the Saint Lucia dataset that are likely to fall within both weathering Grade V (completely decomposed rock) and weathering Grade VI (soil) classes; due to the weathering grade classification being based largely upon visual inspection of samples (GCO Geotechnical Control Office, 1982), and the fact that the Saint Lucia soil sample descriptions did not explicitly record the weathering grade. This is also in keeping with Fookes (1997a, p12) who recognised that using the term 'soil' for weathering Grade VI only is "… somewhat restrictive for engineering purposes as much material normally described as 'soil' occurs below this Grade in the weathered profile".
The three classes of soils, and weathering subsets, identified for this study are as follows: & Class A: Tropical residual soils are the classical 'tropical residual soils'. They are reddish-brown, clayey soils with deep weathering profiles, and are particularly associated with the latosolic soils and andesitic polysols of the 1966 Regional Research Laboratory soil survey. Colluvium derived from these soils is also included as Grade VI material, based on field observations. Two subsets are defined, where 'A1' (n = 17) describes the upper layer of weathering Grades V and VI (soil), and 'A' (n = 23) is the underlying layer of weathered Grade IV material in which some of the relict rock structure and cementation is still present (which, in turn is underlain by Grades III, II and the Grade I parent material). Subset 'A1' contains only data which are known to pertain to the topsoil layer based on knowledge of the original samples. Subset 'A' contains all other 'A' data, excluding 'A1'. & Class B: Agglomerate soils derived from poorly sorted pyroclastic deposits. This soil is recognisable by the presence of boulders and gravels with a weakly cemented (and sometimes clayey) weathered soil matrix. Two subsets are defined, where 'B1' describes the most weathered layer in which the soil matrix is fully weathered and few large particles remain (large particles are manually removed from samples before direct shear testing), and 'B' the underlying less weathered layer. Because the soil matrix of Soil 'B' tested in the DSA is thought to be similar to that of Soil 'A', Soil 'A' is included in the statistical analysis of angle of friction: Subset 'B1' contains all data for Soil 'A', 'A1' and 'B' (n = 47). Subset 'B' contains all data for Soil 'A' (excluding A1) and Soil 'B' (n = 30). & Class C: Volcanic ash soils derived from pumice or tuff. This class contains only data points described recognised as a distinctive grey ash (or tuff) soil, sometimes highly cemented, but otherwise friable (n = 14). It is typically associated with the southern geological series (Fig. 1). Wright et al. 1984 give further commentary on the soils of Saint Lucia and is part of the so called 'Belfond Pumice'.
Of the 91 samples collected, 47 can be classed as residual soils ('A1', 'A', 'B1' and 'B'), according to these geological maps and soil descriptions (see column 2, Table 4). 1 1 Some soils in the database were classified as being comprised of multiple soil classes. Fig. 7 Weibull probability density function (PDF) fitted to (a) ϕ' peak , (b) c' for the SL database

Classification of soils based on soil index parameters
The data from Saint Lucia (Table 4) have been classified using the Casagrande chart in Fig. 6. The number of data points shown in Fig. 6 is below that in the whole database, as for many samples both a liquid limit and plasticity index are not available. The residual soils plot close to the A-line. According to the chart, 32 soils were classified as silts and 24 as clays. All data plot below the U-line.

Statistical Models
Probability density functions can be used to describe the variability of soil parameters and account for parametric uncertainty in reliability-based design, the development of decisionsupport tools, and stochastic physics-based modelling. When data for a particular site or region are limited, a geotechnical database and derived probability distributions, can provide useful information. Vardanega and Bolton (2016) suggest that reliability-based design analysis is best used for assessing the performance of geo-structures, as opposed to failure analysis. For instance, analysis of various predictors of undrained shear strength using a database and probability density function has been recently reported Phoon, 2014a, 2014b). Pham et al. (2016) and Pham et al. (2017) have used decision support trees for landslide susceptibility assessments in Vietnam and India. Chen et al. (2018aChen et al. ( , 2018b used datamining approaches to develop landslide susceptibility maps for the Shangnan County and Shangzhou district in China. Numerical models can be coupled with probability distributions when studying the effects of different sources of uncertainty on computed slope factors of safety, and the sensitivity of the system to different model inputs (e.g., Singh et al., 2014;Pianosi et al., 2016;Almeida et al., 2017). Recent work by Almeida et al. (2017) used Classification and Regression Tree (CART) Analysis (Breiman et al. 1984) in conjunction with the CHASM model (see Anderson and Howes, 1985;Wilkinson et al., 2002) to evaluate thresholds for slope failure due to parameter variation. Almeida et al. (2017) assigned uniform, normal and lognormal distributions to the parameter inputs for CHASM simulations in order to investigate the effects of various parameters on the failure of a modelled slope in Saint Lucia. 2 Lumb (1966) suggested the normal distribution gives a good fit for the w L and w P (and thus also PI), and the strength parameters ϕ' and c' for some Hong Kong soils. A review by Scott et al. (2003) noted that while the normal distribution may be the "least biased" model, many soil parameters cannot have negative valuestherefore the lognormal distribution may be more appropriate. Similarly, Lumb (1970) used the beta distribution to adapt the normal distribution to skewed data. However, to fit a beta distribution the data must be transformed into the range [0,1], thus increasing the number of computational steps required.

Goodness of fit tests
The goodness of fit of a probability distribution can be evaluated using the Anderson-Darling test statistic (Anderson and Darling, 1954), for which a smaller test statistic indicates a better fit (Minitab 17 Support, 2017). For the parameters in the SL database, Anderson-Darling tests have been calculated using Minitab (Minitab 17 Support, 2017), and the results quoted in Table 1. The ranking of different SL database candidate distributions can be done by directly comparing the Anderson-Darling tests statistics because they are all two-parameter distributions. Note that this comparison excludes the exponential distribution, which is a one-parameter distribution and shows relatively poor performance with respect to the other distributions. Table 1 shows that based on the Anderson and Darling criterion the Weibull is the preferred model for ϕ' peak , c' and PI. The lognormal distribution is favoured for w nat.
Another suitable approach is to rank distributions according to the maximum likelihood (e.g., Faber, 2012). Such ranking approaches are effective for the cases in which the candidate family distributions have the same number of parameters. For the cases in which the number of parameters is different, criteria employed in the information theory are more suitable as they "weight" the number of parameters considered allowing a trade-off between the goodness-of-fit of the model and its simplicity (i.e., the number of parameters). The corrected Akaike information criterion (AICc) is employed herein to rank the different fits (Akaike, 1974), to take into account the size of the sample and also to include a three-parameter distributionin this case the Generalised Extreme value (GEV). In general, Fig. 8 Probability plots for ϕ' peak data considering Normal, Lognormal, Weibull, Exponential and Generalised Extreme Value (GEV) fits the GEV distribution is the one with maximum likelihood, but the Akaike criterion accounts for the risk of 'over-fitting' the data (De Risi et al. 2017). The sample of data considered is not large enough, so the corrected Akaike formulation is considered as shown in Eq. (5), where k is the number of parameters, L is the likelihood and n is the size of the data sample considered.
As previously mentioned, from the comparison of Anderson-Darling statistics for each candidate distribution, it appears quite clear that, of the distributions assessed, the twoparameter Weibull (Eq. 6) (Montgomery et al., 2007, p.80) is the best-fitted for most of the data and, even when not ranked as the most suitable, it has test statistics that are relatively close to those of the best distribution (i.e., normal and lognormal), see Table 2. The same results are obtained if the fits are ranked using the corrected AICc: the Weibull always represents the best fit with the exception of the PI in which the GEV distribution provides a slightly better result.
The probability density function of Weibull is presented in Eq. 6, where β is a shape factor and δ is a scale factor, (e.g., Murthy et al., 2004, p.10).
These fitting parameters make the Weibull distribution highly adaptable to data. The Weibull model has been used to model aggregate size (e.g., Nimmo and Perkins, 2002). Fig. 7 shows that the shape factor β estimated for the variables are higher than 1.
The Weibull distribution appears to be a sensible choice of fit for the soils in the SL database. Like the lognormal distribution, the Weibull does not allow for values below zero; while the shape parameter in Eq. 6 allows the skewness to be represented (Weibull, 1951;Murthy et al., 2004). Figure 8 shows the comparison of the probability plots for the effective friction angle data fitted with normal, lognormal, Weibull, exponential and GEV. From the plots it is can be seen that the Weibull and GEV fit the data best. However, the Weibull is preferred, as it gets a result similar to the GEV, but with just two parameters, and is more efficient with respect to the size of the database. Table 3 shows the best fit distributions for the ϕ' peak and c' parameters when the database is divided by soil type (i.e. A1, A, B1, B, C and U), and again the Weibull distribution is either the most favoured or second most favoured based on the Anderson-Darling test statistics.

Conclusions
A database comprising information on 91 soil samples from Saint Lucia has been presented and classified according to soil type and formation. Simple regression analysis has been performed for some soil parameters, and it is seen that the strongest correlation for the friction angle is found with the natural water content. Additionally, a variety of probability distributions have been fitted to key parameters from the SL database. According to two ranking criteria (i.e. Anderson-Darling and Akaike), the Weibull distribution is preferred for ϕ' peak and c'.
These results are directly applicable for slope stability assessments in Saint Lucia. For local engineers, the database, soil property correlations and statistical distributions provide a basis for estimating soil properties in preliminary geotechnical analyses and for prioritising data acquisition (from basic soil sampling to potentially costly geotechnical investigations). For slope stability modellers the identification of appropriate parameter ranges and probability distributions can inform parametric studies at specific sites (e.g. Holcombe et al., 2016) or stochastic physics-based modelling of slope stability or over wide areas (Almeida et al., 2017), thus accounting for the effect on slope factor of safety of uncertainties in geotechnical and other slope properties.
In broader terms, this paper demonstrates that even in locations where traditional, well-curated geotechnical data are relatively scarce, it can be possible to compile useful databases of soil properties. Where statistically significant correlations and probability distributions fitted, this can provide a basis for geotechnical analysis and reliability-based design. In this paper, based on the available data, the Weibull distribution is shown to be more appropriate for certain geotechnical parameters in Saint Lucia. This contrasts with the usual assumption by engineers and environmental modellers that the lognormal distribution is often the best representation for soil parameters (Kulhawy, 2010, Hamm et al., 2006. It is suggested that while the selection of a lognormal distribution has a strong precedent and is the default choice for "operational reasons" (Rackwitz, 2000, p201), if sufficient local soil data can be compiled, it is worthwhile exploring the possibility that another distribution may give a better fit.
Acknowledgements The authors acknowledge the support of 'Landslide risk assessment of lifeline roads for public asset management and rainfallbased index insurance' which formed part of EP/P510920/1 'EPSRC