Inversion of Polar Motion Data: Chandler Wobble, Phase Jumps, and Geomagnetic Jerks Dominique Gibert, Jean-Louis Le Mouël
Total Page:16
File Type:pdf, Size:1020Kb
Inversion of polar motion data: Chandler wobble, phase jumps, and geomagnetic jerks Dominique Gibert, Jean-Louis Le Mouël To cite this version: Dominique Gibert, Jean-Louis Le Mouël. Inversion of polar motion data: Chandler wobble, phase jumps, and geomagnetic jerks. Journal of Geophysical Research : Solid Earth, American Geophysical Union, 2008, 113 (B10), pp.B10405. 10.1029/2008JB005700. insu-00348004 HAL Id: insu-00348004 https://hal-insu.archives-ouvertes.fr/insu-00348004 Submitted on 1 Apr 2016 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113, B10405, doi:10.1029/2008JB005700, 2008 Inversion of polar motion data: Chandler wobble, phase jumps, and geomagnetic jerks Dominique Gibert1 and Jean-Louis Le Moue¨l1 Received 18 March 2008; revised 11 July 2008; accepted 24 July 2008; published 23 October 2008. [1] We reconsider the analysis of the polar motion with a method totally different from the wavelet analysis used in previous papers. The total polar motion is represented as the sum of oscillating annual and Chandler terms whose amplitude and phase perturbations are inverted with a nonlinear simulated annealing method. The phase variations found in previous papers are confirmed with the huge phase change (3p/2) occurring in the 1926–1942 period and another less important one (p/3) in the 1970–1980 epoch. The best Chandler period, Tc, is found to be 434 ± 0.5 mean solar days. The epochs of significant phase changes found in the Chandler wobble are also those where important geomagnetic jerks occur in the secular variation of the geomagnetic field. Turbulent viscous friction produced by small-scale topography at the core-mantle boundary might be responsible for the observed phase jumps. Citation: Gibert, D., and J.-L. Le Moue¨l (2008), Inversion of polar motion data: Chandler wobble, phase jumps, and geomagnetic jerks, J. Geophys. Res., 113, B10405, doi:10.1029/2008JB005700. 1. Introduction and constitutes a new clue to the role of core-mantle coupling onto the Earth’s polar motion. [2] The movement of the Earth’s pole has been observed [4] In the present paper, we reconsider the analysis of the for more than 150 years, and, at the end of the 19th century, Earth’s polar motion with a new method totally different Chandler [1891a, 1891b] showed that this motion is the from the wavelet analysis used in our earlier papers [Gibert sum of two periodic oscillations, the forced annual oscilla- et al., 1998; Bellanger et al., 2002]. Also, the inverse tion and the Chandler wobble with a period close to 434 approach followed in the present study is more direct in mean solar days [e.g., Guinot, 1972; Lambeck, 1980; Wahr, the sense that the detrended total polar motion is analyzed 1988; Cazenave and Feigl, 1994]. Both the amplitude and without the complicated preprocessing stages of the wavelet the phase of the Chandler wobble vary with time [e.g., analysis. Moreover, we hope this new and simpler method Lambeck, 1980], and Danjon and Guinot [1954] first will be more easily used by independent teams in order to reported a phase change of about 150° which occurred perform analysis similar to the one presented here. around 1926 and was achieved in less than 10 years [5] In the remaining of the present paper, we first present [Guinot, 1972]. This phase jump remained poorly docu- the polar motion series and describe the model used to mented until other jumps were reported by Gibert et al. compute realistic synthetic series to be compared with the [1998] and suggested to be in temporal coincidence with data. Then, we detail the parameterization of the inverse geomagnetic jerks [Gibert et al., 1998; Bellanger et al., problem where the amplitudes of the annual and Chandler 2002]. Similar phase jumps have been detected in the Free components of the polar motion are to be retrieved together Core Nutation time series and also found in coincidence with the phase curve of the Chandler term. Next, the with geomagnetic jerks [Shirai et al., 2005]. nonlinear inversion is performed with the simulated anneal- [3] Different hypotheses have been proposed as the ing method and a Bayesian assessment of the uncertainty excitation mechanism of the Chandler wobble (see Dehant attached to the estimated parameters is done with the and de Viron [2002] for a short review), the main ones being Metropolis algorithm. the atmosphere [e.g., Aoyama and Naito, 2001], the pres- sure at the bottom of the oceans [Gross, 2000; Brzezin˜ski and Nastula, 2002; Gross et al., 2003], the earthquakes 2. Modeling of the Polar Motion [Mansinha and Smylie, 1967; Gross, 1986], or core motions 2.1. Polar Motion Data [Wahr,1988;Dickman, 1993]). The proposed temporal [6] The total polar motion, m(t), provided by the Interna- coincidence between geomagnetic jerks and phase jumps tional Earth Rotation Service (IERS) (data series eopc01) is in the Chandler wobble constitutes an exciting challenge shown in Figure 1a. Let us recall that the coordinate system is left-handed with the X axis coinciding with the Green- wich meridian and the Y axis corresponding to the W90 1 meridian. This signal is the superimposition of a smooth Institut de Physique du Globe, CNRS/INSU UMR 7154, Paris, France. secular drift d(t) and of an oscillating term o(t) resulting Copyright 2008 by the American Geophysical Union. from interferences between the annual prograde and retro- 0148-0227/08/2008JB005700$09.00 grade components, ap(t) ar(t), and the prograde Chandler B10405 1of8 B10405 GIBERT AND LE MOUE¨ L: INVERSION OF POLAR MOTION DATA B10405 wobble c(t)[Lambeck, 1980]. In previous studies [Gibert et [10] As can be seen in equation (2), the amplitudes are al., 1998; Bellanger et al., 2002], we showed that the annual linear parameters, while the phase variation is a nonlinear prograde component is dominant with respect to the retro- one. In what follows, the inversion will not deal directly grade one which, however, remains significant (see the inset with these unknown functions, but instead will be per- in Figure 2c for an illustration of the importance of the formed for primary parameters from which the Apa(t), Ara(t), retrograde annual component). In the same studies, it was Ac(t) and Df(t) functions are deduced. More explicitly, the also shown that the Chandler wobble is affected by a small amplitude functions are obtained by linearly interpolating number of phase perturbations Df(t) whose temporal shape coarsely sampled amplitude functions. In practice, this reads may be modeled as sigmoid phase jumps with a short ÂÃÀÁÀÁÀÁ duration of the order of one year. Accounting for these ApaðÞ¼LIt Apa;1; t1 ; ÁÁÁ; Apa;k ; tk ; ÁÁÁ; Apa;N ; tN ; ð4Þ observations the total polar motion may be written as, ÂÃÀÁÀÁÀÁ mtðÞ¼otðÞþdtðÞþntðÞ AraðÞ¼LIt Ara;1; t1 ; ÁÁÁ; Ara;k ; tk ; ÁÁÁ; Ara;N ; tN ; ð5Þ ¼ apðÞþt arðÞþt ctðÞþdtðÞþntðÞ; ð1Þ ÂÃÀÁÀÁÀÁ where n(t) represents the noise presents in the data AcðÞ¼LIt Ac;1; t1 ; ÁÁÁ; Ac;k ; tk ; ÁÁÁ; Ac;N ; tN ; ð6Þ (Figure 1c). [7] In the present study we remove the secular drift before where LI represents the linear interpolation operator, and inverting the data. Owing to the smooth nature of d(t), we where the dates tk are fixed and regularly distributed with a represent it by the least squares fitted polynomials of degree sampling interval of 3 months. Consequently, for ampli- 5 shown in Figure 1a. No significant differences were tudes, the parameters to be actually determined are the observed when using either fourth and sixth degree poly- elements of the three sets {Apa,k, k =1,N}, {Ara,k, k =1,N}, nomials. The noisy oscillating part, o(t)+n(t), is obtained and {Ac,k, k =1,N}. by subtracting the fitted polynomials from the total polar [11] An identical parameterization is adopted for the motion and is shown in Figure 1b. phase variation function Df(t), [8] We remark that IERS data at face value, and the polar ÂÃÀÁÀÁÀÁÀÁ motion series shown in Figure 1a results from the merging 0 0 0 0 DfðÞ¼LIt Df1; t1 ; Df2; t2 ; ÁÁÁ; Dfk ; tk ; ÁÁÁ; DfM ; tM of series with data of different nature (for more details, see 7 IERS Web site http://www.iers.org). For instance, it cannot ð Þ be eliminated that some inconsistencies affect the data in the 0 1970–1980 period due to the merging of classical astrome- with the important difference that the dates t k are not fixed tric and space geodetic measurements [Gross and Vondra´k, but, instead, are to be determined together with the phase 1999]. variations Dfk. This parameterization makes the inversion highly non linear but, as will be seen later, it permits us to 2.2. Parameterization of the Polar Motion Model eventually represent the phase function with a very small [9] The signal to be inverted is the o(t) part of the total number M of parameters. polar motion and, in the present study, the model of signal [12] To summarize, the model to be determined is the set used in the forward problem reads, of parameters, ÈÉ M¼ A ; ÁÁÁ; A ; ÁÁÁ; A [paÈÉ;1 pa;k pa;N m 2ipt o ðÞ¼t ApaðÞt exp þ iDfpa Ara;1; ÁÁÁ; Ara;k ; ÁÁÁ; Ara;N Ta [ ÈÉ A ; ÁÁÁ; A ; ÁÁÁ; A 2ipt [ c;1 c;k c;N þ AraðÞt exp À À iDfra Ta fgDf1; ÁÁÁ; Dfk ; ÁÁÁ; DfM [ ÈÉ 2ipt 0 0 0 þ AcðÞt exp þ iDfðÞt ; ð2Þ t1; ÁÁÁ; tk ; ÁÁÁ; tM Tc [ ÈÉ Dfpa; Dfra ; ð8Þ where the real amplitudes Apa(t), Ara(t), Ac(t), the phase from which the functions Apa(t), Ara(t), Ac(t) and Df(t)are variation Df(t), and the constant phase shifts Dfpa, Dfra of the annual components are to be determined.