Inner core rotation from event-pair analysis

Size: px
Start display at page:

Download "Inner core rotation from event-pair analysis"

Transcription

1 Earth and Planetary Science Letters 261 (2007) Inner core rotation from event-pair analysis Xiaodong Song a,, Georges Poupinet b a Department of Geology, University of Illinois, Urbana, IL 61801, United States b LGIT-CNRS, Universite Joseph Fourier, Grenoble, France Received 21 April 2007; received in revised form 18 June 2007; accepted 29 June 2007 Available online 19 July 2007 Editor: R.D. van der Hilst Abstract The last decade has witnessed an animated debate on whether the inner core rotation is a fact or an artifact. Here we examine the temporal change of inner core waves using a technique that compares differential travel times at the same station but between two events. The method does not require precise knowledge of earthquake locations and earth models. The pairing of the events creates a large data set for the application of statistical tools. Using measurements from 87 events in the South Sandwich Islands recorded at College, Alaska station, we conclude the temporal change is robust. The estimates of the temporal change range from about 0.07 to 0.10 s/decade over the past 50 yr. If we used only pairs with small inter-event distances, which reduce the influence of mantle heterogeneity, the rates range from to s/decade, nearly identical to the rate inferred by Zhang et al. [Zhang, J., Song, X.D., Li, Y.C., Richards, P.G., Sun, X.L., Waldhauser, F., Inner core differential motion confirmed by earthquake waveform doublets, Science 309 (5739) (2005) ] from waveform doublets. The rate of the DF change seems to change with time, which may be explained by lateral variation of the inner core structure or the change in rotation rate on decadal time scale Elsevier B.V. All rights reserved. Keywords: inner core rotation; PKP differential times; bootstrap 1. Introduction The electromagnetic torque from the geodynamo is likely to drive the solid inner core to rotate relative to the mantle (e.g., Glatzmaier and Roberts, 1995). Since the first observation of the inner core rotation from seismology (Song and Richards, 1996), many studies have given support for an inner core motion, including most recent ones (Zhang et al., 2005; Wen, 2006; Cao et al., 2007). The last decade has also witnessed an animated debate on whether the inner core rotation is real, because Corresponding author. address: xsong@uiuc.edu (X. Song). of potential biases such as earthquake mislocation (e.g., Poupinet and Souriau, 2001), failure to detect the motion (e.g., Souriau, 1998), and discrepancy in the inferred rotation rate (e.g., Creager, 1997; Vidale et al., 2000; Laske and Masters, 2003). Recent effort attempted to use earthquake waveform doublets to study the inner core rotation, which are pairs of earthquakes that have identical waveforms recorded the same station, indicating close spatial locations of the two events (Li and Richards, 1999; Poupinet et al., 2000). Zhang et al. (2005) has presented perhaps the strongest evidence for inner core rotation so far with a set of 18 high-quality waveform doublets from South Sandwich Islands (SSI) earthquakes. These doublets show systematic shifts in travel times and changes in the X/$ - see front matter 2007 Elsevier B.V. All rights reserved. doi: /j.epsl

2 260 X. Song, G. Poupinet / Earth and Planetary Science Letters 261 (2007) waveforms through the inner core. The use of waveform doublets avoids artifacts of earthquake mislocation and contamination from small-scale mantle heterogeneities, and allows them to pin down precisely that the temporal changes must have occurred when the waves travel through the inner core. However, waveform doublets that show clear signals through the inner core are relatively rare. In this study, we examine the temporal change of inner core waves using a technique that compares differential travel times between two events recorded at the same station. Similar to the waveform doublet approach, the method does not require precise knowledge of earthquake locations and earth models. The method avoids the stringent requirement of high waveform similarity for waveform doublets, so the data quality is not comparable to that of waveform doublets and the measurements of the relative-time difference are not as precise. However, the pairing of the events creates a large number of event pairs that we can work with using standard statistics. 2. Data and method Our data are differential PKP travel-time measurements of SSI events at College, Alaska (COL) station (Fig. 1). Part of the data has been used extensively for studying inner core anisotropy (Creager, 1992; Song and Helmberger, 1998) and inner core rotation (Song and Richards, 1996; Song, 2000; Poupinet et al., 2000; Zhang et al., 2005). The differential travel times involve three branches of the PKP waves, which traverses the inner core (DF), turns at the bottom of the outer core (BC), and turns in mid-outer core (AB), respectively. We select events for which both differential AB BC time and BC DF time can be measured at the same trace (Fig. 1D). We obtain a total of 87 events, spanning 50 yr Fig. 1. SSI events used in this study. (A) Spatial distribution of all SSI events used (a total of 87). Both BC DF and AB BC times at COL have been measured for these events. The tight cluster (box) contains a total of 30 events. (B, C) Latitudes and longitudes plotted as a function of the origin time of the events. The spatial distribution of the events is uniform over time. (D) Example PKP seismograms recorded at COL station. The two events are Mar 15, 1986 (depth 90 km, mb 5.9) and Mar, 26, 1986 (depth 111 km, mb 5.7). The epicentral distances (Δ) to COL are and , respectively. All three branches of PKP are clear in these traces. Differential times between BC and DF and between AB and BC are measured for each trace. In this study, we form the differential travel time differences between a pair of events, i.e. d(bc DF) and d(ab BC). The difference in Δ affects both d(ab BC) and d(bc DF) (more on the former than on the latter), thus we can calibrate the distance difference using d(ab BC) to correct d(bc DF).

3 X. Song, G. Poupinet / Earth and Planetary Science Letters 261 (2007) from 1958 to 2007 (Fig. 1A C). The historical data are from various sources, which are described in detail in Song (2000). Differential PKP travel times are measured using waveform cross-correlations for each event. All the seismograms (vertical component) have been converted to the short-period WWSSN instrument and, if helpful, are filtered at frequency band of about 0.6 to 3 Hz to increase signal to noise ratio. The DF and BC phases are correlated directly; the BC and AB phases are correlated after correcting for the Hilbert transform of the AB phase. The cross-correlation coefficients are larger than about The median value for the DF and BC correlation is 0.73 and the median value for the BC and AB correlation is Our basic method is to compare the difference in the differential BC DF times, d(bc DF), and the difference in the differential AB BC times, d(ab BC), between a pair of events (Fig. 1D). Denoting the subscript 1 for the earlier event, and subscript 2 for the later event, we define dðbc DFÞ ¼ðBC DFÞ 2 ðbc DFÞ 1 ; and dðab BCÞ ¼ðAB BCÞ 2 ðab BCÞ 1 : The difference in the locations between the two events affects both the d(ab BC) and the d(bc DF) values. In fact, it changes the d(ab DF) value more than the d(bc DF) value. Thus we can correct the influence of location difference on d(bc DF) using the d(ab BC) value. The differential time difference ddt (either d(bc DF) or d(ab BC)) is affected by the location (epicentral distance to the station Δ and depth h) difference, heterogeneities of Earth structure, and measurement errors: ddt ¼ f 1 ðddþþf 2 ðdhþþf 3 ðstructureþ þ f 4 ðmeasurement errorsþ: Because these rays are nearly vertical, the depth effect f 2 is much smaller than the distance effect f 1.For the same amount (50 km) of difference in epicentral distance (dδ) and in depth (dh), d(bc DF) would be about 0.36 s and 0.06 s, respectively, and d(ab BC) would be about 0.82 s and 0.20 s, respectively. The SSI earthquakes are relatively shallow. Almost all events are less than 150 km in depth, with the average depth of 80 km and the standard deviation of 45 km. Thus we ignore the f 2 term in our following analysis. We regard the measurement error term f 4 as random noise. The structure term f 3 comes from mantle for both d(bc DF) and d(ab BC) and inner core heterogeneities for d(bc DF). Our earthquake distributions are quite uniform for different time periods (Fig. 1B,C). Our data sample small patches of the mantle and the inner core. The coverage of these patches is quite uniform for event pairs with different time lapses. Thus, we also regard f 3 as random noise. We will discuss further the depth ( f 1 )and mantle ( f 3 ) effects below. When there are enough pairs, the scaling factor from d(ab BC) to d(bc DF) can be calibrated, which is primarily affected by the distance term f 1.Ifthe DF travel time changes with time, the temporal change would affect d(bc DF) but not d(ab BC). Thus we can detect the temporal change by analyzing how d(bc DF), after being corrected by d(ab BC), changes with the time lapse between the two members of the pair. To summarize, our technique involves the following steps. (1) Measure BC DF and AB BC times for each event. (2) Form d(bc DF) and d(ab BC) for pairs of events. (3) Find the scaling factor between d(bc DF) and d(ab BC) using pairs of short time lapses to limit the influence of temporal changes. (4) Correct the d(bc DF) measurement using the d(ab BC) measurement and the scaling factor, i.e., the corrected d(bc DF)=d(BC DF) factor d(ab BC). As discussed above, this step removes the main influence from the location difference (epicentral distance difference) on d(bc DF). (5) Analyze how the corrected d(bc DF) value changes with the time lapse. In our analysis, we use heavily the standard statistical method of bootstrapping (e.g., Efron and Tibshirani, 1993) to estimate the parameters and its errors. When the differential travel times are measured, wrong cycles may be picked in the cross-correlation functions. We corrected or eliminated a few events with apparent cycle-skipping, which would form characteristic bands in the pairing plots. We conduct synthetic tests to demonstrate the methodology and to examine the resolution of the method (Fig. 2). Using the catalog locations of the same events, we calculate the PKP travel times for 1D reference model AK135 (Kennettetal.,1995). We then add a temporal change of s/yr to the DF travel times and calculate the differential BC DF and AB BC times for each event. In the first test, we don't include any noise. Using the data with time lapses less than 5 yr, we find the scaling factor to be / (1σ) by bootstrapping. Bootstrapping the corrected d(bc DF) values, we find a temporal change that is virtually the same as the input model with a very small error ( / s/yr) (Fig. 2A). In the above test (noise free), the full recovery of the temporal change is achieved without considering the

4 262 X. Song, G. Poupinet / Earth and Planetary Science Letters 261 (2007) Results Fig. 2. Synthetic tests of the recovery of the imposed temporal trend (0.008 s/yr), with (B) and without (A) noise. See text for description. depth difference of the events. As discussed previously, the depth influence (the f 2 term) is much smaller than the distance influence (the f 1 term). In addition, the random pairing of the events makes the f 2 term to be random as well, as evidenced in the figure. If there were no depth difference for all the pairs, all the synthetic data would fall into the straight line. The scatter (about 0.1 s) can be considered as the limit of the depth term. Thus, the term is expected to contribute to the scatter of the observations but does not affect the rate of the temporal change. In the second test, we add Gaussian noise to both BC DF and AB BC measurements with standard deviations of 0.1 s and 0.2 s, respectively. Typical measurement errors of the differential travel times are less than 0.1 s. The larger errors we impose take into account some influences from the mantle and inner core heterogeneities. The scaling factor in this case differs slightly: / The corrected d(bc DF) values show the scatter similar to or perhaps larger than that of our actual data (see below). Despite the large scatter, we can recover the temporal change quite well with a relatively small error bar ( / s/yr) (Fig. 2B). We first determine the scaling factor between d(bc DF) and d(ab BC) using all the pairs with time lapses less than 5 yr (total 733 pairs) (Fig. 3). The optimal value is / from bootstrapping. The scaling factor is significantly larger than that of 0.43 for model AK135 from our synthetic test. This may be caused by measurement errors but more likely by the departure of the average velocity of the Earth for this path from the global reference model. Possible sources include the inner core anisotropy and the velocity gradients in the lowermost mantle and the lowermost outer core. Using the scaling factor of 0.53, Fig. 4 shows all the corrected d(bc DF) values as a function of time lapse (a total of 3741 pairs). The corrected values show a clear temporal trend with a slope of about s/yr. We explore the robustness of the temporal change of d(bc DF) using bootstrapping (Table 1). We test various groups and subgroups of data sets. For each data set, we run 1000 iterations of bootstrapping by selecting events from the original event list at random, obtain the regression of the corrected d(bc DF) values with respect to the time lapse for each iteration, and finally derive the average and the standard deviation of the regressions. For each data set, we use the scaling factor of 0.53, our optimal value. We also test the scaling factor of 0.43, which is appropriate for AK135 for our data set, to examine the sensitivity of the regression results to the scaling factor. In all cases we have tested, the results do not change significantly with the scaling factor between the two values (well within 1σ) (Table 1). We consider the whole data set as well as various subgroups. (1) The distance (d) between the two events calculated from the catalog locations varies from within 0.2 to within 2. The smaller the distance, the more Fig. 3. All pairs with time lapse less than 5 yr used to find the scaling factor between d(bc DF) and d(ab BC). The optimal value (0.528+/ 0.013) is found from bootstrap regression with 1000 iterations.

5 X. Song, G. Poupinet / Earth and Planetary Science Letters 261 (2007) Fig. 4. Corrected d(bc DF) as a function of time lapse (dyr) between the two events. The correction is made using d(ab BC) and the scaling factor of The slope of line, obtained from the bootstrap regression of all the data (Table 1, group 6), is the measure of the DF travel time change. likely that the two events are in fact closer, and thus the less bias the Earth's structure heterogeneity has on d(bc DF) and d(ab BC). However, the number of observations decreases with d. (2) We test groups with various d (AB BC) values, essentially using the d(ab BC) as a proxy for the difference in epicentral distances to the station. (3) We test data after 1964, and after 1981, respectively. Between 1958 and 1965, we only have a few events (less than 10%). After 1981, the records became digital, allowing us to obtain differential times more accurately. (4) The epicentral distance (Δ) of our events ranges from to 153.3, sampling the top 300 km of the inner core. The inner core rotation may cause slightly different time shifts for waves traveling at different depths and spending different times in the inner core. We consider the data set with Δ between 150 and 152, which contains more than 90% of the total event pairs. (5) More than one third of the events are clustered in a very small area of the study region (Fig. 1A), providing a natural data set without the need to impose other restrictions. Considering that the number of observations varies by an order of magnitude (Table 1), the temporal change in the d(bc DF) value, ranging from to s/decade (excluding the subgroups using shorter period ), is remarkably robust. The standard deviation of the mean depends on the number of observations and the temporal coverage of the data. For pairs with smaller d or d(ab BC) values or in the cluster, which have smaller spatial separation and are less affected by mantle heterogeneity, the inferred rates (0.084 to s/decade) are nearly identical to the rate inferred by Zhang et al. (2005) from waveform doublets (0.092+/ s/decade). We have assumed that the f 3 term (from mantle heterogeneity) to be random as discussed above. Precise knowledge of the term is not possible as it requires the determination of very fine-scale structure (km to hundredkm scale) of the mantle, particularly the lowermost mantle. However, we can now estimate the influence of this term from our tests using different parameters (Table 1). The range of the values (0.71 to 0.98 s/decade) we obtained for pairs in Groups 1, 2, 5, and 6 (Table 1)canbe regarded as the maximum variation caused by mantle heterogeneity, as other factors (measurement errors and differences in the number of pairs) also contribute. The smaller the spatial separation between the two events of the pair is, the less influence the mantle heterogeneity would have. It is encouraging that for pairs with small d or d(ab BC) values or in the cluster, the inferred rates are very similar to the rate inferred from waveform doublets (Zhang et al., 2005). All rate estimates using larger data sets (looser requirement on inter-event distance) are smaller, by 0.05 to s/decade or by up to 30%. Estimates of the temporal change using the shorter time period are consistently smaller than the estimates using the whole data set (Fig. 5). The difference is significant at 68% confidence (1σ) but not significant at 95% confidence (2σ). We have also examined the data of even shorter time period ( ). The average rate using all the data in this time period appears even smaller: / s/decade. However, the uncertainty is very large and the issue requires further investigation. Using various groups and Table 1 Estimates of temporal change of d(bc DF) Group Data selection f=0.53 f=0.43 N 1 Restriction on inter-event distance (d) db (0.026) (0.025) 371 db (0.017) (0.016) 979 db (0.013) (0.013) 1927 db (0.012) (0.013) Restriction on d(ab BC) value (dab) dabb0.2 s (0.019) (0.019) 403 dabb0.5 s (0.014) (0.014) 916 dabb1.0 s (0.012) (0.012) 1630 dabb3.0 s (0.012) (0.013) Restriction on origin times of both events (0.016) (0.022) (0.014) (0.026) Restriction on epicentral distances (Δ) of both events 150 bδb (0.018) (0.017) Pairs with both events in the cluster shown in Fig. 1A Clustered pairs (0.024) (0.022) All data (0.012) (0.015) 3741 The temporal change and its standard deviation (in parentheses) are in s/decade. The parameter f is the scaling factor from d(ab BC) to d(bc DF), and N is the average number of observations over 1000 iterations of bootstrapping.

6 264 X. Song, G. Poupinet / Earth and Planetary Science Letters 261 (2007) Fig. 5. Estimates of the temporal change of d(bc DF) using various parameters. The solid dots are for data of all time periods ( ) and the open circles are for events in with modern digital records. The error bar is +/ one standard deviation. Models 1 through 4 are from Group 1 (Table 1) with d from 0.2 to 2, respectively; Models 5 through 8 are from Group 2 with dab from 0.2 s to 3.0 s, respectively; Models 9, 10, 11 correspond to Group 4, 5, and 6, respectively. The scaling factor f=0.53 is used for these calculations. subgroups of the data as in Table 1, the estimated rates range from to The cause of the change of the temporal rate is not clear. The possibilities include the following. (1) Errors in the estimation of the rate. Both the time lapse and the number of observations decrease. As the time lapse decreases, the amount of the temporal shift also decreases. The total number of pairs (1431) for is only 38% of that of the whole data set (3741). Thus the estimate of the rate is more subject to errors from measurements and mantle heterogeneity. If we consider only the time period in the synthetic test (Fig. 2), the rate is also considerably different with a large error ( / s/yr). However, the rate estimate errors for are comparable to those of the whole data set (Fig. 5). Thus this may not be the main cause. (2) Changes of seismic instrumentation over time. Although the instrument responses to different frequency contents of DF and BC and AB have not been well tested, this may not be the cause either. All the seismograms have been converted to the same narrowband instrument (short-period WWSSN) before the differential travel times are measured using waveform correlation. (3) Sampling of laterally-varying small-scale inner core structures, which have been suggested (Vidale and Earle, 2000; Cormier and Li, 2002; Koper et al., 2004). If the inner core rotates at 0.3 to 0.5 /yr (Zhang et al., 2005), significant lateral variation over km distance in the inner core would explain the change in the d(bc DF) temporal shift. (4) A change in the inner core rotation rate. If the lateral structure of the inner core at this region varies smoothly (linearly), a Fig. 6. Probability density distribution of raw data and the means. (A) Distribution of the corrected d(bc DF) value, i.e., d(bc DF) 0.53 d(ab BC), for pairs with different time lapses. The time lapses are within 5 yr for each group, except the top group with largest time lapses, which include data with time lapses between 40 and 48 yr. The number of pairs in each group is labeled at the left. (B) Distribution of the means of the corrected d(bc DF) value for groups with different time lapses. A sliding time-lapse window of 5 yr and a sliding step of 2.5 yr are used to form the groups, except for the last group (top), which has the time window of 6.3 yr. The means are calculated using 5000 bootstrapping iterations. The dashed line is bootstrap result using all the data (Table 1) with a slope of / s/yr.

7 X. Song, G. Poupinet / Earth and Planetary Science Letters 261 (2007) change of rotation rate could cause the change in the DF travel time. Cases 3 and 4 cannot be distinguished with the data at hand. The robust temporal change is demonstrated graphically in Fig. 6. It is clearly visible in the distribution of the raw data when they are divided according to the time lapse (Fig. 6A). It is even clearer in the distribution of the means from boot-strapping (Fig. 6B). The means generally follow a linear trend with time lapse along the average rate with some fluctuations, particularly at larger time lapses when the number of observations decreases. 4. Conclusion We propose a technique that compares differential PKP travel time measurements between event pairs to detect the temporal change of the inner core. Using 87 SSI events recorded at COL station, we estimate the average temporal change in the DF times to be about 0.07 to 0.10 s/decade over the past 50 yr. For pairs that are close in location (as indicated by smaller d, d(ab BC) values, or in the cluster), the inferred rates (0.084 to s/decade) are nearly identical to the rate inferred by Zhang et al. (2005) from waveform doublets ( / s/decade). Bootstrapping analyses of various groups and subgroups of the data set suggest the temporal trend is very robust. The rate of the DF change seems to change with time, which may be explained by lateral variation of the inner core structure or the change in rotation rate in decadal time scale. If the rotation rate change is true, it may explain the discrepancy inferred from the data of different time periods (e.g., Song, 2000; Vidale et al., 2000; Laske and Masters, 2003). The most appealing aspect of the new method is the minimum assumptions required. We rely solely on a small initial data set (174 differential time measurements from 87 events) and do not need precise knowledge of earthquake locations and earth models. The calibration for location difference is purely based on the original data set. Another aspect of the method is the large data set created by the pairing of events, which allow us to use statistical tools to find coherent signals and reduce the influence of random noises (from measurement errors or from structure variations). Our synthetic tests suggest that the temporal rate can be well resolved if the noise is well behaved (random) and the signal to noise ratio is high enough. The technique provides an alternative when perfect doublets are hard to find. It can be similarly applied to other regions and other seismic phases and may be applicable to other processes (such as temporal changes in crustal properties and in magma chambers). Acknowledgements Historical data at COL station were acquired with the help of many people (Song, 2000). Recent digital data were obtained from the IRIS DMC. We thank Ping Ma for comments on statistics and John Vidale for helpful review. The research was supported by NSF EAR and UIUC/CNRS Cooperative Research Program. References Cao, A.M., Masson, Y., Romanowicz, B., Short wavelength topography on the inner-core boundary. Proc. Natl. Acad. Sci. U. S. A. 104 (1), Cormier, V.F., Li, X., Frequency-dependent seismic attenuation in the inner core 2. A scattering and fabric interpretation. J. Geophys. Res. Solid Earth 107 (B12), Creager, K.C., Anisotropy of the inner core from differential travel-times of the phases PKP and PKIKP. Nature 356 (6367), Creager, K.C., Inner core rotation rate from small-scale heterogeneity and time-varying travel times. Science 278 (5341), Efron, B., Tibshirani, R.J., An Introduction to the Bootstrap. Chapman & Hall, New York. Glatzmaier, G.A., Roberts, P.H., A 3-dimensional convective dynamo solution with rotating and finitely conducting inner-core and mantle. Phys. Earth Planet. Inter. 91 (1 3), Kennett, B.L.N., Engdahl, E.R., Buland, R., Constraints on seismic velocities in the earth from travel-times. Geophys. J. Int. 122 (1), Koper, K.D., Franks, J.M., Dombrovskaya, M., Evidence for small-scale heterogeneity in earth's inner core from a global study of pkikp coda waves. Earth Planet. Sci. Lett. 228 (3 4), Laske, G., Masters, G., The Earth's free oscillation and differential rotation of the inner core. Vol. 31 of Earth's Core: Geodynamics, Structure, Rotation, Geodyn. Ser.. American Geophysics Union. Li, A.Y., Richards, P.G., On Detection of Inner Core Rotation using Doublets, Fall AGU Meet., San Francisco. Poupinet, G., Souriau, A., Reply to xiadong song's comment on the existence of an inner core super-rotation questioned by teleseismic doublets. Phys. Earth Planet. Inter. 124 (3 4), Poupinet, G., Souriau, A., Coutant, O., The existence of an inner core super-rotation questioned by teleseismic doublets. Phys. Earth Planet. Inter. 118 (1 2), Song, X.D., Joint inversion for inner core rotation, inner core anisotropy, and mantle heterogeneity. J. Geophys. Res. Solid Earth 105 (B4), Song, X.D., Helmberger, D.V., Seismic evidence for an inner core transition zone. Science 282 (5390), Song, X.D., Richards, P.G., Seismological evidence for differential rotation of the earth's inner core. Nature 382 (6588), Souriau, A., New seismological constraints on differential rotation of the inner core from novaya zemlya events recorded at drv, Antarctica. Geophys. J. Int. 134 (2), F1 F5.

8 266 X. Song, G. Poupinet / Earth and Planetary Science Letters 261 (2007) Vidale, J.E., Earle, P.S., Fine-scale heterogeneity in the Earth's inner core. Nature 404 (6775), Vidale, J.E., Dodge, D.A., Earle, P.S., Slow differential rotation of the earth's inner core indicated by temporal changes in scattering. Nature 405 (6785), Wen, L.X., Localized temporal change of the earth's inner core boundary. Science 314 (5801), Zhang, J., Song, X.D., Li, Y.C., Richards, P.G., Sun, X.L., Waldhauser, F., Inner core differential motion confirmed by earthquake waveform doublets. Science 309 (5739),

Time dependence of PKP(BC) PKP(DF) times: could this be an artifact of systematic earthquake mislocations?

Time dependence of PKP(BC) PKP(DF) times: could this be an artifact of systematic earthquake mislocations? Physics of the Earth and Planetary Interiors 122 (2000) 221 228 Time dependence of PKP(BC) PKP(DF) times: could this be an artifact of systematic earthquake mislocations? Xiaodong Song Department of Geology,

More information

Examination of systematic mislocation of South Sandwich Islands earthquakes using station pairs: Implications for inner core rotation

Examination of systematic mislocation of South Sandwich Islands earthquakes using station pairs: Implications for inner core rotation JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:10.1029/2005jb004175, 2006 Examination of systematic mislocation of South Sandwich Islands earthquakes using station pairs: Implications for inner core rotation

More information

Localized temporal variation of Earth s inner-core boundary from high-quality waveform doublets

Localized temporal variation of Earth s inner-core boundary from high-quality waveform doublets Earthq Sci (2015) 28(3):175 185 DOI 10.1007/s11589-015-0125-0 RESEARCH PAPER Localized temporal variation of Earth s inner-core boundary from high-quality waveform doublets Danhua Xin. Xiaodong Song. Tao

More information

Supporting Online Material for

Supporting Online Material for www.sciencemag.org/cgi/content/full/1131692/dc1 Supporting Online Material for Localized Temporal Change of the Earth s Inner Core Boundary This PDF file includes: Materials and Methods Figs. S1 to S3

More information

Detection of motion and heterogeneity in Earth s liquid outer core

Detection of motion and heterogeneity in Earth s liquid outer core GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L16311, doi:10.1029/2008gl034895, 2008 Detection of motion and heterogeneity in Earth s liquid outer core Wei Dai 1 and Xiaodong Song 1 Received 5 June 2008; revised

More information

Probing Mid-Mantle Heterogeneity Using PKP Coda Waves

Probing Mid-Mantle Heterogeneity Using PKP Coda Waves Probing Mid-Mantle Heterogeneity Using PKP Coda Waves Michael A.H. Hedlin and Peter M. Shearer Cecil H. and Ida M. Green Institute of Geophysics and Planetary Physics Scripps Institution of Oceanography,

More information

High-precision location of North Korea s 2009 nuclear test

High-precision location of North Korea s 2009 nuclear test Copyright, Seismological Research Letters, Seismological Society of America 1 High-precision location of North Korea s 2009 nuclear test Lianxing Wen & Hui Long Department of Geosciences State University

More information

Inner core rotation and its variability from nonparametric modeling

Inner core rotation and its variability from nonparametric modeling Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2009jb006294, 2010 Inner core rotation and its variability from nonparametric modeling Daniela Lindner, 1 Xiaodong Song,

More information

G 3. AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society

G 3. AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society Geosystems G 3 AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society Article Volume 4, Number 9 9 September 2003 1072, doi:10.1029/2002gc000379 ISSN: 1525-2027 Using

More information

Inner core attenuation anisotropy

Inner core attenuation anisotropy Earth and Planetary Science Letters 245 (2006) 581 594 www.elsevier.com/locate/epsl Inner core attenuation anisotropy Wen-che Yu, Lianxing Wen Department of Geosciences, State University of New York at

More information

Geophysical Journal International

Geophysical Journal International Geophysical Journal International Geophys. J. Int. (2011) 187, 355 370 doi: 10.1111/j.1365-246X.2011.05146.x Global seismic body-wave observations of temporal variations in the Earth s inner core, and

More information

ARTICLE IN PRESS UNCORRECTED PROOF

ARTICLE IN PRESS UNCORRECTED PROOF Earth and Planetary Science Letters 6022 (2001) 1^15 www.elsevier.com/locate/epsl Estimate of inner core rotation rate from United Kingdom regional seismic network data and consequences for inner core

More information

Constraints on density and shear velocity contrasts at the inner core boundary

Constraints on density and shear velocity contrasts at the inner core boundary Geophys. J. Int. (00) 57, 6 5 doi: 0./j.65-6X.00.00.x FAST TRACK PAPER Constraints on density and shear velocity contrasts at the inner core boundary Aimin Cao and Barbara Romanowicz Seismological Laboratory,

More information

Wide-scale detection of earthquake waveform doublets and further evidence for inner core super-rotation

Wide-scale detection of earthquake waveform doublets and further evidence for inner core super-rotation Geophys. J. Int. (2008) 174, 993 1006 doi: 10.1111/j.1365-246X.2008.03856.x Wide-scale detection of earthquake waveform doublets and further evidence for inner core super-rotation Jian Zhang, 1, Paul G.

More information

Waveform search for the innermost inner core

Waveform search for the innermost inner core Waveform search for the innermost inner core Vernon F. Cormier 1 and Anastasia Stroujkova 1,2 University of Connecticut Storrs, CT 06269-3046 Abstract Waveforms of the PKIKP seismic phase in the distance

More information

A Local, Crossing-path Study of Attenuation and Anisotropy of the Inner Core

A Local, Crossing-path Study of Attenuation and Anisotropy of the Inner Core 1 A Local, Crossing-path Study of Attenuation and Anisotropy of the Inner Core George Helffrich 1, Satoshi Kaneshima Earth and Planetary Sciences, Tokyo Institute of Technology J.-M. Kendall Earth Sciences,

More information

On the observation of high frequency PKiKP and its coda in Australia

On the observation of high frequency PKiKP and its coda in Australia Physics of the Earth and Planetary Interiors 146 (2004) 497 511 On the observation of high frequency PKiKP and its coda in Australia G. Poupinet a,b,, B.L.N. Kennett a a RSES, Australian National University,

More information

FOCAL MECHANISM DETERMINATION USING WAVEFORM DATA FROM A BROADBAND STATION IN THE PHILIPPINES

FOCAL MECHANISM DETERMINATION USING WAVEFORM DATA FROM A BROADBAND STATION IN THE PHILIPPINES FOCAL MECHANISM DETERMINATION USING WAVEFORM DATA FROM A BROADBAND STATION IN THE PHILIPPINES Vilma Castillejos Hernandez Supervisor: Tatsuhiko Hara MEE10508 ABSTRACT We performed time domain moment tensor

More information

Existence of finite rigidity layer at the base of the Earth s liquid outer core inferred from anomalous splitting of normal modes

Existence of finite rigidity layer at the base of the Earth s liquid outer core inferred from anomalous splitting of normal modes LETTER Earth Planets Space, 54, 67 7, 22 Existence of finite rigidity layer at the base of the Earth s liquid outer core inferred from anomalous splitting of normal modes Seiji Tsuboi,2 and Masanori Saito

More information

PKP travel times at near antipodal distances: implications for inner core anisotropy and lowermost mantle structure

PKP travel times at near antipodal distances: implications for inner core anisotropy and lowermost mantle structure Earth and Planetary Science Letters 199 (2002) 429^445 www.elsevier.com/locate/epsl PKP travel times at near antipodal distances: implications for inner core anisotropy and lowermost mantle structure Xinlei

More information

On the Nature of the P-Wave Velocity Gradient in the Inner Core beneath Central America

On the Nature of the P-Wave Velocity Gradient in the Inner Core beneath Central America Journal of Earth Science, Vol. 24, No. 5, p. 699 705, October 2013 ISSN 1674-487X Printed in China DOI: 10.1007/s12583-013-0365-7 On the Nature of the P-Wave Velocity Gradient in the Inner Core beneath

More information

Nonparametric Methods for Detecting Structure and Dynamics of Earth s Deep Interior. The Wahba Conference June 6, 2014

Nonparametric Methods for Detecting Structure and Dynamics of Earth s Deep Interior. The Wahba Conference June 6, 2014 Nonparametric Methods for Detecting Structure and Dynamics of Earth s Deep Interior The Wahba Conference June 6, 2014 Connection with Grace Ph.D with Chong Gu (First generation) Postdoc with Jun S Liu

More information

One-Dimensional Modeling of Multiple Scattering in the Upper Inner Core: Depth Extent of a Scattering Region in the Eastern Hemisphere

One-Dimensional Modeling of Multiple Scattering in the Upper Inner Core: Depth Extent of a Scattering Region in the Eastern Hemisphere Journal of Earth Science, Vol. 24, No. 5, p. 706 715, October 2013 ISSN 1674-487X Printed in China DOI: 10.1007/s12583-013-0366-6 One-Dimensional Modeling of Multiple Scattering in the Upper Inner Core:

More information

29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies

29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies TRANSITION ZONE WAVE PROPAGATION: CHARACTERIZING TRAVEL-TIME AND AMPLITUDE INFORMATION Peter M. Shearer and Jesse F. Lawrence University of California San Diego, Institute of Geophysics and Planetary Physics

More information

Seismic interferometry with antipodal station pairs

Seismic interferometry with antipodal station pairs GEOPHYSICAL RESEARCH LETTERS, VOL. 4, 1 5, doi:1.12/grl.597, 213 Seismic interferometry with antipodal station pairs Fan-Chi Lin 1 and Victor C. Tsai 1 Received 25 June 213; revised 19 August 213; accepted

More information

Distinct layering in the hemispherical seismic velocity structure of Earth s upper inner core

Distinct layering in the hemispherical seismic velocity structure of Earth s upper inner core JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:10.1029/2011jb008650, 2011 Distinct layering in the hemispherical seismic velocity structure of Earth s upper inner core Lauren Waszek 1 and Arwen Deuss

More information

Geophysical Journal International

Geophysical Journal International Geophysical Journal International Geophys. J. Int. (2013) 195, 2005 2015 Advance Access publication 2013 October 1 doi: 10.1093/gji/ggt368 A low attenuation layer in the Earth s uppermost inner core Lauren

More information

Imaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges

Imaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges 2017 SCEC Proposal Report #17133 Imaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges Principal Investigator Zhongwen Zhan Seismological Laboratory, California

More information

Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival

Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 24, 2248, doi:10.1029/2003gl018413, 2003 Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival

More information

The inner inner core of the Earth: Texturing of iron crystals from three-dimensional seismic anisotropy

The inner inner core of the Earth: Texturing of iron crystals from three-dimensional seismic anisotropy Available online at www.sciencedirect.com Earth and Planetary Science Letters 269 (2008) 56 65 www.elsevier.com/locate/epsl The inner inner core of the Earth: Texturing of iron crystals from three-dimensional

More information

Crust and upper mantle P- and S-wave delay times at Eurasian seismic stations

Crust and upper mantle P- and S-wave delay times at Eurasian seismic stations Physics of the Earth and Planetary Interiors 123 (2001) 205 219 Crust and upper mantle P- and S-wave delay times at Eurasian seismic stations E.R. Engdahl, M.H. Ritzwoller Center for Imaging the Earth

More information

Complex seismic anisotropy in the top of the Earth s inner core beneath Africa

Complex seismic anisotropy in the top of the Earth s inner core beneath Africa Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jb004868, 2007 Complex seismic anisotropy in the top of the Earth s inner core beneath Africa Wen-che Yu 1 and Lianxing

More information

Sharp and seismically transparent inner core boundary region revealed by an entire network observation of near-vertical PKiKP

Sharp and seismically transparent inner core boundary region revealed by an entire network observation of near-vertical PKiKP Earth Planets Space, 58, 855 863, 2006 Sharp and seismically transparent inner core boundary region revealed by an entire network observation of near-vertical PKiKP Hitoshi Kawakatsu Earthquake Research

More information

APPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY

APPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY APPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY Timur TEZEL Supervisor: Takuo SHIBUTANI MEE07169 ABSTRACT In this study I tried to determine the shear wave velocity structure in the crust

More information

Evidence for P P asymmetrical scattering at near podal distances

Evidence for P P asymmetrical scattering at near podal distances GEOPHYSICAL RESEARCH LETTERS, VOL. 39,, doi:10.1029/2012gl052179, 2012 Evidence for P P asymmetrical scattering at near podal distances Wenbo Wu, 1,2 Sidao Ni, 2 and Xiangfang Zeng 1 Received 3 May 2012;

More information

29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies

29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies IMPROVING MAGNITUDE DETECTION THRESHOLDS USING MULTI-STATION, MULTI-EVENT, AND MULTI-PHASE METHODS David Schaff and Felix Waldhauser Lamont-Doherty Earth Observatory, Columbia University Sponsored by Air

More information

Some aspects of seismic tomography

Some aspects of seismic tomography Some aspects of seismic tomography Peter Shearer IGPP/SIO/U.C. San Diego September 7, 2009 Earthquake Research Institute Part 1: Global Tomography P velocity perturbations 200 km 1000 km 2700 km MIT 2006

More information

Physics of the Earth and Planetary Interiors

Physics of the Earth and Planetary Interiors hysics of the Earth and lanetary Interiors 167 (2008) 155 160 Contents lists available at ScienceDirect hysics of the Earth and lanetary Interiors journal homepage: www.elsevier.com/locate/pepi Depth-dependent

More information

Constraints on earthquake epicentres independent of seismic velocity models

Constraints on earthquake epicentres independent of seismic velocity models Geophys. J. Int. (2004) 156, 648 654 doi: 10.1111/j.1365-246X.2004.02123.x RESEARCH NOTE Constraints on earthquake epicentres independent of seismic velocity models T. Nicholson, 1 Ó. Gudmundsson 2 and

More information

29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies MODELING P WAVE MULTIPATHING IN SOUTHEAST ASIA

29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies MODELING P WAVE MULTIPATHING IN SOUTHEAST ASIA MODELING P WAVE MULTIPATHING IN SOUTHEAST ASIA Ali Fatehi and Keith D. Koper Saint Louis University Sponsored by the Air Force Research Laboratory ABSTRACT Contract No. FA8718-06-C-0003 We have used data

More information

Improving Global Seismic Event Locations Using Source-Receiver Reciprocity

Improving Global Seismic Event Locations Using Source-Receiver Reciprocity Bulletin of the Seismological Society of America, 91, 3, pp. 594 603, June 2001 Improving Global Seismic Event Locations Using Source-Receiver Reciprocity by Peter M. Shearer Abstract The leading source

More information

MYRES Seismic Constraints on Boundary Layers. Christine Thomas

MYRES Seismic Constraints on Boundary Layers. Christine Thomas MYRES 2004 Seismic Constraints on Boundary Layers Christine Thomas Outline Seismic constraints: what can we resolve? how can we detect seismic structures? how well can we resolve these structures? (resolution

More information

Earth s inner core: revealing the structures behind the PKP core phase triplication

Earth s inner core: revealing the structures behind the PKP core phase triplication GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:10.1029/, Earth s inner core: revealing the structures behind the PKP core phase triplication N.A. Blom 1 supervisors: A.F. Deuss, 2 H. Paulssen, 1 The

More information

Depth dependence of anisotropy of Earth's inner core

Depth dependence of anisotropy of Earth's inner core JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 100, NO. 87, PAGES 9805-9816, JUNE 10, 1995 Depth dependence of anisotropy of Earth's inner core Xiaodong Song 1 and Don V. Heimberger Seismological La.bora.tory,

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION Seismic evidence for a global low velocity layer within the Earth s upper mantle SUPPLEMENTARY MATERIAL Benoît Tauzin 1, Eric Debayle 2 & Gérard Wittlinger 3 1 Department of Earth

More information

24th Seismic Research Review Nuclear Explosion Monitoring: Innovation and Integration

24th Seismic Research Review Nuclear Explosion Monitoring: Innovation and Integration A REFERENCE DATA SET FOR VALIDATING 3-D MODELS E. R. Engdahl, Eric A. Bergman, Michael H. Ritzwoller, Nikolai M. Shapiro and Anatoli L. Levshin University of Colorado, Boulder Sponsored by Defense Threat

More information

Lateral variation of the D 00 discontinuity beneath the Cocos Plate

Lateral variation of the D 00 discontinuity beneath the Cocos Plate GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L15612, doi:10.1029/2004gl020300, 2004 Lateral variation of the D 00 discontinuity beneath the Cocos Plate T. Lay Earth Sciences Department, University of California,

More information

Scattering objects in the lower mantle beneath northeastern China observed with a short-period seismic array

Scattering objects in the lower mantle beneath northeastern China observed with a short-period seismic array Physics of the Earth and Planetary Interiors 138 (2003) 55 69 Scattering objects in the lower mantle beneath northeastern China observed with a short-period seismic array Tadashi Kito a,, Takuo Shibutani

More information

Short scale heterogeneity in the lowermost mantle: insights from PcP-P and ScS-S data

Short scale heterogeneity in the lowermost mantle: insights from PcP-P and ScS-S data Earth and Planetary Science Letters 201 (2002) 57^68 www.elsevier.com/locate/epsl Short scale heterogeneity in the lowermost mantle: insights from PcP-P and ScS-S data Hrvoje Tkalc ic, Barbara Romanowicz

More information

29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies

29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies MODELING TRAVEL-TIME CORRELATIONS BASED ON SENSITIVITY KERNELS AND CORRELATED VELOCITY ANOMALIES William L. Rodi 1 and Stephen C. Myers 2 Massachusetts Institute of Technology 1 and Lawrence Livermore

More information

Seismogram Interpretation. Seismogram Interpretation

Seismogram Interpretation. Seismogram Interpretation Travel times in the Earth Ray paths, phases and their name Wavefields in the Earth: SH waves, P-SV waves Seismic Tomography Receiver Functions Seismogram Example Long-period transverse displacement for

More information

Continent-sized anomalous zones with low seismic velocity at the base of Earth s mantle

Continent-sized anomalous zones with low seismic velocity at the base of Earth s mantle SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO2733 Continent-sized anomalous zones with low seismic velocity at the base of Earth s mantle Edward J. Garnero 1, Allen K. McNamara 1, and Sang-Heon D. Shim 1

More information

DETERMINATION OF EARTHQUAKE PARAMETERS USING SINGLE STATION BROADBAND DATA IN SRI LANKA

DETERMINATION OF EARTHQUAKE PARAMETERS USING SINGLE STATION BROADBAND DATA IN SRI LANKA DETERMINATION OF EARTHQUAKE PARAMETERS USING SINGLE STATION BROADBAND DATA IN SRI LANKA S.W.M. SENEVIRATNE* MEE71 Supervisors: Yasuhiro YOSHIDA** Tatsuhiko HARA*** ABSTRACT We determined epicenters and

More information

Figures S1 S4 show the measurements taken from the synthetic vespagrams where a)

Figures S1 S4 show the measurements taken from the synthetic vespagrams where a) Figures S1 S4 show the measurements taken from the synthetic vespagrams where a) is the differential travel time versus the Dʺ discontinuity height, b) is the differential travel time versus δv S, c) is

More information

Texture of the Uppermost Inner Core from Forward and Back. Scattered Seismic Waves

Texture of the Uppermost Inner Core from Forward and Back. Scattered Seismic Waves 1 1 2 3 4 5 6 Texture of the Uppermost Inner Core from Forward and Back Scattered Seismic Waves Vernon F. Cormier Physics Department University of Connecticut Storrs, CT 06269-3046 7 8 9 10 11 12 13 14

More information

From Core to Crust: Towards an Integrated Vision of Earth's Interior July Earth's Core: Seismological Perspective

From Core to Crust: Towards an Integrated Vision of Earth's Interior July Earth's Core: Seismological Perspective 2048-6 From Core to Crust: Towards an Integrated Vision of Earth's Interior 20-24 July 2009 Earth's Core: Seismological Perspective H. Tkalcic The Australian National University, Canberra, Australia Strada

More information

Supporting Online Material for

Supporting Online Material for www.sciencemag.org/cgi/content/full/326/5949/112/dc1 Supporting Online Material for Global Surface Wave Tomography Using Seismic Hum Kiwamu Nishida,* Jean-Paul Montagner, Hitoshi Kawakatsu *To whom correspondence

More information

Seismic Scattering in the Deep Earth

Seismic Scattering in the Deep Earth Seismic Scattering in the Deep Earth Peter Shearer IGPP/SIO/U.C. San Diego September 2, 2009 Earthquake Research Institute Mantle mixing calculations Heterogeneity is likely at all scales Davies (2002)

More information

Data Repository: Seismic and Geodetic Evidence For Extensive, Long-Lived Fault Damage Zones

Data Repository: Seismic and Geodetic Evidence For Extensive, Long-Lived Fault Damage Zones DR2009082 Data Repository: Seismic and Geodetic Evidence For Extensive, Long-Lived Fault Damage Zones Fault Zone Trapped Wave Data and Methods Fault zone trapped waves observed for 2 shots and 5 local

More information

Fenglin Niu and Hitoshi Kawakatsu. Earthquake Research Institute, University of Tokyo, Yayoi, Bunkyo-ku, Tokyo , Japan

Fenglin Niu and Hitoshi Kawakatsu. Earthquake Research Institute, University of Tokyo, Yayoi, Bunkyo-ku, Tokyo , Japan Earth Planets Space, 50, 965 975, 1998 Determination of the absolute depths of the mantle transition zone discontinuities beneath China: Effect of stagnant slabs on transition zone discontinuities Fenglin

More information

TOMOGRAPHY S VELOCITY STRUCTURE BETWEEN WASHINGTON S EARTHQUAKE C022801L AND OBSERVATIONAL STATION TUC THROUGH SEISMOGRAM ANALYSIS

TOMOGRAPHY S VELOCITY STRUCTURE BETWEEN WASHINGTON S EARTHQUAKE C022801L AND OBSERVATIONAL STATION TUC THROUGH SEISMOGRAM ANALYSIS 70 TOMOGRAPHY S VELOCITY STRUCTURE BETWEEN WASHINGTON S EARTHQUAKE C022801L AND OBSERVATIONAL STATION TUC THROUGH SEISMOGRAM ANALYSIS Bagus Jaya Santosa Jurusan Fisika, FMIPA, Institut Teknologi Sepuluh

More information

boundaries with additional record sections, as emphasized in Fig. S2. The observations at the

boundaries with additional record sections, as emphasized in Fig. S2. The observations at the Data used to Sample African Anomaly. The great circle paths of the source-receiver combinations used in this study are shown in Fig. S1. The event information is given in Table S1. Abrupt Changes across

More information

Theory. Summary. Introduction

Theory. Summary. Introduction Waveform similarity for quality control of event locations, time picking and moment tensor solutions Fernando Castellanos, University of Alberta. Edmonton, AB. Mirko van der Baan, University of Alberta.

More information

SOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION

SOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION Synopses of Master Papers Bulletin of IISEE, 47, 19-24, 2013 SOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION Fabricio Moquete Everth* Supervisor:

More information

24th Seismic Research Review Nuclear Explosion Monitoring: Innovation and Integration

24th Seismic Research Review Nuclear Explosion Monitoring: Innovation and Integration PROBABILITY DENSITY FUNCTIONS FOR SECONDARY SEISMIC PHASE ARRIVALS Eric A. Bergman 1 and E. R. Engdahl 2 Global Seismological Services 1 and University of Colorado 2 Sponsored by Defense Threat Reduction

More information

LAB 6 SUPPLEMENT. G141 Earthquakes & Volcanoes

LAB 6 SUPPLEMENT. G141 Earthquakes & Volcanoes G141 Earthquakes & Volcanoes Name LAB 6 SUPPLEMENT Using Earthquake Arrival Times to Locate Earthquakes In last week's lab, you used arrival times of P- and S-waves to locate earthquakes from a local seismograph

More information

Supporting Information for An automatically updated S-wave model of the upper mantle and the depth extent of azimuthal anisotropy

Supporting Information for An automatically updated S-wave model of the upper mantle and the depth extent of azimuthal anisotropy GEOPHYSICAL RESEARCH LETTERS Supporting Information for An automatically updated S-wave model of the upper mantle and the depth extent of azimuthal anisotropy Eric Debayle 1, Fabien Dubuffet 1 and Stéphanie

More information

Subduction of the Chile Ridge: Upper Mantle Structure and Flow

Subduction of the Chile Ridge: Upper Mantle Structure and Flow GSA Supplemental Data Item 2010263 1 Subduction of the Chile Ridge: Upper Mantle Structure and Flow R. M. Russo 1, John C. VanDecar 2, Diana Comte 3, Victor I. Mocanu 4, Alejandro Gallego 1, and Ruth E.

More information

LLNL SEISMIC LOCATION: VALIDATING IMPROVEMENT THROUGH INTEGRATION OF REGIONALIZED MODELS AND EMPIRICAL CORRECTIONS

LLNL SEISMIC LOCATION: VALIDATING IMPROVEMENT THROUGH INTEGRATION OF REGIONALIZED MODELS AND EMPIRICAL CORRECTIONS LLNL SEISMIC LOCATION: VALIDATING IMPROVEMENT THROUGH INTEGRATION OF REGIONALIZED MODELS AND EMPIRICAL CORRECTIONS Craig A. Schultz, Megan P. Flanagan, Stephen C. Myers, Mike E. Pasyanos, Jennifer L. Swenson,

More information

Geographic boundary and shear wave velocity structure of the Pacific anomaly near the core mantle boundary beneath western Pacific

Geographic boundary and shear wave velocity structure of the Pacific anomaly near the core mantle boundary beneath western Pacific Earth and Planetary Science Letters 244 (2006) 302 314 www.elsevier.com/locate/epsl Geographic boundary and shear wave velocity structure of the Pacific anomaly near the core mantle boundary beneath western

More information

Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story

Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Petr Kolinsky 1, Leo Eisner 1, Vladimir Grechka 2, Dana Jurick 3, Peter Duncan 1 Summary Shear

More information

Geophysical Journal International

Geophysical Journal International Geophysical Journal International Geophys. J. Int. (2015) 203, 2189 2199 GJI Seismology doi: 10.1093/gji/ggv435 Regional seismic variations in the inner core under the North Pacific J.C.E. Irving 1 and

More information

Global P, PP, and PKP wave microseisms observed from distant storms

Global P, PP, and PKP wave microseisms observed from distant storms GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L23306, doi:10.1029/2008gl036111, 2008 Global P, PP, and PKP wave microseisms observed from distant storms Peter Gerstoft, 1 Peter M. Shearer, 1 Nick Harmon, 1 and

More information

Data Repository Item For: Kinematics and geometry of active detachment faulting beneath the TAG hydrothermal field on the Mid-Atlantic Ridge

Data Repository Item For: Kinematics and geometry of active detachment faulting beneath the TAG hydrothermal field on the Mid-Atlantic Ridge GSA Data Repository Item: 2007183 Data Repository Item For: Kinematics and geometry of active detachment faulting beneath the TAG hydrothermal field on the Mid-Atlantic Ridge Brian J. demartin 1*, Robert

More information

Geophysical Research Letters. Supporting Information for

Geophysical Research Letters. Supporting Information for Geophysical Research Letters Supporting Information for A P-wave based, on-site method for earthquake early warning S. Colombelli(1), A. Caruso(1), A. Zollo(1), G. Festa(1) and H. Kanamori(2) (1) Department

More information

Supplementary Online Material for. Seismic evidence for a chemically distinct thermochemical reservoir in Earth s deep mantle beneath Hawaii

Supplementary Online Material for. Seismic evidence for a chemically distinct thermochemical reservoir in Earth s deep mantle beneath Hawaii Supplementary Online Material for Seismic evidence for a chemically distinct thermochemical reservoir in Earth s deep mantle beneath Hawaii Authors: Chunpeng Zhao 1, Edward J. Garnero 1,*, Allen K. McNamara

More information

Small-scale ultralow-velocity zone structure imaged by ScP

Small-scale ultralow-velocity zone structure imaged by ScP JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. B1, 2056, doi:10.1029/2001jb001627, 2003 Small-scale ultralow-velocity zone structure imaged by ScP Sebastian Rost and Justin Revenaugh Department of Earth

More information

Analysis of the 29th May 2008 Ölfus earthquake and aftershock sequence using three-component t processing on ICEARRAY

Analysis of the 29th May 2008 Ölfus earthquake and aftershock sequence using three-component t processing on ICEARRAY Analysis of the 29th May 2008 Ölfus earthquake and aftershock sequence using three-component t processing on ICEARRAY Benedikt Halldórsson Steven J. Gibbons International Symposium on Strong-motion Earthquake

More information

3D IMAGING OF THE EARTH S MANTLE: FROM SLABS TO PLUMES

3D IMAGING OF THE EARTH S MANTLE: FROM SLABS TO PLUMES 3D IMAGING OF THE EARTH S MANTLE: FROM SLABS TO PLUMES Barbara Romanowicz Department of Earth and Planetary Science, U. C. Berkeley Dr. Barbara Romanowicz, UC Berkeley (KITP Colloquium 9/11/02) 1 Cartoon

More information

A glassy lowermost outer core

A glassy lowermost outer core Geophys. J. Int. (2009) 179, 374 380 doi: 10.1111/j.1365-246X.2009.04283.x A glassy lowermost outer core Vernon F. Cormier University of Connecticut, Physics Department, 2152 Hillside Road, Storrs, CT

More information

Online material: Updated 3-D tomography model and plotting script with examples

Online material: Updated 3-D tomography model and plotting script with examples Model Update January 2010: Upper Mantle Heterogeneity beneath North merica from Traveltime Tomography with Global and usarray Transportable rray Data Scott urdick, Robert D. van der Hilst, Frank L. Vernon,

More information

Vertical coherence in mantle heterogeneity from global seismic data

Vertical coherence in mantle heterogeneity from global seismic data GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl049281, 2011 Vertical coherence in mantle heterogeneity from global seismic data L. Boschi 1,2 and T. W. Becker 3 Received 11 August 2011; revised

More information

Robust Normal Mode Constraints on Inner Core Anisotropy From Model Space Search

Robust Normal Mode Constraints on Inner Core Anisotropy From Model Space Search Robust Normal Mode Constraints on Inner Core Anisotropy From Model Space Search Caroline Beghein, Jeannot Trampert, Faculty of Earth Sciences, Utrecht University, Post Office Box 821, 358 TA utrecht, Netherlands

More information

INNER-CORE ANISOTROPY AND ROTATION

INNER-CORE ANISOTROPY AND ROTATION Annu. Rev. Earth Planet. Sci. 2001. 29:47 69 Copyright c 2001 by Annual Reviews. All rights reserved INNER-CORE ANISOTROPY AND ROTATION Jeroen Tromp Seismological Laboratory, California Institute of Technology,

More information

Seismic ray path variations in a 3D global velocity model

Seismic ray path variations in a 3D global velocity model Physics of the Earth and Planetary Interiors 141 (2004) 153 166 Seismic ray path variations in a 3D global velocity model Dapeng Zhao, Jianshe Lei Geodynamics Research Center, Ehime University, Matsuyama

More information

Precise seismic-wave velocity atop Earth s core: No evidence for outer-core stratification. Authors: Catherine Alexandrakis *a and David W.

Precise seismic-wave velocity atop Earth s core: No evidence for outer-core stratification. Authors: Catherine Alexandrakis *a and David W. Precise seismic-wave velocity atop Earth s core: No evidence for outer-core stratification Authors: Catherine Alexandrakis *a and David W. Eaton b University of Calgary, Department of Geoscience, 2500

More information

Negative repeating doublets in an aftershock sequence

Negative repeating doublets in an aftershock sequence LETTER Earth Planets Space, 65, 923 927, 2013 Negative repeating doublets in an aftershock sequence X. J. Ma and Z. L. Wu Institute of Geophysics, China Earthquake Administration, 100081 Beijing, China

More information

Sharpness of the D 00 discontinuity beneath the Cocos Plate: Implications for the perovskite to post-perovskite phase transition

Sharpness of the D 00 discontinuity beneath the Cocos Plate: Implications for the perovskite to post-perovskite phase transition Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L03304, doi:10.1029/2007gl032465, 2008 Sharpness of the D 00 discontinuity beneath the Cocos Plate: Implications for the perovskite to

More information

Synthetic Seismograms for a Synthetic Earth

Synthetic Seismograms for a Synthetic Earth LUDWIG- MAXIMILIANS- UNIVERSITÄT MÜNCHEN Geophysics Department of Earth and Environmental Sciences Synthetic Seismograms for a Synthetic Earth Computing 3-D wavefields in mantle circulation models to test

More information

ANEWJOINTP AND S VELOCITY MODEL OF THE MANTLE PARAMETERIZED IN CUBIC B-SPLINES

ANEWJOINTP AND S VELOCITY MODEL OF THE MANTLE PARAMETERIZED IN CUBIC B-SPLINES ANEWJOINTP AND S VELOCITY MODEL OF THE MANTLE PARAMETERIZED IN CUBIC B-SPLINES Michael Antolik, Göran Ekström, Adam M. Dziewonski, Yu J. Gu, Jian-feng Pan, and Lapo Boschi Department of Earth and Planetary

More information

Local Magnitude Scale for the Philippines: Preliminary Results

Local Magnitude Scale for the Philippines: Preliminary Results Local Magnitude Scale for the Philippines: Preliminary Results Johnlery Pino Deximo 1 and Tatsuhiko Hara 2 1. Corresponding Author. Science Research Analyst, Philippine Institute of Volcanology and Seismology,

More information

Seismic Activity near the Sunda and Andaman Trenches in the Sumatra Subduction Zone

Seismic Activity near the Sunda and Andaman Trenches in the Sumatra Subduction Zone IJMS 2017 vol. 4 (2): 49-54 International Journal of Multidisciplinary Studies (IJMS) Volume 4, Issue 2, 2017 DOI: http://doi.org/10.4038/ijms.v4i2.22 Seismic Activity near the Sunda and Andaman Trenches

More information

ARTICLE IN PRESS. Error bars for the global seismic Q profile

ARTICLE IN PRESS. Error bars for the global seismic Q profile DTD 5 Earth and Planetary Science Letters xx (2004) xxx xxx www.elsevier.com/locate/epsl Error bars for the global seismic Q profile Joseph Resovsky*, Jeannot Trampert, R.D. Van der ilst 1 Department of

More information

A Glassy Lowermost Outer Core. (km/s) V S. Frequency (Hz) = Pa-s = 10 9 Pa-s e-03 1e-02 1e-01 1e+00.

A Glassy Lowermost Outer Core. (km/s) V S. Frequency (Hz) = Pa-s = 10 9 Pa-s e-03 1e-02 1e-01 1e+00. A Glassy Lowermost Outer Core Vernon F. Cormier (University of Connecticut) New theories for the viscosity of metallic melts at core pressures and temperatures, together with observations of translational

More information

Evidence from P-to-S mantle converted waves for a flat b660-kmq discontinuity beneath Iceland

Evidence from P-to-S mantle converted waves for a flat b660-kmq discontinuity beneath Iceland H4.SMR/1775-32 "8th Workshop on Three-Dimensional Modelling of Seismic Waves Generation, Propagation and their Inversion" 25 September - 7 October 2006 Evidence from P-to-S mantle converted waves for a

More information

Physics of the Earth and Planetary Interiors

Physics of the Earth and Planetary Interiors Physics of the Earth and Planetary Interiors 186 (2011) 59 69 Contents lists available at ScienceDirect Physics of the Earth and Planetary Interiors journal homepage: www.elsevier.com/locate/pepi Stratified

More information

SUMMARY INTRODUCTION THEORY

SUMMARY INTRODUCTION THEORY Stabilizing time-shift estimates in coda wave interferometry with the dynamic time warping method T. Dylan Mikesell, Alison Malcolm and Di Yang, Earth Resources Laboratory, MIT; Matt M. Haney, Alaska Volcano

More information

Thorne, Garnero, Jahnke, Igel, McNamara Supplementary Material - 1 -

Thorne, Garnero, Jahnke, Igel, McNamara Supplementary Material - 1 - Supplementary Material S1. Bin Location Map Location of 2.5 2.5 bins for which data was obtained to infer the ULVZ distribution. The preferred ULVZ model is drawn (yellow areas). Gray shaded regions indicate

More information

Three-dimensional structure of the African superplume from waveform modelling

Three-dimensional structure of the African superplume from waveform modelling Geophys. J. Int. (5) 161, 283 294 doi: 1.1111/j.1365-X.5.258.x Three-dimensional structure of the African superplume from waveform modelling Sidao Ni, 1,2 Don V. Helmberger 2 and Jeroen Tromp 2 1 School

More information

Performance of the GSN station KONO-IU,

Performance of the GSN station KONO-IU, Performance of the GSN station KONO-IU, 1991-2009 A report in a series documenting the status of the Global Seismographic Network WQC Report 2010:9 February 28, 2010 Göran Ekström and Meredith Nettles

More information