Interpretation of field tests using geo-statistics and Kriging to assess the deep vibratory compaction of the Dike A21, Diavik Diamond Mine

This article presents an interpretation and statistical analysis of heavy dynamic probing (DPH) tests carried out on the compacted core of Dike A21, of the Diavik Diamond Mine. Due to the inhomogeneous and inherently random nature of soil and the uncertainties associated with even well controlled field tests such as DPH tests, the evaluation of in situ properties (such as density) should take into consideration the natural variability of the soil. Aided by the large amount of relatively concentrated field data, this paper presents a methodology for not only evaluating the field tests statistically through Kriging, but also for judging the compaction effort, taking into consideration the variability. This paper documents the necessary normalization procedure, analyzing a number of external factors such as the temperature, time since compaction (aging) and the presence of nearby works that may affect the field data. It was shown that the appraisal of soil improvement measures should consider the presence of natural and unavoidable fluctuations of the fill material. The Kriging estimate combined with an averaging procedure and subsequent comparison with the expected variability, quantified by the variogram, was applied which enabled the assessment of areas with insufficient compaction. Additionally, the influence of aging of the fill was demonstrated, where a clear increase in the DPH tests between 0 and 14 days since compaction was observed. It was also found that over-compaction leads typically to a lower relative density. The authors note that this statistical analysis was performed following the control and acceptance of the compaction of the dike and as such had no impact on the acceptance of the dike construction.


Introduction
The inherent variability of soil has been widely documented in the literature, particularly within the realm of field testing, where even homogenous soil layer can often show significant variation (e.g., [4,19,31,37,40]). Even in improved soils, significant variations in the results of field tests have been observed (e.g., [1,35]). In the construction industry, the treatment of areas subjected to soil improvements, which have been identified as not meeting a certain conformance criteria through field testing or other testing means can be ambiguously described in the relevant contract documents, e.g., a target value is specified, but the relationship for the interpretation of the target quantity is not. With regard to soil improvement, it is important that a differentiation is made between insufficient compaction and the effects of soil variability. The application of methods of geostatistical analysis affords an opportunity to assist in making such a judgment.
The following study is based on an analysis of field testing with the heavy dynamic probe (DPH), tests which were carried out within the inner core of a section between Stations 0?983 and 1?366 of the Dike A21, a part of the Diavik Diamond Mine. The location of Dike A21, a plan view with chainages, as well as a typical cross-section can be seen in Fig. 1 and Fig. 2. The goal of this article was to assess the success of the vibro-compaction effort based on state of the art statistical methods. To evaluate the relative density from the DPH tests, an interpretation method based on the cavity expansion theory and calibration parameters from a similar soil was used.
In total, the results of 201 DPH tests were made available between Station 0?983 and 1?366, of which 45 tests were excluded due to incompleteness of the data. For each of the tests, the following information was made available: • Number of blow counts per 10 cm penetration; • Date since compaction; • Location of the test-corresponding to a northing and easting; • A bucket count number, which describes the additional fill material dumped into the hole resulting from the dynamic testing; and • Aditional information regarding vibratory compaction works being carried out in the vicinity.
This well-documented and large dataset corresponded to an area of approximately 20 m 2 per DPH (based on a 10-mwide inner core zone) test and is ideally suited to stochastic analysis, a field of science which was pioneered by Vanmarcke [39], and has been gaining significant attention in recent times (e.g., [30] or [28]).
The use of dynamic probing to assess the in situ density of soils has long been established in the geotechnical community, e.g., Melzer [24] demonstrated its use in sand, Biedermann [2] in silt, and Butcher et al. [5] demonstrated its applicability to the interpretation of clay soils. The European Standards, EN 1997-2:2007 and EN 22476-2:2012 allow for the use of such methods for the assessment of the strength and deformation properties of soil, primarily for non-cohesive soils, but also under certain circumstances it allows for its use in the interpretation of cohesive soils. EN 1997-2:2007 provides empirical relationships between the dynamic probing results and the density index, however, these relationships do not account for the mean pressure, which strongly influences the penetration resistance, and thus must be used with caution.
One common limitation of penetration tests such as the cone penetration test (CPT), standard penetration test (SPT) and the dynamic probing (DP) in coarse-grained soils is related to the grain size and more specifically to the ratio between the mean grain size D 50 and the diameter of the tip of the probe d p . In the case of the CPT, it is generally accepted that the influence of the single grain on the penetration resistance can be disregarded for a ratio d p /D 50 of about 10 to 20 (e.g., [3,21,33]). It is however to be noted that recent simulations with the discrete element method (DEM) by Khosravi et al. [17] report an increase in the calculated tip resistance of only 10% for a d p /D 50 ratio of 2.1 compared to the baseline of 3.1 (with however, a 100% increase in the standard deviation, indicating a greater fluctuation), a finding that shows a reasonable agreement with the 20% increase in cone penetration resistances recorded by Bolton et al. [3] for Leighton Buzzard Sand. As far as the authors are aware, little research has been carried out into the dependence of the penetration resistance of dynamic probing in sandy gravels with D 50 around 2-10 mm as present in Dike A21 (e.g., Melzer [24] examined sands with a maximum D 50 of 4 mm and Biedermann [2] does not provided relationships for gravels). A newer study investigating the correlation between CPT and DPH results by Mahler and Szendefy [22] (who investigated soils with a D 50 of up to 10 mm) suggests that such correlations can be applied up to d p / D 50 & 5, with a larger scattering being observed for D 50 [ 2 mm. Furthermore, the work by Karray and Hussien [16] which synthesized a number of studies comparing the cone penetration resistance obtained from the CPT with shear wave velocities indicates an increase in the fluctuation with increasing mean grain size.
According to Ghafghazi et al. [12] the Becker penetration test (BPT) is more suitable than other field testing tools in the characterization of coarse soils due to the significantly larger probe diameter of 168 mm, compared to the DPH with a diameter of 50 mm, thus reducing single grain effects on the dynamic penetration resistance. For this reason and also because experience with the compaction control of previous neighboring and comparable dikes was available, acceptance criteria in the Diavik project were based on BPT results, which ultimately provided the confidence on the compaction of the core of Dike A21. Nevertheless, since the purpose of this contribution is to develop a geostatistical procedure to evaluate the success of deep compaction, it is assumed that the penetration resistance of the dynamic probing correlates with the density of the soils in the dike. Based on the experimental investigations of Melzer [24], who found out that the influence of the rod friction on penetration in granular materials can be disregarded for sounding depths less than around 10 m (the average height of the compacted core was 12 m) it is assumed for this study that the rod friction is negligible.

Site description
The A21 dike and the neighboring dikes A418 and A154 are located within the Lac de Gras, situated approximately 220 km south of the Arctic circle and 300 km north of Yellowknife, in Canada's Northwest Territories and were constructed to facilitate dewatering and the subsequent open pit extraction of diamonds.
Located within the arctic climate region, the temperatures (daily average) at the Lac de Gras range from -32°C in winter to ? 15°C in summer. These extremely low temperatures could account for the very high penetration resistance encountered in the upper layers of the dike fill and meant that boring had to be undertaken before DPH tests could be carried out in some areas. It is noted that that permafrost is present within the shallower regions of the lake, corresponding to approximately Sta. 0?450 m (Fig. 2). The geology at the A21 site consists predominantly of a thin layer of highly compressible silt with low plasticity clay, followed by a non-plastic glacial till and a bedrock of tonalite and quartz diorite.
The cross section of the dike, depicted in Fig. 2b, was constructed out of waste rock from a nearby stockpile consisting of an outer zone of primary and secondary crushed material (Zone 3 and Zone 2, respectively) with grain sizes of 10-1000 mm and an inner zone from tertiary crushing (Zone 1), consisting of sandy gravel with grain sizes between 0.1 and 50 mm with a fine content of \ 10% (see Fig. 2

c).
To prevent water seepage through the dike from the upstream side, a cutoff wall consisting of low strength plastic concrete was to be installed. To reduce the expected vertical and horizontal deformations of the dike, facilitate the construction of the cutoff wall, as well as to reduce the axial and bending forces on the wall and thus the extreme fiber stress through the expected negative skin friction as a result of the fill settling greater than the concrete following the dewatering, the supporting approximately 10-m-wide core area of the dike was to be compacted with the vibro flotation technique (see e.g., [34]). The planned depth of the treatment varied between 5 m and 22.5 m with an average depth of approx. 12 m. The DPH tests were carried out exclusively in the compacted inner core area.

Interpretation of the DPH tests
For the analysis it is assumed that a successful vibrocompaction of the inner core zone corresponds to a dense state and a relative density of 70%, with the success of the vibro-compaction being assessed by the results of DPH tests based on an estimation of the relative density from the number of blow counts for 10 cm of penetration, N 10 .
Several interpretation models for DPH test results can be found in the literature, one of which is the empirical correlation from Kralik [18] between the cone penetration resistance, q c (in MPa), and the number of blow counts for 20 cm penetration, N 20 , which simplifies to: Mahler and Szendefy [22] proposed a relationship based on a normalized CPT-DPH ratio, taking into consideration the mean grain size and the vertical stress, with after considering an average grain size of the inner core of 10 mm, treatment depth of 12 m and buoyant unit weight of c 0 ¼ 10:6 kN=m 3 the following relationship can be obtained: meaning that the relationship by Kralik [18] is conservative, corresponding to an average grain size of approximately 0.3 mm. For the statistical analysis, the more conservative relationship from Kralik [18] was used for the estimation of the q c as it was deemed to be more representative.
The relative density, D R (in %), can be estimated based on the cone penetration resistance q c as a function of the vertical effective overburden stress r 0 v0 by Jamiolkowski et al. [15] with for young, normally consolidated, quartz-based sands derived from calibration chamber tests. However, in a subsequent publication the correlation from Jamiolkowski et al. [14] was modified to take into account the effect of the calibration chamber boundaries on q c resulting in where p 0 is the mean effective stress and along with q c can be calculated in kPa. Assuming an effective peak frictional angle of u 0 of 42°based on site investigations for the densified fill from the stress deformation analysis of the dike, p 0 can be calculated from the well-known relationship K NC ¼ 1 À sinu 0 [13] leading to p 0 % 1:66r 0 v =3. However, to allow for the effects of compaction, for the subsequent statistical analysis a conservative estimate of K 0 ¼ 0:5 was selected, leading to p 0 % 2r 0 v =3 and thus a lower estimation of the relative density. Using the conservative estimate of K 0 this modification caused a reduction in the estimated q c of approximately 20%.
According to Cudmani [9] the following semi-empirical relationship can be used to estimate the density index: where the empirical shape factor (for granular soils) k q , and the limit pressure p LS , obtained through the expansion of a spherical cavity in a corresponding granular material that considers properties such as the grain shape and the grain hardness, are both dependent on the density index I D .
The shape factor can be calculated using the empirical relationship and the limit pressure with where a i and b i are material constants calculated by numerical simulations of the spherical cavity expansion using the hypoplastic model (see [9]), independent of I D and p 0 . For the interpretation of the DPH tests, material constants corresponding to Ticino sand were selected, a poorly graded quartz river sand (SP) according to the unified soil classification (see [9]). For Ticino sand the material constants are: a 1 =3.055, a 2 =-6.686, a 3 =-1.255, b 1 =0.794, b 2 =0.133 and b 3 =-1.379.
Thus, the cone penetration resistance was first determined from the DPH tests based on Eq. (1), after which the density index of the fill material was calculated as the solution to Eq. (5), (6) and (7). Though not completely representative of the fill material in the inner core-Ticino sand typically does not have any significant gravel component-it was chosen for use in this study due to the lack of experimental data for material calibration of the used fill material. According to Cudmani [9] however, no relationship was found between q c and the mean grain size, rather on the grain hardness and the grain shape.

Comparison of the DPH results
The required number of blow counts per 10 cm to achieve an assumed target relative density D R of 70% according to the relationship proposed by Jamiolkowski et al. [14,15] and Cudmani [9] is compared in Fig. 3 with the blow counts from the DPH tests average blow count per depth increment N 10 , the bounds N 10 AE s N 10 , where s N 10 is the standard deviation. Here, the interpretation models are also compared with the average and lower bound of the results of the DPH testing on previously completed dikes A154 and A418 several years following the construction of the dikes, labelled Criteria in Fig. 3.
This comparison reveals several conclusions. Firstly, if only the average blow count for all of the DPH tests is considered, the criteria of D R [ 70% has been fulfilledboth according to the interpretation models and according to experience from the neighboring dikes. Secondly, it would appear that the recorded blow counts from the DPH tests at dike A21 do not appear to show a sharp increase after approximately 14 m, as per the DPH tests at neighboring dikes A154 and A418. One possible cause for this difference is that the DPH tests were undertaken at dike A418 approximately 10 years after construction. It could be conjectured that the sharp increase in the recorded blow count be attributed to either 'aging' of the fill material, a well-documented phenomena-even in granular material (e.g., [26]), or to an increase of the shear strength of the material due to a decrease in the temperature (e.g., [20]). Thirdly, assuming an average compaction of [ 70% (based on the assumption that the compaction procedure was approved), comparison of the actual blow counts shows a greater agreement with the required blow counts with the methods according to Jamiolkowski et al. [14] and Cudmani [9]. These models, along with Jamiolkowski et al. [15], involve a logarithmic increase in the recorded blow counts with the depth (or overburden stress). It should be noted that despite a homogenous fill material and a consistent method of compaction, inspection of N 10 AE s N 10 both reveals a very high degree of fluctuation. Due to the fact that the upper bound of the recorded blow counts N 10 þ s N 10 lies above the theoretical attainable relative density according to Jamiolkowski et al. [14], the relationship proposed by Cudmani [9] appears to provide a better framework for interpreting the DPH data.

Stochastic analysis
Uzielli et al. [37] describe the practical modeling of the soil variability in terms of a geotechnical property as a random field based on the results of cone penetration tests with the following relationship: where t s ð Þ is a deterministic trend function, w s ð Þ is a fluctuating component, Z s ð Þ is the in situ soil property, and s is a spatial position (e.g., the vertical direction for this example of DPH testing). This approach will be followed in the analysis of the DPH tests by first de-trending the data through a normalization procedure.

Normalization of the DPH test results
The relative densities were determined at all data points through the interpretation relationships proposed by Jamiolkowski et al. [14] and Cudmani [9], using Eq. (1) and the inverse solution to Eqs. (5)-(7). In Fig. 4, the average of the achieved density index per depth increment is compared with the project requirement and the number of conforming tests is displayed (right). As can be seen from the figure, the criteria of I D ! 0:7 (D r ! 70%Þ are generally met by the average achieved density, for both interpretation methods. The average density indexes achieved according to Cudmani [9] show some areas below the criteria, e.g., at depths 6-8 m and 17-20 m; however, in Jamiolkowski et al. [14] all averages were above. This can be explained by the fact that the density index of the average blow counts is not equal to the average of the density index calculated from the individual blow counts due to the nonlinear relationship. Furthermore, the interpretation according to Cudmani [9] appears to be more conservative than that of Jamiolkowski et al. [14] with average density indexes being recorded of 0.75 and 0.85, respectively.
Without direct measurements of the actual density, it is impossible to judge from Fig. 4 which interpretation method provides the most accurate evaluation of the density index in the core of the dyke. However, under the assumption that a homogenous compaction of the fill material was achieved, including down to the base of the fill, a constant density index would be expected. The results of the normalization seem to confirm this, as aside from in the upper 5 m and lower 5 m, a relatively constant relationship of the density with the height can be seen. It can be conjectured that the higher values at the surface are caused by a combination of the additional compaction due to the construction traffic and the freezing of the upper soil layers in the winter months. The large variations in the lower depths are attributed to the small number of conforming records, which can be seen in the right of the figure, and the Part of the DPH documentation included recording: 1. The number of days since compaction was carried out at the compaction point(s) adjacent to the DPH tests; 2. Whether or not additional vibro-compaction was carried out since at other locations within the vicinity (50 m); and 3. Whether or not vibro-compaction within the vicinity (50 m) was being carried out during the DPH testing. Additionally for each vibro-compaction location, the following was made available: 4. The number of loader bucket counts filled during the compaction process.
The effect of ''aging'' of the granular fill-quantified by an increase in the density index over time, following placement or dynamic compaction-has been well documented in the literature (e.g., [6,25,27,34]). The physiochemical causes for this phenomenon are not the focus of this article, rather the influence of this phenomenonwhich is explored in Fig. 5b, where the average number of blow counts over the depth corresponding to: 1-7 days, 8-14 days, 15-21 days and more than 21 days since vibrocompaction at the compaction point adjacent to the test was carried out, and in Fig. 5c where the distribution of the estimated density index was plotted as a boxplot (± 25% quantiles, median, maximum and minimum limits as well as outliers) over the time since compaction. Analysis of the data in Fig. 5b shows a large increase in the number of recorded blow counts per 10 cm between tests carried out \ 7 days compared to tests carried out between 15 and 21 days following compaction; however, this increase is significantly less pronounced when compared with the results of the tests carried out after 21 days. This is confirmed in Fig. 5c where a pronounced increase in the density index can be observed within the first 14 days, both in terms of the median value and the ± 25% quantiles. From approximately 21 days, the density index appears to decrease slightly and a significantly higher variance can be observed, evidenced by the scattering of the medians, possibly due to a reduction in the data points or to other unkown factors. The time dependence can also be assessed based on the relationship of the normalized (vertical stress) tip resistance proposed by Charlie et al. [6], however, this relationship was not found to be able to explain the increase of the density index within the first two weeks. For further analysis, only results of DPH tests conducted 14 days or more after compaction were considered.
The influence of vibro-compaction within 50 m of the DPH tests, both prior to the test and during the test, is analyzed in Fig. 5a. Here it is evident that the vibro-compaction conducted concurrently to DPH testing results in a lower number of recorded blow counts, probably caused by excess pore-water pressures and the large vibrations travelling through the soil causing movement of the grains and allowing for an easier penetration of the DPH probe. Alternatively, vibro-compaction being carried out nearby, prior to performing the DPH tests (but not during) and in addition to the compaction point(s) adjacent to the test, resulted in a marked increase in the blow count, as compared to DPH tests where no (additional) vibro-compaction was carried out previously in the vicinity. This observation can be contrasted with the results in Fig. 6a, where the DPH tests corresponding to the 'Double Point' as illustrated in Fig. 6b record higher blow counts than at the 'Tripe Point.' One possible explanation for this observation is that extensive vibro-compaction causes a breakage of the grains in the immediate vicinity of the vibroprobe, thus causing a reduction in the shear strength and resistance to the penetration of the DPH probe. Another explanation is the nature of compacted coarse grained soil to dilate as a consequence of large shearing. To assess this effect, additional analysis of the test compaction results would need to be undertaken, including additional DPH testing at varying distances from the compaction points.
Based on the assumption of a homogenous fill material it is expected that the amount of fill material added during the vibro-compaction process, quantified by the bucket count would be representative of the amount of compaction experienced by the soil grains in the vicinity of the vibrocompaction. This hypothesis is explored by comparing the average density index with the adjacent average amount of fill material added at the compaction points as seen in Fig. 7a. Here no correlation (a negative correlation could also be argued) between the achieved density index and the amount of material added at the compaction point(s) can be observed, indicating a variation in the shear strength of the material-which governs the resistance to penetration from the DPH probe. From Fig. 7b it is evident that there is a no correlation between the amount of additional fill material  Cudmani [9] vs. the number of days since vibro-compaction; with the sample size n fed into the subsided compaction hole and the date of the vibro-compaction, indicating no dependence on the sequence of the vibro-compaction. Following from Fig. 7c there appears to be an incidental relationship between the amount of volume added during the vibro-compaction and when the DPH tests was carried out, possibly indicating that the areas tested with the DPH test first corresponded to an area with a naturally higher variability-which also appears to correlate with Fig. 7b where the vibro-compaction undertaken during Sep-Oct 2016 showed a greater variability.
The spatial correlation of the density index along the length of the dam (which in turn would be influenced by the variability of the source material) can be assessed from Fig. 8a, where there appears to be a slight negative trend. The amount of fill material added during the vibro-compaction process, quantified by the bucket counts, is also evaluated in Fig. 8b. Here also no significant trend was observed. In both plots an area at approximately Station 1?120 to 1?250, can be seen where the amount of fill material added is quite high and the density index appears to vary significantly. From Station 1?250 to 1?400, the fill material appears to be consistent, albeit with less scatter than Station 1?000 to 1?120 (based on the bucket counts). It can be concluded that this variability in the volume of material added is more likely to be the result of a spatial variability of the source material-e.g., more angular fill material leading to a greater propensity of the grains to break during the vibro-compaction than of varying compaction methods. In order to maintain a conservative estimate of the spatial variability of the dam fill and due to the absence of a physical justification, it was decided not to detrend the data to account for any linear correlation-which would have the effect of reducing the variability and increasing the correlation length.
In summary, the interpreted DPH test data were filtered, removing: (1) tests conducted 1-14 days after compaction at the adjacent compaction point; (2) tests conducted while vibro-compaction was being performed in the vicinity; and (3) tests conducted where additional vibro-compaction was undertaken with in the vicinity (in addition to the adjacent compaction point).

Spatial correlation
The typically observed increasing dissimilarity of soil samples with increasing separation distance is caused by the inherent natural variability of natural materials, a phenomena that occurs even in relatively homogenous layers, and has been the focus of significant research (e.g., [10, 11, 31, 32, 38, 29 and 36 ]). The similarity of soils from two locations can be expressed through the variogram 2c h ð Þ, which describes the variance of the difference between two quantities with a separation distance s as where Z s ð Þ is the soil property of interest. The interpreted and filtered DPH data are typically evaluated using the method of moments (MOM) estimator proposed by Matheron [23] of whereĉ h ð Þ is the experimental semi-variogram, and N h ð Þ is the number of distinct sample pairs (or test locations) with lag distance h. According to Cressie [8] the estimator proposed by Matheron [23] can be significantly affected by the presence of outliers and as such the test data were also evaluated using the MOM estimator from Cressie and Hawkins [7] These two estimators were compared, and it was found that the estimator from Cressie and Hawkins [7] was shown to reduce the scatter as well as the nugget variance.
The DPH data were evaluated considering for possible anisotropy of the spatial correlation by evaluating the horizontal separation between the DPH tests as XY ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi and the angle as u ¼ tan XY= DZ j j ð Þwhere u ¼ 0 represents a separation in the vertical direction and u ¼ p=2 a separation in the horizontal direction. The experimental semi-variogram was fitted with an single model exponential variogram (nugget = 0.002, range = 5.6 m, sill = 0.048) with the results being depicted in Fig. 9. From Fig. 9 there appears to be a certain degree of anisotropy; this was however not considered for because of the possibility that external factors (such as time since vibro-compaction), which by definition are not present when considering the vertical direction, could be causing the increased variance.

Interpolation
The Kriging interpolation method was used to interpolate the DPH test data in order to obtain a continuous representation of the density of the core fill over the length of the dike. Kriging, also known as the best linear unbiased predictor-best in that the prediction variance is minimized and unbiased in that clusters of sample points do not dominate the prediction, uses the semi-variogram to calculate the weights k i from samples (or observations) x i for Fig. 8 a Variation in the estimated density index over the dike and b corresponding bucket counts, accounting for number of days since compaction, after which the DPH test was performed the interpolation and can be performed either punctually or as a block average.
Many variations of Kriging exist; however, for the interpolation of the DPH test results Ordinary Kriging will be used, which assumes stationarity of the data, i.e., where Z Ã x ð Þ is the estimate of the random process Z x ð Þ. The theory can be summarized as follows: A value Z can be predicted at a location x 0 from the observation points aŝ with the condition that The weights k in vector form are determined through the relationship where A is the matrix of semi-variograms between the sample points i and j and b is the vector of semi-variograms between the sample points i and the interpolation point. The variance of the Kriging estimate can be calculated aŝ and the lower bound aŝ based on the 95% confidence interval and assuming a normal distribution of the interpolated values. The interpolation points are chosen such as to be along the centerline of the dike core based on a grid with points every 2 m (h orizontally) and 0.5 m (vertically) along the length and depth of the dam, respectively. The interpolated density indexes at the interpolation points according to the relationship proposed by Cudmani [9] can be seen in Fig. 10, wherein (a) the interpolated density index is shown using only the DPH tests subject to the three above-mentioned criterion (n = 83). In the figure blue corresponds to the theoretical maximum density index of 1, green corresponds to the target density index of 0.7, and red corresponds to a density index of 0.4. In Fig. 10b the interpolated density index is shown, this time considering all the DPH tests over 14 days of age (n = 123), in (c) a lower bound of the Kriging estimate [see Eq. (17)] is shown and in (d) the Kriging variance. From this plot it is evident that according to (a) the very large majority of the dike has a density index far exceeding the target. Isolated areas do appear to have a lower target density such as Station 1?150 to 1?200 (corresponding to the large bucket counts noted previously) and at approx. Station 1?350. It can also be seen in a comparison of (a) and (b) that excluding DPH tests carried out while vibro-compaction is being carried out in the area has a large impact in the interpolation of the density index, e.g., large areas at Station 1?000 to 1?050 are identified as not being adequately vibro-compacted. Due to the significantly high Kriging variance (due to relatively small correlation length and relative sparcity of the data) and spatial scattering of DPH tests the lower bound estimate is very low, with only few areas having a density index of higher than a loose to medium dense state. In summary, isolated pockets of areas with a density index lower than the specification can be identified (so called ''soft spots''); however, it is important to note that these areas do not generally extend to the entire depth of the dam fill. Additionally, care should be taken in using the Kriging estimate in areas where there is a high variances, as there are insufficient data to provide an accurate estimation-as such consideration should be given to the lower bound estimate, which as has been shown here, can be excessively conservative.

Assessment of the compaction effort
To assess the adequacy of the compaction, (1) the presence of areas with an apparent lower density must be ascertained, and (2) it must be evaluated whether this lower density is the result of a spatial variability of the shear strength of the fill material. For (1) the authors suggest a  Fig. 10a one can see that there are several isolated areas along the length of the dike which have been identified as probably having a density index lower than the target (e.g., between Station 1?150 and Station 1?350). For (2), without knowledge of the state of the fill material prior to compaction it is difficult to make a distinction between the natural variability of the soil and inadequate compaction. In this case, the authors recommend a procedure involving averaging the interpolated density indexes using a moving window based on the scale of fluctuations d (where d ¼ 2 Á r, and r is equal to the identified range of 5.6 m ; see [32]) and comparing the averaged interpolated density index to a so called ''expected density'' The expected density is determined based on the premise that the density indexes in the vicinity of a given density index I i D (where i is a point on the interpolation grid) will, with increasing distance h from i, approach the mean density index of all the points I D , and that this similarity can be modelled using the variogram. Based on this relationship, the expected density index averaged over a distance R equal to d, given I i D and I D , can be calculated as follows: where the function f I D ð Þ is defined by: (19). Under the assumption that the presence of poor compaction can be defined by the average interpolated density index within a moving window, over a certain distance R from i, defined a s e I R D , being significantly lower than the above described, expected density index I R D . Thus, the points that most likely represent insufficient compaction can judged based on comparing e I R D and I R D . In Fig. 11 the Kriging estimate (a), the average interpolated density over d (b) and the expected density index (c) based on Eq. (18) can be seen. Through a comparison between Fig. 11b and c one can see that there are areas where the moving window averaging results in a higher than expected I D , as well as areas with a lower than expected I D . This is to be expected, as the data cannot be described perfectly by the variogram. However it can be seen from a visual comparison that several areas, e.g., at Station 1?340, can be identified as having a lower averaged density index than expected and as such can be deemed to have inadequate compaction. The implication is that other areas, which appear to show an under-compaction (e.g., Station 1?150, due to their relatively limited size, can rather be attributed to the natural variability of the fill.

Conclusion
In summary, a method for interpreting the DPH tests within the core zone of Dike A21 within a statistical framework was presented; this involved determining the spatial correlation of the compacted fill material, interpolating the interpreted density index and using the spatial correlation to test for the presence of areas of insufficient compaction.
The normalization of the data to account for factors such as the mean pressure (depth) showed that the influence of the ''aging'' of the material, characterized as the time since vibro-densification, was found to be significant, where a linear increase in the interpreted density index was observed within the first 14 days. Thus, it is recommended that the control of densification shall be carried out at least two weeks after compaction. The increase penetration resistance with time was also observed in the recent DPH tests carried out on the Dike A418, which was constructed 10 years prior, where the resistances recorded were significantly larger than the values obtained for Dike A21. Thus, the penetration resistance obtained immediately after vibro-compaction is a conservative estimative, in comparison with the values expected sometime after construction, e.g., during the open pit mining. Additionally, it was shown through a comparison of the ''double point'' and ''triple point'' DPH tests that the over-compacted areas tended to perform worse, where it was hypothesized that this is due to breakage of the soil grains. As demonstrated through the normalization of the data and the analysis of the bucket counts, the fluctuation in the in situ density is more likely to be caused by fluctuations in the homogeneity of the fill material.
Following the normalization of the data, the spatial correlation in the vertical and horizontal was determined, where a single exponential model was fitted to the experimental variogram. The correlation model was then used to interpolate the density indexes along the center-line of the dike and at various depths with Kriging. Some isolated areas with a density lower than the target were identified; however, it was shown that these areas of lower density could mainly be attributed to natural variability of the shear strength of the fill. A significant area was identified at Station 1?340 and possibly, also at Station 1?200, which could be interpreted as being under-compacted. It is the opinion of the authors that re-compaction is not however expected to improve the situation considerably. Rather, this can be counterproductive as the improvement of stiffness and strength due to aging may be erased by vibrations.
In conclusion, the natural variability of the material should be considered when assessing the performance of soil improvement works and a methodology was presented for performing this assessment objectively, should the quantity of data allow for it, and thereby the risk of an incorrect assessment of the compaction effort can be reduced. For large projects validation of the interpretation method (e.g., in the form of chamber calibration tests), considering both the natural and improved state, should ideally be performed which would further reduce the uncertainty.
Acknowledgements The authors wish to thank the Diavik Diamond Mines Inc. (DDMI) for supporting this research by providing the original DPH data and the relevant project information and for permitting their use in this publication. Finally, we would like to make clear that none of the parties involved in the planning, design and construction of the Diavik Diamond Mine project participated in the preparation of this research work.
Funding Open Access funding enabled and organized by Projekt DEAL.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.
Data availability Some or all data, models or code that supports the findings of this study is available from the corresponding author upon reasonable request. Some or all data, models or code generated or used during the study are proprietary or confidential in nature and may only be provided with restrictions. Some or all data, models or code used during the study were provided by a third party. Direct requests for these materials may be made to the provider as indicated in Acknowledgements.