Evaluation of the liquefaction potential of a liquefaction-prone area is important for geotechnical earthquake engineering, both for assessment for site selection and for planning and new constructions. The liquefaction potential index for the city of Duzce in northwestern Turkey using the empirical relationships between the Standard Penetration Test (SPT) and the Shear Wave Velocity Test (VS) was investigated in this study. After, VS values based on SPT blow counts (N) were obtained from the alluvial soils in the city of Duzce. The liquefaction potential indexes of the soils were determined using the empirical relationships between the Standard Penetration Test (SPT) and the Shear Wave Velocity Test (VS) calculating for a probable earthquake of MW=7.2. In the result of the study, the liquefaction potential index (LPI) values were interpreted and compared evaluating the SPT N blow count values obtained from the study area. Based on the empirical relationships assumed for the soils, it was observed that there was not a perfect agreement between the results of the two methods. The liquefaction potential index values using the SPT N blow counts were found to be lower than those of the VS method.
1. Introduction
The liquefaction resistance of soils can be evaluated using laboratory tests such as the cyclic simple shear and cyclic triaxial and cyclic torsional shear tests. Additionally, field methods such as the Standard Penetration Test (SPT), Cone Penetration Test (CPT), and Shear Wave Velocity Test (VS) can be employed. The occurrence of liquefaction in soils is often evaluated using the simplified procedure originally developed and proposed by Seed and Idriss [1] based on the SPT blow counts correlated with the cyclic stress ratio (CSR), a parameter representing the seismic loading on the soil. This procedure has undergone several revisions and updates [2–4]. In addition to these, procedures have been developed based on the Cone Penetration Test (CPT), Becker Penetration Test (BPT), and small-strain Shear Wave Velocity (VS) measurements. Youd et al. [2] provided and enhanced a recent review of the Seed and Idriss simplified procedure and the in situ test methods commonly used to evaluate the liquefaction resistance of soils.
The use of VS to determine the liquefaction resistance is influenced by factors such as confining stress, plasticity, and relative density [5–7]. In situ VS can be measured by several seismic tests, including cross hole, down hole, seismic cone penetrometer (SCPT), suspension logger, and spectral analysis of the surface waves (SASW) [8].
During the past two decades, several procedures have been proposed to estimate liquefaction resistance based on VS [8]. These procedures were developed from laboratory studies [8–15], analytical studies [16, 17], penetration VS equations [18, 19], and in situ VS measurements at earthquake sites [20–22]. Some of these procedures follow the general approach of the Seed-Idriss simplified procedure, in which the VS is corrected with the cyclic stress ratio. This paper presents the results of the comparison between the VS and SPT methods of soil liquefaction potential evaluation carried out in Duzce Province in Turkey. Furthermore, the liquefaction potential indexes (LPI) for both aforementioned methods were calculated using the procedure of Iwasaki et al. [23].
2. Study Area
Duzce Province is located in northwestern Turkey (Figure 1). It is under the effect of the North Anatolian Fault Zone (NAFZ) and is about 30 km distant from the Black Sea. The provincial capital of Duzce is situated on an alluvial soil site.
The study area (Duzce Province).
2.1. Geomorphological and Geological Setting
The study area is situated in an active seismic earthquake zone [24]. Duzce has been affected by the active faults. The 1957 Bolu (M=7) and 1967 Adapazarı (M=7.1) earthquakes occurred on the Bolu-Abant Dokurcun segments of the NAFZ. The active and probable active faults of Duzce, Hendek, and Çilimli are in close proximity to the study area [25] (Figure 2).
View of the Duzce fault [24].
During the (M=7.4) earthquake of 1999, the 30 km eastern segment of the 130 km fault rupture occurred on the western part of the Duzce Fault reaching to Efteni Lake [25]. Duzce plain is an active subsidence and deposition area controlled by lateral strike slip faults surrounded by pre-Quaterney-aged rocks. The oldest is the Yığılca Unit of Eocene-aged caycuma. Formations include volcanic sandstone, rottenstone, andestic and basaltic lavas and volcanic breccia [26]. Quaternary-aged fan, deltaic, and marsh type deposits cover this unit, which consists of gravel, sand, silt, and clay material (Figure 3).
Geological map of the study area and some borehole locations [26].
Due to the elevation of the surrounding rocks and as a result of the basin drainage, deposition occurred mostly in the Duzce area and its surroundings (Figure 4). Little Melen (Kucuk Melen) River discharges into the lake and continues flowing towards the only discharge from the lake, the Big Melen (Buyuk Melen) River. In addition, Aksu, Ugursuyu and Aydınpınar Creeks deposit alluvial fans and join in the lake basin. The surrounding rocks are extremely weathered by eroding, allowing the increase in sedimentation. The thickness of the sedimentation is about 260–300 m. The deposition areas have displaced laterally and the horizontal stratigraphy has been altered [25] (Figure 4).
View of the deposits in the Duzce area.
2.2. Seismotectonics of the Study Area
Duzce plain is a pull apart type basin that is controlled by the lateral strike slip fault system in the NAFZ [28] (Figure 5). Paleo- and neotectonic period active faults exist at the north and south of the plain. There are several faults which are parallel and oblique to these major faults. During the 12 November 1999’s earthquake, the surface rupture ranged through Golyaka on the south towards Kaynaşlı on the east, ending in the Asarsu Valley and the Bolu Tunnel. The city of Duzce is situated in the middle of the plain on a pressure ridge type hill, probably tectonically controlled. Major earthquake records are given in Table 1. Historical earthquakes have been recorded on the Abant-Bayramoren segment in the south. There were 12 earthquakes between 1967 and 1890. The great earthquake of 17 August, 1668 (MS=8), caused a disaster in Anatolia [29], with aftershocks continuing for 6 months [30]. The Bolu-Gerede earthquake (Mw=7.3) on January 2, 1944, was a major one, recorded after the implementation of instruments for scientific measurement of magnitude. It was noted that 2,381 people died and 50,000 houses were damaged [31]. Although the 17 August, 1999, Marmara and 12 November, 1999, Duzce earthquakes occurred on the western segment of the North Anatolian Fault, the measured average value of the horizontal ground acceleration was 0.54 g in Duzce [32].
Magnitude and damage records of earthquakes around the study area [24].
Locations
Date
Epicenter
Mw
Total damage on structures
Death
Injury
Murefte
09.08.1912
40.60–27.20
7.3
5.540
216
466
Hendek
20.06.1943
40.85–30.51
6.6
No
336
No
Gerede
01.02.1944
41.41–32.69
7.2
20.865
3.959
No
Duzce
10.02.1944
41.00–32.30
5.4
900
No
No
Mudurnu
05.04.1944
40.84–31.12
5.6
900
30
No
Yenice
18.03.1953
39.99–27.36
7.4
9.670
265
336
Abant
26.05.1957
40.60–31.20
7.1
4.201
52
100
Çınarcık
18.09.1963
40.77–29.12
6.3
230
1
26
Adapazarı
22.07.1967
40.60–30.89
7.2
5.569
89
235
Gelibolu
27.03.1975
40.45–26.12
6.4
980
7
No
Golcuk
13.09.1999
40.80–30.03
5.7
No
No
Unknown
Duzce
12.11.1999
40.79–31.21
7.2
15.389
845
4.948
Bolu
17.11.1999
40.83–31.51
5.0
No
No
No
Bolu
22.03.2000
40.94–31.58
5.4
No
No
No
Yigilca
26.08.2001
40.93–31.53
5.1
No
No
No
Epicenters of the earthquakes and historical earthquakes in the region [25].
3. Materials and Methods3.1. Field Investigations
Geotechnical bore holes were drilled at 40 locations. The depths of the boreholes ranged from 10 to 30 m, which measured to a total of 296 m. These boreholes were utilized to determine the consistency of fine-grained soils and the stiffness of the coarse soils, to obtain undisturbed and disturbed samples and to measure the groundwater level. A Standard Penetration Test (SPT) [33] was carried out during the drillings and SPT N blow counts were obtained in the boreholes. In situ unit weight and moisture content values were obtained from the trial pits. Then, representative soil samples were obtained in order to determine the geomechanical properties of the soils. The 296 m-thick alluvium was very heterogeneous and included confined and umconfined aquifers. The groundwater level was mostly at the surface and ranged between 1.5 and 3.9 m.
In this study, the shear wave velocity (VS) measurements were based on Andrus et al. [27] process for assessing liquefaction potential; VS values were calculated using empirical equations between shear wave velocity and SPT blow count (N) for all types as follows [34]. The VS values based on SPT blow count (N) were given below (Table 2):
(1)VS=61·N0.5,(2)VS=97·N0.314,(3)VS=76·N0.33,(4)VS=121·N0.27,(5)VS=22·N0.85.
Field data of Duzce Province.
Boreholes number
SPT-NDepth (m)
SPT-N
Water Level
Shear Wave Velocity
Equation (1) (61·N0.5)
Equation (2) (97·N0.314)
Equation (3) (76·N0.33)
Equation (4) (121·N0.27)
Equation (5) (22·N0.85)
BH1
3
32
3.5
345.06
287.99
238.5
308.45
418.7
BH1
9
32
3.5
435.06
287.99
238.5
308.45
418.7
BH2
3
25
3
305
265.5
219.9
288.6
339.4
BH2
6
15
3
236.3
227.02
185.75
251.4
219.9
BH2
15
45
3
409.2
320.5
266.9
338.2
559.3
BH3
4.5
32
4.00
345.06
227.02
185.75
251.4
219.9
BH3
6.00
32
4.00
345.06
227.02
185.75
251.4
219.9
BH3
7.5
28
4.00
322.8
276.17
228.4
297.52
373.7
BH3
12
31
4.00
339.7
285.14
236.02
305.9
407.5
BH4
3.00
31
2.5
339.7
285.14
236.02
305.9
407.5
BH4
6.00
13
2.5
219.95
217.01
177.2
241.9
194.7
BH4
7.5
16
2.5
244
231.7
189.8
255.8
232.3
BH4
9.00
13
2.5
219.94
217.1
177.2
242
194.6
BH4
17.00
13
2.5
219.94
217.1
177.2
242
194.6
BH5
No
No
No
No
No
No
No
No
BH6
3.00
28
2.5
322.8
276.2
228.24
297.52
373.7
BH6
4.5
23
2.5
292.54
260
213.9
282.2
316.2
BH6
7.5
21
2.5
279.6
252.3
207.56
275.3
292.7
BH6
12.00
32
2.5
345.06
252.3
207.56
275.3
292.7
BH7
No
No
No
No
No
No
No
No
BH8
9.00
11
4.0
202.3
206
167.7
231.2
168.9
BH8
10.50
26
4.0
311.1
269.9
222.7
292
350.9
BH8
13.500
28
4.0
322.8
276.2
228.24
297.52
373.7
BH9
9.00
28
4.00
322.8
276.2
228.24
297.52
373.7
BH9
16.5
81
4.00
549
385.6
324.1
396.4
922
BH10
No
No
No
No
No
No
No
No
BH11
7.5
14
3.5
228.3
222.2
182
246.7
207.4
BH11
10.5
40
3.5
386
308.9
256.7
327.6
506.1
BH12
4.5
32
3.25
345.1
252.3
207.56
275.3
292.7
BH12
6.0
36
3.25
366
298.9
248
318.5
463
BH12
10.5
42
3.25
395.4
313.7
260.9
331.9
527.8
BH13
No
No
No
No
No
No
No
No
BH14
3.0
15
3.5
236.3
227.02
185.8
251.4
219.9
BH14
4.5
37
3.5
372
302
250.1
320.8
473.6
BH14
6
26
3.5
132.6
269.9
222.8
291.7
350.9
BH14
7.5
30
3.5
334.5
282.3
233.5
303.2
396.7
BH15
No
No
No
No
No
No
No
No
BH16
4.5
22
2.5
286.2
256.1
210.8
278.8
304.5
BH17
3.0
13
3.5
219.94
217.1
177.2
242
194.6
BH17
4.5
10
3.5
192.9
199.9
162.5
225.4
155.8
BH17
6.0
21
3.5
279.6
252.3
207.56
275.3
292.7
BH17
7.0
26
3.5
311.1
269.9
222.7
292
350.9
BH18
15
34
4.5
355.7
293.6
243.4
313.6
440.1
BH19
3.0
13
3
219.94
217.1
177.2
242
194.6
BH19
4.5
18
3
258.9
240.1
197.3
264.1
256.7
BH19
7.5
25
3
305
266.6
219.9
288.6
339.4
BH19
13.5
27
3
317
273
255.6
294.9
362.4
BH20
No
No
No
No
No
No
No
No
BH21
No
No
No
No
No
No
No
No
BH22
No
No
No
No
No
No
No
No
BH23
4.5
29
2.0
328.5
279.3
230.8
300.35
385.1
BH23
9.0
8
2.0
172.6
186.4
8.96
212.2
128.9
BH23
10.5
10
2.0
192.9
199.9
162.5
225.4
155.8
BH23
13.5
9
2.0
183
194
157
219
142.5
BH23
15.0
12
2.0
211.4
211.7
172.6
236.7
181.9
BH24
7.5
14
3.5
228.3
222.2
182
246.7
207.4
BH24
10.5
40
3.5
385.8
308.9
256.8
327.6
506.1
BH25
3.0
44
4.5
404.7
318.3
264.9
336.2
548.8
BH25
4.5
14
4.5
228.3
222.2
182
246.7
207.4
BH25
7.5
9
4.5
183
194
157
219
142.5
BH25
10.5
19
4.5
265.9
244.6
200.9
267.9
268.9
BH25
12.0
22
4.5
286.2
256.1
210.8
278.8
304.5
BH26
6.0
13
2.0
219.94
217.1
177.2
242
194.6
BH26
9.0
32
2.0
345.1
252.3
207.56
275.3
292.7
BH27
NO
No
No
No
No
No
No
No
BH28
No
No
No
No
No
No
No
No
BH29
7.5
27
3.5
317
273
255.6
294.9
362.4
BH29
10.5
23
3.5
292.54
260
213.9
282.2
316.2
BH29
12.0
50
3.5
431.2
331.2
276.4
347.95
612
BH30
No
No
No
No
No
No
No
No
BH31
No
No
No
No
No
No
No
No
BH32
4.5
7.0
3.0
161.3
178.8
144.5
204.7
115.1
BH32
6.0
32
3.0
345.1
252.3
207.56
275.3
292.7
BH32
9.0
40
3.0
385.8
308.9
256.8
327.6
506.1
BH32
12
54
3.0
448.3
339.5
283.5
355.3
653.1
BH33
No
No
No
No
No
No
No
No
BH34
3.0
20
3.0
272.9
248.5
204.3
271.7
280.74
BH34
4.5
21
3.0
279.6
252.3
207.56
275.3
292.7
BH35
3.0
9
2.5
183
194
157
219
142.5
BH35
4.5
8
2.5
172.6
186.4
150.9
212.2
128.9
BH35
7.5
10
2.5
192.9
199.9
162.5
225.4
155.8
BH35
12
12
2.5
211.4
211.7
172.6
236.7
181.9
BH36
3
16
3.5
244
232
189.8
255.8
232.3
BH36
4.5
16
3.5
244
232
189.8
255.8
232.3
BH37
No
No
No
No
No
No
No
No
BH38
No
No
No
No
No
No
No
No
BH39
No
No
No
No
No
No
No
No
BH40
3
20
2.25
272.9
248.5
204.3
271.7
280.74
BH40
4.5
49
2.25
427
329.3
274.5
346.1
602
BH40
9
16
2.25
244
232
189.8
255.8
232.3
3.2. Calculation of Seismic Hazard Design Parameters
The Duzce Fault Zone is situated 13 km south of the study area, the North Anatolian Fault Zone is located 73 km south of the study area, and the Hendek Fault is found 29 km north-west of the study field (Figure 6). The fault zone with the highest possible acceleration in the study site is the North Anatolian Fault Zone. A circle with a radius of 100 km was drawn around the study area in order to identify the seismic design parameters. Within this circle, active seismic sources thought to affect the study field were vertically connected to calculate the shortest routes to the study field in km (Figure 6). These investigations and measurements showed that there were three main fault zones inside the circle. Then, the horizontal flying distances to the study field were calculated as 13 km for the Duzce Fault, 29 km for the Hendek Fault, and 73 km for the North Anatolian Fault [35].
Identification of seismic sources within 100 km radius of the study area.
The map of Turkey’s active faults published by the Mineral Research and Exploration Institute indicates the total length of the Duzce, Hendek, and North Anatolian Faults as 85 km, 60 km, and 200 km, respectively [35]. The DuzceFault, which is the shortest distance to the study area and has the potential to produce an earthquake, was taken into consideration in the main objective of the study when estimating the next earthquake expected to occur.
According to Mark [39], it is assumed that 1/3 of this fault zone could be ruptured. Therefore, the moment size of the probable seismic design was calculated by using the equation of Wells and Coppersmith [40], as seen below (6):
(6)M=4.86+1.32LogL,
where M is moment magnitude and L is fault length (km).
According to this approach, moment magnitude was calculated to be 7.2 in the case of a rupture of 1/3 of the fault length.
Horizontal earthquake acceleration (peak ground acceleration, PGA) was calculated by using the attenuation relationship (7) developed for faults based on the earthquakes in Turkey [41]:
(7)PGA=2.18e0,0218(33,3Mw-Re+7.8427SA+18.9282SB),
where SA=0 and SB=1 values are used for soft soils, Re is the shortest horizontal flying length to the respective fault zone from the settlement, and Mw is the magnitude of the earthquake. The peak horizontal earthquake acceleration that can be created by the seismic design was found to be 0.28 g.
4. Assessment of Liquefaction Potential
Prediction of the liquefaction potential of soils is based on cyclic laboratory testing on soil samples and use of in-situ tests and empirical methods. However, the use of laboratory testing is complicated due to difficulties associated with sample disturbance during both sampling and reconsolidation. Thus, empirical approaches based on the in-situ penetration test results have gained popularity in engineering practice as well as in engineering codes [42].
In this study, after obtaining the in situ test results, the evaluation of the liquefaction procedure was begun. The evaluation procedures based on the Standard Penetration Test (SPT) [43] and measurement of shear wave velocity (VS) [27] require the measurement of three parameters: (1) the level of cyclic loading on the soil caused by the earthquake, expressed as an acyclic stress ratio (CSR); (2) the stiffness of the soil, expressed as over burden stress (corrected SPT blow count) due to shear wave velocity; and (3) the resistance of the soil to liquefaction, expressed as a cyclic resistance ratio (CRR). Guidelines for calculating each parameter are presented below [44].
4.1. Cyclic Stress Ratio (CSR)
The cyclic stress ratio (CSR) characterizes the seismic demand induced by a given earthquake, and it can be determined from peak ground surface acceleration that depends upon site-specific ground motions [45]. The expression for the CSR induced by earthquake ground motions formulated by Idriss and Boulanger [46] is as follows:
(8)CSR=0.65amaxgσVσV′rd1MSF1Kσ,
where 0.65 is a weighing factor to calculate the equivalent uniform stress cycles required to generate the same pore water pressure during an earthquake; the amax is the peak horizontal ground acceleration; g is the acceleration of gravity; σV and σV′ are total vertical overburden stress and effective vertical overburden stress, respectively, at a given depth below the ground surface; rd is the depth-dependent stress reduction factor; MSF is the magnitude scaling factor; and Kσ is the overburden correction factor.
The stress reduction factor (rd) accounts for the dynamic response of the soil column and represents the variation of shear stress amplitude with depth. Idriss and Boulanger [46] formulated the following expressions to calculate the stress reduction factor (rd) (9)–(11):
(9)rd=exp[α(z)+β(z)Mw],(10)α(z)=-1.012-1.126sin(z11.73+5.133),(11)β(z)=0.106+0.118sin(z11.28+5.142),
where z is the depth (m) and Mw is the moment magnitude. The arguments inside the sine terms in (10) and (11) are in radians. The above expression for rd is valid up to a depth of z≤34 m, and the depths of the boreholes considered in the present analysis were less than 34 m.
The values of CSR that pertain to the equivalent uniform shear stress induced by an earthquake of magnitude, Mw=7.5, were adjusted to an equivalent CSR for an earthquake of magnitude Mw=7.5 through the introduction of the magnitude scaling factor (MSF), which accounts for the duration effect of ground motions. The MSF for Mw<7.5 is expressed as follows (12):
(12)MSF=6.9exp(-Mw4)-0.058≤1.8.
Since the liquefaction resistance increases with increasing confining stress, the overburden correction factor (Kσ) was applied such that the values of CSR were adjusted to an equivalent overburden pressure σV′ of 1 atmosphere equations (13)-(14):
(13)Kσ=1-Cσ(lnσv)≤1.0,
where
(14)Cσ=118.9-2.5507(N1)60≤0.3Pa is the atmospheric pressure (= 00 kPa).
4.2. Corrected SPT Blowcount and Shear Wave Velocity
In this study, the measured SPT N values (N) were corrected for overburden stress, energy ratio, diameter of boreholes, length of sampling rod, and the type of sampler by introducing a series of correction factors. N60 is the corrected Nm value for a 60% energy ratio with an assumption that 60% of the energy was transferred from the falling hammer to the SPT sampler. The corrected (N1)60 values were calculated as follows (15):
(15)(N1)60=NmCNCECBCRCR,
where CN is a factor to normalize Nm to a common reference effective overburden stress; CE is the correction for the hammer energy ratio (ER); CB is the correction factor for borehole diameter; CR is the correction factor for rod length; and CS is the correction for samplers with or without liners. The value of CN was calculated as per (15) and was limited to a maximum value of 1.7. CS, CB, and CE were assumed to be 1.1, 1.0, and 0.6, respectively. Rod length correction with respect to the depth (CR) at each borehole location was corrected as shown in Table 3, as suggested by Youd and Idriss [47].
The rod length correction with respect to the depth.
Depth
Correction for rod length
d
CR
d<3 m
0.75
d = 3-4 m
0.8
d = 4–6 m
0.85
d = 6–10 m
0.95
d = 10–30 m
1.0 m
The overburden correction (CN) factor to normalize (N1)60 to a common reference effective overburden stress is as follows (16)-(17):
(16)CN=(PaσV′)α≤1.7,
where
(17)α=0.784-0.0768√(N1)60.
It can be observed from (16) and (17) that (N1)60 and CN are interdependent. A series of iterations were carried out to determine (N1)60 and CN until the difference between successive iteration values was less than 0.001.
In addition, shear wave velocities had to be corrected. In the procedure of liquefaction potential evaluation proposed by Andrus et al. [27], shear wave velocity was corrected to overburden stress and (18) was suggested:
(18)VS1=VS1(PaσV′)0.25(0.5KO′)0.125,
where VS is the shear wave velocity (m/s); VS1 is the stress-corrected shear wave velocity (m/s); Pa is the atmosphere pressure equal to 100 kPa; σV′ shows the the effective overburden stress; and KO′ is the coefficient of effective earth pressure (in this study assumed equal to 0.5) [44].
4.3. Evaluation of the Cyclic Resistance Ratio (CRR)
Determination of the cyclic resistance ratio (CRR) requires fines content (FC) of the soil to correct updated SPT blow count (N1)60 to an equivalent clean sand standard penetration resistance value (N1)60cs. Idriss and Boulanger [46] determined the CRR value for cohesionless soil with any fines content using the following expression (19)–(21):
(19)CRR=exp{(N1)60cs14.1+((N1)60cs126)2-((N1)60cs23.6)3kkkkkk+((N1)60cs25.4)3-2.81},(20)(N1)60cs=(N1)60+Δ(N1)60,
where Δ(N1)60 is the correction for fines content in percent (FC) present in the soil and is expressed as
(21)Δ(N1)60=exp(1.63+9.7FC+0.1-(15.7FC+0.1)2).
Separately, as regarding the use of VS as an index of liquefaction resistance, it has been illustrated by several authors. The most popular CRR-VS correlation (Figure 7) was proposed by Andrus and Stokoe [22] for uncemented Holocene-age soils, based on a database including 26 earthquakes and more than 70 test sites. The CRR is obtained as a function of an overburden-stress corrected shear wave velocity VS1=VS(Pa/σV0′)0.25, where VS= measured shear wave velocity, Pa= atmospheric pressure (≈100kpa), and σV0′= initial effective vertical stress (same units as Pa). Andrus et al. [27] introduced age correction factors to extend the original correlation of Andrus and Stokoe [22] to soils older than Holocene. Their CRR-VS1 relationship (curves in Figure 7, for various fines contents) is approximated by:
(22)CRR=[0.022(Ka1VS1100)2+2.8(1VS1*-Ka1Vs1-1VS1*)]×Ka2,
where VS1* = the limiting upper value of VS1 for liquefaction occurrence (VS1*= 200 m/s for the curve for fines content ≥35%); VS1*= 215 m/s for the curve for fines content ≤5%; VS1* varies linearly from 200 to 215 m/s for fines content between 35 and 5%; Ka1= factor to correct for high VS1 values caused by aging; and Ka2= Factor to correct for influence of age on CRR. Magnitude scaling factors should be used to scale (22) (for magnitude Mw=7.5 earthquakes) to different magnitudes. Both Ka1 and Ka2 are 1 for uncemented soils of Holocene age. For older soils, suggested Ka1 values (mostly in the range 0.6 to 0.8) are derived from SPT-VS1 relationships (e.g. Ohta and Goto [48], Rollins et al. [49], or site specific). Lower-bound values of Ka2 (1.1 to 1.5) are based on the study by Arango et al. [50]. However, Andrus et al. [27] noted the associated high uncertainty and the need for additional work to quantify the influence of age on CRR, as well as on Vs.
Recommended curves for evaluating CRR from shear wave velocity VS for clean, uncemented soils with liquefaction data from compiled case histories [27].
4.3.1. Determination of the Factor of Safety
The factor of safety against liquefaction (FS) is commonly used to quantify liquefaction potential. The factor of safety against liquefaction (FS) can be defined as follows:
(23)FS=(CRR)Mw=7.5(CSR)Mw=7.5,σV′MSF.
Both CSR and CRR vary with depth and, therefore, the liquefaction potential is evaluated at corresponding depths within the soil profile.
4.3.2. Determination of the Liquefaction Potential Index
The liquefaction potential index (LPI) is a single-valued parameter to evaluate regional liquefaction potential. The LPI at a site is computed by integrating the factors of safety (FS) along the soil column up to a depth of 20 m. A weighting function is added to give more weight to the layers closer to the ground surface. The liquefaction potential index (LPI) proposed by Iwasaki et al. [36, 51] is expressed as follows (24). The criteria of the level of liquefaction severity indexes were given below (Table 4):
(24)LPI=∫020F(z)W(z)dz,
where z is the depth of the midpoint of the soil layer (0 to 20 m) and dz is the differential increment of depth. The weighting factor, W(z), and the severity factor, F(z), are calculated as per the following expressions (25):
(25)F(z)=1-FSforFS<1.0,F(z)=0forFS≥1.0,W(z)=10-0.5zforz<20m,W(z)=0forz>20m.
For the soil profiles with depths of less than 20 m. The LPI was calculated using the following expression [37] (26)-(27):
(26)LPI=∑i=1nWiFiHi
with
(27)Fi=1-FSiforFSi<1.0,Fi=0forFSi≥1.0,
where Hi is the thickness of the discretized soil layers; n is the number of layers; Fi is liquefaction severity for ith layer; FSi is the factor of safety for ith layer; Wi is the the weighting factor (=10-0.5Zi); and Zi is the depth of ith layer (m).
The level of liquefaction severity.
LPI
Iwasaki et al. [23]
Luna and Frost [37]
MERM [38]
LPI = 0
Very low
Little to none
None
0 < LPI < 5
Low
Minor
Low
5 < LPI < 15
High
Moderate
Medium
15 < LPI
Very high
Major
High
5. Assessment of the Liquefaction Potential Index
The city of Duzce has been reconstructed since the 12 November earthquake of 1999. The general form of construction had typically been a 4-5-storey reinforced-concrete frame and masonry structure. After the 12 November Duzce earthquake experience, regulations were changed to limit construction to 2-3-storey buildings. The city is located over deep alluvial deposits. The main soils deposited at this site are comprised of alluvial sand and silt. The boreholes in Duzce were drilled around the Efteni Lake at a depth of 200 m and did not reach bedrock. The shallow soils at approximately 10 m are recent deposits laid down by the Aksu and Melen Rivers that flooded the area.
Turkey is located in the active tectonic region of the Alp Himalayan Earthquake Zone, so this area is an active region seismologically. There are several active tectonic sections in Turkey, such as the North Anatolian Fault Zone, the East Anatolian Fault Zone, the West Anatolian Grabens, the Ecemiş Fault Zone, and the Tuzgolu Fault Zone [52]. Duzce Province is located near the Duzce Fault segment of the North Anatolian Fault Zone which is active in the Western Black Sea Region. Furthermore, this area consists of granular alluvial deposits which are loose to the surface. The groundwater is between 2.5 and 4 m below the surface and changes seasonally. For the analysis of the liquefaction potential index of Duzce Province, a total of 40 geotechnical boreholes were drilled by the General Directorate of Mineral Research and Exploration. The field data of the works were assessed for the liquefaction potential index for Duzce Province. The SPT samples were implemented at depth intervals of 1.5 m from the first to the last of the boreholes, and the disturbed samples were used to describe the grain size distribution and Atterberg limits of the soils. The boundaries of the soil layers, SPT-N values, fines content, and the liquid limit for all layers throughout the boreholes were employed as input parameters to determine the liquefaction potential index.
In addition, the magnitude of the earthquake and the maximum horizontal acceleration of those parameters to be created due to local faults were used here for evaluating the liquefaction potential index. The Duzce Fault Zone of the North Anatolian Fault Zone and surrounding zones were generated and showed an average of 7.2 moment magnitudes. For this reason, the magnitude of the projected earthquake was found by using 7.2 for the calculations. In this context, the typical computation of factors of safety against liquefaction for earthquakes (Mw=7.2) yielded by the Duzce Fault Zone was carried out at the chosen borehole using (2) through (20). The LPI at this particular site was calculated from the FS values based on the expression by Luna and Frost [37]. The LPI values were computed at the study site for magnitudes of Mw=7.2.
Great effort was taken in the analysis of the other input parameter for determining the liquefaction potential, the maximum ground acceleration (amax). However, some researchers have offered empirical equations for the maximum ground acceleration [41, 53, 54]. In particular, the comprehensive study of Ulusay et al. [41] should be mentioned as it relates to the iso-acceleration map of Turkey. In this study, the amax values were calculated as approximately 502 gal for the Duzce Fault Zone segment. The liquefaction potential index indices for 40 boreholes were calculated and are given in Table 5 and Figures 8, 9, 10, and 11. In addition, the distributions of liquefaction potential indexes are presented in Figure 12 as a pie chart.
Liquefaction potential index results.
Boreholes number
The level of liquefaction severity (LPI)
Equation (1)
Equation (2)
Equation (3)
Equation (4)
Equation (5)
BH1
0
0
0
0
0
BH2
0
0
14.14
19.79
0
BH3
0
0
35.63
0
0
BH4
30.32
13.35
39.17
0
28.97
BH5
Not performed
0
Not performed
Not performed
Not performed
BH6
0
0
23.09
0
0
BH7
Not performed
Not performed
Not performed
Not performed
Not performed
BH8
46.83
33.46
26.34
5.82
40.87
BH9
0
0
0
0
0
BH10
Not performed
Not performed
Not performed
Not performed
Not performed
BH11
5.43
0
0.8
0
24.39
BH12
0
0
0
0
0
BH13
Not performed
Not performed
Not performed
Not performed
Not performed
BH14
0
9.14
0
0
0
BH15
Not performed
Not performed
Not performed
Not performed
Not performed
BH16
0
0
0
0
0
BH17
10.5
5.37
0
0
8.90
BH18
0
0
0
0
0
BH19
0
0
24.17
0
0
BH20
Not performed
Not performed
Not performed
Not performed
Not performed
BH21
Not performed
Not performed
Not performed
Not performed
Not performed
BH22
Not performed
Not performed
Not performed
Not performed
Not performed
BH23
23.86
36.96
35.26
42.53
40.84
BH24
1.5
0
19.89
35.34
26.39
BH25
17.85
14.52
14.68
8.48
19.94
BH26
5.68
0
41.81
0
22.44
BH27
Not performed
Not performed
Not performed
Not performed
Not performed
BH28
Not performed
Not performed
Not performed
Not performed
Not performed
BH29
0
8
9.51
0
0
BH30
Not performed
Not performed
Not performed
Not performed
Not performed
BH31
Not performed
Not performed
Not performed
Not performed
Not performed
BH32
33.05
22.94
32.74
0
30.74
BH33
Not performed
Not performed
Not performed
Not performed
Not performed
BH34
0
0
0
0
0
BH35
66.53
31.64
56.95
2.84
59.78
BH36
0
0
9.67
0
0
BH37
Not performed
Not performed
Not performed
Not performed
Not performed
BH38
Not performed
Not performed
Not performed
Not performed
Not performed
BH39
Not performed
Not performed
Not performed
Not performed
Not performed
BH40
0
0
18.88
0
0
Borehole 23 was chosen to check and compare for liquefaction potential as a sample according to shear wave velocity values.
Borehole 23 was chosen to check and compare for liquefaction potential as a sample for shear wave velocity.
Borehole 23 was chosen to check and compare for liquefaction potential as a sample according to the SPT-N values.
Borehole 23 was chosen to check and compare for liquefaction potential as scattering liquefaction points of sample for SPT-N values.
Pie charts showing the areas of the potential zones.
6. Results and Conclusions
The evaluation of the liquefaction potential of a liquefaction-prone area is of vital importance in geotechnical earthquake engineering, both for assessment for site selection and for planning and construction. This study investigated two field methods used to evaluate the liquefaction potential of soils, the Standard Penetration Test (SPT) and the Shear Wave Velocity Test (VS), based on the empirical relationships between them. Attempts were made to evaluate the factors of safety against liquefaction (FS) and corresponding liquefaction potential indices (LPI) for a local fault zone in order to produce the seismic movement for the province using SPT-N-based semiempirical procedures.
The concept of the liquefaction potential index was used in this study for liquefaction susceptibility, as proposed by Iwasaki et al. [36]. The distribution of the LPI was generated in order to predict the occurrence of damaging liquefaction for an earthquake to be yielded by the local fault zone in Duzce Province in the Western Black Sea Region of Turkey. This study area is under the effect of the North Anatolian Fault Zone through its segment, the Duzce Fault Zone, which was evaluated for producing the liquefaction potential indices by calculating for a probable earthquake of Mw=7.2.
The comparison of the safety factors and liquefaction potential indexes reveal that the severity of liquefaction occurrences in the study area based on the VS methods (Equation (1) = 43.86, equation (2) = 40.84, equation (3) = 42.53, equation (4) = 36.96, equation (5) = 43.86) are bigger than the one based on the SPT method (35,36). Moreover, it can be observed that the relationships between the SPT method and the shear wave velocity are not suitable. Because the relationships used in the present study are dependent on soil type, fines content, type of tests, and their accuracy, it might be more valid to perform both methods for the same place and then compare the results in order to evaluate the liquefaction potential.
Finally, a very high susceptibility category of liquefaction was observed for the potential earthquake of Mw=7.2; however, 3.8–10.2% of the study area is in the highly susceptible liquefaction class in five distribution charts according to (1)–(5). The percentage that is moderately susceptible takes up the least area from the other class: 1.2–4.1% for all locations in the distribution charts. The low susceptibility areas are 28.76–65%, respectively.
In conclusion, the areas developed on reclaimed land having large, thick deposits of soft soil and shallow groundwater levels were observed to be more prone to liquefaction. This paper reveals that some of the areas are more highly prone to liquefaction due to the greater thickness of the soft soil deposits and groundwater table at shallow depths. It can be observed from the distribution of the LPI that a high degree of liquefaction would occur at several sites in the Province of Duzce during a seismic event. These LPI distributions will help the structural designers and city planners to check the vulnerability of the area against liquefaction.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to express their gratitude to the Duzce Provincial Governorate and the General Directorate of Mineral Research and Exploration for information and the logs of geotechnical boreholes.
SeedH. B.IdrissI. M.Simplified procedure for evaluating soil liquefaction potential1971979124912732-s2.0-0015128601YoudT. L.IdrissI. M.AndrusR. D.ArangoI.CastroG.ChristianJ. T.DobryR.Liam FinnW. D.Harder L.F.J.HynesM. E.IshiharaK.KoesterJ. P.LiaoS. S. C.MarcusonW. F.IIIMartinG. R.MitchellJ. K.MoriwakiY.PowerM. S.RobertsonP. K.SeedR. B.StokoeK. H.IIILiquefaction resistance of soils: summary report from the 1996 NCEER and 1998 NCEER/NSF workshops on evaluation of liquefaction resistance of soils2001127108178332-s2.0-003547467110.1061/(ASCE)1090-0241(2001)127:10(817)SeedH. B.IdrisI. M.ArangoI.Evaluation of liquefaction potential using field performance data198310934584822-s2.0-002097158810.1061/(ASCE)0733-9410(1983)109:3(458)SeedH. B.TokimatsuK.HarderL. F.ChungR. M.The influence of SPT procedures in soil liquefaction resistance evaluation198511112142514452-s2.0-002236008810.1061/(ASCE)0733-9410(1985)111:12(1425)HardinB. O.DrnevichV. P.Shear modulus and damping in soils: design equation and curves19729876676922-s2.0-0015360542KramerS.1996Upper Saddle River, NJ, USAPrentice HallIshiharaK.1996Oxford, UKOxford University PressOxford Engineering Science SeriesRauchA. F.DuffyM.StokoeK. H.IILaboratory correlation of liquefaction resistance with shear wave velocity110Proceedings of the Sessions of Computer Simulation of Earthquake Effects (Geo-Denver '00)August 2000Denver, Colo, USA66802-s2.0-003372573310.1061/40523(298)4HuangY.-T.HuangA.-B.KuoY.-C.TsaiM.-D.A laboratory study on the undrained strength of a silty sand from Central Western Taiwan2004249-107337432-s2.0-424318274110.1016/j.soildyn.2004.06.013YunminC.HanK.Ren-PengC.Correlation of shear wave velocity with liquefaction resistance based on laboratory tests20052564614692-s2.0-2294448698210.1016/j.soildyn.2005.03.003ZhouY.-G.ChenY.-M.HuangB.Experimental study of seismic cyclic loading effects on small strain shear modulus of saturated sands2005632292362-s2.0-1474430010710.1631/jzus.2005.A0229LiuN.MitchellJ. K.HonM.Influence of nonplastic fines on shear wave velocity-based assessment of liquefaction20061328109110972-s2.0-3374617867810.1061/(ASCE)1090-0241(2006)132:8(1091)ZhouY.-G.ChenY.-M.Laboratory investigation on assessing liquefaction resistance of sandy soils by shear wave velocity200713389599722-s2.0-3454713209210.1061/(ASCE)1090-0241(2007)133:8(959)BaxterC. D. P.BradshawA. S.GreenR. A.WangJ.-H.Correlation between cyclic resistance and shear-wave velocity for providence silts2008134137462-s2.0-3754905507510.1061/(ASCE)1090-0241(2008)134:1(37)AskariF.DabiriR.ShafieeA.JafariM. K.Liquefaction resistance of sand-silt mixtures using laboratory-based shear wave velocity2011921351442-s2.0-79958828624AndrusR. D.1994Austin, Tex, USAUniversity Of TexasDabiriR.AskariF.ShafieeA.JafariM. K.Shear wave velocity-based liquefaction resistance of sand-silt mixtures: deterministic versus probabilistic approach20113521992152-s2.0-80755176866LodgeA. L.1996Berkeley, Calif, USAUniversity of CaliforniaKayabaliK.Soil liquefaction evaluation using shear wave velocity1996441–41211272-s2.0-0030265310AndrusR. D.StokoeK. H.YoudT. L.IdrissI. M.Liquefaction resistance based on shear wave velocityProceedings of the NCEER Workshop on Evaluation of Liquefaction Resistance of SoilsJanuary 1997 1996Salt lake City, Utah, USANational Center for Earthquake Engineering Research, Buffalo, NY, USA89128Technical Report NCEER-97-0022JuangC. H.ChenC. J.TangW. H.RosowskyD. V.CPT-based liquefaction analysis, part 1: determination of limit state function20005055835922-s2.0-0034302748AndrusR. D.StokoeK. H.IILiquefaction resistance of soils from shear-wave velocity200012611101510252-s2.0-003432684210.1061/(ASCE)1090-0241(2000)126:11(1015)IwasakiT.TokidaK.TatsuokaF.WatanabeS.YasudaS.SatoH.Microzonation for soil liquefaction potential using simplified methodsProceedings of the 3rd International Conference on MicrozonationJune-July 1982Seattle, Wash, USA13191330NEMCB.U. Kandıllı Observatory and Earthquake Research Instıtute, Natıonal Earthquake Monıtorıng Centre, 2014, http://www.koeri.boun.edu.tr/sismoTubitak-AU-MTA1999Ankara, TurkeyEmreO.DumanT. Y.DoğanA.1999Ankara, TurkeyGeneral Directorate of Mineral Research and ExplorationNo. 22AndrusR. D.PiratheepanP.EllisB. S.ZhangJ.JuangC. H.Comparing liquefaction evaluation methods using penetration VS relationships2004249-107137212-s2.0-414305577310.1016/j.soildyn.2004.06.001GoktenE.OzaksoyV.DemirtasR.Some neotectonical properties of the North Anatolian fault zone between bayramo-ren-abant (Turkey)Active Tectonic Research Group First Meeting Proceeding1998ITU Mining Faculty-Eurasia Earth Sciences Institute6877DemirtasR.2000Ankara, TurkeyUniversity of AnkaraAmbraseysN. N.Engineering seismology1988171110510.1002/eqe.4290170101TasmanC.The Report About The Gerede-Bolu EarthquakeDAD/At. 46, 1944 (Turkish)UgrasT.ColakogluZ.CekenU.KuruT.KokyerU.CinatC.TuzelB.BeyhanM.CoruhE.17 Agustos_Izmit Korfezi Depremi Kuvvetli Yer Hareketi Kayıtlarının Değerlendirilmesi17 Ağustos_Izmit Korfezi Deprem Raporu, Bayındırlık Ve Iskan Bakanlığı, Afet Isleri Genel Müdürlügü, Deprem Araştırma Dairesi, Ankara, Pp 156–189, 2000, (Turkish)American Society for TestingmaterialsStandard test for penetration test and split-barrel sampling of soils (ASTM D1596)199404. 09Section 41092JafariM. K.ShafieeA.RazmkhahA.Dynamic properties of fine grained soils in south of tehran2002412536ŞaroğluF.EmreO.KuşçuI.1992Ankara, TurkeyMTA Genel MüdürlüğüIwasakiT.TokidaK.TatsukoF.YasudaS.A practical method for assessing soil liquefaction potential based on case studies at various sites in JapanProceedings of the 2nd International Conference on MicrozonationNovember-December 1978San Francisco, Calif, USA885896LunaR.FrostJ. D.Spatial liquefaction analysis system199812148562-s2.0-0031977546Microzonation for Earthquake Risk Mitigation (MERM)2003MarkR. K.Application of linear statistical model of earthquake magnitude versus fault length in estimating maximum expectable earthquakes19775846446610.1130/0091-7613(1977)5LTHEXA464:AOLSMO>2.0.CO;2WellsD. L.CoppersmithK. J.New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement199484497410022-s2.0-0028555013UlusayR.TuncayE.SonmezH.GokceogluC.An attenuation relationship based on Turkish strong motion data and iso-acceleration map of Turkey2004743-42652912-s2.0-2104443689610.1016/j.enggeo.2004.04.002AnsalA.TonukG.PitilakisK.Source and site effects for microzonation2007Chapter 4Geotechnical Geological and Earthquake Engineering7392Geotechnical Geological and Earthquake EngineeringSeedH. B.TokimatsuK.HarderL. F.ChungR. M.The influence of SPT procedures in soil liquefaction resistance evaluations198511112142514452-s2.0-002236008810.1061/(ASCE)0733-9410(1985)111:12(1425))NoutashM. K.DabiriR.BonabM. H.The evaluation of soil liquefaction potential using shear wave velocity based on empirical relationships201264DixitJ.DewaikarD. M.JangidR. S.Assessment of liquefaction potential index for Mumbai city2012129275927682-s2.0-8486745951610.5194/nhess-12-2759-2012IdrissI. M.BoulangerR. W.Semi-empirical procedures for evaluating liquefaction potential during earthquakes2006262–41151302-s2.0-3244443706510.1016/j.soildyn.2004.11.023YoudT. L.IdrissI. M.Liquefaction resistance of soils: summary report from the 1996 NCEER and 1998 NCEER/NSF workshops on evaluation of liquefaction resistance of soils200112742973132-s2.0-003531151110.1061/(ASCE)1090-0241(2001)127:4(297)OhtaY.GotoN.Empirical shear wave velocity equations in terms of characteristic soil indexes1978621671872-s2.0-0017946399RollinsK. M.DiehlN. B.WeaverT. J.Implications of VS-BPT (N1)60 correlations for liquefaction assessment in gravels19981ASCE506517Geotechnical Special Publication No. 75ArangoI.LewisM. R.KramerC.Updated liquefaction potential analysis eliminates foundation retrofitting of two critical structures2000201–417252-s2.0-003435126710.1016/S0267-7261(00)00034-8IwasakiT.TokidaK.TatsuokaF.WatanabeS.YasudaS.SatoH.Microzonation for soil liquefaction potential using simplified methodsProceedings of the 2nd International Conference on Microzonation1982Seattle, Wash, USA13191330YalcinA.GokceogluC.SönmezH.Liquefaction severity map for Aksaray city center (Central Anatolia, Turkey)2008846416492-s2.0-47049093164JoynerW. B.BooreD. M.Peak horizontal acceleration and velocity from strong-motion records from the 1997 imperial valley, California earthquake198171620112038AydanO.SezaiM.YararT.The seismic characteristics of Turkish earthquakesProceedings of the 11th International Conference on Earthquake Engineering1986Acapulco, Mexico18