Layout 6 A case study of detecting the triplet of 3S1 using superconducting gravimeter records with an alternative data preprocessing technique Wen-Bin Shen1,2,3*, Bo Wu1 1 Wuhan University, School of Geodesy and Geomatics, Department of Geophysics, Wuhan, China 2 Wuhan University, State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan, China 3 Wuhan University, Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan, China ANNALS OF GEOPHYSICS, 55, 2, 2012; doi: 10.4401/ag-4944 ABSTRACT Due to their very low noise levels in the low frequency band (<1 mHz), superconducting gravimeters (SGs) are particularly suitable to observe long-period free oscillations of the Earth. This case study is dedicated to the detection of the triplet of the seismic normal mode 3S1 that was excited by the December 26, 2004, Sumatra-Andaman earthquake (Mw = 9.3). Using SG records, the Hilbert-Huang transformation is used as an alternative data preprocessing technique, instead of the traditional detiding method. After removal of atmospheric pressure effects from the original SG records, we applied the Hilbert-Huang transformation to the SG residues, to select the signals that included the frequency band of interest, and to construct a new data series. Then, by applying the multi-station experimental technique to five 273-h-long common new data series recorded at different SG stations, we clearly observed all of the three singlets of the mode 3S1, with the central singlet more evident compared to previous studies. Observations of the low-frequency modes 3S1 (n = 0, 1, 2, ...; l = = 1,2, ...) provide constraints on the inner and outer core structure. This case study provides an alternative data-preprocessing approach to observe the splitting frequencies of the low-frequency mode type 3S1 (n = 0, 1, 2, ...). 1. Introduction After a large earthquake, various seismic normal modes (free oscillations) of the Earth are excited for several days to months (e.g., Alterman et al. 1974, Widmer-Schnidrig 2003). The typical length of the time that a seismic normal mode lasts for depends on its quality factor Q (which is inversely proportional to the ratio of the energy loss per cycle to peak strain energy), whereby the larger the Q, the longer the vibration of the normal mode (e.g., Aki and Richards 1980, Chao and Gilbert 1980). The frequencies of the modes are closely related to the structure of the Earth. For instance, the triplet of the seismic normal mode 3S1 (generally the normal modes below 3 mHz) is sensitive to the density, compressibility, and velocity of P-waves in the mantle and core [e.g., Ritzwoller and Lavely 1995, Resovsky and Ritzwoller 1998, see also http://phys-geophys.colorado.edu/geophysics/nm.dir/3s/ 3s1.kerplot.gif]. Theoretically, if the Earth is taken as an idealized spherically symmetric, non-rotating, purely elastic, isotropic body, modes nS l m with the same n and l have the same eigenfrequency, which is referred to as degeneracy [Dahlen and Tromp 1998]. Here n is the radial overtone number, l the angular degree, and m the azimuthal order. However, for the real Earth, its rotation, ellipticity and lateral heterogeneity remove the degeneracy, which results in the appearances of split peaks [e.g., Alterman et al. 1974, Dahlen and Sailor 1979, Rogister 2003]. The splitting of the modes below 1 mHz is highly sensitive to the three-dimensional density structure of the Earth mantle and core [e.g., Okal 1978, Ritzwoller and Lavely 1995, Widmer-Schnidrig 2003], and therefore, observations of the splitting of modes below 1 mHz might help to constrain Earth models. Due to the Global Geodynamics Project (GGP) [e.g., Crossley et al. 1999, Crossley and Hinderer 2008], superconducting gravimeters (SGs) are extremely sensitive to gravity variations that are related to various geophysical processes (e.g., tides, inner core wobble, tectonic activity, Earth free oscillations). They also have very low noise levels in the low frequency band, and thus they are particularly suitable for the observation of long-period signals of interest. Previous studies [e.g., Rosat et al. 2003a, 2004; Hu et al. 2006a, b, Xu and Sun 2009] have shown that for frequencies below 1 mHz, the SGs can reach a better signal-to-noise ratio than most of the broadband seismometers. Hence, SGs have specific advantages for the study of the gravest normal modes. The 3S1 multiplet was first observed by Chao and Gilbert [1980], using seven records at the International Deployment Accelerometer spring gravimeter stations after the 1977 DATA AND EXPERIMENT DESCRIPTIONS 293 Article history Received December 29, 2010; accepted December 7, 2011. Subject classification: Seismic normal mode, Superconducting gravimeters data, Hilbert-Huang transformation, Multi-station experiment, 3S1 triplet. Indonesia earthquake (Mw = 7.8) and then by Roult et al. (2010) using 247 recordings at 157 Federation of Digital Seismograph Networks broadband seismometer stations after the 2004 Sumatra earthquake (Mw = 9.3). However, following a closer examination of the spectra shown in Figure 10 of Roult et al. (2010), we find that the middle spectrum line (m = 0) of the mode 3S1 is very weak and cannot be easily observed. Here, it will be shown that by combining only five SG records at three stations after the 2004 Sumatra-Andaman earthquake with an alternative analysis technique, instead of traditional detiding process for data preprocessing, we can clearly resolve all three singlets of the seismic normal mode 3S1. 2. Methods The preprocessing of the SG raw data generally involves tides removal and air correction (necessary for better observation of low frequency modes below 1 mHz), in addition to correcting for instrumental errors, such as spikes, gaps, and abrupt offsets. The detiding process is usually done 294 SHEN ET AL. Figure 1. Illustration of the Hilbert-Huang transformation using synthetic data. (a) Eight IMFs labeled by the sequences they are sifted, which are decomposed from the synthetic signal with frequencies IMF1-8 (2.5, 2.1, 1.92, 1.2, 0.31, 0.24, 0.03, 0.02) mHz, and amplitudes (1.5, 1.3, 0.8, 0.55, 0.7, 1.0, 4.0, 3.0) nm/s2, respectively. (b) Instantaneous frequencies corresponding to each of IMF1-8 according to time. (a) (b) 295 by subtracting syntheses computed from certain tide models, although removing the tides by this method is not very thorough, as tidal effects on the records can be different for individual places, and the tidal factors used to compute synthetic tides are just average values obtained over several years. Alternatively, the SG data can be bandpassed to a certain frequency range of interest, although it must be noted that the presence of Gibbs phenomenon and marginal effects might result in some undesired effects in the original observations. While wavelet-based detection for weak signals [Hu et al. 2006a, b, Rosat et al. 2007] has great potential, the uniform frequency resolution will show inferior resolution in time and frequency compared to the Hilbert analysis representation [see Figure 1 of Huang et al. 1999]. Herein, the Hilbert-Huang transformation (HHT) analysis is used to locate our desired frequency band directly, which will generate a new data series that will be naturally free from influences of long-period signals, e.g., tidal effects and some other possible influences. To enhance each singlet peak, the multi-station experiment (MSE) technique [Cummins et al. 1991] will be applied. The following two paragraphs are short descriptions of these two techniques. The HHT analysis technique was first proposed by Huang et al. [1998], to decompose a complicated data series into a finite, and often small number of, intrinsic mode functions (IMFs) that are based on empirical mode decomposition (EMD) (i.e., the sifting process) [see details in Huang et al. 1998]. The essence of the method is to empirically identify the intrinsic oscillatory modes by their characteristic time scales in the data, and then to decompose the data accordingly. When the first IMF is obtained, the sifting process will be repeated again to find the next IMF in the residue. If the last residue becomes a monotonic function from which no more IMFs can be extracted, the sifting process stops. By virtue of the definition of IMF, the earlier sifted IMFs (signals) generally have higher frequencies. Figure 1 illustrates the method with synthetic data. Figure 1a shows the extracted IMFs based on EMD, and the corresponding instantaneous frequency of each IMF is shown in Figure 1b. Using these, by zooming in (horizontal or vertical zoom), we can roughly determine the frequency range of each IMF. By this criterion, with concentration on one or a few frequency bands of interest, we can select certain IMFs for further study. Figure 1c shows the corresponding Fourier amplitude spectrum of each IMF in Figure 1a, and is just a verification of our judgment from Figure 1b. The MSE technique was first proposed by Cummins et al. [1991], and used by Courtier et al. [2000] to detect the triplet of the Slichter mode 1S1 [Slichter 1961; see also e.g. Smylie 1992]. It was also used by Rosat et al. [2003b] to detect the mode 2S1 and by Rosat et al. [2006] and Guo et al. [2006] to search for the 1S1 triplet. This method takes into account the temporal and spatial properties of degree-one spheroidal modes to generate three new time series, each of which contains only one of the prograde equatorial (m=–1), retrograde equatorial (m=1) and axial (m=0) signals. Spectral analysis is then applied to every time series. In the present study, the MSE technique is used to effectively isolate and enhance the three singlets that correspond to the mode nS l m (m=–1, 0, 1). 3. Data and data analysis The datasets are corrected minute data (gaps and disturbances filled with synthetic signals by station operator after decimation to 1 min), downloaded from the GGP data centre (http://www.eas.slu.edu/GGP/ggphome.html), and corrected for the local atmospheric pressure effect using a nominal constant admittance of –3 nms–2/hPa. The air pressure correction is necessary to better detect the frequencies less than 1 mHz [Zürn and Widmer 1995]. As a case study, five records from three SG stations, Canberra (C1), Bad Homburg (H1, H2), Sutherland (S1, S2), are used in the present study. The length of the records is 273 h, about one Q-cycle, as the length suggested by Dahlen [1982] for better frequency estimation of a certain mode. After correction of the local atmospheric pressure effect, the gravity residues are decomposed into a series of IMFs, as shown by Figure 2a, based on the EMD [Huang et al. 1998], and after applying Hilbert transformation to every IMF and by examining the variation of the instantaneous frequencies with time of the IMFs (see Figure 2b), we chose the correct IMFs that included the signal of the mode 3S1. In the present study, we chose the first three IMFs of all of the records, and accordingly constructed five common new data series, which A CASE STUDY OF DETECTING THE TRIPLET OF 3S1 Figure 1. (c) Fourier amplitude spectra for the corresponding IMF1-8, demonstrating that each IMF in Figure 1a indeed contains certain frequencies, as can be seen from Figure 1b. Therefore certain IMFs can be selected to study the signal of interest. included the signal of the mode 3S1, which will be further used to detect the 3S1 triplet based on the MSE technique. As the typical frequency of 3S1 is very close to that of 1S3, and as it also obeys certain selection rules [Alterman et al. 1974], the splitting frequencies of 3S1 might be seriously contaminated by those of 1S3, because the typical frequencies (m = 0) of 1S3 and 3S1 are 0.944139 mHz and 0.944364 mHz, respectively [Resovsky and Ritzwoller 1998]. However, as the quality factor Q of 3S1 (826.9) is almost three times of that of 1S3 (282.7) [Dziwonski and Anderson 1981], we can consider the Earth as a filter and use the records with a later starting point after the earthquake, to weaken the influence of 1S3. As shown by Figure 3, we find that about 56 h after the Sumatra-Andaman event, the amplitude of 1S3 is largely reduced compared to that of five hours after the event. Indeed, based on the decay law A(t)=A(t0)exp[–r(t – t0)f/Q] [Aki and Richards 1980], 56 h after the event, the amplitudes of the 1S3 and 3S1 decay to 14.7% and 51.8%, respectively, of those standing 5 h after the event. 4. Results We applied the MSE technique to the five common new data series (i.e., the new data series constructed by the HHT technique, as indicated above), and we obtain three time-series, each of which contains only one of the SHEN ET AL. 296 (a) Figure 2. The Hilbert-Huang transformation applied to the SG data following correction for pressure, as recorded at Bad Homberg (H2). (a) Extracted IMFs (1-7) of the SG data based on EMD. (b) Instantaneous frequencies corresponding to each of IMF1-7 of panel (a). (b) 297 prograde, retrograde and axial modes that have been successfully enhanced in the corresponding amplitude spectra, as shown in Figure 4a (m=–1), b (m=1), and c (m=0), respectively. Hence, all of the three singlets of 3S1 have been clearly extracted and are close to the frequencies computed from PREM. By fitting a synthetic Lorentzian resonance function [e.g., Dahlen and Tromp 1998] to each singlet of the spectrum, we obtained the three splitting frequencies, as listed in Table 1. The errors given in the present study are standard deviations of the estimated frequencies obtained by least-square fitting with five to eight points. The exact number of points used for estimation of a certain singlet might vary, depending upon how many points around the peak fit the Lorentzian resonance function best. In Table 1, the theoretical predictions are based on the model PREM-re [Roult et al. 2010] and the model PREM- re+SAW12D [Li and Romanowicz 1996, He and Tromp, 1996], respectively. Both models are based on PREM, but the former includes only the Earth rotation and ellipticity, and the A CASE STUDY OF DETECTING THE TRIPLET OF 3S1 Figure 3. Comparison of the Fourier amplitude spectrum of the 273-h SG record from Bad Homburg (H2) starting 5 h (top slot) and 56 h (bottom slot) after the earthquake. The spectra were performed on the first two IMFs, after corrections for pressure. The vertical dotted lines denote degenerate frequencies for the model PREM. (a) (b) (c) Figure 4. Amplitude spectra of the three singlets of the mode 3S1 obtained from the multi-station experimental analysis technique. (a) m=–1. (b) m=1. (c) m=0. latter includes not only the rotation and ellipticity, but also the three-dimensional elastic structure heterogeneity of the mantle. Chao and Gilbert [1980] first observed the triplet frequencies of 3S1 by applying the spherical harmonic stacking technique [Buland et al. 1978] to seven spring gravimeters records, Roult et al. [2010] also observed the triplet frequencies of the mode 3S1 based on 247 seismic records by taking simply the average value of the frequencies corresponding to the mode of interest. We note that, the triplet frequencies of the mode 3S1 were observed by only the two above-mentioned studies, to our present knowledge. For instance, although Resovsky and Ritzwoller [1998] observed the frequency of 3S1 (m = 0) using a dataset of more than 4,500 seismograms, they did not provide the frequencies of 3S1 (m = ±1). As a comment here, we note that in the study of Roult et al. [2010], the middle spectrum line (m = 0) of the mode 3S1 is very weak. In our study, the middle spectrum line (m = 0) of 3S1 was clearly extracted from the records due to the combination of the HHT and MSE techniques. As can be noted from Table 2, the observed 3S1 triplet in the present study are very close to the model PREM- re+SAW12D predictions, except that the central singlet deviates from the theoretical value a little more, as do the results of Chao and Gilbert [1980] and Roult et al. [2010]. The causes for this are not clear. Besides, all of the observed frequencies [given by Chao and Gilbert 1980, Roult et al. 2010, and the present study] corresponding to m=±1 are higher than the model PREM-re predictions, and smaller than the model PREM-re + SAW12D predictions. This implies that the mode 3S1 is very sensitive to the physical parameters of the mantle and core, and the observed triplet of 3S1 might provide significant constraints on the three-dimensional structure of the Earth. We can preliminarily conclude that based on the observations of the 3S1 triplet (by previous studies and the present study, see Tables 1 and 2), both the models PREM-re and PREM-re + SAW12D might need to be adjusted as the predictions based on these models deviate from the «most likely real values» (the observations given by different authors) by slightly systematic shifts of about �–3.0 nHz and about 3.0 nHz in the frequency-decreasing and -increasing directions, respectively. In addition, from Table 1 we note that as model PREM-re + SAW12D predicted, the splitting width of the mode 3S1 is 3.22 nHz, close to the estimate (3.21 nHz) given by He and Tromp (1996), of 3.21 nHz, and the observed splitting widths given by Chao and Gilbert [1980], Roult et al. [2010], and the present study are 2.90 nHz, 3.23 nHz, and 3.26 nHz, respectively. 5. Conclusions As the mode 3S1 is quite sensitive to the mantle and outer core, the observation of the splitting frequencies of the mode 3S1 can provide significant information of the deep SHEN ET AL. 298 Model or author(s) m=–1 (mHz) m=0 (mHz) m=1 (mHz) PREM-re 0.942267 0.944215 0.945472 PREM-re + SAW12D 0.942925 0.944914 0.946144 Chao and Gilbert (1980) 0.94270 ±5.5#10–5 0.94535 ±9.0#10–5 0.94563 ±4.0#10–5 Resovsky and Ritzwoller (1998)a _ 0.944364±1.0#10–4 _ Roult et al. (2010) 0.94256 ±1.241#10–4 0.94419 ±3.444#10–4 0.94579 ±1.493#10–4 Present study 0.942598 ±4.25#10–5 0.944113 ±2.65#10–4 0.945864 ±2.13#10–4 a As a reference, the observation of the frequency (m = 0) of 3S1 given by Resovsky and Ritzwoller [1998], who did not provide the triplet 3S1. Table 1. Comparison of the observed triplet frequencies of 3S1 from previous studies with the present study, according to the model predictions. Table 2. Differences between the observations and the model predictions for the 3S1 triplet. Differences Author(s) m=–1 (mHz) m=0 (mHz) m=1 (mHz) Obs - PREM-re Chao and Gilbert (1980) 4.33#10–4 1.14#10–3 1.58#10–4 Resovsky and Ritzwoller (1998) – 1.49#10–4 – Roult et al. (2010) 2.93#10–4 –2.50#10–5 3.18#10–4 Present study 3.31#10–4 –1.02#10–4 3.92#10–4 Obs - PREM-re + SAW12D Chao and Gilbert (1980) –2.25#10–4 4.36#10–4 –5.14#10–4 Resovsky and Ritzwoller (1998) – –5.50#10–4 – Roult et al. (2010) –3.65#10–4 –7.24#10–4 –3.54#10–4 Present study –3.27#10–4 –8.01#10–4 –2.80#10–4 299 interior of the Earth, and consequently can improve the Earth model. The application of the HHT and MSE analysis techniques to five common SG records leads to a clear observation of the triplet frequencies of the mode 3S1. The observed 3S1 triplet frequencies from the SG records at the Canberra (CB), Sutherland (S1, S2) and Bad Homburg (H1, H2) stations are in close agreements with the theoretical predictions provided by the model PREM-re, and especially by the model PREM-re + SAW12D. The present study shows that the approaches (i.e., a combination of the HHT and MSE analysis techniques) proposed here are effective for the detection of the splitting frequencies of the mode 3S1, as only five SG records were used, and consequently, it can potentially be applied to the detection of the Slichter triplet 1Sl m (m=–1, 0, 1), of which the claimed observations are still controversial. Acknowledgements. The authors are grateful to the GGP station managers who provided the SG datasets used in this study. The authors also thank the Reviewers and the Associated Editor for their valuable comments and suggestions, which have greatly improved the manuscript. This study was supported by NSFC (grant No. 41174011, No. 40974015, and No. 40637034), and the Special Scientific Research Fund of State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing (China). References Aki, K. and P.G. Richards (1980). Quantitative seismology: theory and methods. W. H. Freeman, San Francisco. Alterman, Z.S., Y. Eyal and A.M. Merzeron (1974). On free oscillations of the Earth, Geophys. Survey, 1, 409-428. Buland, R., J. Berger and F. Gilbert (1978). Improved resolu- tion of complex eigenfrequencies in analytically continued seismic spectra, Geophys. J. R. Astron. Soc., 52, 457-470. Chao, B.F. and F. Gilbert (1980). Autoregressive estimation of complex eigenfrequencies in low frequency seismic spectra, Geophys. J. R. Astron. Soc., 63, 641-657. Courtier, N., B. Ducarme, J. Goodkind, J. Hinderer, Y. Ima- nishi, N. Seama, H. Sun, J. Merriam, B. Bengert and D.E. Smylie (2000). Global superconducting gravimeter ob- servations and the search for the translational modes of the inner core, Phys. Earth Planet. Int., 117, 3-20. Crossley, D. and J. Hinderer (2008). Report of GGP activities to Commission 3, completing 10 years for the worldwide net- work of superconducting gravimeters, In: Observing our changing Earth, IAG Symposia, 133, 511-521. Crossley, D.G., J. Hinderer, G. Sasula, O. Francis, H.T. Hsu, Y. Imanishi, G. Jentzsch, J. Kaarianen, J. Merriam, B. Meurers, J. Neumeyer, B. Richter, K. Shibuya, T. Sato and T. van Dam (1999). Network of superconducting gravimeters benefits a number of disciplines, EOS Trans. AGU, 80, 121/125-126. Cummins, P., J.M. Wahr, D.C. Agnew and Y. Tamura (1991). Constraining core undertones using stacked In- ternational Deployment Accelerometer gravity records, Geophys. J. Int., 106, 189-198. Dahlen, F.A. (1982). The effect of data windows on the esti- mation of free oscillations parameters, Geophys. J. R. Astron. Soc., 69, 537-549. Dahlen, F.A. and J. Tromp (1998). Theoretical Global Sei- smology, Princeton Univ. Press, Princeton, pp. 234-235. Dahlen, F.A. and R. Sailor (1979). Rotational and elliptical splitting of the free oscillations of the Earth, Geophys. J. R. Astr. Soc., 58, 609-623. Dziewonski, A.M. and D.L. Anderson (1981). Preliminary re- ference Earth model, Phys. Earth Planet. Int., 25, 297-356. Guo, J.Y., O. Dierks, J. Neumeyer and C.K. Shum (2006). Weighting algorithms to stack superconducting gravi- meter data for the potential detection of the Slichter modes, J. Geodyn., 41, 326-333. He, X. and J. Tromp (1996). Normal-mode constraints on the structure of the Earth, J. Geophys. Res., 110(B9), 20053-20082. Hu, X.-G, L.-T. Liu, J. Hinderer, H.T. Hsu and H.-P. Sun (2006a). Wavelet filter analysis of atmospheric pressure effects in the long-period seismic mode band, Phys. Earth Planet. Int., 154, 70-84. Hu, X, L. Liu, H. Sun, H. Xu, J. Hinderer and X. Ke (2006b). Wavelet filter analysis of splitting and coupling of sei- smic normal modes below 1.5 mHz with superconduc- ting gravimeter records after the December 26, 2004 Sumatra earthquake, Science in China Series D: Earth Sciences, 49, 1259-1269. Huang, N.E., Z. Shen, S.R. Long, M.C. Wu, H.H. Shih, Q. Zheng, N.-C. Yen, C.C. Tung and H.H. Liu (1998). The empirical mode decomposition and the Hilbert spec- trum for nonlinear and non-stationary time series ana- lysis, Proc. Roy. Soc. Lond., A454, 903-995. Huang, N.E., Z. Shen and S.R. Long (1999). A new view of nonlinear water waves: the Hilbert spectrum, Annu. Rev. Fluid Mech., 31, 417-457. Li, X.D. and B. Romanowicz (1996). Global shear velocity model developed using nonlinear asymptotic coupling theory, J. Geophys. Res., 101, 22245-22272. Okal, E.A. (1978). A physical classification of the Earth’s spheroidal modes, J. Geophys. Earth, 26, 75-103. Resovsky, J.S. and M.H. Ritzwoller (1998). New and refi- ned constrains on three-dimensional Earth structure from normal modes below 3 mHz, J. Geophys. Res., 103, 783-810. Ritzwoller, M.H. and E.M. Lavely (1995). Three-dimensio- nal seismic models of the Earth’s mantle, Rev. Geo- phys., 33, 1-66. Rogister, Y. (2003). Splitting of seismic free oscillations and of the Slichter triplet using the normal mode theory of a rotating, ellipsoidal Earth, Phys. Earth Planet. Int., 140, 169-182. A CASE STUDY OF DETECTING THE TRIPLET OF 3S1 Rosat, S., J. Hinderer, D. Crossley and L. Rivera (2003a). The search for the Slichter mode: comparison of noise levels of superconducting gravimeters and investigation of a stacking method, Phys. Earth Planet. Int., 140, 183-202. Rosat, S., J. Hinderer and L. Rivera (2003b). First observa- tion of 2S1 and study of the splitting of the football mode 0S2 after the June 2001 Peru earthquake of ma- gnitude 8.4, Geophys. Res. Lett., 30, 10-11. Rosat, S., J. Hinderer, D. Crossley and J.P. Boy (2004). Per- formance of superconducting gravimeters from long- period seismology to tides, J. Geodyn., 38, 461-476. Rosat, S., P. Sailhac and P. Gegout (2007). A wavelet-based detection and characterization of damped transient waves occurring in geophysical time-series: theory and application to the search for the translational oscilla- tions of the inner core, Geophys. J. Int., 171, 55-70. Rosat, S., Y. Rogister, D. Crossley and J. Hinderer (2006). A search for the Slichter triplet with superconducting gra- vimeters: impact of the density jump at the inner core boundary, J. Geodyn., 41, 296-306. Roult, G., J. Roch and E. Clévédé (2010). Observation of split modes from the 26th December 2004 Sumatra-An- daman mega-event, Phys. Earth Planet. Int., 179, 45-59. Slichter, L.B. (1961). The fundamental free mode of the Earth’s inner core, Proc. Nat. Acad. Sci., 47(2), 186-190. Smylie, D.E. (1992). The inner core translational triplet and the density near Earth’s center, Science, 255, 1678-1682. Widmer-Schnidrig, R. (2003). What can superconducting gravimeters contribute to normal mode seismology?, Bull. Seism. Soc. Am., 93(3), 1370-1380. Xu, J. and H. Sun (2009). Temporal variations in free core nutation period, Earthq. Sci., 22, 331−336. Zürn, W. and R. Widmer (1995). On noise reduction in ver- tical seismic records below 2mHz using local barome- tric pressure. Geophys. Res. Lett., 22, 3537-3540. *Corresponding author: Wen-Bin Shen, Wuhan University, School of Geodesy and Geomatics, Department of Geophysics, Wuhan, China; e-mail: wbshen@sgg.whu.edu.cn. © 2012 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. SHEN ET AL. 300