The Effect of Regional Variation of Seismic Wave Attenuation on the Strong Ground Motion from Earthquakes

Total Page:16

File Type:pdf, Size:1020Kb

The Effect of Regional Variation of Seismic Wave Attenuation on the Strong Ground Motion from Earthquakes NUREG/CR-1655 UCRL-53004 XA04NO347 fNIS-XA-N--048 The Effect of Regional Variation of Seismic Wave Attenuation on the Strong Ground Motion from Earthquakes D. H Chung and D. L. Bernreuter Prepared for U.S. Nuclear Regulatory Commission I-AMIRENCE UVERMORE LABORATORY NOTICE This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, or any of their employees, makes any warranty, expressed or implied, or assumes any legal liability or responsibility for any third party's use, or the results of such use, of any information, apparatus product or process disclosed in this report, or represents that its use by such third party would not infringe privately owned rights. Available from GPO Sales Program Division of Technical Information and Document Control U. S. Nuclear Regulatory Commission Washington, D. C. 20555 Printed copy price: 425 and National Technical Information Service Springfield, Virginia 22161 NUREG/CR-1655 UCRL-53004 RM, RD, RA, R6 The Effect of Regional Variation of Seismic Wave Attenuation on the Strong Ground Motion from Earthquakes Manuscript Completed: July 1981 Date Published: October 1981 Prepared by D. H. Chung and D. L. Bernreuter Lawrence Livermore Laboratory 7000 East Avenue Livermore, CA 94550 Prepared for Office of Nuclear Regulatory Research U.S. Nuclear Regulatory Commission Washington, D.C. 20555 NRC FIN No. A0139 Availability of Reference Materials Cited in NRC Publications Most documents cited in NRC publications will be available from one of the following sources: 1 . The NRC Public Document Room, 1717 H Street., N.W. Washington, DC 20555 2 The NRC/GPO Sales Program, U.S. Nuclear Regulatory Commission, Washington, DC 20555 3. The National Technical Information Service, Springfield, VA 22161 Although the listing that follows represents the majority of documents cited in NRC publications, it is not intended to be exhaustive. Referenced documents available for inspection and copying for a fee from the NRC Public Document Room include NRC correspondence and internal NRC memoranda-, NRC Office of Inspection and Enforce- ment bulletins, circulars, information notices, inspection and investigation notices, Licensee Event Reports; vendor reports and correspondence; Commission papers; and applicant and licensee documents and correspondence. The following documents in the NUREG series are available for purchase from the NRC/GPO Sales Pro- gram: formal NRC staff and contractor reports, NRC-sponsored conference proceedings, and NRC booklets and brochures. Also available are Regulatory Guides, NRC regulations in the Code of Federal Regulations, and NuclearRegulatory Commission Issuances. Documents available from the National Technical Information Service include NUREG series reports and technical reports prepared by other federal agencies and reports prepared by the Atomic Energy Commis- sion, forerunner agency to the Nuclear Regulatory Commission. Documents available from public and special technical libraries include all open literature items, such as books, journal and periodical articles, transactions, and codes and standards. FederalRegister notices, federal and state legislation, and congressional reports can usually be obtained from these libraries. Documents such as theses, dissertations, foreign reports and translations, and non-NRC conference pro- ceedings are available for purchase from the organization sponsoring the publication cited. Single copies of NRC draft reports are available free upon written request to the Division of Technical Infor- mation and Document Control, U.S. Nuclear Regulatory Commission, Washington, DC 20555. ABSTRACT Attenuation is caused by geometric spreading and absorption. Geometric spreading is almost independent of crustal geology and physiographic region, but absorption depends strongly on crustal geology and the state of the earth's upper mantle. Except for very high frequency waves, absorption does not affect ground motion at distances less than about 25 to 50 km. Thus, in the near-field zone, the attenuation in the eastern United States is similar to that in the western United States. Beyond the near field, differences in ground motion can best be accounted for by differences in attenuation caused hy differences in absorption. The stress drop of eastern earthquakes may be higher than for western earthquakes of the same seismic moment, which would affect the high-frequency spectral content. But we believe this factor is of much less significance than differences in absorption in explaining the differences in ground motion between the East and the West. The characteristics of strong ground motion in the conterminous United States are discussed in light of these considerations, and estimates are made of the epicentral qround motions in the central and eastern United States. iii CONTENTS Abstract . Figures . vi Tables . vii Foreword . ix Glossary of Symbols . xi Executive Summary . 1 Introduction . 3 Evidence for Regionalization . 7 Evidence from P- and S-Wave Velocities . 7 Evidence from Surface-Wave Studies . 14 Implications of Regionalization: Discussion . 23 Strong Ground Motion Characteristics in the Conterminous United States . 32 Inferring EUS Attenuations: Four Possible Methods . 33 Method 1 . 36 Method 2 . 36 Method 3 . 37 Method 4 . 38 Examples and Discussion . 39 Estimates of the Epicentral Ground Motion in the Central and Eastern United States . 46 Approach and Results . 47 Comparison with Nuttli's Analysis . 55 Appendix: Published Attenuation Functions . 57 References . 65 v FIGURES 1. Typical ground motion attenuation . 4 2. The Pn-velocity distribution in the conterminous United States . 8 3. Contour map of mean teleseismic P-wave station corrections for the conterminous United States . 9 4. Seismic-wave transmission characteristics in North America . 12 5. Seismograms along three profiles from the Salmon nuclear explosion . 15 6. Comparison of the upper-mantle Q models for the eastern and western United States . 21 7. Contours of upper mantle a from Pn paths in the conterminous United States . 22 8. Peak particle velocity as a function of distance from a salvo-type underground nuclear explosion source (Buggy) . 25 9. Attenuation of vertical ground velocity from the 1971 San Fernando earthquake . 26 10. Comparison of vertical components of velocity from earthquakes in the western United States with magnitudes between 4.5 and 5.8 . 27 11. Variation of amplitude ratio with epicentral distance for three frequencies . 31 12. The intensity attenuation relationships for five historical events in the EUS . 40 13. A comparison between the results of the four attenuation models, for an event with b = 65 . 42 Vi 14. Measured peak accelerations from the Rio Blanco event as a function of distance . 49 15. Estimated peak horizontal acceleration for the 1811-1812 New Madrid earthquake series, based on the 1968 southern Illinois event . 53 TABLES 1. Values for the constants A. in Eq. 14) . 39 1 2. Peak ground accelerations (in cm/s2) calculated for three body-wave magnitudes, using methods 3 and 4 . 44 3. Attenuation coefficients for selected events . 54 vii FOREWORD This report was prepared as part of Phase I of the Seismic Safety Margins Research Program at the Lawrence Livermore National Laboratory. It was supported by the U.S. Nuclear Regulatory Commission under a memorandum of understanding with the U.S. Department of Energy. ix GLOSSARY OF SYMBOLS a Wave amplitude; exponential constant in Eq. 12) A Wave amplitude in Eq. 12) ah Peak horizontal ground acceleration C(z) Factor to account for the variation of Q with depth f Wave frequency GMP Ground motion parameter I0 Epicentral intensity Is Site intensity L Length of salvo explosion array M Generalized earthquake magnitude mb Body-wave magnitude ML Local magnitude PGA Peak ground acceleration PGV Peak ground velocity PSA Peak spectral amplitude Generalized seismic quality factor a Seismic quality factor for P-waves Seismic quality factor for S-waves r Hypocentral distance R Epicentral distance; hypocentral or epicentral distance in Eq. 4) s Ray path coordinate (p. 11) S Factor to account for geometric spreading V Group velocity z Depth a Compressional wave velocity Shear-wave velocity Coefficient of anelastic attenuation Epicentral distance Exponent in attenuation equation on page 5 W Angular frequency xi EXECUTIVE SUMMARY Attenuation is the decay in the amplitude of earthquake motion--motion generated by sudden release of the earth's strain about the earthquake source. The key factors influencing attenuation are the source conditions, the transmission path characteristics, and the local site conditions. The source conditions influence the initial level and frequency content of the motion. The transmission path characteristics influence the path length and the rate of amplitude decay. The local site conditions further influence the level and frequency content of the motion. The attenuation of ground motion from an earthquake is complex in nature, owing in part to the complexity of the ground motion, which consists of a combination of harmonic motions of varying frequency. The frequency dependency of typical ground motion attenuation shows that higher frequency motions attenuate faster than lower frequency motions. The composition of frequencies and the percentage of strong motion caused by different waves i.e., P waves, SH and SV waves, surface waves, etc.) undergo changes along the transmission path. During
Recommended publications
  • Seismic Wavefield Imaging of Earth's Interior Across Scales
    TECHNICAL REVIEWS Seismic wavefield imaging of Earth’s interior across scales Jeroen Tromp Abstract | Seismic full- waveform inversion (FWI) for imaging Earth’s interior was introduced in the late 1970s. Its ultimate goal is to use all of the information in a seismogram to understand the structure and dynamics of Earth, such as hydrocarbon reservoirs, the nature of hotspots and the forces behind plate motions and earthquakes. Thanks to developments in high- performance computing and advances in modern numerical methods in the past 10 years, 3D FWI has become feasible for a wide range of applications and is currently used across nine orders of magnitude in frequency and wavelength. A typical FWI workflow includes selecting seismic sources and a starting model, conducting forward simulations, calculating and evaluating the misfit, and optimizing the simulated model until the observed and modelled seismograms converge on a single model. This method has revealed Pleistocene ice scrapes beneath a gas cloud in the Valhall oil field, overthrusted Iberian crust in the western Pyrenees mountains, deep slabs in subduction zones throughout the world and the shape of the African superplume. The increased use of multi- parameter inversions, improved computational and algorithmic efficiency , and the inclusion of Bayesian statistics in the optimization process all stand to substantially improve FWI, overcoming current computational or data- quality constraints. In this Technical Review, FWI methods and applications in controlled- source and earthquake seismology are discussed, followed by a perspective on the future of FWI, which will ultimately result in increased insight into the physics and chemistry of Earth’s interior.
    [Show full text]
  • 2D Seismic Survey in Block AD- 10, Offshore Myanmar
    2D Seismic Survey in Block AD- 10, Offshore Myanmar Initial Environmental Examination 02 December 2015 Environmental Resources Management www.erm.com The world’s leading sustainability consultancy 2D Seismic Survey in Block AD-10, Environmental Resources Management Offshore Myanmar ERM-Hong Kong, Limited 16/F, Berkshire House 25 Westlands Road Initial Environmental Examination Quarry Bay Hong Kong Telephone: (852) 2271 3000 Facsimile: (852) 2723 5660 Document Code: 0267094_IEE_Cover_AD10_EN.docx http://www.erm.com Client: Project No: Statoil Myanmar Private Limited 0267094 Summary: Date: 02 December 2015 Approved by: This document presents the Initial Environmental Examination (IEE) for 2D Seismic Survey in Block AD-10, as required under current Draft Environmental Impact Assessment Procedures Craig A. Reid Partner 1 Addressing MOECAF Comments, Final for MOGE RS CAR CAR 02/12/2015 0 Draft Final RS JNG CAR 31/08/2015 Revision Description By Checked Approved Date Distribution Internal Public Confidential CONTENTS 1 EXECUTIVE SUMMARY 1-1 1.1 PURPOSE AND EXTENT OF THE IEE REPORT 1-1 1.2 SUMMARY OF THE ACTIVITIES UNDERTAKEN DURING THE IEE STUDY 1-2 1.3 PROJECT ALTERNATIVES 1-2 1.4 DESCRIPTION OF THE ENVIRONMENT TO BE AFFECTED BY THE PROJECT 1-4 1.5 SIGNIFICANT ENVIRONMENTAL IMPACTS 1-5 1.6 THE PUBLIC CONSULTATION AND PARTICIPATION PROCESS 1-6 1.7 SUMMARY OF THE EMP 1-7 1.8 CONCLUSIONS AND RECOMMENDATIONS OF THE IEE REPORT 1-8 2 INTRODUCTION 2-1 2.1 PROJECT OVERVIEW 2-1 2.2 PROJECT PROPONENT 2-1 2.3 THIS INITIAL ENVIRONMENTAL EVALUATION (IEE)
    [Show full text]
  • Chapter 2 the Evolution of Seismic Monitoring Systems at the Hawaiian Volcano Observatory
    Characteristics of Hawaiian Volcanoes Editors: Michael P. Poland, Taeko Jane Takahashi, and Claire M. Landowski U.S. Geological Survey Professional Paper 1801, 2014 Chapter 2 The Evolution of Seismic Monitoring Systems at the Hawaiian Volcano Observatory By Paul G. Okubo1, Jennifer S. Nakata1, and Robert Y. Koyanagi1 Abstract the Island of Hawai‘i. Over the past century, thousands of sci- entific reports and articles have been published in connection In the century since the Hawaiian Volcano Observatory with Hawaiian volcanism, and an extensive bibliography has (HVO) put its first seismographs into operation at the edge of accumulated, including numerous discussions of the history of Kīlauea Volcano’s summit caldera, seismic monitoring at HVO HVO and its seismic monitoring operations, as well as research (now administered by the U.S. Geological Survey [USGS]) has results. From among these references, we point to Klein and evolved considerably. The HVO seismic network extends across Koyanagi (1980), Apple (1987), Eaton (1996), and Klein and the entire Island of Hawai‘i and is complemented by stations Wright (2000) for details of the early growth of HVO’s seismic installed and operated by monitoring partners in both the USGS network. In particular, the work of Klein and Wright stands and the National Oceanic and Atmospheric Administration. The out because their compilation uses newspaper accounts and seismic data stream that is available to HVO for its monitoring other reports of the effects of historical earthquakes to extend of volcanic and seismic activity in Hawai‘i, therefore, is built Hawai‘i’s detailed seismic history to nearly a century before from hundreds of data channels from a diverse collection of instrumental monitoring began at HVO.
    [Show full text]
  • ADUA Azerbaijan 2D-3D Seismic Survey Environmental Impact
    Environmental Impact Assessment (EIA) for 2D-3D Doc. No. seismic survey in the Ashrafi-Dan Ulduzu-Aypara (ADUA) Exploration area, Azerbaijan Valid from Rev. no. 0 01.03.2019 Environmental Impact Assessment (EIA) for 2D-3D seismic survey in the Ashrafi-Dan Ulduzu-Aypara (ADUA) Exploration area, Azerbaijan March 2019 Environmental Impact Assessment (EIA) for 2D-3D seismic survey in the Ashrafi-Dan Ulduzu-Aypara (ADUA) Exploration area, Azerbaijan Valid from 01.03.2019 Rev. no. 0 Table of contents Acronyms ...................................................................................................................................................... 10 Executive Summary ...................................................................................................................................... 13 Regulatory Framework .................................................................................................................................... 14 Project Description .......................................................................................................................................... 17 Description of the Environmental and Social Baseline ................................................................................... 18 Impact Assessment and Mitigation ................................................................................................................. 22 Environmental Management Plan ..................................................................................................................
    [Show full text]
  • Types of Seismic Energy Sources for Petroleum Exploration in Desert, Dry-Land, Swamp and Marine Environments in Nigeria and Other Sub- Saharan Africa
    International Journal of Science and Research (IJSR) ISSN (Online): 2319-7064 Index Copernicus Value (2015): 78.96 | Impact Factor (2015): 6.391 Types of Seismic Energy Sources for Petroleum Exploration in Desert, Dry-Land, Swamp and Marine Environments in Nigeria and Other Sub- Saharan Africa Madu Anthony Joseph Chinenyeze Department of Geology, College of Physical And Applied Sciences, Michael Okpara University of Agriculture Umudike, Abia State, Nigeria Abstract: The various seismic energy sources used for petroleum exploration in Nigeria (Niger Delta, Benue Trough, Chad Basin) and sub-Saharan Africa (Niger Republic, Sudan, Ethiopia inclusive) have been undergoing improvement over the years. Explosive source usually Dynamite, Airgun, Thumper/Weight Drop, and Vibrators are remarkable seismic energy sources today with frequencies for primary seismic production, bandwidth 10Hz – 75Hz. These sources contribute to higher S/N ratio and greater bandwidth. The spectral analysis of shot records from the afore-mentioned energy sources confirm their corresponding suitability, effectiveness, and adequacy for wave propagation into the ground, until it encountered an impedance in the subsurface. The latter condition results in reflections and refractions which are picked by the receivers or sensors (geophones or hydrophones). The spectral analyses from these selected sources reveal comparable range of signal frequencies. The Thumpers or vibrator shots at any given shot point location SP, is repeated severally for the purpose of enhancement and then summed or stacked. This is also called vibrator sweep, as it vibrates repeatedly according to specification of program issue. In the repetition of frequency sweep, amplitude increases are optimized to a specific level within the same duration.
    [Show full text]
  • PEAT8002 - SEISMOLOGY Lecture 13: Earthquake Magnitudes and Moment
    PEAT8002 - SEISMOLOGY Lecture 13: Earthquake magnitudes and moment Nick Rawlinson Research School of Earth Sciences Australian National University Earthquake magnitudes and moment Introduction In the last two lectures, the effects of the source rupture process on the pattern of radiated seismic energy was discussed. However, even before earthquake mechanisms were studied, the priority of seismologists, after locating an earthquake, was to quantify their size, both for scientific purposes and hazard assessment. The first measure introduced was the magnitude, which is based on the amplitude of the emanating waves recorded on a seismogram. The idea is that the wave amplitude reflects the earthquake size once the amplitudes are corrected for the decrease with distance due to geometric spreading and attenuation. Earthquake magnitudes and moment Introduction Magnitude scales thus have the general form: A M = log + F(h, ∆) + C T where A is the amplitude of the signal, T is its dominant period, F is a correction for the variation of amplitude with the earthquake’s depth h and angular distance ∆ from the seismometer, and C is a regional scaling factor. Magnitude scales are logarithmic, so an increase in one unit e.g. from 5 to 6, indicates a ten-fold increase in seismic wave amplitude. Note that since a log10 scale is used, magnitudes can be negative for very small displacements. For example, a magnitude -1 earthquake might correspond to a hammer blow. Earthquake magnitudes and moment Richter magnitude The concept of earthquake magnitude was introduced by Charles Richter in 1935 for southern California earthquakes. He originally defined earthquake magnitude as the logarithm (to the base 10) of maximum amplitude measured in microns on the record of a standard torsion seismograph with a pendulum period of 0.8 s, magnification of 2800, and damping factor 0.8, located at a distance of 100 km from the epicenter.
    [Show full text]
  • UNITED STATES DEPARTMENT of the INTERIOR BUREAU of OCEAN ENERGY MANAGEMENT SITE-SPECIFIC ENVIRONMENTAL ASSESSMENT GEOLOGICAL &Am
    GEOPHYSICAL EXPLORATION FOR MINERAL RESOURCES SEA NO. L19-009 UNITED STATES DEPARTMENT OF THE INTERIOR BUREAU OF OCEAN ENERGY MANAGEMENT GULF OF MEXICO OCS REGION NEW ORLEANS, LOUISIANA SITE-SPECIFIC ENVIRONMENTAL ASSESSMENT OF GEOLOGICAL & GEOPHYSICAL SURVEY APPLICATION NO. L19-009 FOR SHELL OFFSHORE INC. April 2, 2019 RELATED ENVIRONMENTAL DOCUMENTS GulfofMexico OCS Proposed Geological and Geophysical Activities Westem, Central, and Eastem Planning Areas Final Programmatic Environmental Impact Statement (OCS EIS/EA BOEM 2017-051) GulfofMexico OCS Oil and Gas Lease Sales: 2017-2022 GulfofMexico Lease Sales 249, 250, 251, 252, 253, 254, 256, 257, 259, and 261 Final Environmental Impact Statement (OCS EIS/EA BOEM 2017-009) Gulf ofMexico OCS Lease Sale Final Supplemental Environmental Impact Statement 2018 (OCS EIS/EA BOEM 2017-074) FINDING OF NO SIGNIFICANT IMPACT (FONSI) The Bureau of Ocean Energy Management (BOEM) has prepared a Site-Specific Environmental Assessment (SEA) (No. L19-009) complying with the National Environmental Policy Act (NEPA). NEPA regulations under the Council on Environmental Quality (CEQ) (40 CFR § 150 E3 and § 1508.9), die United States Department of the Interior (DOI) NEPA implementing regulations (43 CFR § 46), and BOEM policy require an evaluation of proposed major federal actions, which under BOEM jurisdiction includes approving a plan for oil and gas exploration or development activity on the Outer Continental Shelf (OCS). NEPA regulation 40 CFR § 1508.27(b) requires significance to be evaluated in terms of
    [Show full text]
  • Paper 4: Seismic Methods for Determining Earthquake Source
    Paper 4: Seismic methods for determining earthquake source parameters and lithospheric structure WALTER D. MOONEY U.S. Geological Survey, MS 977, 345 Middlefield Road, Menlo Park, California 94025 For referring to this paper: Mooney, W. D., 1989, Seismic methods for determining earthquake source parameters and lithospheric structure, in Pakiser, L. C., and Mooney, W. D., Geophysical framework of the continental United States: Boulder, Colorado, Geological Society of America Memoir 172. ABSTRACT The seismologic methods most commonly used in studies of earthquakes and the structure of the continental lithosphere are reviewed in three main sections: earthquake source parameter determinations, the determination of earth structure using natural sources, and controlled-source seismology. The emphasis in each section is on a description of data, the principles behind the analysis techniques, and the assumptions and uncertainties in interpretation. Rather than focusing on future directions in seismology, the goal here is to summarize past and current practice as a companion to the review papers in this volume. Reliable earthquake hypocenters and focal mechanisms require seismograph locations with a broad distribution in azimuth and distance from the earthquakes; a recording within one focal depth of the epicenter provides excellent hypocentral depth control. For earthquakes of magnitude greater than 4.5, waveform modeling methods may be used to determine source parameters. The seismic moment tensor provides the most complete and accurate measure of earthquake source parameters, and offers a dynamic picture of the faulting process. Methods for determining the Earth's structure from natural sources exist for local, regional, and teleseismic sources. One-dimensional models of structure are obtained from body and surface waves using both forward and inverse modeling.
    [Show full text]
  • Source, Scattering and Attenuation Effects on High
    SOURCE, SCATTERING AND ATTENUATION EFFECTS ON HIGH FREQUENCY SEISMIC vlAVES by BERNARD ALFRED CHOUET Electrical Engineer Federal Institute of Technology Lausanne, Switzerland (1968) S.M., Aeronautics and Astronautics, M.I.T. (1972) S.M., Earth and Planetary Sciences, M.I.T. (1973) SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY at the MASSACHUSETTS INSTITUTE OF TECHNOLOGY Sept.ember, 1976 Signature of the Author Department: of Earth and Planetary Sciences, September 1976 Certified by ................................... Thesis Supervisor Accepted by ..............-.................... Chairman, Departmental Committee on Graduate Students ii SOURCE, SCATTERING AND ATTENUATION EFFECTS ON HIGH FREQUENCY SEISMIC WAVES by BERNARD ALFRED CHOUET Submitted to the Department of Earth and Planetary Sciences on August 12, 1976 in partial fulfillment of the requirements for the degree of Doctor of Philosophy ABSTRACT High frequency coda waves from small local earthquakes are interpreted as backscattering waves from numerous hete- rogeneities distributed uniformly in the earth's crust. Two extreme models of the wave medium that account for the ob- servations on the coda are proposed. In the first model we use a simple statistical approach and consider the back- scattering waves as a superposition of independent singly scattered wavelets. The basic assumption underlying this single backscattering model is that the scattering is a weak process so that the loss of seismic energy as well as the multiple scattering can be neglected. In the second model the seismic energy transfer is considered as a diffusion process. These two approaches lead to similar formulas that allow an accurate separation of the effect of earth- quake source from the effects of scattering and attenuation on the coda spectra.
    [Show full text]
  • Modelling Seismic Wave Propagation for Geophysical Imaging J
    Modelling Seismic Wave Propagation for Geophysical Imaging J. Virieux, V. Etienne, V. Cruz-Atienza, R. Brossier, E. Chaljub, O. Coutant, S. Garambois, D. Mercerat, V. Prieux, S. Operto, et al. To cite this version: J. Virieux, V. Etienne, V. Cruz-Atienza, R. Brossier, E. Chaljub, et al.. Modelling Seismic Wave Propagation for Geophysical Imaging. Seismic Waves - Research and Analysis, Masaki Kanao, 253- 304, Chap.13, 2012, 978-953-307-944-8. hal-00682707 HAL Id: hal-00682707 https://hal.archives-ouvertes.fr/hal-00682707 Submitted on 5 Apr 2012 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. 130 Modelling Seismic Wave Propagation for Geophysical Imaging Jean Virieux et al.1,* Vincent Etienne et al.2†and Victor Cruz-Atienza et al.3‡ 1ISTerre, Université Joseph Fourier, Grenoble 2GeoAzur, Centre National de la Recherche Scientifique, Institut de Recherche pour le développement 3Instituto de Geofisica, Departamento de Sismologia, Universidad Nacional Autónoma de México 1,2France 3Mexico 1. Introduction − The Earth is an heterogeneous complex media from
    [Show full text]
  • Three New Fundamental Problems in Physics of the Aftershocks
    arXiv 2020 Three New Fundamental Problems in Physics of the Aftershocks Anatol V. Guglielmi, Alexey D. Zavyalov, Oleg D. Zotov Institute of Physics of the Earth RAS, Moscow, Russia [email protected], [email protected], [email protected] Abstract Recently, the physics of aftershocks has been enriched by three new problems. We will conditionally call them dynamic, inverse, and morphological problems. They were clearly formulated, partially resolved and are fundamental. Dynamic problem is to search for the cumulative effect of a round-the-world seismic echo appearing after the mainshock of earthquake. According to the theory, a converging surface seismic wave excited by the main shock returns to the epicenter about 3 hours after the main shock and stimulates the excitation of a strong aftershock. Our observations confirm the theoretical expectation. The second problem is an adequate description of the averaged evolution of the aftershock flow. We introduced a new concept about the coefficient of deactivation of the earthquake source, “cooling down” after the main mainshock, and proposed an equation that describes the evolution of aftershocks. Based on the equation of evolution, we posed and solved the inverse problem of the physics of the earthquake source and created an Atlas of aftershocks demonstrating the diversity of variations of the deactivation coefficient. The third task is to model the spatial and spatio-temporal distribution of aftershocks. Its solution refines our understanding of the structure of the earthquake source. We also discuss in detail interesting unsolved problems in the physics of aftershocks. Keywords: round-the-world echo, deactivation coefficient, evolution equation, inverse problem, spatio-temporal distribution of aftershocks 1 1.
    [Show full text]
  • The Story of Wolfspar®, BP's Low-Frequency Marine Seismic Source
    The story of Wolfspar®, BP’s low-frequency marine seismic source; or, “Sometimes the ‘brute force’ approach to a problem can work”! Joe Dellinger, Drew Brenders, John Etgen, Scott Michell In exploration seismology we image the Earth’s interior using reflected sound waves. We do this using an algorithm called “migration”, which converts the reflected sound waves recorded at the Earth’s surface into an image of subsurface structures. There is a catch-22 problem in doing this, however. For migration to work we need a reasonably good numerical model of the sound speeds in the Earth. We often find ourselves trying to guess what those sound speeds may be. If our guess is not too far wrong, it allows us to produce a degraded image, which we can then use to update our subsurface sound speed model and try migrating again, etc. Sometimes, however, we just can’t get this process off the ground. The standard response when confronted with such difficulties is to construct ever more sophisticated imaging / inversion algorithms. The holy grail would be an algorithm that could produce a good image with minimal to no human intervention or a-priori information required. Some of the best minds in the industry have been working on developing such algorithms for decades now, and while there has been some success, large areas of the US Gulf of Mexico have remained intractable. That’s a prize of many billions of dollars to the US economy that we still cannot access. It is well known that “the lower the frequency, the easier the problem becomes”.
    [Show full text]