Abstract

This study presents a case study on ground response analysis of one of the important cultural heritages in Hanoi, Vietnam. One-dimensional nonlinear and equivalent linear site response analyses which are commonly applied to solve the problem of seismic stress wave propagation are performed at the Ba Dinh square area. A measured in-situ shear wave velocity profile and corresponding geotechnical site investigation and laboratory test data are utilized to develop the site model for site-specific ground response analysis. A suite of earthquake records compatible with Vietnamese Design Code TCVN 9386: 2012 rock design spectrum is used as input ground motions at the bedrock. A few concerns associated with site-specific ground response evaluation are analyzed for both nonlinear and equivalent linear procedures, including shear strains, mobilized shear strength, and peak ground acceleration along with the depth. The results show that the mean maximum shear strains at any soil layer are less than 0.2% in the study area. A deamplification portion within the soil profile is observed at the layer interface with shear wave velocity reversal. The maximum peak ground acceleration (PGA) at the surface is about 0.2 g for equivalent linear analysis and 0.16 g for nonlinear analysis. The ground motions are amplified near the site natural period 0.72 s. The soil factors calculated in this study are 1.95 and 2.07 for nonlinear and equivalent linear analyses, respectively. These values are much different from the current value of 1.15 for site class C in TCVN 9386: 2012. A comparison of calculated response spectra and amplification factors with the local standard code of practice revealed significant discrepancies. It is demonstrated that the TCVN 9386: 2012 soil design spectrum is unable to capture the calculated site amplification in the study area.

1. Introduction

Seismic hazard assessment is needed for seismically active regions and plays a significant role in the sustainable development of urban infrastructure. Ground response analyses are commonly employed for this purpose [1].

Vietnam is located in the center of Southeast Asia and is strongly affected by relative movements between the Mediterranean and Himalaya belts [2]. Hanoi is the capital of Vietnam and lies on two main faults, which are the Red River fault and the Chay River fault. These faults belong to the Red River Shear Zone that comes from Tibet and runs more than 1000 km to the Gulf of Tonkin (located off the coasts of North Vietnam and South China). The Red River fault is also considered as the border of the Indochina Plate and the South China Plate. Recently, the geodynamic research in the Hanoi basin shows that the Red River fault is reactivated [35]. Besides, Hanoi is also affected by nearby active faults such as the Lai Chau-Dien Bien (LC-DB), Ma River, Son La, Lo River, and Chay River faults as shown in Figure 1 [6, 7].

From 114 to 2003, a total of 1645 earthquakes, either by measurement or by studying the historical archives, were recorded in Vietnam with a local magnitude () of 3 or higher [8]. The typical earthquakes that occurred in the northern part of Vietnam nearby the Hanoi are the Dienbien earthquake (1935,  = 6.75, Ma River fault) and the Tuangiao earthquake (1983,  = 6.75, Son La fault) [9]. Recently, the earthquake activity was recorded in the adjacent provinces of the capital, Hanoi. These include Dien Bien in 2001, Ha Giang in 2005, Ninh Binh in 2005, and Hoa Binh in 2010 with a of 5.3, 4.7, 3.1, and 4.0, respectively. They caused vibration and damage to buildings in Hanoi. Furthermore, in the past, three earthquakes were recorded at Hanoi in 1276, 1278, and 1285 with  = 5.5–6.0 [7].

The study of Nhung and Phuong [10] presented that the geology of Hanoi is mainly site classes D, E, F according to EC-8 [11]. The surface is mainly sandy and clay sediments in Holocene or Pleistocene. The bottom of Holocene deposits is at a depth of 10–45 m. The bottom of Pleistocene deposits is at a depth of 45 m (in the north of Hanoi) to 110 m (in the south of Hanoi) [5]. With the seismic and geological conditions, it revealed that a high amplification can occur in Hanoi. A site-specific study is needed to carry out (i.e. seismic hazard analysis), especially at important locations such as the Ba Dinh square area which includes the President’s place and One Pillar Pagoda.

One-dimensional (1D) ground response analysis is commonly performed to evaluate the seismic hazard of the site-specific area [1215]. Shylamoni et al. [12] evaluated ground response at the Kashiwazaki-Kariwa Nuclear Power Plant site (Japan). Mahmood et al. [13] implemented equivalent linear and nonlinear site response analysis of Pashto Cultural Museum Peshawar, Pakistan. Satyam and Towhata [14] assessed the ground response and liquefaction of the Vijayawada city (India). Aaqib et al. [15] Aaqib et al. [16] investigated the strength correction effect of sites in Korea. Recently, several studies have examined seismic hazard assessment of the Hanoi region ground [2, 7, 17, 18]. Alexandr and Chi [2] and Nhung et al. [17] accessed the liquefaction possibility. Wen et al. [18] and Nguyen et al. [7] reported that sites in the Hanoi area may experience a high amplification and site-specific site response analysis was recommended. This is because the Eurocode-8 (EC-8) [11] forms the basis of the current Vietnamese design code, TCVN 9386: 2012 [19]. The soil factors and the design spectrum used from EC-8 may not be suitable for the geological conditions of Vietnam. In addition, a range of ground motion characteristics such as frequency content and the magnitude has not been investigated yet in the previous studies.

This study performs a site-specific analysis and aims to develop a design spectrum for the Ba Dinh square area where the most important historical, cultural, and government buildings in Vietnam are located. A representative soil profile with a comprehensive site and laboratory investigation is utilized. A suite of input ground motions compatible with the site class A design spectrum in TCVN 9386: 2012 are used. The amplification and shear strain characteristics of the area are studied. Discrepancies with the local standard design of practice are also discussed.

2. The Importance of the Ba Dinh Square Area

Located in the heart of Hanoi, the Ba Dinh square is the largest square in the Ba Dinh district, Vietnam, with a dimension of 550 m (length) and 300 m (width). The square is the historical destination where the Ho Chi Minh Mausoleum was built. Ho Chi Minh is considered the founder of the Democratic Republic of Vietnam. This square also locates many important buildings around it, including the President’s Palace, the Ministry of Foreign Affairs, the Ministry of Planning and Investment, and the National Assembly Building. Cultural-historical buildings such as One Pillar Pagoda, the Imperial Citadel of Thang Long (UNESCO World Heritage Site), and the Ho Chi Minh Museum are set nearby. The location of the Ba Dinh square area is shown in Figure 2. It should be added that One Pillar Pagoda is a historic Buddhist temple in Hanoi. It is regarded, alongside the Perfume Temple, as one of Vietnam’s two most iconic temples. Thang Long Citadel is a complex of historic imperial buildings and was made a UNESCO World Heritage Site in 2010. It is also known as Hanoi Citadel.

The ground in the Ba Dinh district is mainly site class E according to EC-8 [10]. The Ba Dinh district soil strata typically include four layers [20]. The first layer is river, lake or marsh sediment. The thickness of this layer can be reached to 10 m near the river or lake and 0–0.5 m at other locations. It is followed by white stiff clay, clay mixed with pebbles, and gritstone or sandstone. The thickness of these layers is 10–35 m, 10–20 m, and 5–30 m, respectively. Furthermore, this square is very close to the Red River fault and is only 2.3 km away. The Red River has been shown to have the potential to cause earthquakes with  = 7-8 [21]. Therefore, it is important to clarify the significance of ground responses in this area.

3. Soil Profile

A site investigation was performed to characterize the soil profile. The location of the borehole is shown in Figure 2(b). Boring logs show that the depth to the bedrock is approximately 35 m. The site consists of fill, sandy clay, clayey sand, and sand. Shear wave velocity profile and soil properties are presented in Figure 3 and Table 1, respectively. The site natural period () is 0.72 s; the time-averaged shear wave velocity to 30 m depth () is 181 m/s and considered as site class C according to EC-8 and TCVN 9386: 2012. It should be noted that the site class of the Ba Dinh square is different than the typical site class in Ba Dinh district. The groundwater level has a significant effect on the ground response, and hence it is taken into account in this study. The water table is located at a depth of 9 m below the ground surface.

4. Input Ground Motions

Owing to a lack of recorded motions in the Ba Dinh square area, a total of 15 input motions from 12 earthquake events used in the study of Nguyen et al. [7] are adopted in this study. Notably, two acceleration-time history components (i.e. North-South and East-West) of earthquake events 5th, 8th, and 10th were used. The NGA-west2 database (http://ngawest2.berkeley.edu/) was utilized to select the ground motions. The target characteristics of the selected motions are compatible with those of the Hanoi city such as moment magnitude () of 5.5–7.0, rupture distance () from 0.90 to 85 km, and of 760–1500 m/s. The selected input ground motions were scaled to 0.0976 g which is the representative shaking intensity level for an earthquake return period of 475 years for the Ba Dinh district [19]. The ground motions were selected such that the mean response spectrum matches the target spectrum of site class A in TCVN 9386: 2012 [19], as shown in Figure 4. Table 2 presents in detail the characteristics of the selected earthquake events.

5. One-Dimensional Ground Response Analysis

The pattern of damage in urban areas during an earthquake depends on the characteristics of the event and interaction between site response and the vulnerability of the built structure. Most of the urban settlements have occurred along river valleys with soft and young soil deposits which were prone to severe damage during the earthquake [14]. The Ba Dinh square area is located nearby the Red River; the site-specific response analysis will help predict the local site effects and approximating the dynamic behavior of soil during seismic loading.

One-dimensional equivalent linear (EQL) and nonlinear (NL) response analysis is performed using DEEPSOIL v7.0 [22] in this study. This program has been widely used in previous studies [7, 16, 2326]. One-dimensional site response analysis is the process of propagating shear waves numerically from the reference rock through the overlying soil layers to the surface. One-dimensional site response allows vertical propagation of shear waves, and soil deposits are assumed as horizontally homogeneous layers. Conventionally, all the site response procedures require the specification of acceleration-time history, and a suite of input motions is usually applied to develop a statistically stable response. EQL is a method in which stiffness and damping of soil layers are calculated for a given shear strain based on EQL iterative procedure [1]. In the case of NL, the dynamic soil properties are determined by the stress–strain relationship at each time step of the loading. A recently developed Generalized Quadratic/Hyperbolic (GQ/H) constitutive model [27] implemented in DEEPSOIL v7.0 was used in this study. This is a widely used model in site response analysis [15, 24, 25, 2830] and can capture an initial shear modulus at zero shear strain as well as a limiting shear strength at large shear strains. The pressure-dependent nonlinear dynamic curves proposed by Darendeli [31] is used to generate dynamic properties of soil layers. The required parameters are over the consolidation ratio () of 1.0, the coefficient of lateral Earth pressure () of 0.5, the number of the cycle of loading () of 10, and the frequency of the loading () of 1.0. The plasticity index () of sand and clay were 0 and 15, respectively. These chosen parameters of Darendeli are widely used in site response analysis as reported in Aaqib et al. [16], Hashash et al. [28], Nguyen et al. [23], and Nguyen et al. [7]. Figure 5 presents a flowchart of the site response procedure implemented in this study. The flowchart illustrates that a suite of equivalent linear and nonlinear site response analyses are used to determine the site-specific spectral shapes, peak ground acceleration, and maximum shear strain profiles of the Ba Dinh square area. Moreover, the amplification factors resulting from this study will provide insights into the discrepancies between the design and site-specific spectra.

6. Results and Discussion

The analysis of induced shear strains is an important factor in any site-specific ground response study. On one hand, in case of the EQL procedure, the initial shear modulus () and damping () of the soil layer are used to calculate the initial shear strain time history. The degraded shear modulus and damping are then calculated using a shear strain of 65%. These soil properties are then subsequently used in the next iterations until the convergence with the target degradation curves is achieved. On the other hand, in the case of the NL procedure, shear stresses are induced at different layers during the propagation of input ground motion from bedrock to the surface, which, in turn, causes shear strains in these soil layers. It is clear that both NL and EQL procedures use distinct ways to calculate the dynamic properties of soil.

Figure 6 presents the variation of shear strains produced by a suite of ground motions with depth for both NL and EQL analyses. Notably, the curves in the gray color in Figures 6 to 8 are the results corresponding to the different ground motions. The solid red line and dashed black line are average results from NL and EQL analyses, respectively. It is demonstrated that the maximum strains calculated by both the procedures differ for a set of ground motions. The NL procedure results in maximum strains up to 0.47%, whereas the EQL procedure results in strains up to 0.4%. However, there are insignificant discrepancies in the mean shear strains calculated by both procedures. The average maximum shear strains are encountered at a depth of 11 m in both the procedures and go up to 0.2%. The high strains at 11 m are attributed to the presence of shear wave velocity reversal in the soil profile.

Figure 7 presents the variation of peak ground acceleration (PGA) as a function of depth for all the used ground motions and the two procedures. It is demonstrated that PGA increases as the surface is approached, and the peak PGA is observed at the surface by both NL and EQL procedures. The NL procedure calculates surface PGA up to 0.21 g, whereas the EQL results in a peak surface PGA of 0.24 g. The mean surface PGA calculated by the NL procedure is lower than that by the EQL procedure. A deamplification portion is observed at 11 m which is attributed to the presence of shear wave velocity reversal in soil profile at that depth.

Shear strains and shear modulus are also utilized in the calculation of mobilized shear strength. Figure 8 illustrates the mobilized shear strength profiles as a function of depth for all the ground motions. Mobilized shear strength is defined as yield stress developed due to the applied load to resist the deformation of the soil. As a general trend, mobilized shear strength is shown to increase with depth. Some discrepancies are observed where mobilized shear strength reversals are produced by a small set of input motions. In a site-specific study, this behavior is attributed to the use of different shear strains and shear modulus used by each ground motion.

The evaluation of a site for design purposes is most often performed using the acceleration response spectrum. It reflects the frequency content of the ground motion at the surface. Figure 9 depicts the mean surface response spectrum calculated using both the site response procedures. The input response spectrum and the rock code design spectrum () are also shown. EQL procedure results in a higher response than the NL procedure. A spectral bump is observed at the site natural period of 0.72 s, demonstrating the effect of the site period on amplification.

The amplification factor at a particular spectral period is calculated as the ratio of spectral acceleration at the ground surface and the spectral acceleration at the bedrock. The propagation of seismic waves results in the development of a standing wave at each natural frequency of the soil deposit which causes the amplification at the ground surface. Figure 10 compares the amplification factors calculated from NL and EQL procedures. Peak amplification is observed at the site natural period of 0.9 s in the case of NL analysis, whereas it shifts to 0.72 s in the case of EQL analysis. The soil factor, , is also calculated and shown in Figure 10. As TCVN 9386: 2012 is based on the Eurocode-8, is calculated as the average of amplification factors within the interval of 0.05–2.5 s, defined as follows [32]:where and are the spectral acceleration indexes of soil and rock, respectively, defined as follows: is calculated as 1.95 by NL and 2.07 by EQL analysis. The calculated soil factor is much higher than that defined in TCVN 9386: 2012 ( = 1.15 for ground type C).

Figure 11 compares the calculated mean response spectra and proposed design spectrum with the TCVN 9386: 2012 design spectrum (). The rock design spectrum () is also shown in this figure. The proposed design spectrum was obtained as a product of the rock design spectrum () and amplification factor (). The calculated spectra are higher than the code design spectrum across the whole period range. These questions the applicability of code provisions for Vietnam site conditions, which are primarily based on EC-8. These discrepancies were also reported by Nguyen et al. [7]. The calculated response spectra and proposed design spectrum match favorably except at the spectral periods of 0.65–1.26 s. This is consistent with the results of both NL and EQL approaches. This mismatch is attributed to the use of a single amplification factor for both short-period and long-period range in TCVN 9386: 2012 as well as EC-8. This was also highlighted in the study of Pitilakis et al. [33], where it was recommended that separate amplification factors should be used for short and long periods, as used in NEHRP 2015 provisions [34].

7. Conclusions

This study evaluates ground response at a crucial location in Vietnam which is the Ba Dinh square area. A representative soil profile and 15 compatible earthquake ground motions scaled to the representative design shaking level in the region are used. The one-dimensional nonlinear and equivalent linear site-specific ground response analyses are performed. The results are also compared to the current seismic design code TCVN 9386: 2012. The following conclusions are drawn:

A set of in-situ and laboratory test measurements were utilized. The site under investigation has a of 181 m/s and is classified as site C according to TCVN 9386: 2012. It is revealed that the site class in the region of study is different from the typical class in the Ba Dinh district (site E).

The NL procedure is usually preferred over the EQL procedure because of its ability to appropriately capture the actual response. Meanwhile, the EQL procedure is cost-effective and easier to apply. However, the EQL procedure was reported to be effective for shear strains <0.4% [35]. In this study, the maximum shear strains calculated by the EQL procedure are less than 0.4%. Therefore, both NL and EQL procedures can be used for the study area. Maximum shear strains in both cases were produced at the interface of the soil layer 2 (sandy clay) and layer 3 (clayey sand) where a shear wave velocity reversal is encountered. The average maximum shear strain reaches up to 2%. Mobilized shear strength is shown to increase as a function of depth.

The PGA profiles were shown to increase as the surface is approached except for the shear wave reversal portion. The average maximum PGA at surface is about 0.2 g and 0.16 g for NL and EQL analyses, respectively. The response spectral comparison indicated that the ground motions are clearly amplified near the site natural period of 0.72 s.

Comparisons of the calculated amplification factor and response spectra with the local design code of practice revealed significant discrepancies across the whole period range. The computed soil factors are 1.95 and 2.07 for NL and EQL analyses, respectively, and they are much higher than that defined in TCVN 9386: 2012 ( = 1.15 for ground type C). The calculated response spectra and design spectrum developed using the soil factor calculated in this study also shown discrepancies at the period range 0.65–1.2 s. It is recommended to use different soil factors for both short-period and long-period ranges rather than using a single factor in TCVN 9386: 2012.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.