Institute for Coastal Marine Environment, Italian National Research Council, Naples, Italy 2

Transcription

1 Microtremor HVSR analysis to highlight the effects of potential liquefaction: experimental results at Magoodhoo Island (Faafu Atoll, Maldives, Indian Ocean) M. Punzo 1, G. Cavuoto 1, D. Tarallo 1, N. Pelosi 1, P. Scotto Di Vettimo 1, M. Iavarone 1, S. Mazzola 2, P. Galli 3,4, V. Di Fiore 1 1 Institute for Coastal Marine Environment, Italian National Research Council, Naples, Italy 2 Institute for Coastal Marine Environment, Italian National Research Council, Torretta Granitola-Trapani, Italy 3 Department of Biotechnology and Biosciences, University of Milano-Bicocca, Milan, Italy 4 MaRHE Center (Marine Research and High Education Center), Magoodhoo Island, Faafu Atoll, Maldives Introduction. The Magoodhoo island (Maldives) is a part of coral atoll which presents several hazard linked to the effects of sea tides and waves, earthquakes and tsunamis. In particular, the island coralline coarse-grained sands lies on a coral rubble base; this sandy non-cohesive could experience liquefaction and therefore local subsidence following earthquake loads, severe storm wave loads, or a combination of the two. Cyclic load may cause an excess pore pressure (due to ground water mounding) and the shallow sands may be liquefied, and the island structures will become unstable. Today the effects induced by these phenomena were still visible differently on Magoodhoo island. For this purposes we performed a study by microtremor analysis (Horizontal-to-Vertical Spectral Ratio - HVSR: Nakamura, 1989) to identify possible site effects and to investigate the characteristics of soil liquefaction. Microtremor method is commonly used to detect the rigid basement which present shear wave velocity 750 m/s (Monge et al., 1999; Guéguen et al., 2000; Régnier et al., 2000). The HVSR method, in fact, is considered by many authors as giving a good estimation of the site fundamental frequency resonance (Lachet and Bard, 1995; Goula et al., 1997; Mucciarelli, 1998; Lebrun et al., 2001) and a very effective tool to decipher the local soil conditions (Singh et al., 2017). The passive recording of ambient vibrations may provide a mapping tool of site features where geotechnical information is usually difficult to obtain. Microtremors were recorded at 20 sites with single stations. Study area. The Maldives archipelago comprises a chain of coral reefs and reef islands situated 700 km SW of Sri Lanka (Fig. 1A). The archipelago extends 868 km from Ihavandiffulu Atoll in the north (7 10 N) to Addu Atoll (0 43 S) just south of the equator, and is 130 km wide. Its central sector consists of a double chain of atolls rising from a submerged plateau, whose Fig. 1 - A: The Maldives archipelago south-west of India, in the central equatorial Indian Ocean, is an isolated tropical carbonate platform. The archipelago comprises about 1200 smaller atolls. B: Faafu Atoll. C: Magoodhoo Island; the red circles indicate free-field measurements for which the resonance frequency was determined. 421

2 depth varies from over 500 m in the north to m in the south. The archipelago is made up of 22 atolls of circular or elongated shape, the sizes varying from some kilometers to many tens of kilometers, and contains more than 1000 islands which cover a surface of 298 km 2. Every atoll is formed by a marginal rim surrounding a lagoon commonly less than m deep, though some reach depths of more than 80 m. The atoll lagoons are characterized mainly by sandy sediments; a large amount of the sand (more than 50%) derives from mechanical erosion and from bioerosion of the reef. A minor amount is due to benthonic organisms living in the sand with fewer contributions from planktonic foraminifers. Magoodhoo, is a remote small island on the Faafu Atoll (Fig. 1B), at a distance of 134 km from Malé, Republic of Maldives capital, with a population of 683 inhabitants distributed on a surface of less than 0.36 km 2. Field measurements. The basic goal of single-station ambient vibration measurements is the detection of seismic impedance contrasts in the subsoil, that are responsible for seismic resonance phenomena (e.g., Kramer, 1996). The technique consists of recording and analyzing the seismic ambient noise to retrieve information about subsoil stratigraphy from seismic tremor. In particular, we use HVSR technique, which consists of studying the ratio between the Vertical and Horizontal spectral components of the tremor recorded at a single station. The method, provides a direct estimate of the resonance frequencies of subsoil. This well known capability allows to translate the H/V curves into stratigraphy (Castellaro et al., 2005). The theoretical bases of this technique are simple only in a 2 layer Fig. 2 - Results of HVSRs using microtremor measurements. 422

3 1-D model in which parameters are constant within each layer. In the latter case, the depth h of the seismic discontinuity is immediately derived from the free-surface resonance equation h=v/ (4f 0 ), where V is the seismic wave velocity of the upper layer and f 0 is the resonance frequency which appears as a peak in the H/V curve. Another important parameter is the vulnerability index (Kg) which can be determined using the empirical relation Kg=(A f0 ) 2 /f 0, in which A f0 is the amplitude at the fundamental frequency, which determines the damage level due to local site effects (Nakamura, 1996). It may be considered for predicting the potential for soil liquefaction. This value, in fact, is clearly site-specific and value might be useful in demarcating weak regions, which are susceptible to liquefaction (Huang and Tseng, 2002; Beroya et al., 2009). Tab. 1 - HVSR stations: Coordinates, peaks frequency and vulnerability index. NAME COORDINATE RESONANCE FREQUENCY VULNERABILITY INDEX STATION 1 E N STATION 2 E N STATION 3 E N STATION 4 E N STATION 5 E N STATION 6 E N STATION 7 E N STATION 8 E N STATION 9 E N STATION 10 E N STATION 11 E N STATION 12 E N STATION 13 E N STATION 15 E N STATION 16 E N STATION 17 E N STATION 18 E N STATION 19 E N STATION 20 E N

4 The tremor recordings were acquired with a three-component short period sensors (Lennartz 3D-5s) specifically designed to record seismic noise. Seismic noise was recorded at least 60 minutes at each site with a sampling frequency of 50 Hz. Microtremor data were analyzed using the open source software Geopsy which contains a dedicated tool to obtain Horizontal to Vertical spectral ratio starting from 3 signals registration. The seismic data were processed using the following parameters: short-term averaging equals 2 s, long-term averaging equals 30 s, with low and high thresholds of 0.2 and 3.0, respectively. Before averaging, the individual spectra were smoothed using the Konno and Ohmachi (1998) method, using a bandwidth coefficient (b) value of 40. Finally, a mean horizontal-to-vertical (H/V) curve was estimated for each site (Fig. 2) following the standard criteria adopted globally (SESAME, 2004). The location of 20 site tests (Tab. 1) are shown in Fig. 1C. Results. HVSR analyses of 20 free-field measurements taken on an approximate 0.16 Km 2 area. The frequency distribution is in a relatively large range: Hz. The maximum Kg values are found at stations 18 and 19, in the central south-eastern part of the island. In the south-western part the stations 1, 2,3, 4, 5, 6, 7, 8, 9, 10, 11 show peaks in the frequency range Hz, which have medium amplitudes (about 2). This indicates a relatively verythin sedimentary cover, implying that the Kg values are lowest at these sites ( ). Points 13, 15, 18 which are located in the north-eastern part of the Island show very low frequencies ( Hz), which have relatively high amplitudes (4-6). This indicates a thick layer of sediments and a strong impedance contrast with the bedrock; the estimated vulnerability index range 5.91 to In the north-eastern part, points 12, 16, 17, 20 are characterized by slightly higher frequencies ( Hz) and slightly better expressed peaks. The frequencies and the amplitude of the HVSR peaks determined at 20 points were used to contour the iso-frequency and the ground vulnerability index contour maps (Figs. 3A and 3B), showing the resonance frequency of the sediments. Conclusions. Since the geotechnical characteristics of sediments and their thickness are poorly known in the Magoodhoo Island, due to lack of borehole or geophysical data, we performed single microtremor measurements to determine the characteristics of site response of the shallow subsurface. The entire Magoodhoo island area was surveyed with a dense grid of free-field measurements, in a reasonable time and with relatively low costs in order to detect the potential danger of soil-structure resonance. A map of the fundamental frequencies of the sediments shows a distribution in a wide range Hz (Fig. 3A). There is a good correlation of different frequency ranges with the supposed lithology and thickness of the Quaternary sediments. The high-frequency range is characteristic of the south-western part of the Island, where the sediments are thin. The northeastern part of the Island is characterized by medium-low sediment frequencies, indicating thick Plio-Quaternary sediments. The transition between higher and lower frequencies zone is clean. According to our study there is a strong difference in impedance contrast between the south-western and north-eastern area of the Island. The area located at south-western is characterized by a high resonance frequency due to the presence of very shallow beach-rock. It is evident because there aren t structures building. On the other side, in the north-eastern part of Magoodho, the strong influence of incoherent and large thickness sediments (very low resonance frequency) and the presence of water table variation in shallow sand deposit, produce partial or total liquefaction effects with serious building damages. In particular, the weakest zone is in the central south-eastern part of the Island (Fig. 3B); larger Kg values are found in the stations 18 and 19 characterized by high amplitude at lower fundamental frequency. Therefore the extensive damage would be expected in this area of the island. Acknowledgements The authors would like to acknowledge the help and cooperation of the Bicocca University staff: Paolo Galli (Responsible Center of Milano-Bicocca University - Magoodhoo Island), Simone Montano, Andrea Di Pietro, Laura Bernasconi de Luca, Angelica Cajiao, Davide Seveso. We wish to thank Bicocca Milan University 424

5 GNGTS 2017 Sessione 2.2 Fig. 3 - Contour maps of (A) sediments resonance frequency and (B) ground vulnerability index for support the realization of this research activities. We also wish to thank all Maldivian Authorities and Magoodhoo People for patience keep during the field research activities. References Beroya M.A.A., Aydin A., Tiglao R., Lasala M.; 2009: Use of microtremor in liquefaction hazard mapping. Eng. Geol. 107, Castellaro S., Mulargia F., Bianconi L.; 2005: Passive seismic stratigraphy: a new efficient, fast and economic technique. J. Geotech. Environ. Geol. 3, Geopsy; Geophysical Signal database for noise array processing. Version 2.9.0, Copyright by Marc Wathelet (marc@geopsy.org), available at Goula X., Susagna T., Figueras S., Farres P., Cid X., Alfaro A., Barchiesi A.; 1997: Analysis of Site Effect in the City of Barcelona (Spain). XIX General Assembly of the IASPEI, Thessaloniki, Abstract Book, p

6 Guéguen Ph., Chatelain J-L., Guillier B., Yepes H.; 2000: An indication of the topmost layer response in Quito (Ecuador) using noise H/V spectral ratio. Soil Dyn Earth Eng 19: Huang H.C., Tseng Y.S.; 2002: Characteristics of soil liquefaction using H/V of microtremors in Yuan-Lin area, Taiwan. Terr. Atmos. Ocean. Sci. 13, no. 3, Konno K., Ohmachi T.; 1998: Ground-motion characteristics estimated from spectral ratio between horizontal and vertical components of microtremor. Bull Seismological Soc Am 88(1): Kramer S.L.; 1996: Geotechnical Earthquake Engineering. Prentice Hall, Inc., 653 pp. Lachet C., Bard P-Y.; 1995: Theoretical investigation on the Nakamura s technique. Proceeding 3rd International Conference On Recent Advances in Geot. Earthq. Eng. & Soil Dyn., April 2 7, II, St Louis, MI, Paper no Lebrun B., Hatzfeld D., Bard P-Y.; 2001: A site effect study in urban area: experimental results in Grenoble (France). Pure Appl Geophys 158: Monge O., Chassagneux D., Martin C., Sedan O., Vermeersch F.; 1999: Evaluation de l aléa sismique local: partie 1, microzonages de Fort-de-France et Pointe-à-Pitre. Génie parasismique et réponse dynamique des ouvrages, proc. 5ème Coll. Nat. AFPS, Cachan, vol. 1, p , octobre Mucciarelli M.; 1998: Reliability and applicability of Nakamura s technique using microtremors: an experimental approach. J Earthq Eng 2(4): Nakamura Y.; 1989: A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface. Q Rep Railway Tech Res Inst 30(1): Nogoshi M., Igarashi T.; 1970: On the propagation characteristics of microtremors. J. Seism. Soc. Jpn. 23, Régnier M., Moris S., Shapira A., Malitzky A., Shorten G.; 2000: Microzonation of the expected seismic site effects across Port Vila, Vanuatu. J Earthq Eng 4(2): SESAME; 2004: Guidelines for the implementation of the H/V spectral ratio technique on ambient vibrations: measurements, processing and interpretation. Delivrables/Del-D23 HV_ User_Guidelines.pdf. Singh A.P., Shukla A., Ravi Kumar M., Thakkar M. G.; 2017: Characterizing Surface Geology, Liquefaction Potential, and Maximum Intensity in the Kachchh Seismic Zone, Western India, through Microtremor Analysis. Bulletin of the Seismological Society of America, Vol. 107, No. 3, pp