Three-Dimensional P- and S-Wave Velocity Profiling of Geotechnical Sites Using Full-Waveform Inversion Driven by Field Data

Total Page:16

File Type:pdf, Size:1020Kb

Three-Dimensional P- and S-Wave Velocity Profiling of Geotechnical Sites Using Full-Waveform Inversion Driven by Field Data Soil Dynamics and Earthquake Engineering 87 (2016) 63–81 Contents lists available at ScienceDirect Soil Dynamics and Earthquake Engineering journal homepage: www.elsevier.com/locate/soildyn Three-dimensional P- and S-wave velocity profiling of geotechnical sites using full-waveform inversion driven by field data Arash Fathi a, Babak Poursartip b, Kenneth H. Stokoe II b, Loukas F. Kallivokas a,b,n a Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX, USA b Department of Civil, Architectural and Environmental Engineering, The University of Texas at Austin, Austin, TX, USA article info abstract Article history: We discuss the application of a recently developed full-waveform-inversion-based technique to the Received 22 December 2015 imaging of geotechnical sites using field-collected data. Specifically, we address the profiling of arbitrarily Received in revised form heterogeneous sites in terms of P- and S-wave velocities in three dimensions, using elastic waves as 18 April 2016 probing agents. We cast the problem of finding the spacial distribution of the elastic soil properties as an Accepted 19 April 2016 inverse medium problem, directly in the time domain, and use perfectly-matched-layers (PMLs) to ac- Available online 7 May 2016 count for the semi-infinite extent of the site under investigation. Keywords: After briefly reviewing the theoretical and computational aspects of the employed technique, we Site characterization focus on the characterization of the George E. Brown Jr. Network for Earthquake Engineering Simulation Full-waveform inversion (FWI) site in Garner Valley, California (NEES@UCSB). We compare the profiles obtained from our full-wave- Elastic waves form-inversion-based methodology against the profiling obtained from the Spectral-Analysis-of-Surface- perfectly-matched-layers (PMLs) Spectral-Analysis-of-Surface-Waves (SASW) Waves (SASW) method, and report agreement. In an attempt to validate our methodology, we also Vibroseis compare the recorded field data at select control sensors that were not used for the full-waveform in- Field data processing version, against the response at the same sensors, computed based on the full-waveform-inverted profiles. We report very good agreement at the control sensors, which is a strong indicator of the cor- rectness of the inverted profiles. Overall, the systematic framework discussed herein seems robust, general, practical, and promising for three-dimensional site characterization purposes. & 2016 Elsevier Ltd. All rights reserved. 1. Introduction profiling of the soil shares common elements with medical ima- ging [17]. In both cases, elastic waves are used to interrogate the Reliable three-dimensional P- and S-wave velocity profiles of unseeable medium, and the medium's response to this probing is the near-surface soil have great value in geotechnical engineering subsequently used for driving the imaging. In the mathematical practice, and can play a vital role in the safe design of critical sciences, this process is referred to as solving an inverse medium components of the civil infrastructure, such as nuclear power problem. fi plants, bridges, and hospitals, among others. Wave velocity pro- Our goal is to nd the distribution of the P-wave velocity cp(x), 1 filing of the soil is referred to as seismic site characterization. and S-wave velocity cs(x) of the soil in three-space dimensions. We use a seismic vibrator to apply loads on the ground surface, Seismic site characterization techniques can be broadly classified which results in the propagation of elastic waves in the soil. Due to into two categories: invasive and non-invasive. Invasive techni- the possibly heterogeneous character of the site, multiple reflec- ques, such as borehole-based methods, are costly and can typically tions and refractions occur. The time-history response of the soil provide only localized information of the site under investigation, medium to this loading is recorded by using geophones, spatially while, non-invasive techniques are capable of providing global arranged over the formation's surface, and is, hereafter, referred to fi imaging, and are generally more cost-ef cient. Non-invasive as measured response (Fig. 1). Arriving at a material profile can then be achieved by minimizing the difference between the measured n Corresponding author at: Department of Civil, Architectural and Environmental response at sensor locations, and a computed response corre- Engineering, The University of Texas at Austin, Austin, TX, USA. sponding to a trial distribution of the material parameters, under E-mail addresses: [email protected] (A. Fathi), [email protected] (B. Poursartip), [email protected] (K.H. Stokoe II), [email protected] (L.F. Kallivokas). 1 x =(xyz,,) is the position vector. http://dx.doi.org/10.1016/j.soildyn.2016.04.010 0267-7261/& 2016 Elsevier Ltd. All rights reserved. 64 A. Fathi et al. / Soil Dynamics and Earthquake Engineering 87 (2016) 63–81 Fig. 1. Problem definition: (a) interrogation of a heterogeneous semi-infinite domain by an active source; and (b) computational model truncated from the semi-infinite medium via the introduction of perfectly-matched-layers (PMLs). the same loading, using a computational model. The computa- data, which are inherently three-dimensional, into two-dimensional tional model should mimic the physics of the problem as closely as full-waveform-inversion-based software, whenever the site condi- possible. A key challenge is proper termination of the semi-infinite tions would support plane strain assumptions. We demonstrated the extent of the domain of interest.2 We use perfectly matched-layers method's viability by characterizing a site in Austin, Texas [19].Inthis (PMLs) for domain truncation (Fig. 1) [5,15,39], since they are communication, we abandon the plane strain assumption, and con- among the best available tools in the presence of heterogeneity sider arbitrary heterogeneity in three space dimensions. To be able to [15]. A PML is a (practically reflectionless) buffer zone that sur- extend the framework described in [19] to the fully three-dimen- rounds the domain of interest, where outgoing waves are forced to sional case, various algorithmic modifications are needed, including decay. modifications to the forward three-dimensional wave simulator, per Obtaining the P- and S-wave velocity profiles of the soil in [15], to the inversion methodology, per [14], and to the manner the three-dimensional, arbitrarily heterogeneous domains, remains a field experiment is designed and the associated field data are pro- challenging problem. The inverse medium problem at hand is cessed post-recording. We use the mathematical apparatus that we notoriously ill-posed, computationally expensive, posed over a described recently in [13-15], and excerpt herein, but focus on the semi-infinite medium, and has challenging physics. Due to these three-dimensional characterization of the NEES@UCSB site in Garner complexities, and in order to render the problem tractable, current Valley, California, as a case study. To this end, we integrate recent techniques rely on simplifying assumptions. We classify these advances in several areas. Specifically, we use: (a) a parallel, three- simplifications as follows: (a) limiting the spatial variability of the dimensional, state-of-the-art wave simulation tool for domains ter- soil properties, whereby it is assumed that the soil is horizontally minated by PMLs [15]; (b) a partial-differential-equation (PDE)-con- layered (one-dimensional) [25,29,33], or has properties varying strained optimization framework, through which the misfit mini- only within a plane (two-dimensional) [3,12,19,21]; (b) assuming mization between the measured response and the computed response that the measured response of the soil at sensor locations is due is enforced [14]; (c) regularization schemes to alleviate the ill-po- only to Rayleigh waves, thus neglecting other wave types, such as sedness and solution multiplicity inherent in inverse problems, compressional and shear waves, as is the case in the SASW [33],or where the regularization factor is computed adaptively at each in- its close variant, the MASW method [30]; (c) idealizing the soil version iteration [14,20]; (d) continuation schemes that provide al- medium, which is porous, and, generally, partially or fully satu- gorithmic robustness [20]; and (e) a biasing scheme for the robust rated, as an elastic solid; (d) imaging only one elastic property, fi simultaneous inversion of both cp(x) and cs(x) pro les [14,21]. such as the shear wave velocity or an equivalent counterpart Most inversion approaches result in profile reconstructions that [2,11,27,28]; and (e) grossly simplifying the boundary conditions are physically acceptable. Validation of the inverted profiles, in the fi associated with the semi-in nite extent of the medium, due to the absence of borehole data, becomes challenging, if not impossible. complexity and computational cost that a rigorous treatment Herein, to aid in the validation of the inverted profiles, we use a set would require [11,36]. In recent years, the ubiquity of parallel of control sensors. Specifically, out of the 35 sensors that we de- computers, and significant advances in computational geosciences, ployed at the NEES@UCSB site, we use only 30 of them to feed our has created the opportunity of developing a toolkit that is capable full-waveform-inversion algorithms. The remaining 5 sensors of robust, accurate, and three-dimensional characterization
Recommended publications
  • Effect of Fabric Anisotropy on the Dynamic Mechanical Behavior Of
    EFFECT OF FABRIC ANISOTROPY ON THE DYNAMIC MECHANICAL BEHAVIOR OF GRANULAR MATERIALS by Bo Li Submitted in partial fulfillment of the requirements For the degree of Doctor of Philosophy Dissertation Advisor: Prof. Xiangwu Zeng Department of Civil Engineering Case Western Reserve University January, 2011 CASE WESTERN RESERVE UNIVERSITY SCHOOL OF GRADUATE STUDIES We hereby approve the thesis/dissertation of Bo Li Candidate for Ph. D. degree*. (signed) Xiangwu Zeng Chair of Commitee Adel Saada Wojbor A. Woyczynski Xiong Yu (date) 11/15/2010 *We also certify that written approval has been obtained for any proprietary material contained therein Table of Contents Table of Contents I List of Tables IV List of Figures VI Acknowledgement XI Abstract XII Notation XIV Chapter 1 Introduction ........................................................................................ 1 1.1 Effect of Anisotropy on Granular Materials ........................................................................... 1 1.2 Research Objectives ............................................................................................................... 1 1.3 Outline of the Dissertation ...................................................................................................... 2 1.4 Review of the Effect of Fabric Anisotropy on Clay ............................................................... 3 1.4.1 Experimental Study on Effect of Fabric Anisotropy on Clay ......................................... 3 1.4.2 Analytical Method on Effect of Fabric Anisotropy
    [Show full text]
  • Confidential Manuscript Submitted to Engineering Geology
    Confidential manuscript submitted to Engineering Geology 1 Landslide monitoring using seismic refraction tomography – The 2 importance of incorporating topographic variations 3 J S Whiteley1,2, J E Chambers1, S Uhlemann1,3, J Boyd1,4, M O Cimpoiasu1,5, J L 4 Holmes1,6, C M Inauen1, A Watlet1, L R Hawley-Sibbett1,5, C Sujitapan2, R T Swift1,7 5 and J M Kendall2 6 1 British Geological Survey, Environmental Science Centre, Nicker Hill, Keyworth, 7 Nottingham, NG12 5GG, United Kingdom. 2 School of Earth Sciences, University 8 of Bristol, Wills Memorial Building, Queens Road, Bristol, BS8 1RJ, United 9 Kingdom. 3 Lawrence Berkeley National Laboratory (LBNL), Earth and 10 Environmental Sciences Area, 1 Cyclotron Road, Berkeley, CA 94720, United 11 States of America. 4 Lancaster Environment Center (LEC), Lancaster University, 12 Lancaster, LA1 4YQ, United Kingdom 5 Division of Agriculture and Environmental 13 Science, School of Bioscience, University of Nottingham, Sutton Bonington, 14 Leicestershire, LE12 5RD, United Kingdom 6 Queen’s University Belfast, School of 15 Natural and Built Environment, Stranmillis Road, Belfast, BT9 5AG, United 16 Kingdom 7 University of Liege, Applied Geophysics, Department ArGEnCo, 17 Engineering Faculty, B52, 4000 Liege, Belgium 18 19 Corresponding author: Jim Whiteley ([email protected]) 20 21 Copyright British Geological Survey © UKRI 2020/ University of Bristol 2020 22 1 Confidential manuscript submitted to Engineering Geology 23 Abstract 24 Seismic refraction tomography provides images of the elastic properties of 25 subsurface materials in landslide settings. Seismic velocities are sensitive to 26 changes in moisture content, which is a triggering factor in the initiation of many 27 landslides.
    [Show full text]
  • On the Wave Bottom Shear Stress in Shallow Depths: the Role of Wave Period and Bed Roughness
    water Article On the Wave Bottom Shear Stress in Shallow Depths: The Role of Wave Period and Bed Roughness Sara Pascolo *, Marco Petti and Silvia Bosa Dipartimento Politecnico di Ingegneria e Architettura, University of Udine, 33100 Udine, Italy; [email protected] (M.P.); [email protected] (S.B.) * Correspondence: [email protected]; Tel.: +39-0432-558-713 Received: 6 September 2018; Accepted: 25 September 2018; Published: 28 September 2018 Abstract: Lagoons and coastal semi-enclosed basins morphologically evolve depending on local waves, currents, and tidal conditions. In very shallow water depths, typical of tidal flats and mudflats, the bed shear stress due to the wind waves is a key factor governing sediment resuspension. A current line of research focuses on the distribution of wave shear stress with depth, this being a very important aspect related to the dynamic equilibrium of transitional areas. In this work a relevant contribution to this study is provided, by means of the comparison between experimental growth curves which predict the finite depth wave characteristics and the numerical results obtained by means a spectral model. In particular, the dominant role of the bottom friction dissipation is underlined, especially in the presence of irregular and heterogeneous sea beds. The effects of this energy loss on the wave field is investigated, highlighting that both the variability of the wave period and the relative bottom roughness can change the bed shear stress trend substantially. Keywords: finite depth wind waves; bottom friction dissipation; friction factor; tidal flats; bottom shear stress 1. Introduction The morphological evolution of estuarine and coastal lagoon environments is a widely discussed topic in literature, because of its complexity and the great variability of morphological patterns.
    [Show full text]
  • Evaluation of Liquefaction Potential of Soil by Down Borehole Method
    Missouri University of Science and Technology Scholars' Mine International Conferences on Recent Advances 2010 - Fifth International Conference on Recent in Geotechnical Earthquake Engineering and Advances in Geotechnical Earthquake Soil Dynamics Engineering and Soil Dynamics 26 May 2010, 4:45 pm - 6:45 pm Evaluation of Liquefaction Potential of Soil by Down Borehole Method Sumit Ghose Bengal Engineering and Science University, India Ambarish Ghosh Bengal Engineering and Science University, India Follow this and additional works at: https://scholarsmine.mst.edu/icrageesd Part of the Geotechnical Engineering Commons Recommended Citation Ghose, Sumit and Ghosh, Ambarish, "Evaluation of Liquefaction Potential of Soil by Down Borehole Method" (2010). International Conferences on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics. 15. https://scholarsmine.mst.edu/icrageesd/05icrageesd/session01/15 This work is licensed under a Creative Commons Attribution-Noncommercial-No Derivative Works 4.0 License. This Article - Conference proceedings is brought to you for free and open access by Scholars' Mine. It has been accepted for inclusion in International Conferences on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics by an authorized administrator of Scholars' Mine. This work is protected by U. S. Copyright Law. Unauthorized use including reproduction for redistribution requires the permission of the copyright holder. For more information, please contact [email protected]. EVALUATION OF LIQUEFACTION
    [Show full text]
  • Defining P-Wave and S-Wave Stratal Surfaces with Nine-Component Vsps
    INTERPRETER’S CORNER Coordinated by Rebecca B. Latimer Defining P-wave and S-wave stratal surfaces with nine-component VSPs BOB A. HARDAGE, MICHAEL DEANGELO, and PAUL MURRAY, Bureau of Economic Geology, Austin, Texas, U.S. Nine-component vertical seismic profile (9C VSP) data were acquired across a three-state area (Texas, Kansas, Colorado) to evaluate the relative merits of imaging Morrow and post-Morrow stratigraphy with compressional (P-wave) seismic data and shear (S-wave) seismic data. 9C VSP data generated using three orthogonal vector sources were used in this study rather than 3C data generated by only a ver- tical-displacement source because the important SH shear mode would not have been available if only the latter had been acquired. The popular SV converted mode (or C wave) utilized in 3C seismic technology is created by P-to-SV mode conversions at nonvertical angles of incidence at subsurface interfaces rather than propagating directly from an SV source as in 9C data acquisition. In contrast, 9C acquisition generates a P-wave mode and both fundamental S-wave modes (SH and SV) directly at the source station and allows these three basic compo- nents of the elastic seismic wavefield to be isolated from each other for imaging purposes. The purpose of these field tests was to compare SH, SV, and P reflectivity at targeted inter- faces. Because C waves were not an imaging objective, we used zero-offset acquisition geometry, which caused C-wave reflectivity to be quite small. We consider zero offset a source- to-receiver offset that is less than one-tenth of the depth to the shallowest receiver station.
    [Show full text]
  • Seismic Site Classifications for the St. Louis Urban Area
    Missouri University of Science and Technology Scholars' Mine Geosciences and Geological and Petroleum Geosciences and Geological and Petroleum Engineering Faculty Research & Creative Works Engineering 01 Jun 2012 Seismic Site Classifications for the St. Louis Urban Area Jaewon Chung J. David Rogers Missouri University of Science and Technology, [email protected] Follow this and additional works at: https://scholarsmine.mst.edu/geosci_geo_peteng_facwork Part of the Geology Commons Recommended Citation J. Chung and J. D. Rogers, "Seismic Site Classifications for the St. Louis Urban Area," Bulletin of the Seismological Society of America, vol. 102, no. 3, pp. 980-990, Seismological Society of America, Jun 2012. The definitive version is available at https://doi.org/10.1785/0120110275 This Article - Journal is brought to you for free and open access by Scholars' Mine. It has been accepted for inclusion in Geosciences and Geological and Petroleum Engineering Faculty Research & Creative Works by an authorized administrator of Scholars' Mine. This work is protected by U. S. Copyright Law. Unauthorized use including reproduction for redistribution requires the permission of the copyright holder. For more information, please contact [email protected]. Bulletin of the Seismological Society of America, Vol. 102, No. 3, pp. 980–990, June 2012, doi: 10.1785/0120110275 Seismic Site Classifications for the St. Louis Urban Area by Jae-won Chung and J. David Rogers Abstract Regional National Earthquake Hazards Reduction Program (NEHRP) soil class maps have become important input parameters for seismic site character- ization and hazard studies. The broad range of shallow shear-wave velocity (VS30, the average shear-wave velocity in the upper 30 m) measurements in the St.
    [Show full text]
  • Shear Wave Velocity and Its Effect on Seismic Design Forces and Liquefaction Assessment
    Missouri University of Science and Technology Scholars' Mine International Conference on Case Histories in (2004) - Fifth International Conference on Case Geotechnical Engineering Histories in Geotechnical Engineering 17 Apr 2004, 10:30am - 12:30pm Shear Wave Velocity and Its Effect on Seismic Design Forces and Liquefaction Assessment Thomas G. Thomann URS Corporation, New York, New York Khaled Chowdhury URS Corporation, New York, New York Follow this and additional works at: https://scholarsmine.mst.edu/icchge Part of the Geotechnical Engineering Commons Recommended Citation Thomann, Thomas G. and Chowdhury, Khaled, "Shear Wave Velocity and Its Effect on Seismic Design Forces and Liquefaction Assessment" (2004). International Conference on Case Histories in Geotechnical Engineering. 2. https://scholarsmine.mst.edu/icchge/5icchge/session03/2 This work is licensed under a Creative Commons Attribution-Noncommercial-No Derivative Works 4.0 License. This Article - Conference proceedings is brought to you for free and open access by Scholars' Mine. It has been accepted for inclusion in International Conference on Case Histories in Geotechnical Engineering by an authorized administrator of Scholars' Mine. This work is protected by U. S. Copyright Law. Unauthorized use including reproduction for redistribution requires the permission of the copyright holder. For more information, please contact [email protected]. SHEAR WAVE VELOCITY AND ITS EFFECT ON SEISMIC DESIGN FORCES AND LIQUEFACTION ASSESSMENT Thomas G. Thomann Khaled Chowdhury URS Corporation URS Corporation New York, NY (USA) New York, NY (USA) ABSTRACT The shear wave velocity of soil and rock is one of the key components in establishing the design response spectra, and therefore the seismic design forces, for a building, bridge, or other structure.
    [Show full text]
  • Relationship Between Shear-Wave Velocity and Geotechnical Parameters for Norwegian Clays
    Relationship between Shear-Wave Velocity and Geotechnical Parameters for Norwegian Clays Jean-Sebastien L’Heureux, Ph.D.1; and Michael Long2 Abstract: A database of shear-wave velocity (Vs) measurements using a variety of techniques and soil properties measured on high-quality samples for 28 Norwegian sites has been established. The purpose was to evaluate the different methods of measuring Vs, to present guide- lines and correlations to assist in estimating Vs profiles in these clays in the absence of site-specific data, and to outline relationships that can be used to give first-order estimates of soil properties. It was found that consistent measurements of Vs can be obtained from a variety of techniques and that for practical engineering purposes the Vs values obtained from the different methods are similar. Surface wave techniques can be particularly useful but careful survey design is necessary and in particular the inversion process needs to be carefully controlled. Differences of about 15–20% can be obtained in the Vs values depending on the algorithm used. Vs values for Norwegian clays are consistent with well-established frameworks for other materials, based on relationships between effective stress and index parameters. Piezocone pen- etration testing (CPTU) can be used to give acceptable estimates of Vs and this includes techniques which utilize the CPTU data only and are independent of any index property. Vs correlates well with triaxial compression and direct simple shear derived undrained shear strength (su) 0 values. There appears to be a particularly good link between Vs and preconsolidation stress (pc). Satisfactory relationships also exist 0 between Vs and the tangent moduli of the clays at in situ stress (M0) and at pc (ML).
    [Show full text]
  • Liquefaction Potential Assesment of Some Coastal Sands of Lagos, South-Western Nigeria, Using Geophysical and Geotechnical Methods
    LIQUEFACTION POTENTIAL ASSESMENT OF SOME COASTAL SANDS OF LAGOS, SOUTH-WESTERN NIGERIA, USING GEOPHYSICAL AND GEOTECHNICAL METHODS Abstract This paper presents the results of the investigation conducted in some parts of wetland areas of Lagos, Nigeria, using Multi-channel analysis of surface waves (MASW), Cone penetration test (CPT) and Standard penetration test (SPT). The primary aim was to delineate probable areas that were prone to induced seismicity. The MASW of the seismic method was used to generate the shear wave velocity (Vs) of the near surface soil while the CPT and SPT were employed to infer the penetration resistance and the blow count in the assessment of the stress-based liquefaction potential of these soils respectively. SeisImager and liquefaction assessment software (known as Cliq) were used to process both the MASW and CPT data respectively. The values of shear wave velocity generated for most sands in the study areas range from 120m/s -200m/s. This value fell within potentially liquefiable sands. From the CPT results, 41.67% of Ikoyi data showed a very high risk to liquefaction while 37% of Badore data indicated severe liquefaction potential. The factors of safety (F.S) against liquefaction potential obtained from the plot of cyclic stress ratio (CSR) against s-wave velocity (Vs) were less than 1. This also corroborated the presence of liquefiable sands within the study areas. Background of the Study Liquefaction phenomenon is one of the geological hazards that resulted from sudden ground shaken (Earthquake). It leads to the loss of strength/stiffness of saturated or partially saturated cohesionless soil which makes it to lose their bearing capacity as they behave like liquid.
    [Show full text]
  • Download/Menu.Php), Reid, M.E., T.W
    Zhang and Wang Geoenvironmental Disasters (2019) 6:9 Geoenvironmental Disasters https://doi.org/10.1186/s40677-019-0125-9 RESEARCH Open Access Three-dimensional seismic slope stability assessment with the application of Scoops3D and GIS: a case study in Atsuma, Hokkaido Shuai Zhang* and Fawu Wang Abstract Background: The 2018 Hokkaido Eastern Iburi earthquake (Mj 6.7, Mw 6.6), occurred on the 6th of September 2018 and triggered thousands of landslides in pyroclastic fall deposits. The triggered landslides caused destructive damages to the structures and resulted in serious causalities. Hundreds of earthquakes persisted after the mainshock and the Iburi region is of high possibility to suffer major earthquakes in the future. An effective method to assess the seismic slope stability is of great importance for the disaster prevention and mitigation in the Iburi region. Results: Based on the isopachs of different pyroclastic fall deposits mantled in the study area, GIS was employed to process the input soil layers and construct the 3D domain. By applying different horizontal pseudo-acceleration coefficients in the Scoops3D program, the factor-of-safety maps of eight cases were obtained. After validating by the coseismic landslide inventory, the performance of the computed results was evaluated. Conclusions: A horizontal pseudo-acceleration coefficient between 1/2 and 2/3 of PHGA is suitable for seismic slope stability assessment in pyroclastic fall deposits. The catastrophic Tomisato landslide and Yoshinoya landslide were correctly predicted. Scoops3D proves to be an effective and efficient method for guiding disaster mitigation and management. Keywords: Scoops3D, 3D seismic slope stability analysis, Iburi earthquake, Coseismic landslides, Pyroclastic fall deposits Introduction equilibrium theory.
    [Show full text]
  • Detectability of Seismic Waves from the Submarine Landslide That Caused the 1998 Papua New Guinea Tsunami
    Detectability of seismic waves from the submarine landslide that caused the 1998 Papua New Guinea tsunami Akio Katsumata1, Yasuhiro Yoshida2, Kenji Nakata1, Kenichi Fujita1, Masayuki Tanaka1, Koji Tamaribuchi1, Takahito Nishimiya1, and Akio Kobayashi1 1Meteorological Research Institute, Japan Meteorological Agency, 1-1 Nagamine, Tsukuba, Ibaraki Prefecture, Japan 2Meteorological College, Japan Meteorological Agency, 7-4-81 Asahi-cho, Kashiwa, Chiba Prefecture, Japan Correspondence: Akio Katsumata ([email protected]) Abstract. On 17 July 1998, a tsunami caused serious damage on the northern coast of Papua New Guinea about 20 min after the mainshock of an Mw 7.0 earthquake. The tsunami has been attributed to a submarine landslide that occurred about 13 min after the mainshock because its arrival at the coast was too late and its height too great to be the direct result of the fault slip of the earthquake. Bathymetric data recorded after the tsunami revealed an amphitheater-like structure that was consistent with a 5 recent submarine landslide. Most current tsunami warning systems are based on analysis of the early arrivals of seismic waves generated by an earthquake. In this study we investigated whether evidence of the landslide could be identified in the coda waves recorded after the mainshock. Based on previous studies of the tsunami source, we constructed synthetic seismograms to represent the submarine landslide and compared them to the observed coda waves of the preceding earthquake, with particular attention to the period around 13 min after the mainshock ::in ::::::::frequency::::::ranges::::close::to:::the::::::::landslide:::::::duration. We found phases 10 possibly corresponding to the landslide event.
    [Show full text]
  • Activity–Seismic Slinky©
    Activity–Seismic Slinky© Slinkies prove to be a good tool for modeling the behavior of compressional P waves and shearing S waves. We recommend reading about the behavior of seismic waves and watching the variety of animations below to understand how they travel and how the P, S, and surface waves differ from each other. Science Standards (NGSS; pg. 287) • From Molecules to Organisms—Structures and Processes: MS-LS1-8 Master Teacher Bonnie Magura demonstrates S waves using • Motion and Stability—Forces two metal slinkies taped together. and Interactions: HS-PS2-1, MS-PS2-2 to the taped joint to emphasize the movement directions. • Energy: MS-PS3-1, MS-PS3-2, HS-PS3-2, NOTE that when a metal slinky is taped to a plastic slinky, the P MS-PS3-5 • Waves and Their Applications in boundary due to the change in material response. Technologies for Information Transfer: SEE question 2 on Page 4 on the worksheet for this activity. MS-PS4-1, HS-PS4-1, MS-PS4-2, HS-PS4-5 Resources on this DVD & Internet relevant to Build a Better Wall VIDEOS —Demonstrations using the Slinky© are in the 2. Earthquakes & Tsunamis folder on this DVD: 3. VIDEOS_Earthquake & Tsunami > DEMO_SlinkySeismicWaves_Groom.mov Roger Groom’s class demo: http://www.youtube.com/watch?v=KZaI4MEWdc4 ANIMATIONS —Watch the following groups of animations in the 2. Earthquakes & Tsunamis folder : 2. ANIMATIONS_Earthquake & Tsunami > Seismic Wave Behavior_1station Seismic Wave Behavior-Travel time Seismic Wave Motion-Brailes INTERNET —Other descriptions of this activity: http://www.exo.net/~pauld/summer_institute/summer_day10waves/wavetypes.html
    [Show full text]