Joint Mineral Physics and Seismic Wave Traveltime Analysis of Upper Mantle Temperature
Total Page:16
File Type:pdf, Size:1020Kb
Joint mineral physics and seismic wave traveltime analysis of upper mantle temperature Jeroen Ritsema1, Paul Cupillard2, Benoit Tauzin3, Wenbo Xu1, Lars Stixrude4, Carolina Lithgow-Bertelloni4 1Department of Geological Sciences, University of Michigan, Ann Arbor, Michigan 48109, USA 2Seismological Laboratory, University of California–Berkeley, Berkeley, California 94720, USA 3Ecole et Observatoire des Sciences de la Terre, Institut de Physique du Globe de Strasbourg, 67084 Strasbourg, France 4Department of Earth Sciences, University College London, London WC1E 6BT, UK ABSTRACT CONVERTED WAVE TRAVELTIMES We employ a new thermodynamic method for self-consistent Teleseismic P wave conversions, denoted as Pds, propagate as com- computation of compositional and thermal effects on phase transition pressional (P) waves through the deep mantle and convert to shear (S) depths, density, and seismic velocities. Using these profi les, we compare waves at an upper mantle discontinuity at depth d beneath the seismom- theoretical and observed differential traveltimes between P410s and eter. The phases P410s and P660s, for example, are waves that have P (T410) and between P600s and P410s (T660–410) that are affected only by converted from P to S at the 410 km and the 660 km transitions, respec- seismic structure in the upper mantle. The anticorrelation between T410 tively (Fig. 1). and T660–410 suggests that variations in T410 and T660–410 of ~8 s are due to Our data are recorded traveltime differences between the phases lateral temperature variations in the upper mantle transition zone of P410s and P (T410) and between P660s and P410s (T660–410). T410 and ~400 K. If the mantle is a mechanical mixture of basaltic and harzbur- T660–410 are unaffected by imprecise earthquake locations and lower mantle gitic components, our traveltime data suggest that the mantle has an hetero geneity. Although P410s and P660s are weak signals, they have pre- average temperature of 1600 ± 50 K, in agreement with temperature dictable slowness (i.e., angle of incidence at the surface). It is therefore estimates from magma compositions of mid-ocean ridge basalts. We straight forward to enhance P410s and P660s amplitudes by summation infer a 100 K hotter mantle if we assume the mantle to have a homoge- (stacking) of a large (>50–100) number of teleseismic (epicentral dis- neous pyrolitic composition. The transition-zone temperature beneath tances >4000 km) receiver functions from regional (Li et al., 2000; Shen hotspots and within subduction zones is relatively high and low, respec- et al., 1998; Simmons and Gurrola, 2000) or global (Chevrot et al., 1999; tively. However, the largest variability in T410 and T660–410 is recorded by Lawrence and Shearer, 2006; Tauzin et al., 2008) networks. global stations far from subduction zones and hotspots. This indicates We consider two sets of T410 and T660–410 measurements. The fi rst set that the 400 K variation in upper mantle temperature is complicated by is derived from ~7000 shallow (depth <50 km) earthquake recordings at tilted upwellings, slab fl attening and accumulation, ancient subduction, broadband seismic instruments in California. These recordings feature and processes unrelated to present-day subduction and plume ascent. simple fi rst-arriving P waves that have amplitudes at least three times larger than the coda. We measure the variation in T410 and T660–410 using INTRODUCTION multiple stacks with common P410s conversion points within circular Constraints on the variation of temperature in the Earth’s mantle are caps with radii of 1.5° across California and portions of adjacent states critical for our understanding of Earth’s convective behavior. The best (Fig. 2). The second set of T410 and T660–410 measurements is from the study constraints on temperature in the deeper parts of the mantle are obtained of Tauzin et al. (2008). T410 and T660–410 from this data set are measured from seismic studies. In particular, analyses of traveltimes of refl ected using single-station stacks for 167 global network stations. It is the most (Shearer, 1991) and converted (Vinnik, 1977) waves that have interacted with wave speed and density discontinuities due to the olivine to wadsley- ite (the 410 km) and the ringwoodite to perovskite + ferropericlase (the 660 km) phase transitions, yield localized constraints on temperature in Surface the transition zone (the region between the 410 km and the 660 km transi- tions). Since their Clapeyron slopes have opposite signs (Katsura and Ito, 1989), the 410 km transition is raised in a relatively cold upper mantle 410 km while the 660 km transition is lowered, thereby thickening the transition zone. Opposite defl ections of the 410 km and 660 km and, hence, thin- ner transition zones, characterize the relatively warm regions of the upper Transition zone 660 km mantle (Bina and Helffrich, 1994). We introduce a new thermometer for the upper mantle based on a P wave combination of mineral physics constraints and the analysis of traveltimes S wave of compressional waves that have converted to shear waves in the tran- P660s P410s P sition zone. Contrary to previous global studies (Chevrot et al., 1999; Lawrence and Shearer, 2006; Tauzin et al., 2008), we do not infer the temperature from estimates of transition zone thickness that rely on upper mantle velocity corrections. Instead, our analysis employs theoretical one- Figure 1. Ray paths of teleseismic P, P410s, and P660s through upper dimensional profi les of seismic velocity based on realistic upper mantle mantle. Phase P410s converts from P to S wave at 410 km discontinu- bulk compositions and geotherms. These profi les account for the full com- ity, while P660s converts from P to S at 660 km discontinuity. P wave plexity of the phase diagram and include all the mineralogical and thermal and S wave segments are indicated by solid and dashed lines, re- spectively. Note that P, P410s, and P660s have nearly parallel P wave effects on phase transition depths, density, and wave speeds. Hence, they propagation paths in the lower mantle and that P410s and P660s differ enable us to constrain the potential temperature of the mantle from travel- only in the transition zone. Differential traveltimes T410 and T660–410 are times directly and self-consistently. affected only by structure in upper 660 km of mantle. © 2009 The Geological Society of America. For permission to copy, contact Copyright Permissions, GSA, or [email protected]. GEOLOGY,Geology, April April 2009; 2009 v. 37; no. 4; p. 363–366; doi: 10.1130/G25428A.1; 4 fi gures. 363 predicted value) in California is due to reduced shear velocities in the rela- A tively warm uppermost mantle beneath western North America. The set of global T410 measurements exhibits a broad variation due to the relatively strong (thermal and compositional) heterogeneity in the 40 uppermost 200–300 km of mantle. T410 peaks near 43.5 s, a value that is nearly identical to the ak135 predicted value. T410 measurements are skewed to low values because the majority of global seismic stations are located on relatively cool continental lithosphere. The global set of T660–410 has a narrower distribution than T410. T660–410 is sensitive only to transition zone thickness and velocity variations. 35 ′ Q Q TRAVELTIME PREDICTIONS We determine theoretical values of T410 and T660–410 from the velocity profi les of Xu et al. (2008), computed using the self-consistent thermo- Latitude (°N) dynamic method of Stixrude and Lithgow-Bertelloni (2005a, 2005b) and the updated parameters in Xu et al. (2008). For a given mantle composition 30 and geotherm, velocity profi les are determined by computing equilibrium phase assemblages, the physical properties of the individual phases via stress and temperature derivatives of the Gibbs free energy, and the elastic T 660–410 properties of the assemblage from the Voigt-Reuss-Hill average. T410 and T660–410 are determined using refl ectivity seismograms (Fuchs and Muller, 22.25 22.65 23.05 23.50 23.95 24.35 24.75 1971) that fully incorporate the effects on P410s and P660s waveforms 25 due to the abrupt 410 km and 660 km and the gradual transformations in 125 120 115 110 the non-olivine components. Longitude (°W) Xu et al. (2008) considered the mantle to be a mixture of basalt (XB) and harzburgite (XH) with a basalt fraction of 18%. They considered two Q Q′ models of identical bulk composition that represent contrasting views of 0 the reequilibration of subducted heterogeneity: (1) The equilibrium assem- B blage (EA) represents perfect re-equilibration of basalt and harzburgite 10 into a homogeneous pyrolite. (2) The mechanical mixture (MM) main- tains perfect disequilibrium. The MM value is motivated by the suggestion 20 that the compositionally stratifi ed oceanic lithosphere is responsible for hetero geneity in the mantle (Allègre and Turcotte, 1986). 30 RESULTS 40 Figure 4 compares T410 and T660–410 for both the California and global data sets to theoretical values for the MM and EA and a range of adia- 50 batic geotherms with potential temperature that vary between 1400 K and 1800 K. Low values for T410 and high values for T660–410 are recorded at Time since P (s) 60 stations in subduction zones. High and low values for T410 and T660–410 are recorded at hotspot stations above putative mantle plumes, as expected if 70 the traveltime variations refl ect temperature heterogeneity in the mantle. However, the total global variation in T410 and T660–410 among the normal 80 stations is as large as the differences in T410 and T660–410 between hotspot 122 120 118 116 114 and subduction zone stations. Longitude (°W) Our computations indicate that the range in T410 and T660–410 values of ~8 s can be explained by MM for geotherms with potential temperatures Figure 2.