Review

Design basis for Arctic infrastructure facilities

  • A discussion related to the selection of a proper design basis for Arctic infrastructure facilities is presented. The design basis, which must be sufficiently robust to ensure safe operations of the facilities during their planned lifetime, should include relevant information available from the local and indigenous people of the area. The need to properly document all data and assumptions made when preparing the design basis is highlighted, and it is emphasized that, in the case of the upgrading of the facilities, all new data collected must be included in the database. Documentation of all updating is necessary and must be available to all involved during maintenance activities and for possible later upgrading of the facilities. This database must be available digitally, as access to the basics for the design may be of particular concern in the sparsely populated Arctic region, where the distribution of paper copies takes a long time, particularly in the case of emergency situations. Therefore, this database must be protected from cyberattacks. A "Plan B" is needed to ensure that a backup of the fully updated design basis documents, as well as documentation of the "as built" facilities, is available at any time.

    Citation: Liv Brita Hætta Myrmel, Ove T Gudmestad. Design basis for Arctic infrastructure facilities[J]. AIMS Geosciences, 2023, 9(1): 86-94. doi: 10.3934/geosci.2023006

    Related Papers:

    [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
  • A discussion related to the selection of a proper design basis for Arctic infrastructure facilities is presented. The design basis, which must be sufficiently robust to ensure safe operations of the facilities during their planned lifetime, should include relevant information available from the local and indigenous people of the area. The need to properly document all data and assumptions made when preparing the design basis is highlighted, and it is emphasized that, in the case of the upgrading of the facilities, all new data collected must be included in the database. Documentation of all updating is necessary and must be available to all involved during maintenance activities and for possible later upgrading of the facilities. This database must be available digitally, as access to the basics for the design may be of particular concern in the sparsely populated Arctic region, where the distribution of paper copies takes a long time, particularly in the case of emergency situations. Therefore, this database must be protected from cyberattacks. A "Plan B" is needed to ensure that a backup of the fully updated design basis documents, as well as documentation of the "as built" facilities, is available at any time.



    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 1.  Schematic-longitudinal-cross-section-of-a-Gilbert-type-delta (top) [12]; typical stratigraphy at the Øysand research site (middle, not to scale); results of two resistivity cone soundings (ERT01 and ERT05) at the Øysand site (bottom) [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.

    Figure 2.  Example of borehole log for Borehole OYSB09 (B for borehole) at the Øysand site [15].

    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.

    Figure 3.  Robertson’s [13] soil behavior classification for Øysand CPTU OYSC42 (CPT 8).
    Figure 4.  Robertson’s [13] soil behavior type (SBT) chart vs depth for CPTU OYSC40 (CPT6).

    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.

    Figure 5.  CPTU locations at Øysand site: (a) Left: Map of all CPTUs (OYSC09 is highlighted in red); (b) Right: Map of 13 CPTUs selected for the spatial variability analysis.

    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].

    Table 1.  Overview of CPTU tests at Øysand site and those selected for spatial variability analysis.
    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.

     | Show Table
    DownLoad: CSV

    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.

    Table 2.  Filter characteristics of the probes used for the 13 tests 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

     | Show Table
    DownLoad: CSV

    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.

    Figure 6.  Corrected cone resistance qt for the 13 CPTUs in the spatial variability analysis.
    Figure 7.  Superposed corrected cone resistance qt, CPTUs in the spatial variability analysis.

    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 8.  Corrected cone resistance qt once the 13 CPTU’s are placed “in-phase” depth-wise.

    Figure 9, Comparison of qt-values for 13 CPTU brought in phase depth-wise.

    Figure 9.  Superposed qt “in phase” depth-wise, using 0 to 4 m as relative depth.

    Figure 10, Polynomial surface of depth adjustments bringing the 13 CPTU’s “in phase”.

    Figure 10.  Interpolation of the “corresponding” depths for the 13 CPTUs.

    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.

    Table 3.  Coordinates and elevation of CPTs included in spatial variability analysis.
    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

     | Show Table
    DownLoad: CSV

    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.

    Table 4.  Autocorrelation functions (after [21]).
    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.

     | Show Table
    DownLoad: CSV
    Figure 11.  Four autocorrelation functions normalized with separation distance ratio Δh/δh.

    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).

    Table 5.  Scale of fluctuation for clays and sands from the literature.
    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]

     | Show Table
    DownLoad: CSV
    Table 6.  Horizontal (δh) and vertical (δv) scales of fluctuation of geotechnical parameters [38].
    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.

     | Show Table
    DownLoad: CSV

    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.

    Figure 12.  Examples of the fitted autocorrelation functions for the vertical SOF, using the ACF method and the single exponential function (SNX).
    Table 7.  Scales of fluctuation (SOF) obtained for qt inferred from CPTU measurements.
    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]
    ---

     | Show Table
    DownLoad: CSV

    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].

    Figure 13.  Example of the horizontal SOF at 8m depth below ground surface, using simplified Vanmarcke method (circles indicate relative location of the CPT’s in the Northing direction).

    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.

    Table 8.  Scales of fluctuation (SOF) obtained for depth-adjusted qt-values.
    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]
    ---

     | Show Table
    DownLoad: CSV

    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] Gudmestad OT (2015) Establishment of design basis in the initial phase of critical infrastructure projects. J Crit Infrastruct 11: 183–193. https://doi.org/10.1504/IJCIS.2015.068616 doi: 10.1504/IJCIS.2015.068616
    [2] Rice E (1996) Building in the North. Alaska Science and Technology Foundation. Available from: https://www.uaa.alaska.edu/academics/college-of-engineering/departments/civil-engineering/arctic/_documents/building_in_the_north.pdf.
    [3] Gudmestad OT (2019) Rationale for Development of Design Basis for Barents Sea field developments. IOP Conf Ser Mater Sci Eng 700: 0120042. https://doi.org/10.1088/1757-899X/700/1/012042 doi: 10.1088/1757-899X/700/1/012042
    [4] Friedrich D, Hirnsberger M, Bauer S (2022) More than "Nature", Research on Infrastructure and Settlements in the North. Vienna: LIT-Verlag. Available from: https://www.thearcticinstitute.org/past-present-future-themes-arctic-infrastructure-settlements/.
    [5] Lennox C (2012) Natural resource development and the rights of minorities and indigenous peoples. In: State of the World's Minorities and Indigenous Peoples. Available from: https://minorityrights.org/wp-content/uploads/old-site-downloads/download-1114-Natural-resource-development-and-the-rights-of-minorities-and-indigenous-peoples.pdf.
    [6] Saxinger G (2018) Community Based Participatory Research as a Long-Term Process: Reflections on Becoming Partners in Understanding Social Dimensions of Mining in the Yukon. North Rev 47: a187–a206. https://doi.org/10.22584/nr47.2018.009
    [7] Orimolade AP, Furevik BR, Noer G, et al. (2016) Waves in polar lows. J Geophys Res Oceans 121: 6470–6481. https://doi.org/10.1002/2016JC012086 doi: 10.1002/2016JC012086
    [8] Buslaev G, Tsvetkov P, Lavrik A, et al. (2021) Ensuring the sustainability of Arctic industrial facilities under conditions of global climate change. Resources 10: 128. https://doi.org/10.3390/resources10120128 doi: 10.3390/resources10120128
    [9] Hjort J, Karjalainen O, Aalto J, et al. (2018) Degrading permafrost puts arctic infrastructure at risk by mid-century. Nat Commun 9: 5147. https://doi.org/10.1038/s41467-018-07557-4 doi: 10.1038/s41467-018-07557-4
    [10] Hjort J, Streletskiy D, Doré G, et al. (2022) Impacts of permafrost degradation on infrastructure. Nat Rev Earth Environ 3: 24–38. https://doi.org/10.1038/s43017-021-00247-8 doi: 10.1038/s43017-021-00247-8
    [11] The Council of the European Union (2008) On the identification and designation of European critical infrastructures and the assessment of the need to improve their protection. Available from: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32008L0114&from=EN
    [12] Wallace I (2020) Territorial Claims in the Arctic Circle: An Explainer. Obs. Available from: https://theobserver-qiaa.org/territorial-claims-in-the-arctic-circle-an-explainer.
    [13] Wei Z, Chen H, Lei R, et al. (2020) Overview of the 9th Chinese National Arctic Research Expedition. Atmos Oceanic Sci Lett 13: 1–7. https://doi.org/10.1080/16742834.2020.1675137 doi: 10.1080/16742834.2020.1675137
    [14] Svalbard Science Forum, Ny-Ålesund. Available from: https://www.forskningsradet.no/en/svalbard-science-forum/svalbard-research-planning/research-communities/ny-alesund/.
    [15] Research in Svalbard, The Research in Svalbard database (RiS). Available from: https://www.researchinsvalbard.no/
    [16] Myrmel LB, Hætta Gudmestad OT (2021) Cyber security for cities and rural areas in the Arctic Region, Rhodes, Greece.
    [17] Ten CW, Manimaran G, Liu CC (2010) Cybersecurity for critical infrastructures: Attack and defense modeling. IEEE T Syst Man Cy Part A Syst Humans 40: 853–865. https://doi.org/10.1109/TSMCA.2010.2048028 doi: 10.1109/TSMCA.2010.2048028
    [18] Guardian (2022) Pope in Canada to apologize for abuse of Indigenous children in church schools. The Guardian. Available from: https://www.theguardian.com/world/2022/jul/24/all
    [19] Rød B, Barabadi A, Gudmestad OT (2016) Characteristics of Arctic Infrastructure Resilience: Application of expert judgment. ISOPE Int Ocean Polar Eng Conf. Rhodes, Greece.
    [20] Mottahedi A, Sereshki F, Ataei M, et al. (2021) The Resilience of Critical Infrastructure Systems: A Systematic Literature Review. Energies 14: 1571. https://doi.org/10.3390/en14061571 doi: 10.3390/en14061571
    [21] Gudmestad OT (2020) Technical and economic challenges for Arctic Coastal settlements due to the melting of the Arctic ice and the permafrost. IOP Conf Ser Earth Environ Sci 612: 012049.
    [22] Orimolade AP, Haver S, Gudmestad OT (2016) Estimation of extreme significant wave heights and the associated uncertainties: A case study using NORA10 hindcast data for the Barents Sea. Mar Struct 49: 1–17. https://doi.org/10.1016/j.marstruc.2016.05.004 doi: 10.1016/j.marstruc.2016.05.004
    [23] Vanem E (2015) Uncertainties in extreme value modelling of wave data in a climate change perspective. J Ocean Eng Mar Energy 1: 339–359. https://doi.org/10.1007/s40722-015-0025-3 doi: 10.1007/s40722-015-0025-3
    [24] von Dimling TS, Lee H, Ingeman-Nielsen T, et al. (2021) Consequences of permafrost degradation for Arctic infrastructure—bridging the model gap between regional and engineering scales. Cryosphere 15: 2451–2471. https://doi.org/10.5194/tc-15-2451-2021 doi: 10.5194/tc-15-2451-2021
    [25] Cheng G, Sun Z, Niu F (2008) Application of the roadbed cooling approach in Qinghai—Tibet railway engineering. Cold Reg Sci Technol 53: 241–258. https://doi.org/10.1016/j.coldregions.2007.02.006 doi: 10.1016/j.coldregions.2007.02.006
    [26] Taarup J, Gudmestad OT (2022) Arctic supply chain reliability in Baffin Bay and Greenland. Polar Geogr 45: 77–100. https://doi.org/10.1080/1088937X.2022.2032447 doi: 10.1080/1088937X.2022.2032447
    [27] Shyu WH, Ding JF (2016). Key factors influencing the building of Arctic shipping routes. J Navigation 69: 1261–1277. https://doi.org/10.1017/S0373463316000254 doi: 10.1017/S0373463316000254
    [28] Aase JG, Jabour J (2015) Can monitoring maritime activities in the European High Arctic by satellite-based Automatic Identification System enhance polar search and rescue? Polar J 5: 386–402. https://doi.org/10.1080/2154896X.2015.1068534 doi: 10.1080/2154896X.2015.1068534
    [29] Andreassen N, Borch OJ (2020) Crisis and Emergency Management in the Arctic: Navigating Complex Environments. Routledge.
    [30] Sheehan R, Dalaklis D, Christodoulou A, et al. (2021) The Northwest Passage in the Arctic: A Brief Assessment of the Relevant Marine Transportation System and Current Availability of Search and Rescue Services. Logistics 5: 23. https://doi.org/10.3390/logistics5020023 doi: 10.3390/logistics5020023
    [31] Solberg KE, Jensen JE, Barane E, et al. (2020). Time to Rescue for Different Paths to Survival Following a Marine Incident. J Mar Sci Eng 8: 997. https://doi.org/10.3390/jmse8120997 doi: 10.3390/jmse8120997
    [32] Hætta LB (2017) The use of Remotely Piloted Aircrafts in reindeer herding, Master thesis in Technology and Safety in The High North. The Arctic University of Norway.
    [33] Kawagley A, Norris-Tull D, Norris-Tull R (1998) The indigenous worldview of Yupiaq culture: Its scientific nature and relevance to the practice and teaching of science. J Res Sci Teach 35: 133–144. https://doi.org/10.1002/(SICI)1098-2736(199802)35:2%3C133::AID-TEA4%3E3.0.CO;2-T doi: 10.1002/(SICI)1098-2736(199802)35:2%3C133::AID-TEA4%3E3.0.CO;2-T
    [34] Riseth J, Tømmervik T, Helander-Renvall E, et al. (2011) Sámi traditional ecological knowledge as a guide to science: Snow, ice and reindeer pasture facing climate change. Polar Rec 47: 202–217. https://doi.org/10.1017/S0032247410000434 doi: 10.1017/S0032247410000434
    [35] World Intellectual Property Organization, World Intellectual Property Organization. Available from: https://www.wipo.int/tk/en/tk/.
    [36] Yevseiev S, Melenti Y, Voitko O, et al. (2021) Development of a concept for building a critical infrastructure facilities security system. East-Eur J Enterp Technol 3: 63–83. https://doi.org/10.15587/1729-4061.2021.233533 doi: 10.15587/1729-4061.2021.233533
    [37] International Maritime Organization, Guidelines on maritime cyber risk management, 2017. Available from: https://wwwcdn.imo.org/localresources/en/OurWork/Security/Documents/MSC-FAL.1-Circ.3%20-%20Guidelines%20On%20Maritime%20Cyber%20Risk%20Management%20(Secretariat).pdf.
    [38] Det Norske Veritas, DNV-RP-0496. Recommended Practice: Cyber security resilience management for ships and mobile offshore units in operation, 2021. Available from: https://rules.dnv.com/docs/pdf/DNV/RP/2021-10/DNV-RP-0496.pdf.
    [39] International Standard Organization, ISO/IEC 2700 suite of standards for overall management of information security management, 2019. Available from: https://cyberwatching.eu/relevant-standards-cybersecurity-risk-management.
  • This article has been cited by:

    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
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(11705) PDF downloads(107) Cited by(1)

Figures and Tables

Figures(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog