
The integrals of diffusion processes are of significant importance in the field of finance, particularly in relation to stochastic volatility models, which are frequently employed to represent the temporal variability of stock prices. In this paper, we consider the strong consistency of the nonparametric kernel estimator of the transition density for second-order diffusion processes, using the moment inequalities of ρ-mixing sequences to demonstrate the strong consistency under some regularity conditions. Furthermore, the asymptotic mean square error is provided based on the deviation and variance of the transition density kernel estimator. The optimal bandwidth is found and thus the convergence rate of the kernel estimator is obtained. At the same time, our results are compared with the conclusions of the univariate density function.
Citation: Yue Li, Yunyan Wang. Strong consistency of the nonparametric kernel estimator of the transition density for the second-order diffusion process[J]. AIMS Mathematics, 2024, 9(7): 19015-19030. doi: 10.3934/math.2024925
[1] | Zhongqiang Liu, Åse Marit Wist Amdal, Jean-Sébastien L'Heureux, Suzanne Lacasse, Farrokh Nadim, Xin Xie . Correction: Spatial variability of medium dense sand deposit. AIMS Geosciences, 2020, 6(2): 149-150. doi: 10.3934/geosci.2020010 |
[2] | Elvis Ishimwe, Richard A. Coffman, Kyle M. Rollins . Predicting blast-induced liquefaction within the New Madrid Seismic Zone. AIMS Geosciences, 2020, 6(1): 71-91. doi: 10.3934/geosci.2020006 |
[3] | Laura Tonni, Guido Gottard . Assessing compressibility characteristics of silty soils from CPTU: lessons learnt from the Treporti Test Site, Venetian Lagoon (Italy). AIMS Geosciences, 2019, 5(2): 117-144. doi: 10.3934/geosci.2019.2.117 |
[4] | David Igoe, Kenneth Gavin . Characterization of the Blessington sand geotechnical test site. AIMS Geosciences, 2019, 5(2): 145-162. doi: 10.3934/geosci.2019.2.145 |
[5] | Aleksander S. Gundersen, Pasquale Carotenuto, Tom Lunne, Axel Walta, Per M. Sparrevik . Field verification tests of the newly developed flow cone tool—In-situ measurements of hydraulic soil properties. AIMS Geosciences, 2019, 5(4): 784-803. doi: 10.3934/geosci.2019.4.784 |
[6] | Paul W. Mayne, Ethan Cargill, Bruce Miller . Geotechnical characteristics of sensitive Leda clay at Canada test site in Gloucester, Ontario. AIMS Geosciences, 2019, 5(3): 390-411. doi: 10.3934/geosci.2019.3.390 |
[7] | Shehab S Agaiby, Paul W Mayne . CPTU identification of regular, sensitive, and organic clays towards evaluating preconsolidation stress profiles. AIMS Geosciences, 2021, 7(4): 553-573. doi: 10.3934/geosci.2021032 |
[8] | Priscilla Paniagua, Marco D'Ignazio, Jean-Sébastien L'Heureux, Tom Lunne, Kjell Karlsrud . CPTU correlations for Norwegian clays: an update. AIMS Geosciences, 2019, 5(2): 82-103. doi: 10.3934/geosci.2019.2.82 |
[9] | Marco D'Ignazio, Tom Lunne, Knut. H. Andersen, Shaoli Yang, Bruno Di Buò, Tim Länsivaara . Estimation of preconsolidation stress of clays from piezocone by means of high-quality calibration data. AIMS Geosciences, 2019, 5(2): 104-116. doi: 10.3934/geosci.2019.2.104 |
[10] | Elmar Buchner, Volker J Sach, Martin Schmieder . Ries impact deposits in the North Alpine Foreland Basin of Germany as a terrestrial analog site for impact-produced seismites and sand spikes on planet Mars. AIMS Geosciences, 2025, 11(1): 68-90. doi: 10.3934/geosci.2025005 |
The integrals of diffusion processes are of significant importance in the field of finance, particularly in relation to stochastic volatility models, which are frequently employed to represent the temporal variability of stock prices. In this paper, we consider the strong consistency of the nonparametric kernel estimator of the transition density for second-order diffusion processes, using the moment inequalities of ρ-mixing sequences to demonstrate the strong consistency under some regularity conditions. Furthermore, the asymptotic mean square error is provided based on the deviation and variance of the transition density kernel estimator. The optimal bandwidth is found and thus the convergence rate of the kernel estimator is obtained. At the same time, our results are compared with the conclusions of the univariate density function.
The engineering properties of a soil always exhibit variation from point to point within a soil volume. It is of interest to quantify the spatial variability in order to assess its importance, help plan future soil investigations and, if possible, reduce the uncertainty in a property by spatial averaging. This paper quantifies the spatial variability of the cone resistance in a sand deposit. The spatial variability of clays has been studied earlier. However, little information exists on the spatial variability of sands and silts.
The piezocone test has become one of the most frequently used methods to obtain the characteristics of a soil in situ. Piezocone tests in sand, even when closely spaced, can show large variability. Quantifying the distance over which the soil parameters in a sand are correlated can help plan site investigations and reduce the uncertainty in the soil parameters derived from the piezocone test. Reducing the uncertainty can have a significant impact on the reliability of geotechnical structures (e.g. [1,2,3,4]). In this paper, the variability in the corrected cone resistance is analyzed through first a visual inspection of the data, and second a calculation of the spatial correlation structure.
The spatial correlation structure of a geotechnical property is characterized by random field modelling [5]. The scale of fluctuation (SOF), describing the random field, defines the distance over which there is a significant correlation of a geotechnical parameter. It is recognized that soil properties are anisotropic by nature and the horizontal and vertical SOFs will probably differ (e.g. [1,5,6]). However, 3D spatial variability characterization remains a challenging task, due to limited site-specific investigation data, especially in the horizontal direction, and the lack of suitable characterization method [4,7,8].
The Norwegian GeoTest Sites (NGTS) project established five research sites in Norway in 2016. The sites are referred to as the sand, soft clay, quick clay, silt and permafrost sites. The establishment of the sites was funded by the Research Council of Norway and the aim of the project is to establish, characterize, share digital data and manage the test sites over the next 20 years [9]. Each of the sites have been carefully and extensively characterized to offer well documented data sets in five typical geo-materials for further testing, physical model testing, constitutive model development and modelling.
The objective of this study is to do a characterization of the inherent spatial variability of the sand in Øysand1. The results of piezocone tests (CPTU, but also denoted CPT in some of the figures and tables in this paper) from the Øysand benchmark site, near Trondheim in Norway, were utilized to establish a correlation model. The paper describes the soil profile at the test site (Section 2) and provides an overview of the results of the cone penetration tests selected for spatial variability analysis (Sections 3 and 4). In Sections 3 and 4, no calculation of the spatial variability is done. Section 5 summarizes briefly concepts from random field theory and the analysis approaches used in the study. The result of the analyses for Øysand sand in terms of scale of fluctuations in the vertical and horizontal directions are compared and discussed in Section 6. Conclusions are drawn in Section 7.
1 Øysand (or Oeysand) means “sand island” in English. The Norwegian letter Ø is often written Oe or simply O.
The Øysand research site is located in central Norway, approximately 15 km south of Trondheim. The locality sits on the south side of the Gaula River, at the head of the Gaulosen, an arm off the main Trondheim Fjord. Over the past thousand years, the river has mostly prograded westwards in the fjord. The ground surface at the site is at an elevation of 2.7 m above mean sea level.
The fluvial and deltaic deposit at Øysand consists of a 20–25 m fine silty sand with occasional high gravel content [10,11]. Figure 1 (top) presents a schematic longitudinal cross-section of a deltaic deposit, depicting its characteristic tripartite architecture (topset, foreset and bottomset) [12]. At Øysand, the stratigraphy features a general coarsening upward sequence as typically observed in deltaic deposits with topset, foreset and bottomset units (Figure 1 (middle)). The layers in these units can have different properties, geometry, fall and dip that can be linked to the depositional history at the mount of the river delta. Figure 1 (bottom) illustrates the layering at Øysand, based on the results of two resistivity cone tests [10].
Figure 2 presents a snapshot of the stratigraphy and index properties of the soils at the site, as obtained from in situ and laboratory tests2. The borehole log is for Borehole OYSB09, which is located very close to CPTU OYSC09.
2 The symbols used in Figure 1 are defined either in the figure text above each profile or at the bottom of the table. Other symbols not directly defined are γs, the density of solids and D10 and D60, the particle diameter (in mm) on the grain size distribution curve with 10 and 60% of the particles by weight passing.
Figures 3 and 4 illustrate two of the piezocone tests on the soil behavior classification charts [13], as interpreted from the pore pressure ratio Bq, and from the friction ratio Fr [14]. Figure 3, with data from CPTU OYSC42 (called CPT8 later in the spatial variability analysis), indicates that for this borehole, there is a fairly clean sand between 4 and about 10 m. Below 10 or 11m, the sand becomes more silty. Figure 4, with data from CPTU OYSC40 (called CPT6 later in the spatial variability analysis), indicates a clean sand to silty sand between 5 and 20 m. Note that the data points are so close to each other that the colors seem darker than in reality. Many of the CPT showed a change in gradient in the CPTU measured values from depth 13 or 14 m. The chart indicates a clean sand to silty sand. As the depth increases, the sand contains more silt and clay.
Based on all the boreholes and CPTUs (NGI [14,15]), a “more uniform” layer of clean sand to lightly silty sand was found between depths of 7 and 12 m was identified. Above this unit, the fine to coarse sand contained some gravel. Below 12 m, the sand contained an increasing portion of silty sand and sandy silt with traces of organic material. For this reason, and because the study was to address specifically medium dense sands, the spatial variability analyses in this paper concentrated on the sand between depths 7 to 12 m. This depth interval is shown with the zone between two pink horizontal lines in Figure 2.
Figure 5 (left) illustrates the locations of the piezocone tests at the Øysand site. There are a large quantity of tests in a cluster. These are shown in more detail in Figure 5 (right). The original testing plan was to have a minimum distance of 1.5 m to 2 m between two adjacent tests, but several tests were run at closer intervals.
The CPTU tests in or close to the cluster were used for the spatial variability analysis. Figure 5 (right) shows the 13 piezocone tests considered in this study: the piezocones from the cluster plus three CPTU tests 10 to 15 m away. For the spatial variability analyses, the tests were renumbered CPT1 to CPT13 (Figure 5 (right) and as will be listed in the second column of Table 1). All the tests except piezocone tests OYSC50, OYSC51 and OYSC52 (CPT11, 12 and 13) were run within an area of 18 m by 15 m (270 m2). Table 1 gives for each test the cone type, the zero drifts on the cone resistance qc measured, the sleeve friction fs and the pore pressure just behind the cone u2, the measurement interval in mm and the air temperature during testing. The air temperature was taken from meteorological records for all tests except tests with cone types 7 and 8 for which the air temperature was measured on site [16].
Test ID1) | CPT No.2) | Cone Type | Zero drifts | Measure interval | Temp.3) | Comment | ||
qc, kPa | fs, kPa | u2, kPa | mm | ℃ | ||||
OYSC21 | -- | 12 | −71.4 | 0.2 | −5.8 | 2 to 4 | 12 | |
OYSC22 | -- | 5 | 9.8 | 0.4 | −11.1 | 2 to 4 | 12 | |
OYSC23 | -- | 3 | 5.1 | −3.9 | −20.2 | 2 to 4 | 12 | |
OYSC24 | -- | 11 | −48.5 | 0.9 | 14.9 | 2 to 4 | 12 | |
OYSC25 | -- | 5 | −44.5 | 0.0 | −0.7 | 2 to 4 | 12 | |
OYSC26 | -- | 4 | −28.5 | −0.5 | −33.2 | 2 to 4 | 12 | Large zero drift u24) |
OYSC27 | -- | 12 | −41.6 | −0.5 | 2.8 | 2 to 4 | 12 | |
OYSC28 | -- | 5 | −4.5 | 0.0 | −4.5 | 2 to 4 | 12 | |
OYSC29 | -- | 11 | −127.3 | −6.4 | 15.6 | 2 to 4 | 12 | Large zero drift q. |
OYSC30 | -- | 3 | −19.4 | −0.9 | −31.1 | 2 to 4 | 12 | Large zero drift u24) |
OYSC31 | -- | 5 | −21.7 | −0.1 | −0.9 | 2 to 4 | 12 | |
OYSC32 | -- | 4 | −62.2 | −1.1 | 7.5 | 2 to 4 | 12 | |
OYSC34 | 1 | 1 | 5.4 | 0.6 | 0.0 | 10 | 18 | |
OYSC35 | 2 | 1 | 21.7 | 0.1 | 0.2 | 10 | 17 | Seismic |
OYSC37 | 3 | 1 | 21.7 | 0.2 | 0.2 | 10 | 18 | |
OYSC38 | 4 | 1 | 16.3 | 0.2 | 0.1 | 10 | 17 | |
OYSC39 | 5 | 6 | 168.0 | 80.3 | 2.0 | 10 | 17 | Large zero drifts. |
OYSC40 | 6 | 6 | 64.0 | −1.1 | −19.5 | 10 | 17 | |
OYSC41 | 7 | 6 | 56.0 | −0.9 | −12.8 | 10 | 17 | |
OYSC42 | 8 | 6 | 26.0 | 62.6 | −23.1 | 10 | 17 | Results above 4 m not incl. |
OYSC43 | 9 | 7 | 6.1 | 0.5 | −1.9 | 20 | 10 | Three result files combined |
OYSC44 | -- | 7 | −13.4 | −1.5 | 0.4 | 20 | 10 | Seismic cone |
OYSC45 | 10 | 7 | −20.2 | 0.0 | 13.1 | 20 | 10 | |
OYSC50 | 11 | 2 | 4.2 | 0.3 | 1.0 | 10 | 16 | |
OYSC51 | 12 | 2 | −18.9 | 0.3 | −1.0 | 10 | 16 | |
OYSC52 | 13 | 2 | −26.3 | 0.5 | −1.1 | 10 | 16 | |
1) Original number of CPTU, in site reports and in Figure 5, left. 2) Numbering of 13 CPTU selected for spatial variability analysis. Numbering, in Figure 5, right, is used in rest of paper. 3) Representative air temperature used to correct measured results. 4) It is assumed that the zero shift occurred at the end of test when the probe hit the harder layer. |
The selection of the CPTU tests was based on proximity and tests run under as similar conditions as possible. The air temperatures were between 16 and 18 ℃, except for three tests where the temperature was 10 ℃. To avoid introducing a statistical uncertainty, only the CPTUs with 10 and 20 mm measurement interval were used to evaluate the scale of fluctuation. Several types of cones were used during the characterization of the sand at the Øysand site. Tests with the seismic cone were excluded. The selected CPTU tests for spatial variability analysis were run with the cone types 1, 2, 6 and 7. Tests with cone type 1 were performed approximately 1.5 m away from other boreholes. The tests with cone type 7 were performed 2 m away from other boreholes.
Table 2 lists the characteristics of the piezocone at the Øysand site, including the zero drifts, the measurement intervals and some of the observations made. Tests CPT 1 to 13 are those that are used in the spatial variability analysis. As mentioned, only the sand between depths 7 and 12 m was included in the spatial variability analysis.
Cone type | Filter type | Saturation fluid |
1 | Bronze | Silicone ISOVG 100 |
2 | Bronze | Glycerine |
6 | Slot | Grease/Oil |
7 | Stainless steel, S/S 10 μ | Silicone oil, DC200, 50 cST |
Before the spatial variability analyses were started, the cone resistance data were looked at in detail. This section examines the cone data without any calculation of the spatial structure. The corrected cone resistance qt was judged to be a more representative parameter to compare the CPTU results than the measured cone resistance, because the derived qt accounts for the pore pressure generated during cone penetration and acting behind the cone tip.
Figure 6 presents the corrected cone resistance, qt, versus depth for each of the selected 13 CPTU tests and Figure 7 superposes the 13 tests in a single graph. There is very wide variability in the corrected cone resistance derived from the cone penetration measurements and the measured pore pressure u2, as qt varies between 1.5 and 4.5 MPa or more in the depth interval 7 to 12 m. This is not unexpected and links to the depositional history at the site. The layers found in the delta foreset (Figure 1) can have varying geometry and varying content of clay, silt and sand.
While examining the data, it was observed that the 13 CPTU seemed slightly out of phase depth-wise because of the structure (strike and dip) in the foreset of the deltaic deposits. A simple depth adjustment was made by shifting the corrected cone resistance qt data in each CPTU sounding up or down in order to match the peaks and troughs in the corrected cone resistances. Three figures illustrate this depth-wise phase difference and adjustment:
Figure 8, Cone resistance qt versus depth once they are “in phase”.
Figure 9, Comparison of qt-values for 13 CPTU brought in phase depth-wise.
Figure 10, Polynomial surface of depth adjustments bringing the 13 CPTU’s “in phase”.
A comparison of Figures 7 and 9 shows how much better the data from the 13 CPTU tests fit together once the depths have been adjusted. Table 3 lists the coordinates and the ground elevation at the location of each of the CPTs in the spatial variability analysis, as well the depths from the ground to the peak values of qt for each CPT. These depths were used to do the depth-wise adjustment of the cone resistance values.
Test ID | CPT No. | Northing | Easting | Elevation of ground surface | Depth from ground to peak values of qt |
OYSC34 | 1 | 7022898.60 | 562574.46 | 2.70 m | 10.9 m |
OYSC35 | 2 | 7022898.96 | 562575.89 | 2.69 m | 10.2 m |
OYSC37 | 3 | 7022897.82 | 562579.20 | 2.68 m | 8.6 m |
OYSC38 | 4 | 7022897.97 | 562580.71 | 2.68 m | 8.0 m |
OYSC39 | 5 | 7022898.50 | 562572.97 | 2.66 m | 12.0 m |
OYSC40 | 6 | 7022898.60 | 562574.46 | 2.70 m | 11.4 m |
OYSC41 | 7 | 7022898.96 | 562575.89 | 2.69 m | 10.5 m |
OYSC42 | 8 | 7022899.18 | 562577.38 | 2.70 m | 9.7 m |
OYSC43 | 9 | 7022900.64 | 562573.44 | 2.78 m | 10.0 m |
OYSC45 | 10 | 7022901.22 | 562577.41 | 2.82 m | 9.4 m |
OYSC50 | 11 | 7022911.04 | 562566.05 | 2.74 m | 9.1 m |
OYSC51 | 12 | 7022913.13 | 562566.82 | 2.72 m | 11.6 m |
OYSC52 | 13 | 702291238 | 562568.75 | 2.72 m | 9.2 m |
Figure 10 presents an interpolation of the relative depths. The three-dimensional figure shows the relative easting and northing coordinates on the horizontal plane and the depth to the peak value of the corrected cone resistance (in m). The bi-harmonic spline interpolation in Matlab was used for the interpolation of the depth to the peak qt value based on the known location of peak values in each CPT. One could also have done such interpolation analysis with kriging techniques, e.g. [17], and obtain an even better spatial representation pf the spatial representation of the variation of the depth-wise adjustment that had to be made.
The spatial variability analyses will be done with the corrected cone resistance qt from the measured cone resistance and pore pressure vs depth values (Figure 6), and for the corrected cone resistance qt adjusted depth-wise (Figure 8).
The uncertainty in soil properties estimates results from many different sources associated with data scatter and systematic errors [1,18]. The spatial correlation of a soil property, g(z), can be modelled as the sum of a deterministic trend component, t(z), and a randomly varying fluctuation term w(z) [2,19]:
g(z)=t(z)+w(z) | (1) |
Vanmarcke [19] proposed that the spatial variability can be modeled as a random field. Among random field models, stationary random fields are widely used due to their simplicity and because stationary random fields are possibly the only practical type of random field that can be characterized from limited data [21]. The removal of trends is the most common means for achieving stationary fields [2]. Linearly de-trended CPT data were used to characterize spatial variability. The scale of fluctuation δ is used to quantify the spatial variability. Small values of the scale of fluctuation δ indicate rapid fluctuations about the mean trend and large values of δ are indicative of long-distance correlations. The scale of fluctuation is often denoted as SOF. Details on the approach and examples of its use for cone resistance profiles in mixed materials can be found in [3] and [20].
Four autocorrelation functions are commonly used in geotechnical applications: the single-exponential (SNX), the cosine exponential (CSX), the second-order Markov (SOM), and the squared exponential (SQX) models [19,21,22]. Table 4 lists the formulations of the autocorrelation models and the relations to scale of fluctuation and Figure 11 illustrates the four autocorrelation functions as a function of normalized separation distance. In Figure 11, △h is the horizontal separation distance and δh is the scale of fluctuation in the horizontal direction. The functions in Figure 11 illustrate that for two CPTUs with a horizontal separation distance of 2δh, the auto-correlations are small.
Autocorrelation function | Autocorrelation model ρ(τ)* |
Single-exponential (SNX) | ρ(τ)=exp(−2|τ|/δ) |
Cosine exponential (CSX) | ρ(τ)=exp(−|τ|/δ)cos(τ/δ) |
Second-order Markov (SMK) | ρ(τ)=exp(−4|τ|/δ)[1+4|τ|δ] |
Squared exponential (SQX) | ρ(τ)=exp(−π(τ/δ)2) |
* where τ is the separation distance; and δ the scale of fluctuation. |
The four function models were tested with the Øysand CPTU data, and the single exponential function (SNX) was found to give the most consistent results. All four functions gave good convergence with the ACF method (Section 5.3). With the MLE method (Section 5.3), the SNX model showed the best convergence. In the rest of the paper, the results with the SNX model are shown.
Tables 5 and 6 list the scales of fluctuation obtained earlier. Table 5 lists earlier spatial variability results per soil type, Table 6 list earlier results per soil property and testing method. There are far fewer evaluations of horizontal correlation distance than evaluations of vertical correlation distance, because of the lack of data horizontally. The earlier data indicate that the horizontal autocorrelation distance is much shorter for sand than for clay. The vertical autocorrelation distance is usually shorter than the horizontal value. This reflects the dependence of the shear strength on the effective stress (or thickness of the overburden).
Soil type | Direction | SOF (m) | Reference (s) |
Offshore clay and sand | Horizontal | 30 | [23,24] |
Offshore sand | Horizontal | 14–38 | [25] |
Silty clay | Horizontal | 5–12 | [26] |
Clean sand | Vertical | 3 | [27] |
Mexico clay | Vertical | 1 | [27] |
Clay | Vertical | 1 | [18] |
Copper tailings | Vertical | 1 | [28] |
Sand | Vertical | 0.13–0.71 | [29] |
Clean sand | Vertical | 3.2 | [30] |
Sensitive clay | Vertical | 2 | [31] |
Silty clay | Vertical | 1 | [26] |
Sensitive clay | Vertical | 2 | [32] |
Sandy clay | Vertical | 0.30–1.22 | [33] |
Sand, clay | Vertical | 0.9 | [5] |
Clay | Vertical | 0.4 | [34] |
Soft clay | Vertical | 0.2–0.5 | [35] |
Stiff clay | Vertical | 0.1–0.55 | [36] |
Sand and clay | Vertical | 0.37–0.80 | [37] |
Sand, clay | Vertical | 0.13–1.11 | [38] |
Clay | Vertical | 2.53 | [42] |
Clay | Vertical | 0.16–1.17 | [7] |
Offshore clay | Horizontal | 317 | [40] |
Offshore clay | Horizontal | 4500 | [6] |
Clay | Vertical | 0.35–0.6 | [41] |
Property* | Soil type | Testing method** | δh(m) | δv(m) | |
su | Clay | Laboratory | - | 0.8–8.6 | |
su | Clay | VST | 46–60 | 2.0–6.2 | |
qc | Sand, clay | CPT | 3–80 | 0.1–3.0 | |
qc | Offshore soils | CPT | 14–38 | 0.3–0.4 | |
1/qc | Alluvial deposits | CPT | - | 0.1–2.6 | |
qt | Clay | CPTU | 23–66 | 0.2–0.5 | |
qc1N | Cohesive-behaviour soils | CPT | - | 0.1–0.6 | |
qc1N | Intermediate-behaviour soils | CPT | - | 0.3–1.0 | |
qc1N | Cohesionless-behaviour soils | CPT | - | 0.4–1.1 | |
fs | Sand | CPT | - | 1.3 | |
fs | Deltaic soils | CPT | - | 0.3–0.4 | |
FR | Cohesive-behaviour soils | CPT | - | 0.1–0.5 | |
FR | Intermediate-behaviour soils | CPT | - | 0.1–0.6 | |
FR | Cohesionless-behaviour soils | CPT | - | 0.2–0.6 | |
Ic | Cohesive-behaviour soils | CPT | - | 0.2–0.5 | |
Ic | Intermediate-behaviour soils | CPT | - | 0.6 | |
Ic | Cohesionless-behaviour soils | CPT | - | 0.3–1.2 | |
N | Sand | SPT | - | 2.4 | |
W | Clay, loam | Laboratory | 170 | 1.6–12.7 | |
wL | Clay, loam | Laboratory | - | 1.6–8.7 | |
γ' | Clay | Laboratory | - | 1.6 | |
γ | Clay, loam | Laboratory | - | 2.4–7.9 | |
E | Organic silty clay | Laboratory | - | 3.0 | |
σ'p | Organic silty clay | Laboratory | 180 | 0.6 | |
KS | Dry sand fill | PLT | 0.3 | - | |
ln(DR) | Sand | SPT | 67.0 | 3.7 | |
N | Sand | Laboratory- | 3.3 | 6.5 | |
* su = undrained shear strength; qc = cone resistance; qt = corrected cone tip resistance; qc1N = dimensionless, stress−normalised cone resistance; fs = sleeve friction; FR = stress−normalised friction ratio; Ic = CPT soil behaviour classification index; N = SPT blow count; w = water content; wL = liquid limit; γ' = submerged unit weight; γ = unit weight; e = void ratio; σ'p = preconsolidation stress; KS = subgrade modulus; DR = relative density; n = porosity. ** VST = vane shear testing; CPT = cone penetration testing; CPTU = piezocone testing; SPT = standard penetration testing; PLT = plate load testing. |
Several statistical techniques are available in the geotechnical literature for the estimation of random field parameters in spatial variability characterization, such as autocorrelation model fitting (AMF) (e.g. [1,7]) and maximum likelihood estimation (MLE) [8,42]. Three approaches were used to characterize the vertical and horizontal variability of the cone resistance qt of Øysand sand.
Auto-correlation model fitting (AMF) method for vertical SOF characterization: The basic principle of the autocorrelation function (ACF) approach is to fit a plausible theoretical autocorrelation function and evaluate the scale of fluctuation based on the parameters in the model (Table 4). There are two key parameters required to characterize a zero-mean stationary random field model using the ACF method, namely (1) the scale of fluctuation (δ) and (2) the auto-correlation function (ACF) (Table 4).
Maximum Likelihood Estimation (MLE) method for vertical and horizontal SOF characterization: The Maximum Likelihood Estimation (MLE) method has the objective to quantify the soil parameters by maximize the likelihood function. However, this method has limitations computationally when large number of data are to be analyzed, because a very large matrix needs to be solved. There is on-going research to improve the computation efficiency of the calculation process of MLE [8,41].
Simplified Vanmarcke method for horizontal SOF characterization: Vanmarcke [19] proposed an approximate relationship between the scale of fluctuation and the average distances between the intersections of the soil property curves and the trend line. This method is very simple to estimate the scale of fluctuation and is widely used for the reliability analysis of geo-structures [3].
The spatial variability analyses were run in the vertical and horizontal directions, and for two sets of data: the corrected cone resistance qt inferred from the CPTU measurements of cone resistance and pore pressure, and the corrected cone resistance qt after adjusting the depths to phase the data depth-wise.
Figure 12 presents four examples of the autocorrelation functions for vertical scale of fluctuation, δv, obtained with the ACF method. In all cases, the vertical scale of fluctuation (SOF) was one meter or less. The average of vertical scale of fluctuation, δv, from the 13 CPTU tests was 0.32 m, with a standard deviation (SD) of 0.07 m (Table 7). The coefficient of variation3 (COV) was 22%. For the fitting of the ACF function, the least square method was used on the linearly detrended data.
3 Coefficient of variation (COV) is the ratio of the Standard deviation (SD) to the Mean and is often expressed as a percentage. It is an indication of the uncertainty in the parameter characterized by a mean and standard deviation.
Analysis | Horizontal SOF, mean and [SD] (m) | Vertical SOF, mean and [SD] (m) |
AMF | --- | 0.3 [0.07] |
MLE | 3.5 | 1.6 |
Simplified Vanmarcke | Northing: 2.3 [1.0] Easting: 1.9 [0.8] |
--- |
Figure 13 presents one example of the autocorrelation functions for the horizontal scale of fluctuation, δh, obtained with the simplified Vanmarcke approach. Table 7 presents the results of the spatial variability analyses. The horizontal scale of fluctuation is between 1.9 and 2.3 m, indicating high variability of the soil over short distance. The COV of the horizontal scale of fluctuation was 42 to 50%, which is very high. The values of scale of fluctuation were very close in the northing and easting directions. It is expected that there is link between spatial variability and the direction of delta progradation which mostly was towards the west [43].
The vertical scale of fluctuation was close to 0.5 m for the AMF method, but as much as 1.6 m for the maximum likelihood evaluation (MLE) method. The MLE method tended to give higher values than the other methods, and was not always stable, when the matrices involved in the calculations became very large.
Table 8 presents the results of the spatial variability analyses for the phased data over an adjusted “common” relative depth of 4 meters (Figure 9). The average vertical scale of fluctuation, δv, from the 13 CPTU tests from the AMF method was 0.27 m with a standard deviation (SD) of 0.03 m (Table 8). The coefficient of variation (COV) was only 11%. The reduction in the COV is a direct result of the phasing with depth of the observed cone resistance data. For the fitting of the ACF function, the least square method was used on the linearly detrended data.
Analysis | Horizontal SOF, mean and [SD] (m) | Vertical SOF, mean and [SD] (m) |
AMF | --- | 0.3 [0.03] |
MLE | 2.5 | 0.8 |
Simplified Vanmarcke | Northing: 3.2 [1.7] Easting: 1.9 [1.2] |
--- |
The horizontal scale of fluctuation is between 1.9 and 3.2 m, indicating high variability of the soil over short distance. The COV of the horizontal scale of fluctuation was 32% in the northing direction and as high as 83% in the easting direction. The value of scale of fluctuation varied by more than one meter in the northing and easting directions. For this dataset, there is a link between spatial variability and the direction of delta progradation, with less autocorrelation in the direction of delta progradation and more correlation perpendicular to the progradation.
The vertical scale of fluctuation was 0.8 m for the maximum likelihood evaluation (MLE) method.
The CPTU data at the Øysand presented some very interesting spatial variability aspects. The study of vertical and horizontal spatial variability of the corrected cone resistance qt leads to the following observations. Only Table 8 with the depth-adjusted CPT-values was considered as it represents a more reliable result for scale of fluctuation in the sand layer than Table 7:
‒The vertical scale of fluctuation in Øysand sand was between 0.3 and 0.8 m (or 1 m rounded off for practice). The vertical scale of fluctuation for Øysand is similar to values obtained earlier for both sand and clay (Tables 5 and 6).
‒The horizontal scale of fluctuation was very short, between 1.9 and 3.2 m (depending on the direction), indicating short autocorrelation distances and high variability of the Øysand sand deposit. These values can be rounded off to 2 to 3 m in practice. Very high standard deviation of the horizontal scale of fluctuation (40 to 50%) were obtained, thus reinforcing the indication of variability of the Øysand sand. Such variability is typical for deltaic deposits, especially in the foreset unit where the layers are dipping at steep angles (Figure1).
‒Based on the results of the analyses and the calculation time it took to obtain the results, the autocorrelation function fitting (ACF) and the Simplified Vanmarcke methods seem to be more reliable than the maximum likelihood estimate (MLE). However, the MLE approach has a stronger and more fundamental mathematical basis. The challenge lies in how to solve large matrix properly and efficiently, not with the method itself.
‒Considering only 10 CPTU values in the spatial variability analysis in the horizontal direction (without the three outer CPTU’s (CPT 11, CPT12 and CPT 13 in Figure 5 right) did not change the results obtained.
‒Neglecting doing the depth-wise “phasing” leads to a larger uncertainty in the calculated scale of fluctuation.
To do reliable analyses of spatial variability, one needs to have access to large datasets comprising closely spaced data [44]. Abundance of data is rarely the case in practice. The paucity of data represents a major limitation for the use of the spatial variability statistical tools in practice.
The scale of fluctuation SOF is a convenient measure for describing the spatial variability of a soil property in a random field. It is a measure of the distance within which points are significantly correlated. Points separated by a larger distance than SOF will show little correlation, and practically no correlation will be observed when points are separated by a significantly larger distance than SOF.
This paper presented a study of the spatial variability of the cone resistance in the medium dense Øysand sand in Norway. The corrected cone resistance qt derived from the cone resistance and pore pressure measurements was used to establish the autocorrelation structure. The qt data were first analyzed visually. A depth interval of 7 to 12 m was selected for further analysis. To study spatial variability numerically the scale of fluctuation was identified. Several autocorrelations functions were considered, and the single exponential function was found to be the one offering the best autocorrelation. The spatial variability in the vertical and horizontal directions was analyzed statistically, using three different approaches, i.e. auto-correlation fitting, maximum likelihood estimation and simplified Vanmarcke method. The results indicate short autocorrelation distances of 3.5 m or less in the horizontal direction, suggesting a very variable sand at the Øysand site. In the vertical direction the scale of fluctuation was less than one meter, as obtained for other soils.
The CPTU data at the Øysand presented some very interesting spatial variability aspects. The study of vertical and horizontal spatial variability of the corrected cone resistance qt led to the following observations.
‒It was possible to establish a pattern of the corrected cone resistance with depth where the peaks and troughs of each CPTU could be “brought into phase” with a small adjustment of the depths (as illustrated by the surface in Figure 10). That aspect should probably be also considered at other sites with complex sand profiles where the data show a wide range of cone resistances (as in Figure 7). Neglecting to do such phasing will result in even greater spatial variability (and larger uncertainty). Such phasing further highlights the importance of understanding the geological history at a given site when assessing the spatial variability in soil properties.
‒The vertical scale of fluctuation was between 0.3 and 0.8 m, when one considers the depths-adjusted values of cone resistance. The calculated values of the vertical scale of fluctuation of Øysand sand compare well with the values obtained in earlier studies of the vertical scale of fluctuation of sand and clay.
‒The horizontal scale of fluctuation of Øysand sand was very short, between 1.9 and 3.2 m, again considering the depths-adjusted values of cone resistance. The results indicate very short autocorrelation distances, and high variability in the Øysand deposit. Very high standard deviations of the horizontal scale of fluctuation (COV of 40 to 50%) were obtained, reinforcing the variability of the measured values at the site. It is expected that such variability is typical for sands found in the foresets of fjord deltas.
‒The analysis results also illustrate the importance of understanding the geological setting, processes and history of a sand deposit to help interpret the result of the spatial variability analyses.
‒Based on the results of the analyses and the calculation time required to obtain the results, the ACF and the Simplified Vanmarcke methods seem to be more reliable at this time than the maximum likelihood estimate (MLE). However, improvement to this approach, mainly solving large matrices in a reasonable time, should gradually lead to good estimates.
To study spatial variability, one needs large datasets, rather evenly and closely spaced. This is rarely the case in practice. This is a major limitation for the use of the elegant statistical approaches. It was not possible from this study to establish general trends that are applicable to other sand sites, except perhaps for the vertical scale of fluctuation one meter or less, which is a similar result to that obtained for other soils earlier.
The authors wish to thank The Research Council of Norway for their generous infrastructure grant to establish the Norwegian GeoTest Sites infrastructure (NGTS) (Grant No. 245650/F50). The authors would like to acknowledge a large number of colleagues at NGI and NTNU who have contributed significantly to the soil characterization presented in this paper.
All authors declare no conflicts of interest in this paper.
[1] |
F. M. Bandi, P. C. B. Phillips, Fully nonparametric estimation of scalar diffusion models, Econometrica, 71 (2003), 241–283. https://doi.org/10.1111/1468-0262.00395 doi: 10.1111/1468-0262.00395
![]() |
[2] |
J. Fan, C. Zhang, A reexamination of diffusion estimators with applications to financial model validation, J. Am. Stat. Assoc., 98 (2003), 118–134. https://doi.org/10.1198/016214503388619157 doi: 10.1198/016214503388619157
![]() |
[3] |
F. Comte, V. G. Catalot, Y. Rozenholc, Penalized nonparametric mean square estimation of the coefficients of diffusion processes, Bernoulli, 13 (2007), 514–543. https://doi.org/10.3150/07-BEJ5173 doi: 10.3150/07-BEJ5173
![]() |
[4] |
J. Nicolau, Nonparametric estimation of second-order stochastic differential equations, Economet. Theor., 23 (2007), 880–898. https://doi.org/10.1017/S0266466607070375 doi: 10.1017/S0266466607070375
![]() |
[5] |
A. Gloter, Discrete sampling of an integrated diffusion process and parameter estimation of the diffusion coefficient, ESAIM-Probab. Stat., 4 (2000), 205–227. https://doi.org/10.1051/ps:2000105 doi: 10.1051/ps:2000105
![]() |
[6] |
A. Gloter, Parameter estimation for a discretely observed integrated diffusion process, Scand. J. Stat., 33 (2006), 83–104. https://doi.org/10.1111/j.1467-9469.2006.00465.x doi: 10.1111/j.1467-9469.2006.00465.x
![]() |
[7] |
J. Nicolau, Modeling financial time series through second-order stochastic differential equations, Stat. Probabil. Lett., 78 (2008), 2700–2704. https://doi.org/10.1016/j.spl.2008.03.024 doi: 10.1016/j.spl.2008.03.024
![]() |
[8] |
A. Carriero, T. E. Clark, M. Marcellino, Large Bayesian vector autoregressions with stochastic volatility and non-conjugate priors, J. Econometrics, 212 (2019), 137–154. https://doi.org/10.1016/j.jeconom.2019.04.024 doi: 10.1016/j.jeconom.2019.04.024
![]() |
[9] |
X. J. He, W. Chen, A closed-form pricing formula for European options under a new stochastic volatility model with a stochastic long-term mean, Math. Financ. Econ., 15 (2021), 381–396. https://doi.org/10.1007/s11579-020-00281-y doi: 10.1007/s11579-020-00281-y
![]() |
[10] | L. Rogers, W. David, Diffusions, Markov processes and martingales: Volume 2, Itô calculus, Cambridge University Press, 2 (2020). https://doi.org/10.1112/blms/21.1.106 |
[11] |
P. D. Ditlevsen, S. Ditlevsen, K. K. Andersen, The fast climate fluctuations during the stadial and interstadial climate states, Ann. Glaciol., 35 (2002), 457–462. https://doi.org/10.3189/172756402781816870 doi: 10.3189/172756402781816870
![]() |
[12] |
H. C. Wang, Z. Y. Lin, Local linear estimation of second-order diffusion models, Commun. Stat.-Theor. M., 40 (2011), 394–407. https://doi.org/10.1080/03610920903391345 doi: 10.1080/03610920903391345
![]() |
[13] |
Y. Y. Wang, L. X. Zhang, M. T. Tang, Re-weighted functional estimation of second-order diffusion processes, Metrika, 75 (2012), 1129–1151. https://doi.org/10.1007/s00184-011-0372-6 doi: 10.1007/s00184-011-0372-6
![]() |
[14] |
M. T. Tang, Y. Y. Wang, Q. Q. Zhan, Non parametric bias reduction of diffusion coefficient in integrated diffusion processes, Commun. Stat.-Theor. M., 51 (2022), 6435–6446. https://doi.org/10.1080/03610926.2020.1861298 doi: 10.1080/03610926.2020.1861298
![]() |
[15] |
S. C. Yang, S. Zhang, G. D. Xing, X. Yang, Strong consistency of nonparametric kernel estimators for integrated diffusion process, Commun. Stat.-Theor. M., 53 (2024), 2792–2815. https://doi.org/10.1080/03610926.2022.2148540 doi: 10.1080/03610926.2022.2148540
![]() |
[16] |
L. C. Zhao, Z. J. Liu, Strong consistency of the kernel estimators of conditional density function, Acta Math., 1 (1985), 314–318. https://doi.org/10.1007/BF02564838 doi: 10.1007/BF02564838
![]() |
[17] | L. G. Xue, M. Z. Fan, Strong convergence rates of the double kernel estimates of conditional density (in Chinese), Chinese J. Eng. Math., 1 (2001), 89–93. |
[18] |
S. Khardani, S. Semmar, Nonparametric conditional density estimation for censored data based on a recursive kernel, Electron. J. Stat., 8 (2015), 2541–2556. https://doi.org/10.1214/14-EJS960 doi: 10.1214/14-EJS960
![]() |
[19] |
A. Benkhaled, F. Madani, S. Khardani, Strong consistency of local linear estimation of a conditional density function under random censorship, Arab. J. Math, 9 (2020), 513–529. https://doi.org/10.1007/s40065-020-00282-1 doi: 10.1007/s40065-020-00282-1
![]() |
[20] |
N. Kadiri, M. Meghnafi, A. Rabhi, Strong uniform consistency rates of conditional density estimation in the single functional index model for functional data under random censorship, REVSTAT-Stat. J., 20 (2022), 221–249. https://doi.org/10.57805/revstat.v20i2.374 doi: 10.57805/revstat.v20i2.374
![]() |
[21] | H. Wang, W. X. Chen, B. J. Li, Large sample properties of maximum likelihood estimator using moving extremes ranked set sampling, J. Korean. Stat. Soc., 2024. https://doi.org/10.1007/s42952-023-00251-2 |
[22] | B. Sayedana, M. Afshari, P. E. Caines, A. Mahajan, Strong consistency and rate of convergence of switched least squares system identification for autonomous markov jump linear systems, IEEE. T. Automat. Contr., 2024. https://doi.org/10.1109/TAC.2024.3351806 |
[23] | Y. Li, Y. Y. Wang, M. T. Tang, Non parametric estimation of transition density for second-order diffusion processes, Commun. Stat.-Theor. M., 2023, 1–13. https://doi.org/10.1080/03610926.2023.2234521 |
[24] | I. Karatzas, E. Shreve, Brownian motion and stochastic calculus, Springer, 2014. https://doi.org/10.1007/978-3-319-57511-7_2 |
[25] | B. W. Silverman, Density estimation for statistics and data analysis, London: Chapman and Hall, 1986. https://doi.org/10.2307/2348849 |
[26] |
R. J. Hyndman, D. M. Bashtannyk, G. K. Grunwald, Estimating and visualizing conditional densities, J. Comput. Graph. Stat., 5 (1996), 315–336. https://doi.org/10.1080/10618600.1996.10474715 doi: 10.1080/10618600.1996.10474715
![]() |
[27] | D. W. Scott, Multivariate density estimation and visualization, Handbook Comput. Stat., 2012,549–569. https://doi.org/10.1007/978-3-642-21551-3_19 |
[28] |
E. Parzen, On the estimation of a probability density function and mode, Ann. Math. Stat., 33 (1962), 1065–1076. https://doi.org/10.2307/2237880 doi: 10.2307/2237880
![]() |
[29] |
S. C. Yang, Moment inequality for mixing sequences and nonparametric estimation, Acta Math. Sin., 40 (1997), 271–279. https://doi.org/10.12386/A1997sxxb0034 doi: 10.12386/A1997sxxb0034
![]() |
1. | Jean-Sebastien L’Heureux, Tom Lunne, Characterization and Engineering properties of Natural Soils used for Geotesting, 2020, 6, 2471-2132, 35, 10.3934/geosci.2020004 | |
2. | Zhongqiang Liu, Åse Marit Wist Amdal, Jean-Sébastien L’Heureux, Suzanne Lacasse, Farrokh Nadim, Xin Xie, Correction: Spatial variability of medium dense sand deposit, 2020, 6, 2471-2132, 149, 10.3934/geosci.2020010 | |
3. | Chao Liu, Hongrui Zhang, Ying Yuan, Aihong Zhou, Weiwen Liu, Wanying Guo, Influence of Correlation Distance of Soil Parameters on Pile Foundation Failure Probability, 2023, 15, 2071-1050, 4298, 10.3390/su15054298 | |
4. | Sheng Zhang, Hongjian Fu, Xueqian Ni, Xiaofeng Chen, Zhao Zhang, Influence of pre-seismic history on liquefaction resistance of seasonal frozen sand, 2024, 433, 09500618, 136638, 10.1016/j.conbuildmat.2024.136638 | |
5. | Shuairong Wang, Te Xiao, Guangyao Li, Yunhong Lv, Cong Dai, Liangtong Zhan, Yunmin Chen, Shuai Zhang, Spatial variability characterization of clayey waste soils and its impact on probabilistic stability assessment of a landfill slope, 2024, 19, 1861-1125, 5157, 10.1007/s11440-024-02244-9 |
Test ID1) | CPT No.2) | Cone Type | Zero drifts | Measure interval | Temp.3) | Comment | ||
qc, kPa | fs, kPa | u2, kPa | mm | ℃ | ||||
OYSC21 | -- | 12 | −71.4 | 0.2 | −5.8 | 2 to 4 | 12 | |
OYSC22 | -- | 5 | 9.8 | 0.4 | −11.1 | 2 to 4 | 12 | |
OYSC23 | -- | 3 | 5.1 | −3.9 | −20.2 | 2 to 4 | 12 | |
OYSC24 | -- | 11 | −48.5 | 0.9 | 14.9 | 2 to 4 | 12 | |
OYSC25 | -- | 5 | −44.5 | 0.0 | −0.7 | 2 to 4 | 12 | |
OYSC26 | -- | 4 | −28.5 | −0.5 | −33.2 | 2 to 4 | 12 | Large zero drift u24) |
OYSC27 | -- | 12 | −41.6 | −0.5 | 2.8 | 2 to 4 | 12 | |
OYSC28 | -- | 5 | −4.5 | 0.0 | −4.5 | 2 to 4 | 12 | |
OYSC29 | -- | 11 | −127.3 | −6.4 | 15.6 | 2 to 4 | 12 | Large zero drift q. |
OYSC30 | -- | 3 | −19.4 | −0.9 | −31.1 | 2 to 4 | 12 | Large zero drift u24) |
OYSC31 | -- | 5 | −21.7 | −0.1 | −0.9 | 2 to 4 | 12 | |
OYSC32 | -- | 4 | −62.2 | −1.1 | 7.5 | 2 to 4 | 12 | |
OYSC34 | 1 | 1 | 5.4 | 0.6 | 0.0 | 10 | 18 | |
OYSC35 | 2 | 1 | 21.7 | 0.1 | 0.2 | 10 | 17 | Seismic |
OYSC37 | 3 | 1 | 21.7 | 0.2 | 0.2 | 10 | 18 | |
OYSC38 | 4 | 1 | 16.3 | 0.2 | 0.1 | 10 | 17 | |
OYSC39 | 5 | 6 | 168.0 | 80.3 | 2.0 | 10 | 17 | Large zero drifts. |
OYSC40 | 6 | 6 | 64.0 | −1.1 | −19.5 | 10 | 17 | |
OYSC41 | 7 | 6 | 56.0 | −0.9 | −12.8 | 10 | 17 | |
OYSC42 | 8 | 6 | 26.0 | 62.6 | −23.1 | 10 | 17 | Results above 4 m not incl. |
OYSC43 | 9 | 7 | 6.1 | 0.5 | −1.9 | 20 | 10 | Three result files combined |
OYSC44 | -- | 7 | −13.4 | −1.5 | 0.4 | 20 | 10 | Seismic cone |
OYSC45 | 10 | 7 | −20.2 | 0.0 | 13.1 | 20 | 10 | |
OYSC50 | 11 | 2 | 4.2 | 0.3 | 1.0 | 10 | 16 | |
OYSC51 | 12 | 2 | −18.9 | 0.3 | −1.0 | 10 | 16 | |
OYSC52 | 13 | 2 | −26.3 | 0.5 | −1.1 | 10 | 16 | |
1) Original number of CPTU, in site reports and in Figure 5, left. 2) Numbering of 13 CPTU selected for spatial variability analysis. Numbering, in Figure 5, right, is used in rest of paper. 3) Representative air temperature used to correct measured results. 4) It is assumed that the zero shift occurred at the end of test when the probe hit the harder layer. |
Cone type | Filter type | Saturation fluid |
1 | Bronze | Silicone ISOVG 100 |
2 | Bronze | Glycerine |
6 | Slot | Grease/Oil |
7 | Stainless steel, S/S 10 μ | Silicone oil, DC200, 50 cST |
Test ID | CPT No. | Northing | Easting | Elevation of ground surface | Depth from ground to peak values of qt |
OYSC34 | 1 | 7022898.60 | 562574.46 | 2.70 m | 10.9 m |
OYSC35 | 2 | 7022898.96 | 562575.89 | 2.69 m | 10.2 m |
OYSC37 | 3 | 7022897.82 | 562579.20 | 2.68 m | 8.6 m |
OYSC38 | 4 | 7022897.97 | 562580.71 | 2.68 m | 8.0 m |
OYSC39 | 5 | 7022898.50 | 562572.97 | 2.66 m | 12.0 m |
OYSC40 | 6 | 7022898.60 | 562574.46 | 2.70 m | 11.4 m |
OYSC41 | 7 | 7022898.96 | 562575.89 | 2.69 m | 10.5 m |
OYSC42 | 8 | 7022899.18 | 562577.38 | 2.70 m | 9.7 m |
OYSC43 | 9 | 7022900.64 | 562573.44 | 2.78 m | 10.0 m |
OYSC45 | 10 | 7022901.22 | 562577.41 | 2.82 m | 9.4 m |
OYSC50 | 11 | 7022911.04 | 562566.05 | 2.74 m | 9.1 m |
OYSC51 | 12 | 7022913.13 | 562566.82 | 2.72 m | 11.6 m |
OYSC52 | 13 | 702291238 | 562568.75 | 2.72 m | 9.2 m |
Autocorrelation function | Autocorrelation model ρ(τ)* |
Single-exponential (SNX) | ρ(τ)=exp(−2|τ|/δ) |
Cosine exponential (CSX) | ρ(τ)=exp(−|τ|/δ)cos(τ/δ) |
Second-order Markov (SMK) | ρ(τ)=exp(−4|τ|/δ)[1+4|τ|δ] |
Squared exponential (SQX) | ρ(τ)=exp(−π(τ/δ)2) |
* where τ is the separation distance; and δ the scale of fluctuation. |
Soil type | Direction | SOF (m) | Reference (s) |
Offshore clay and sand | Horizontal | 30 | [23,24] |
Offshore sand | Horizontal | 14–38 | [25] |
Silty clay | Horizontal | 5–12 | [26] |
Clean sand | Vertical | 3 | [27] |
Mexico clay | Vertical | 1 | [27] |
Clay | Vertical | 1 | [18] |
Copper tailings | Vertical | 1 | [28] |
Sand | Vertical | 0.13–0.71 | [29] |
Clean sand | Vertical | 3.2 | [30] |
Sensitive clay | Vertical | 2 | [31] |
Silty clay | Vertical | 1 | [26] |
Sensitive clay | Vertical | 2 | [32] |
Sandy clay | Vertical | 0.30–1.22 | [33] |
Sand, clay | Vertical | 0.9 | [5] |
Clay | Vertical | 0.4 | [34] |
Soft clay | Vertical | 0.2–0.5 | [35] |
Stiff clay | Vertical | 0.1–0.55 | [36] |
Sand and clay | Vertical | 0.37–0.80 | [37] |
Sand, clay | Vertical | 0.13–1.11 | [38] |
Clay | Vertical | 2.53 | [42] |
Clay | Vertical | 0.16–1.17 | [7] |
Offshore clay | Horizontal | 317 | [40] |
Offshore clay | Horizontal | 4500 | [6] |
Clay | Vertical | 0.35–0.6 | [41] |
Property* | Soil type | Testing method** | δh(m) | δv(m) | |
su | Clay | Laboratory | - | 0.8–8.6 | |
su | Clay | VST | 46–60 | 2.0–6.2 | |
qc | Sand, clay | CPT | 3–80 | 0.1–3.0 | |
qc | Offshore soils | CPT | 14–38 | 0.3–0.4 | |
1/qc | Alluvial deposits | CPT | - | 0.1–2.6 | |
qt | Clay | CPTU | 23–66 | 0.2–0.5 | |
qc1N | Cohesive-behaviour soils | CPT | - | 0.1–0.6 | |
qc1N | Intermediate-behaviour soils | CPT | - | 0.3–1.0 | |
qc1N | Cohesionless-behaviour soils | CPT | - | 0.4–1.1 | |
fs | Sand | CPT | - | 1.3 | |
fs | Deltaic soils | CPT | - | 0.3–0.4 | |
FR | Cohesive-behaviour soils | CPT | - | 0.1–0.5 | |
FR | Intermediate-behaviour soils | CPT | - | 0.1–0.6 | |
FR | Cohesionless-behaviour soils | CPT | - | 0.2–0.6 | |
Ic | Cohesive-behaviour soils | CPT | - | 0.2–0.5 | |
Ic | Intermediate-behaviour soils | CPT | - | 0.6 | |
Ic | Cohesionless-behaviour soils | CPT | - | 0.3–1.2 | |
N | Sand | SPT | - | 2.4 | |
W | Clay, loam | Laboratory | 170 | 1.6–12.7 | |
wL | Clay, loam | Laboratory | - | 1.6–8.7 | |
γ' | Clay | Laboratory | - | 1.6 | |
γ | Clay, loam | Laboratory | - | 2.4–7.9 | |
E | Organic silty clay | Laboratory | - | 3.0 | |
σ'p | Organic silty clay | Laboratory | 180 | 0.6 | |
KS | Dry sand fill | PLT | 0.3 | - | |
ln(DR) | Sand | SPT | 67.0 | 3.7 | |
N | Sand | Laboratory- | 3.3 | 6.5 | |
* su = undrained shear strength; qc = cone resistance; qt = corrected cone tip resistance; qc1N = dimensionless, stress−normalised cone resistance; fs = sleeve friction; FR = stress−normalised friction ratio; Ic = CPT soil behaviour classification index; N = SPT blow count; w = water content; wL = liquid limit; γ' = submerged unit weight; γ = unit weight; e = void ratio; σ'p = preconsolidation stress; KS = subgrade modulus; DR = relative density; n = porosity. ** VST = vane shear testing; CPT = cone penetration testing; CPTU = piezocone testing; SPT = standard penetration testing; PLT = plate load testing. |
Analysis | Horizontal SOF, mean and [SD] (m) | Vertical SOF, mean and [SD] (m) |
AMF | --- | 0.3 [0.07] |
MLE | 3.5 | 1.6 |
Simplified Vanmarcke | Northing: 2.3 [1.0] Easting: 1.9 [0.8] |
--- |
Analysis | Horizontal SOF, mean and [SD] (m) | Vertical SOF, mean and [SD] (m) |
AMF | --- | 0.3 [0.03] |
MLE | 2.5 | 0.8 |
Simplified Vanmarcke | Northing: 3.2 [1.7] Easting: 1.9 [1.2] |
--- |
Test ID1) | CPT No.2) | Cone Type | Zero drifts | Measure interval | Temp.3) | Comment | ||
qc, kPa | fs, kPa | u2, kPa | mm | ℃ | ||||
OYSC21 | -- | 12 | −71.4 | 0.2 | −5.8 | 2 to 4 | 12 | |
OYSC22 | -- | 5 | 9.8 | 0.4 | −11.1 | 2 to 4 | 12 | |
OYSC23 | -- | 3 | 5.1 | −3.9 | −20.2 | 2 to 4 | 12 | |
OYSC24 | -- | 11 | −48.5 | 0.9 | 14.9 | 2 to 4 | 12 | |
OYSC25 | -- | 5 | −44.5 | 0.0 | −0.7 | 2 to 4 | 12 | |
OYSC26 | -- | 4 | −28.5 | −0.5 | −33.2 | 2 to 4 | 12 | Large zero drift u24) |
OYSC27 | -- | 12 | −41.6 | −0.5 | 2.8 | 2 to 4 | 12 | |
OYSC28 | -- | 5 | −4.5 | 0.0 | −4.5 | 2 to 4 | 12 | |
OYSC29 | -- | 11 | −127.3 | −6.4 | 15.6 | 2 to 4 | 12 | Large zero drift q. |
OYSC30 | -- | 3 | −19.4 | −0.9 | −31.1 | 2 to 4 | 12 | Large zero drift u24) |
OYSC31 | -- | 5 | −21.7 | −0.1 | −0.9 | 2 to 4 | 12 | |
OYSC32 | -- | 4 | −62.2 | −1.1 | 7.5 | 2 to 4 | 12 | |
OYSC34 | 1 | 1 | 5.4 | 0.6 | 0.0 | 10 | 18 | |
OYSC35 | 2 | 1 | 21.7 | 0.1 | 0.2 | 10 | 17 | Seismic |
OYSC37 | 3 | 1 | 21.7 | 0.2 | 0.2 | 10 | 18 | |
OYSC38 | 4 | 1 | 16.3 | 0.2 | 0.1 | 10 | 17 | |
OYSC39 | 5 | 6 | 168.0 | 80.3 | 2.0 | 10 | 17 | Large zero drifts. |
OYSC40 | 6 | 6 | 64.0 | −1.1 | −19.5 | 10 | 17 | |
OYSC41 | 7 | 6 | 56.0 | −0.9 | −12.8 | 10 | 17 | |
OYSC42 | 8 | 6 | 26.0 | 62.6 | −23.1 | 10 | 17 | Results above 4 m not incl. |
OYSC43 | 9 | 7 | 6.1 | 0.5 | −1.9 | 20 | 10 | Three result files combined |
OYSC44 | -- | 7 | −13.4 | −1.5 | 0.4 | 20 | 10 | Seismic cone |
OYSC45 | 10 | 7 | −20.2 | 0.0 | 13.1 | 20 | 10 | |
OYSC50 | 11 | 2 | 4.2 | 0.3 | 1.0 | 10 | 16 | |
OYSC51 | 12 | 2 | −18.9 | 0.3 | −1.0 | 10 | 16 | |
OYSC52 | 13 | 2 | −26.3 | 0.5 | −1.1 | 10 | 16 | |
1) Original number of CPTU, in site reports and in Figure 5, left. 2) Numbering of 13 CPTU selected for spatial variability analysis. Numbering, in Figure 5, right, is used in rest of paper. 3) Representative air temperature used to correct measured results. 4) It is assumed that the zero shift occurred at the end of test when the probe hit the harder layer. |
Cone type | Filter type | Saturation fluid |
1 | Bronze | Silicone ISOVG 100 |
2 | Bronze | Glycerine |
6 | Slot | Grease/Oil |
7 | Stainless steel, S/S 10 μ | Silicone oil, DC200, 50 cST |
Test ID | CPT No. | Northing | Easting | Elevation of ground surface | Depth from ground to peak values of qt |
OYSC34 | 1 | 7022898.60 | 562574.46 | 2.70 m | 10.9 m |
OYSC35 | 2 | 7022898.96 | 562575.89 | 2.69 m | 10.2 m |
OYSC37 | 3 | 7022897.82 | 562579.20 | 2.68 m | 8.6 m |
OYSC38 | 4 | 7022897.97 | 562580.71 | 2.68 m | 8.0 m |
OYSC39 | 5 | 7022898.50 | 562572.97 | 2.66 m | 12.0 m |
OYSC40 | 6 | 7022898.60 | 562574.46 | 2.70 m | 11.4 m |
OYSC41 | 7 | 7022898.96 | 562575.89 | 2.69 m | 10.5 m |
OYSC42 | 8 | 7022899.18 | 562577.38 | 2.70 m | 9.7 m |
OYSC43 | 9 | 7022900.64 | 562573.44 | 2.78 m | 10.0 m |
OYSC45 | 10 | 7022901.22 | 562577.41 | 2.82 m | 9.4 m |
OYSC50 | 11 | 7022911.04 | 562566.05 | 2.74 m | 9.1 m |
OYSC51 | 12 | 7022913.13 | 562566.82 | 2.72 m | 11.6 m |
OYSC52 | 13 | 702291238 | 562568.75 | 2.72 m | 9.2 m |
Autocorrelation function | Autocorrelation model ρ(τ)* |
Single-exponential (SNX) | ρ(τ)=exp(−2|τ|/δ) |
Cosine exponential (CSX) | ρ(τ)=exp(−|τ|/δ)cos(τ/δ) |
Second-order Markov (SMK) | ρ(τ)=exp(−4|τ|/δ)[1+4|τ|δ] |
Squared exponential (SQX) | ρ(τ)=exp(−π(τ/δ)2) |
* where τ is the separation distance; and δ the scale of fluctuation. |
Soil type | Direction | SOF (m) | Reference (s) |
Offshore clay and sand | Horizontal | 30 | [23,24] |
Offshore sand | Horizontal | 14–38 | [25] |
Silty clay | Horizontal | 5–12 | [26] |
Clean sand | Vertical | 3 | [27] |
Mexico clay | Vertical | 1 | [27] |
Clay | Vertical | 1 | [18] |
Copper tailings | Vertical | 1 | [28] |
Sand | Vertical | 0.13–0.71 | [29] |
Clean sand | Vertical | 3.2 | [30] |
Sensitive clay | Vertical | 2 | [31] |
Silty clay | Vertical | 1 | [26] |
Sensitive clay | Vertical | 2 | [32] |
Sandy clay | Vertical | 0.30–1.22 | [33] |
Sand, clay | Vertical | 0.9 | [5] |
Clay | Vertical | 0.4 | [34] |
Soft clay | Vertical | 0.2–0.5 | [35] |
Stiff clay | Vertical | 0.1–0.55 | [36] |
Sand and clay | Vertical | 0.37–0.80 | [37] |
Sand, clay | Vertical | 0.13–1.11 | [38] |
Clay | Vertical | 2.53 | [42] |
Clay | Vertical | 0.16–1.17 | [7] |
Offshore clay | Horizontal | 317 | [40] |
Offshore clay | Horizontal | 4500 | [6] |
Clay | Vertical | 0.35–0.6 | [41] |
Property* | Soil type | Testing method** | δh(m) | δv(m) | |
su | Clay | Laboratory | - | 0.8–8.6 | |
su | Clay | VST | 46–60 | 2.0–6.2 | |
qc | Sand, clay | CPT | 3–80 | 0.1–3.0 | |
qc | Offshore soils | CPT | 14–38 | 0.3–0.4 | |
1/qc | Alluvial deposits | CPT | - | 0.1–2.6 | |
qt | Clay | CPTU | 23–66 | 0.2–0.5 | |
qc1N | Cohesive-behaviour soils | CPT | - | 0.1–0.6 | |
qc1N | Intermediate-behaviour soils | CPT | - | 0.3–1.0 | |
qc1N | Cohesionless-behaviour soils | CPT | - | 0.4–1.1 | |
fs | Sand | CPT | - | 1.3 | |
fs | Deltaic soils | CPT | - | 0.3–0.4 | |
FR | Cohesive-behaviour soils | CPT | - | 0.1–0.5 | |
FR | Intermediate-behaviour soils | CPT | - | 0.1–0.6 | |
FR | Cohesionless-behaviour soils | CPT | - | 0.2–0.6 | |
Ic | Cohesive-behaviour soils | CPT | - | 0.2–0.5 | |
Ic | Intermediate-behaviour soils | CPT | - | 0.6 | |
Ic | Cohesionless-behaviour soils | CPT | - | 0.3–1.2 | |
N | Sand | SPT | - | 2.4 | |
W | Clay, loam | Laboratory | 170 | 1.6–12.7 | |
wL | Clay, loam | Laboratory | - | 1.6–8.7 | |
γ' | Clay | Laboratory | - | 1.6 | |
γ | Clay, loam | Laboratory | - | 2.4–7.9 | |
E | Organic silty clay | Laboratory | - | 3.0 | |
σ'p | Organic silty clay | Laboratory | 180 | 0.6 | |
KS | Dry sand fill | PLT | 0.3 | - | |
ln(DR) | Sand | SPT | 67.0 | 3.7 | |
N | Sand | Laboratory- | 3.3 | 6.5 | |
* su = undrained shear strength; qc = cone resistance; qt = corrected cone tip resistance; qc1N = dimensionless, stress−normalised cone resistance; fs = sleeve friction; FR = stress−normalised friction ratio; Ic = CPT soil behaviour classification index; N = SPT blow count; w = water content; wL = liquid limit; γ' = submerged unit weight; γ = unit weight; e = void ratio; σ'p = preconsolidation stress; KS = subgrade modulus; DR = relative density; n = porosity. ** VST = vane shear testing; CPT = cone penetration testing; CPTU = piezocone testing; SPT = standard penetration testing; PLT = plate load testing. |
Analysis | Horizontal SOF, mean and [SD] (m) | Vertical SOF, mean and [SD] (m) |
AMF | --- | 0.3 [0.07] |
MLE | 3.5 | 1.6 |
Simplified Vanmarcke | Northing: 2.3 [1.0] Easting: 1.9 [0.8] |
--- |
Analysis | Horizontal SOF, mean and [SD] (m) | Vertical SOF, mean and [SD] (m) |
AMF | --- | 0.3 [0.03] |
MLE | 2.5 | 0.8 |
Simplified Vanmarcke | Northing: 3.2 [1.7] Easting: 1.9 [1.2] |
--- |