Joint Inversion of Normal Mode and Body Wave Data for Inner Core Anisotropy 2

Total Page:16

File Type:pdf, Size:1020Kb

Joint Inversion of Normal Mode and Body Wave Data for Inner Core Anisotropy 2 JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. B12, 2380, doi:10.1029/2001JB000713, 2002 Joint inversion of normal mode and body wave data for inner core anisotropy 2. Possible complexities Miaki Ishii and Adam M. Dziewon´ski Department of Earth and Planetary Sciences, Harvard University, Cambridge, Massachusetts, USA Jeroen Tromp Seismological Laboratory, California Institute of Technology, Pasadena, California, USA Go¨ran Ekstro¨m Department of Earth and Planetary Sciences, Harvard University, Cambridge, Massachusetts, USA Received 19 June 2001; revised 10 June 2002; accepted 8 July 2002; published 31 December 2002. [1] In this study we focus on inner core sensitive body wave data to investigate lateral and radial heterogeneity within the inner core. Normal mode data are used to constrain a global model and to monitor if complexities introduced by body wave data are compatible with mode measurements. In particular, we investigate the possibilities of an isotropic layer near the inner core boundary and large-scale variations in anisotropy such as hemispheric dependence proposed in studies based upon differential travel time data. Travel time data from distances between 130° and 140° require anisotropy near the inner core boundary. This evidence is supported by differential travel time data based upon diffracted waves, contrasting the previous inferences of isotropy at the surface of the inner core. Our experiments also show that variations at a hemispheric scale are not necessary and that the sources of apparent hemispheric differences can be localized. A comparison of differential and absolute travel time data suggests that differences in inferred inner core anisotropy models arise mainly from a few anomalous paths. These paths are responsible for the strong anisotropy which is characteristic of models based upon differential data. Assuming constant anisotropy in the inner core, we investigate the distribution of residual differential travel times at the entry and exit points of rays turning within the outer core. The results are consistent with existing models of structure near the core-mantle boundary. Placing the source of the difference between differential and absolute travel time data in the lowermost mantle gives a more satisfactory result than attempting to model it with a complex inner core. INDEX TERMS: 7207 Seismology: Core and mantle; 7260 Seismology: Theory and modeling; 7203 Seismology: Body wave propagation; 7255 Seismology: Surface waves and free oscillations; KEYWORDS: inner core, anisotropy, PKIKP, differential travel times, hemispheric variations Citation: Ishii, M., A. M. Dziewon´ski, J. Tromp, and G. Ekstro¨m, Joint inversion of normal mode and body wave data for inner core anisotropy, 2, Possible complexities, J. Geophys. Res., 107(B12), 2380, doi:10.1029/2001JB000713, 2002. 1. Introduction differential travel time data is the similarity of the mantle ray path of waves turning in the outer core (BC or AB) and the [2] The properties of the Earth’s inner core have been inner core (DF) (Figure 1). This similarity allows one to investigated using three different types of data: normal modes assume that differential measurements do not include effects [e.g., Woodhouse et al., 1986; Tromp,1993;Durek and due to the source, receiver, or mantle, and hence that the Romanowicz, 1999], absolute travel times of PKIKP (or travel time anomalies are solely due to the path within the PKP ) [e.g., Morelli et al., 1986; Shearer et al., 1988; DF inner core. Many models of the inner core have been Su and Dziewon´ski, 1995], and the differential travel proposed using differential travel time data, some of which times BC DF (or PKP PKP )andAB DF (or À BCÀ DF À contain structures not seen in models based upon other types PKP PKP ) [e.g., Creager, 1992; Song and Helmberger, ABÀ DF of data, such as an isotropic layer near the inner core 1993; Vinnik et al., 1994; McSweeney et al., 1997; Tanaka boundary [e.g., Song and Helmberger, 1998], a transition and Hamaguchi, 1997]. The advantage of the high-quality zone within the inner core [e.g., Song and Helmberger, 1998], and hemispherically dependent anisotropy [e.g., Copyright 2002 by the American Geophysical Union. Tanaka and Hamaguchi, 1997; Creager, 1999; Niu and 0148-0227/02/2001JB000713$09.00 Wen, 2001]. ESE 21 - 1 ESE 21 - 2 ISHII ET AL.: JOINT INVERSION FOR IC ANISOTROPY, 2 Figure 1. (left) Ray paths for AB, BC, and DF at an epicentral distance range of 150°, and (right) a travel time table for PKP branch using PREM [Dziewon´ski and Anderson, 1981]. For the travel time table we use an earthquake with a source depth of 300 km. [3] In the companion paper (Ishii et al. [2002], hereinafter Vinnik et al., 1994; Song, 1996; McSweeney et al., 1997], referred to as paper 1), we demonstrate that it is relatively whereas absolute travel times tend to result in models with easy to obtain a simple anisotropic inner core model which weaker anisotropy [e.g., Su and Dziewon´ski, 1995]. To fits all three types of data. Normal mode and mantle- investigate this discrepancy, we compare absolute and corrected DF data are compatible with one another and differential travel time data. Contamination of differential good fits can be achieved using a constant anisotropy data due to mantle structure near the core-mantle boundary model. Differential travel time data are harder to fit. In this (CMB) has also been suggested [e.g., Breger et al., 1999, paper we explore the source of this difficulty by (1) 2000], especially for the AB rays. We follow this suggestion considering in detail paths sampled by differential data; and investigate if the residual differential travel time signal (2) introducing structure within the inner core; and (3) based upon a constant anisotropy model can be placed in the investigating the possible effects of D" heterogeneity. Our mantle. focus is on large-scale structure, and therefore the existence of lateral variations [e.g., Su and Dziewon´ski, 1995], except for hemispheric variations or small-scale structure [e.g., 2. Comparison of Body Wave Data Creager, 1997; Dziewon´ski, 2000; Vidale and Earle, [6] We use normal mode splitting functions [He and 2000], is addressed only briefly. We assume that the Tromp, 1996; Resovsky and Ritzwoller, 1998; Durek and symmetry axis of the inner core is aligned with the rotation Romanowicz, 1999], ISC DF residuals, and high-quality axis, and we do not consider differential rotation [e.g., differential travel time data [Creager, 1992; Vinnik et al., Shearer and Toy, 1991; Creager, 1992; Song and Richards, 1994; McSweeney, 1995; Song, 1996; Creager, 1997; 1996; Su et al., 1996; Souriau and Poupinet, 2000] in our McSweeney et al., 1997; Tanaka and Hamaguchi, 1997; analysis. A detailed discussion of the theory, data, and Creager, 1999] to investigate the structure of the inner core. inversion method used in this paper can be found in paper 1. Before using these data, we explore the differences between [4] Using all available data, we investigate if the intro- mantle-corrected absolute and differential travel time data. duction of an isotropic layer near the inner core boundary For details of data processing, see paper 1. If the mantle (ICB) can improve the fit, especially to differential travel correction to absolute and differential travel times is appro- time data [e.g., Song and Helmberger, 1995; Creager, priate, these two sets of data should agree well. The fact that 2000]. We also explore the possibility of hemispherically models obtained from independent inversions of these two distinct anisotropy, previously proposed on the basis of data sets differ significantly suggest that this is not the case. body wave measurements [e.g., Tanaka and Hamaguchi, [7] The BC arrivals are observed in the 145° to 153° 1997; Creager, 1999; Niu and Wen, 2001]. Normal modes epicentral distance range (Figure 1), and the corresponding do not constrain this type of structure, because they are only BCÀDF data are sensitive to inner core structure down to sensitive to even-degree anomalies. Focusing on absolute roughly 300 km beneath the ICB. There are 851 BCÀDF body wave data, we investigate if large-scale variations, measurements in our data set, spanning the distance from including hemispheric variation, in inner core anisotropy 145° to 160°. Data with an epicentral distance greater than can be observed consistently. We use both normal mode and 153° involve diffracted BC traveling along the ICB. We body wave data to constrain the global average and consider discard these data when performing inversions but return to hemispheric deviations as perturbations from this average. them in section 4, where we consider structure near the ICB. [5] In general, inner core models derived using absolute The remaining 512 BCÀDF data are binned into two and differential travel time data differ most in their level of distance ranges, 145° to 150° and 150° to 153°,and anisotropy. Models based upon differential travel times have averaged for each 0.1 increment in cos2 x, where x is the strong anisotropy [e.g., Creager, 1992; Song and Helm- angle between the direction of wave propagation and the berger, 1993], especially near the center of the Earth [e.g., symmetry axis of transverse isotropy. The averaging process ISHII ET AL.: JOINT INVERSION FOR IC ANISOTROPY, 2 ESE 21 - 3 reduces the effects of mantle heterogeneity by cancelling out the isotropic signal if the outer core entry and exit positions of BC or AB cover the globe uniformly. In this paper we show that differential travel time data coverage is uneven compared to absolute travel time data, and suggest biased sampling as the reason for difference in inner core models derived from these two data sets. [8] We compare BCÀDF data in the range from 150° to 153° with DF from the same distance range (Figure 2a).
Recommended publications
  • Adam Marian Dziewoński Was Spread Upon the Permanent Records of the Faculty
    At a meeting of the FACULTY OF ARTS AND SCIENCES on November 7, 2017, the following tribute to the life and service of the late Adam Marian Dziewoński was spread upon the permanent records of the Faculty. ADAM MARIAN DZIEWOŃSKI BORN: November 15, 1936 DIED: March 1, 2016 Adam Marian Dziewoński was born in Lwów, Poland (now Lviv, Ukraine) and received his M.S. in 1960 from the University of Warsaw and his Doctor of Technical Sciences in 1965 from the Academy of Mines and Metallurgy in Krakow. He moved to the United States in 1965 as a research associate at the Graduate Research Center of the Southwest in Dallas and became an assistant professor at University of Texas at Dallas in 1969. Adam moved to Harvard University as an associate professor in 1972, became a full professor in 1976, and was named the Frank B. Baird, Jr., Professor of Science in 1994. He was the chair of the Department of Geological Sciences between 1982 and 1986 (Department of Earth and Planetary Sciences thereafter). He retired in 2009 but remained actively engaged in research until his death. Adam was born in a year of seismological significance, the year in which the Earth’s inner core was discovered. As if fated, Adam was involved in all subsequent major discoveries about this deepest layer of Earth. With James Freeman Gilbert, Adam provided the first direct evidence that the inner core is solid, settling the debate over its nature. He continued his work on the inner core with his colleagues at Harvard and went on to show that it is anisotropic, rotates at a slightly different rate than the mantle, and has another layer at its center.
    [Show full text]
  • DOI: 10.1126/Science.285.5431.1231 , 1231 (1999); 285 Science Miaki
    Normal-Mode and Free-Air Gravity Constraints on Lateral Variations in Velocity and Density of Earth's Mantle Miaki Ishii and Jeroen Tromp Science 285, 1231 (1999); DOI: 10.1126/science.285.5431.1231 This copy is for your personal, non-commercial use only. If you wish to distribute this article to others, you can order high-quality copies for your colleagues, clients, or customers by clicking here. Permission to republish or repurpose articles or portions of articles can be obtained by following the guidelines here. The following resources related to this article are available online at www.sciencemag.org (this information is current as of August 27, 2012 ): Updated information and services, including high-resolution figures, can be found in the online version of this article at: on August 27, 2012 http://www.sciencemag.org/content/285/5431/1231.full.html This article cites 23 articles, 1 of which can be accessed free: http://www.sciencemag.org/content/285/5431/1231.full.html#ref-list-1 This article has been cited by 155 article(s) on the ISI Web of Science This article has been cited by 13 articles hosted by HighWire Press; see: http://www.sciencemag.org/content/285/5431/1231.full.html#related-urls www.sciencemag.org This article appears in the following subject collections: Geochemistry, Geophysics http://www.sciencemag.org/cgi/collection/geochem_phys Downloaded from Science (print ISSN 0036-8075; online ISSN 1095-9203) is published weekly, except the last week in December, by the American Association for the Advancement of Science, 1200 New York Avenue NW, Washington, DC 20005.
    [Show full text]
  • Normal-Mode and Free-Air Gravity Constraints on Lateral Variations In
    R ESEARCH A RTICLES ficient quality. There are density models based on current and past plate motions in which the Normal-Mode and Free-Air sources of heterogeneity are cold subducted slabs (5), and an upper mantle density model, Gravity Constraints on Lateral containing strong high-degree components, has been determined with surface wave data (6). Variations in Velocity and A density model of the mantle is im- portant to geophysics because mantle-flow calculations cannot be performed without Density of Earth’s Mantle it (7, 8). However, a whole-mantle den- Miaki Ishii* and Jeroen Tromp sity model is usually obtained by scaling an S velocity model, which is justifiable With the use of a large collection of free-oscillation data and additional con- if the heterogeneity has a single cause straints imposed by the free-air gravity anomaly, lateral variations in shear that affects both, for example, variations in velocity, compressional velocity, and density within the mantle; dynamic to- temperature. pography on the free surface; and topography on the 660-km discontinuity and In recent years, mainly as a result of the 9 the core-mantle boundary were determined. The velocity models are consistent June 1994 Bolivia earthquake, the quantity and with existing models based on travel-time and waveform inversions. In the quality of normal-mode, or free-oscillation, lowermost mantle, near the core-mantle boundary, denser than average ma- data have improved substantially (9), allowing terial is found beneath regions of upwellings centered on the Pacific Ocean and us to constrain the long-wavelength structure of Africa that are characterized by slow shear velocities.
    [Show full text]
  • For Immediate Release Contact: Nan Broadbent 408-431-9885
    For Immediate Release Contact: Nan Broadbent 408-431-9885, [email protected] Miaki Ishii honored with Richter Early Career Award Harvard professor recognized by Seismological Society of America SAN FRANCISCO -- Since her first forays into studying earthquakes, Miaki Ishii has successfully tackled seismological problems that would likely intimidate even the most seasoned researchers. In her young career, Ishii, an assistant professor of earth and planetary science at Harvard University, has made two groundbreaking discoveries in geophysics that have fostered intense debate and subsequent research that has changed the understanding of deep Earth seismology. For her work, the Seismological Society of America (SSA) will honor the 33-year-old Ishii with its Charles F. Richter Early Career Award, which honors outstanding contributions to the goals of the Society by a member early in her or his career. She will receive the award at a special ceremony at the SSA annual meeting on April 8 in Monterey, CA. Since she entered Harvard to begin her doctoral work in the late 1990s, Ishii has shown a knack for answering big questions. Shortly after her arrival at Harvard, she researched the driving force behind plate tectonics -- lateral variations in mantle density. The research inferred that the slowest parts of the lowermost mantle are denser than average, rather than lighter as most had assumed. The findings flew in the face of the long-held theory of a homogenized mantle and generated significant subsequent research and debate. Recent research is beginning to confirm Ishii’s observations. Her second groundbreaking find built on her previous study with Harvard’s Adam Dziewonski.
    [Show full text]
  • 15 March 2016 Magazine
    VOL. 97 NO. 6 15 MAR 2016 WHERE CURIOSITY HAS TAKEN US Building Better Climate Models No Tolerance for Sexual Harassment Cheers, Yawns for Planet Nine Have an Idea or Topic to Present at the Fall Meeting? Submit a Session Proposal Deadline: 20 April, 11:59 P.M. EDT You must be current in your 2016 membership dues in order to propose a session. fallmeeting.agu.org Earth & Space Science News Contents 15 MARCH 2016 PROJECT UPDATE VOLUME 97, ISSUE 6 22 Better Tools to Build Better Climate Models A Department of Energy collaboration aims to make climate model development faster and more efficient by creating a prototype of a system for testing model components. OPINION 9 Steps to Building a No-Tolerance Culture for Sexual Harassment AGU can lead the way in building an environment where sexual and other types of harassment have no place. 16 AGU NEWS COVER 2017 Fall Meeting 26 to Go to New Orleans, 2018 to Washington, D. C. Where Curiosity Has Taken Us AGU opts to relocate its 26,000-plus- attendee meeting from San Francisco for The Curiosity rover, one of NASA’s flagship missions, analyzes 2 years to avoid detrimental impacts from Martian geology, geochemistry, climatology, and radiation to assess extensive construction at the convention whether Mars could have supported microbial life. center. Earth & Space Science News Eos.org // 1 Contents DEPARTMENTS Editor in Chief Barbara T. Richman: AGU, Washington, D. C., USA; eos_ [email protected] Editors Christina M. S. Cohen Wendy S. Gordon Carol A. Stein California Institute Ecologia Consulting, Department of Earth and of Technology, Pasadena, Austin, Texas, USA; Environmental Sciences, Calif., USA; wendy@ecologiaconsulting University of Illinois at cohen@srl .caltech.edu .com Chicago, Chicago, Ill., USA; [email protected] José D.
    [Show full text]
  • SPRING 2016 TIMES Earthharvard UNIVERSITY DEPARTMENT&Planetary of EARTH and PLANETARY SCIENCES
    SPRING 2016 TIMES EarthHARVARD UNIVERSITY DEPARTMENT&Planetary OF EARTH AND PLANETARY SCIENCES R/V Knorr at sunset in the Atlantic. This vessel was used during a recent five-week MARPEX cruise to collect samples from the Mid-Atlantic Ridge. Earthquakes, movies, sailing: Freshman Seminars not CEDRIC HAMELIN typical Harvard fare 4 Exploring Earth’s Final Frontier What mid-ocean ridges reveal about tectonics and magmatism BY JOCELYN FUENTES, THIRD-YEAR GRADUATE STUDENT Three siblings reflect on hen most people think of the ocean, tonic boundaries, ocean ridges are an essential part 13family traditions—and EPS they imagine a vast expanse of calm of the Earth system where the seafloor is constantly blue water with little thought for what being created through magmatism and faulting. liesW below, or they think of the water that fills the Since the first map of the ocean floor was cre- ocean basins. For many ocean scientists, however, ated, in 1977 by Marie Tharp and Bruce Heezen, what is most interesting lies beneath the water on many advances have been made in the field of the seafloor. Two-thirds of the geology of the marine geology and geophysics. In particular, our planet is there, obscured by the water above. The ability to map and image the seafloor has improved Rock Hard Bodies play most prominent feature of the seafloor is the greatly. Unfortunately, these techniques generally Meanhouse Effect in first mid-ocean ridge system, one of Earth’s most active require a ship, which is costly and time-consuming. department Earth19 Bowl tectonic, volcanic, and biologic features.
    [Show full text]
  • Final Technical Report-Cover Page
    Final Technical Report Award Number G10AP00086 Sponsors Name: Seismological Society of America Conference Title: Seismological Society of America’s 2010 Annual Meeting in Portland, Oregon Conference Dates: April 21‐23, 2010 Summary of Program: This year’s SSA Annual Meeting included 282 oral presentations and 258 poster presentations distributed among 31 sessions. The meeting was attended by 557 participants from 26 countries. The meeting was preceded by the annual board meeting and committee meetings. Selected results of those meetings were shared with the membership at the annual luncheon, the first day of the meeting. Special presentations were given at the two other luncheons. The Joyner Lecture, a special lecture focusing on the boundary of seismology and engineering, was presented on Thursday night. The attached PDF contains the following: ‐ The program schedule (pages 2‐32) ‐ The meeting report including the annual luncheon minutes, the board minutes, awards presentations, and treasurer’s report (pages 33‐62) ‐ A list of registrants (Pages 63‐73) ‐ All the meeting abstracts (Pages 74‐187) SSA Annual meeting program SSA 2010 Annual Meeting Announcement Seismological Society of America Technical Sessions 21–23 April 2010 (Wednesday–Friday) Portland, Oregon IMPORTANT DATES Meeting Pre-registration Deadline 19 March 2010 Hotel Reservation Cut-Off 29 March 2010 Online Registration Cut-Off 9 April 2010 PROGRAM COMMITTEE TECHNICAL PROGRAM This year’s technical program committee is composed of co- Advances in Seismic Hazard Mapping chairs Seth Moran (USGS–Cascades Volcano Observatory) Recently much progress has been made in developing infor- and Nick Beeler (USGS–Vancouver), Ivan Wong (Seismic mative regional seismic hazard maps that incorporate local Hazards Group, Oakland CA), Ray Weldon (University of as well as regional geologic and geotechnical data.
    [Show full text]
  • Earthscope Usarray Facility Final Report
    EarthScope USArray Facility Final Report Award # EAR-0733069 Award Title Collaborative Research: EarthScope-Facility Operations and Maintenance (USArray) Award Term October 1, 2008 – September 30, 2013 Report Prepared on December 30, 2014 Report Prepared for National Science Foundation 4201 Wilson Boulevard Arlington, VA 22230 Report Prepared by Incorporated Research Institutions for Seismology 1200 New York Avenue, NW, Suite 400 Washington, DC 20005 Introduction The USArray component of the National Science Foundation’s EarthScope program facility consists of four major observatory components: (1) a Reference Network of permanent seismic stations; (2) a Transportable Array (TA) of ~400 seismic stations; (3) a Flexible Array (FA) pool of seismic instruments available to Principal Investigators (PIs); and (4) a Magnetotelluric (MT) observatory with permanently and temporarily deployable instruments. USArray also includes comprehensive data management and siting outreach efforts. USArray efforts during FY09–FY13 continue and expand facility activities established under the previous MREFC and O&M awards (EAR-0323309 and EAR-0323311, respectively) and continued under the “Seismological Facilities for the Advancement of Geoscience and EarthScope” award (EAR-1261681). Quarterly and annual reports submitted to the National Science Foundation throughout the award document USArray progress and accomplishments. This present report provides an overview of programmatic activities conducted under EAR- 0733069. Facility Information The Incorporated Research Institutions for Seismology (IRIS) operates USArray. Under the governance structure in place during the O&M award, USArray management received community input via the IRIS Board of Directors, the USArray Advisory Committee, and the Transportable Array and Electromagnetic Working Groups (Figure 1). The organizational structure of the USArray has been designed to leverage the experience and facilities that already exist within IRIS.
    [Show full text]