pith. sign in

arxiv: 2412.07868 · v3 · submitted 2024-12-10 · ⚛️ physics.geo-ph · astro-ph.IM· physics.data-an

Filling the gap in the IERS C01 polar motion series in 1858.9-1860.9

Pith reviewed 2026-05-23 07:29 UTC · model grok-4.3

classification ⚛️ physics.geo-ph astro-ph.IMphysics.data-an
keywords polar motionIERS C01gap fillingsingular spectrum analysisEarth orientation parametersChandler wobbleannual wobble
0
0 comments X

The pith

The 2-year gap in the IERS C01 polar motion series can be filled using either a parametric model or singular spectrum analysis.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The IERS C01 series is the longest reliable record of Earth's polar motion. It has a gap in the Y_p coordinate from 1858.9-1860.9 that prevents even spacing. Two methods were tested to fill it: a parametric astronomical model with bias and wobbles of varying amplitudes, and singular spectrum analysis. The methods agree within errors, but SSA is preferred for using a fuller model of polar motion.

Core claim

The gap in the Y_p series of the IERS C01 polar motion data can be filled reliably by both a parametric model consisting of bias plus Chandler and annual wobbles with linear amplitude changes and singular spectrum analysis. The results from both methods agree within the uncertainties of the original series, but the SSA-based filling is preferable because it rests on a more complete description of polar motion behavior.

What carries the argument

Singular Spectrum Analysis (SSA) applied as a data-driven reconstruction of polar motion signals across the data gap

If this is right

  • The completed series provides an evenly spaced record suitable for studies of long-term polar motion variations.
  • Both filling methods produce values consistent with each other within the reported errors of the IERS C01 series.
  • SSA incorporates additional periodic components of polar motion beyond the basic Chandler and annual terms used in the parametric model.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • SSA reconstruction could be tested on other gaps or incomplete geophysical time series to check consistency.
  • Comparison against any newly uncovered astronomical observations from the gap years would directly test the filled values.
  • Extending SSA to the entire C01 series might expose additional low-amplitude signals not captured by standard parametric models.

Load-bearing premise

That the parametric components and SSA decomposition capture the true polar motion behavior inside the gap without systematic bias from unmodeled signals or data selection effects.

What would settle it

Independent historical measurements of pole position from 1858.9-1860.9 that lie outside the uncertainty range of either set of filled values.

Figures

Figures reproduced from arXiv: 2412.07868 by Nina Golyandina, Roman Olenev, Zinovy Malkin.

Figure 1
Figure 1. Figure 1: The Xp component (top panel) and Y p component (bottom panel) of the IERS C01 series. in the framework of several programs for determination of absolute declinations of stars at three observatories: Greenwich, Pulkovo, and Washington. The Greenwich latitude variations were taken from Chandler (1894b) (years 1836.0–1851.0) and Thackeray (1893) (years 1851.0–1891.5). The Pulkovo latitude variations were take… view at source ↗
Figure 2
Figure 2. Figure 2: Generalized Lomb-Scargle periodogram of the IERS C01 PM [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 1
Figure 1. Figure 1: 3 Parametric model It is not possible to use a straightforward model consisting of two components, CW and AW, because their amplitudes do not remain constant over time in the date range around the gap (Ko laczek and Kosek, 1999; Yatskiv, 2000; Guo et al., 2005; Miller, 2011; Lopes et al., 2022; Zotov et al., 2022). However, these two inde￾pendent analyses of the IERS C01 series showed that the CW amplitude… view at source ↗
Figure 3
Figure 3. Figure 3: IERS C01 Yp series with filled missed data between 1858.9 and 1860.9 using the BCA model (upper panel) and BCA2 (bottom panel). Plot titles show date intervals used to fit the Yp data to the model. To evaluate the real accuracy of the PM series approximation by the parametric models, we computed the differences between the modeled and observed PM values for each of seven variants for the epochs before and … view at source ↗
Figure 4
Figure 4. Figure 4: IERS C01 Yp series with filled missed data between 1858.9 and 1860.9 computed using the CSSA model with the seven models providing the minimal approximation error. Comparison of the data presented in [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: The Y components of CSSA solution with r=6 (in arcseconds) and their Lomb-Scargle periodograms (except the trend component, arbitrary units). The period of the most powerful oscillation and the contribution of the component to the total signal are shown on the right. 10 [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: The Y components of CSSA solution with r=9 (in arcseconds) and their Lomb-Scargle periodograms (except the trend component, arbitrary units). The period of the most powerful oscillation and the contribution of the component to the total signal are shown on the right. 11 [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: IERS C01 Yp series with filled missed data between 1858.9 and 1860.9 using the average values for both BCA and SSA models. wavelet technique, etc., have modifications for processing unevenly spaced data. Many examples of application of such analyses to study PM series can be found in the literature, e. g., Gaposchkin (1972); Ko laczek and Kosek (1999); Gambis (2000); Guo et al. (2005); Zotov and Bizouard (… view at source ↗
read the original abstract

The C01 Earth orientation parameters (EOP) series provided by the International Earth Rotation and Reference Systems Service (IERS) is the longest reliable record of the Earth's rotation. In particular, the polar motion (PM) series beginning from 1846 provides a basis for investigation of the long-term PM variations. However, the pole coordinate $Y_p$ in the IERS C01 PM series has a 2-year gap, which makes this series not completely evenly spaced. This paper presents the results of the first attempt to overcome this problem and discusses some ways to fill this gap. Two novel approaches were considered for this purpose: parametric astronomical model consisting of the bias and the Chandler and annual wobbles with linearly changing amplitudes, and data-driven model based on Singular Spectrum Analysis (SSA). Both methods were tested with various options to ensure robust and reliable results. The results obtained by the two methods generally agree within the $Y_p$ errors in the IERS C01 series, but the results obtained by the SSA approach can be considered preferable because it is based on a more complete PM model.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper presents the first attempt to fill a 2-year gap (1858.9-1860.9) in the Y_p component of the IERS C01 polar motion series using two approaches: a parametric model consisting of bias plus Chandler and annual wobbles with linearly changing amplitudes, and a data-driven Singular Spectrum Analysis (SSA) decomposition. Both methods produce filled values that lie within the published C01 uncertainties, with the authors concluding that the SSA results are preferable because they rest on a more complete PM model.

Significance. A reliably filled gap would yield a fully even-spaced long-term PM record useful for studies of secular variations. The significance is reduced by the absence of any quantitative validation, error propagation, or external checks, so it remains unclear whether the filled values improve or degrade the series for downstream analyses.

major comments (3)
  1. [Abstract] Abstract: the claim that SSA results 'can be considered preferable because it is based on a more complete PM model' is not accompanied by any quantitative comparison (RMS difference, cross-validation score, or recovery metric) or falsification test; without such evidence the preference cannot be justified over the deliberately restricted parametric model.
  2. [Methods / Results] No section describes a synthetic-gap recovery experiment, hold-out interval test, or comparison against independent historical proxies for 1858.9-1860.9; such a test is required to show that the additional SSA modes are unbiased physical signals rather than edge or embedding artifacts.
  3. [Abstract / Results] The manuscript states that 'multiple options were tested and results agree within errors' but supplies no numerical validation metrics, error-propagation formulas, or uncertainty estimates for the filled Y_p values themselves.
minor comments (2)
  1. [Methods] Notation for the parametric model components (bias, Chandler amplitude trend, annual amplitude trend) should be defined explicitly with equations rather than described only in prose.
  2. [Methods] The embedding dimension and window length chosen for SSA should be stated together with a brief justification or sensitivity check.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for the constructive comments highlighting the need for stronger quantitative validation. We respond to each major comment below and will revise the manuscript accordingly where feasible.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the claim that SSA results 'can be considered preferable because it is based on a more complete PM model' is not accompanied by any quantitative comparison (RMS difference, cross-validation score, or recovery metric) or falsification test; without such evidence the preference cannot be justified over the deliberately restricted parametric model.

    Authors: We accept the point that the preference requires quantitative support. The SSA approach is based on a data-driven decomposition that captures the full observed PM spectrum (including any additional low-frequency or irregular components) rather than restricting to pre-specified Chandler and annual terms with linear amplitude trends. In revision we will add explicit RMS differences between the two filled series, the fraction of variance explained by extra SSA modes, and a synthetic recovery test to provide the requested metrics. revision: yes

  2. Referee: [Methods / Results] No section describes a synthetic-gap recovery experiment, hold-out interval test, or comparison against independent historical proxies for 1858.9-1860.9; such a test is required to show that the additional SSA modes are unbiased physical signals rather than edge or embedding artifacts.

    Authors: We agree a synthetic-gap recovery test is needed and will add one by creating artificial gaps in adjacent well-sampled intervals, recovering them with both methods, and reporting recovery statistics and mode stability. No independent historical proxies for polar motion exist for the specific 1858.9-1860.9 gap interval; this is the documented observational gap in the primary series. revision: partial

  3. Referee: [Abstract / Results] The manuscript states that 'multiple options were tested and results agree within errors' but supplies no numerical validation metrics, error-propagation formulas, or uncertainty estimates for the filled Y_p values themselves.

    Authors: We agree that numerical metrics and uncertainty estimates should be supplied. The revised manuscript will include tables of RMS agreement values between methods and options, explicit formulas for propagating the original C01 uncertainties into the filled values, and per-epoch uncertainty estimates for the SSA-filled Y_p series. revision: yes

standing simulated objections not resolved
  • Comparison against independent historical proxies for 1858.9-1860.9, as no such proxies are known to exist.

Circularity Check

0 steps flagged

No circularity: gap fill derived from external data only

full rationale

The paper applies both the parametric model (bias plus linear-amplitude Chandler and annual terms) and SSA to the C01 series outside the 1858.9-1860.9 gap, then uses the fitted models to generate values inside the gap. The filled values are therefore model outputs, not inputs that define the models. The stated preference for SSA rests on the modeling assumption of greater completeness rather than any equation that reduces the result to its own fitted parameters or to a self-citation chain. No self-definitional, fitted-input-called-prediction, or ansatz-smuggling steps are present in the described derivation.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The work rests on the domain assumption that polar motion is adequately described by bias plus Chandler and annual terms (parametric) or by the dominant singular vectors extracted from surrounding data (SSA). No new entities are postulated.

free parameters (1)
  • linear amplitude trends
    Amplitudes of Chandler and annual wobbles are allowed to vary linearly and are adjusted to surrounding data.
axioms (2)
  • domain assumption Polar motion consists of bias plus Chandler and annual wobbles whose amplitudes may change linearly
    Invoked to construct the parametric model for gap filling.
  • domain assumption Singular spectrum analysis extracts the dominant components of polar motion from the observed series
    Basis for the data-driven SSA filling method.

pith-pipeline@v0.9.0 · 5742 in / 1160 out tokens · 19433 ms · 2026-05-23T07:29:58.198638+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

51 extracted references · 51 canonical work pages

  1. [1]

    Advances in Space Research 66(11):2487--2515

    Beutler G, Villiger A, Dach R, Verdun A, J \"a ggi A (2020) Long polar motion series: Facts and insights . Advances in Space Research 66(11):2487--2515. https://doi.org/10.1016/j.asr.2020.08.033

  2. [2]

    Journal of Geodesy 93(5):621--633

    Bizouard C, Lambert S, Gattano C, Becker O, Richard JY (2019) The IERS EOP 14C04 solution for Earth orientation parameters consistent with ITRF 2014 . Journal of Geodesy 93(5):621--633. https://doi.org/10.1007/s00190-018-1186-3

  3. [3]

    14:73--75

    Chandler SC (1894 a ) On the inequalities in the coefficients of the law of latitude-variation . 14:73--75. https://doi.org/10.1086/102074

  4. [4]

    14:57--60

    Chandler SC (1894 b ) Variation of latitude from the Greenwich muralcircle Observations, 1836-51 . 14:57--60. https://doi.org/10.1086/102062

  5. [5]

    178: Polar Motion: Historical and Scientific Problems , Astronomical Society of the Pacific Conference Series, vol 208

    Dick S, McCarthy D, Luzum B (eds) (2000) IAU Colloq. 178: Polar Motion: Historical and Scientific Problems , Astronomical Society of the Pacific Conference Series, vol 208. Astronomical Society of the Pacific

  6. [6]

    Springer

    Elsner JB, Tsonis AA (1996) Singular Spectrum Analysis: A New Tool in Time Series Analysis . Springer

  7. [7]

    Kyiv, Naukova Dumka (in Russian)

    Fedorov EP, Korsun' AA, Major SP, Panchenko NI, Taradij VK, Yatskiv YS (1972) Motion of the Earth's pole from 1890.0 to 1969.0 . Kyiv, Naukova Dumka (in Russian)

  8. [8]

    In: Dick S, McCarthy D, Luzum B (eds) IAU Colloq

    Gambis D (2000) Long-term Earth Orientation Monitoring Using Various Techniques . In: Dick S, McCarthy D, Luzum B (eds) IAU Colloq. 178: Polar Motion: Historical and Scientific Problems, Astronomical Society of the Pacific, Astronomical Society of the Pacific Conference Series, vol 208, pp 337--344

  9. [9]

    In: Melchior P, Yumi S (eds) IAU Symp

    Gaposchkin EM (1972) Analysis of Pole Position from 1846 to 1970 . In: Melchior P, Yumi S (eds) IAU Symp. 48: Rotation of the Earth, Cambridge University Press, vol 48, p 19–32, https://doi.org/10.1017/S0074180900098016

  10. [10]

    Reviews of Geophysics 40(1):1003

    Ghil M, Allen MR, Dettinger MD, Ide K, Kondrashov D, Mann ME, Robertson AW, Saunders A, Tian Y, Varadi F, Yiou P (2002) Advanced Spectral Methods for Climatic Time Series . Reviews of Geophysics 40(1):1003. https://doi.org/10.1029/2000RG000092

  11. [11]

    J Stat Plan Inference 137(8):2642--2653

    Golyandina N, Osipov E (2007) The `` Caterpillar ''-- SSA method for analysis of time series with missing values. J Stat Plan Inference 137(8):2642--2653. https://doi.org/10.1016/j.jspi.2006.05.014

  12. [12]

    Second Edition

    Golyandina N, Zhigljavsky A (2020) Singular Spectrum Analysis for Time Series. Second Edition . Springer, https://doi.org/10.1007/978-3-662-62436-4

  13. [13]

    Chapman and Hall/CRC

    Golyandina N, Nekrutkin V, Zhigljavsky A (2001) Analysis of Time Series Structure: SSA and related techniques . Chapman and Hall/CRC

  14. [14]

    J Stat Softw 67(2):1--78

    Golyandina N, Korobeynikov A, Shlemov A, Usevich K (2015) Multivariate and 2D Extensions of Singular Spectrum Analysis with the Rssa Package. J Stat Softw 67(2):1--78. https://doi.org/10.18637/jss.v067.i02

  15. [15]

    Geophysical Journal International 109:162--170

    Gross RS (1992) Correspondence between theory and observations of polar motion. Geophysical Journal International 109:162--170. https://doi.org/10.1111/j.1365-246X.1992.tb00086.x

  16. [16]

    Journal of Geodesy 78(11-12):654--659

    Guo JY, Greiner-Mai H, Ballani L, Jochmann H, Shum CK (2005) On the double-peak spectrum of the Chandler wobble . Journal of Geodesy 78(11-12):654--659. https://doi.org/10.1007/s00190-004-0431-0

  17. [17]

    Surveys in Geophysics 25(1):1--54

    H \"o pfner J (2004) Low-Frequency Variations, Chandler and Annual Wobbles of Polar Motion as Observed Over One Century . Surveys in Geophysics 25(1):1--54. https://doi.org/10.1023/B:GEOP.0000015345.88410.36

  18. [18]

    Journal of Geodesy 97(8):74

    Ji K, Shen Y, Chen Q, Wang F (2023) Extended singular spectrum analysis for processing incomplete heterogeneous geodetic time series . Journal of Geodesy 97(8):74. https://doi.org/10.1007/s00190-023-01764-8

  19. [19]

    Experimental Long-Lead Bulletin 5(3):39--40

    Keppenne C, Lall U (1996) Complex singular spectrum analysis and multivariate adaptive regression splines applied to forecasting the southern oscillation. Experimental Long-Lead Bulletin 5(3):39--40

  20. [20]

    In: Journ\'ees 1998 ``Syst\`emes de R\'ef\'erence Spatio-Temporels: Conceptual, Conventional and Practical Studies Related to Earth Rotation'', pp 215--220

    Ko aczek B, Kosek W (1999) Variations of the amplitude of the Chandler wobble. In: Journ\'ees 1998 ``Syst\`emes de R\'ef\'erence Spatio-Temporels: Conceptual, Conventional and Practical Studies Related to Earth Rotation'', pp 215--220

  21. [21]

    Nonlinear Processes in Geophysics 13(2):151--159

    Kondrashov D, Ghil M (2006) Spatio-temporal filling of missing points in geophysical data sets . Nonlinear Processes in Geophysics 13(2):151--159. https://doi.org/10.5194/npg-13-151-2006

  22. [22]

    IEEE Trans Acoust 30(6):833--840

    Kumaresan R, Tufts D (1982) Estimating the parameters of exponentially damped sinusoids and pole-zero modeling in noise. IEEE Trans Acoust 30(6):833--840

  23. [23]

    Acta Astronomica Sinica 15:86--92 (in Chinese)

    Li Z, Xu H, Lin B (1974) Coordinates of the earth's instataneous pole from 1825.0 to 1897.9. Acta Astronomica Sinica 15:86--92 (in Chinese)

  24. [24]

    Geosciences 12(12):448

    Lopes F, Courtillot V, Gibert D, Le Mou \"e l JL (2022) Extending the Range of Milankovic Cycles and Resulting Global Temperature Variations to Shorter Periods (1 100 Year Range) . Geosciences 12(12):448. https://doi.org/10.3390/geosciences12120448

  25. [25]

    Earth, Planets and Space 62(12):943--947

    Malkin Z, Miller N (2010) Chandler wobble: two more large phase jumps revealed . Earth, Planets and Space 62(12):943--947. https://doi.org/10.5047/eps.2010.11.002

  26. [26]

    Solar System Research 45(4):342--353

    Miller NO (2011) Chandler wobble in variations of the Pulkovo latitude for 170 years . Solar System Research 45(4):342--353. https://doi.org/10.1134/S0038094611040058

  27. [27]

    Okhotnikov G, Golyandina N (2019) EOP time series prediction using singular spectrum analysis. In: Corpetti T, Ienco D, Interdonato R, et al (eds) Proceedings of MACLEAN : MAChine Learning for EArth ObservatioN Workshop , RWTH Aahen University, Germany, CEUR Workshop Proceedings, Vol. 2466

  28. [28]

    In: Aksent'eva ZN (ed) Orlov A

    Orlov AY (1961) Free nutation based on observations in Pulkovo in 1842--1912 . In: Aksent'eva ZN (ed) Orlov A. Ya., Selected works, Vol. 1 , Kyiv: Academy of Sciences of USSR, pp 95--113 (in Russian)

  29. [29]

    ur Kartographie und Geod\

    Petit G, Luzum B (eds) (2010) IERS Conventions (2010). IERS Technical Note No. 36, Verlag des Bundesamts f\"ur Kartographie und Geod\"asie, Frankfurt am Main

  30. [30]

    Soobshcheniya Gosudarstvennogo Astronomicheskogo Instituta, Moscow State University 163:3--10 (in Russian)

    Rykhlova LV (1970 a ) The coordinates of the Earth's pole for the years 1846.0-1891.5. Soobshcheniya Gosudarstvennogo Astronomicheskogo Instituta, Moscow State University 163:3--10 (in Russian)

  31. [31]

    Trudy Gosudarstvennogo Astronomicheskogo Instituta, Moscow State University 39:200--212 (in Russian)

    Rykhlova LV (1970 b ) Variations of the Washington latitude in 1845--1891 . Trudy Gosudarstvennogo Astronomicheskogo Instituta, Moscow State University 39:200--212 (in Russian)

  32. [32]

    In: Lieske JH, Abalakin VK (eds) Inertial Coordinate System on the Sky, vol 141, p 157

    Rykhlova LV, Kurbasova GS (1990) The Study of the Structure of 142-Year Series of Pole Coordinates . In: Lieske JH, Abalakin VK (eds) Inertial Coordinate System on the Sky, vol 141, p 157

  33. [33]

    Journal of Geodesy 74(10):701--710

    Schuh H, Nagel S, Seitz T (2001) Linear drift and periodic variations observed in long time series of polar motion . Journal of Geodesy 74(10):701--710. https://doi.org/10.1007/s001900000133

  34. [34]

    Journal of Geophysical Research (Solid Earth) 117(B9):B09403

    Seitz F, Kirschner S, Neubersch D (2012) Determination of the Earth's pole tide Love number k _ 2 from observations of polar motion using an adaptive Kalman filter approach . Journal of Geophysical Research (Solid Earth) 117(B9):B09403. https://doi.org/10.1029/2012JB009296

  35. [35]

    Journal of the Geodetic Society of Japan 21:131--141

    Sekiguchi N (1975) On the latitude variations of the interval between 1830 and 1860. Journal of the Geodetic Society of Japan 21:131--141

  36. [36]

    Nonlinear Processes in Geophysics 22(4):371--376

    Shen Y, Peng F, Li B (2015) Improved singular spectrum analysis for time series with missing data . Nonlinear Processes in Geophysics 22(4):371--376. https://doi.org/10.5194/npg-22-371-201510.5194/npgd-1-1947-2014

  37. [37]

    Thackeray WG (1893) Latitude variation and Greenwich observations, 1851-91 . 53:120. https://doi.org/10.1093/mnras/53.3.120

  38. [38]

    Physica D Nonlinear Phenomena 35:395--424

    Vautard R, Ghil M (1989) Singular spectrum analysis in nonlinear dynamics, with applications to paleoclimatic time series . Physica D Nonlinear Phenomena 35:395--424. https://doi.org/10.1016/0167-2789(89)90077-8

  39. [39]

    Vityazev VV, Miller NO, Prudnikova EJ (2010) Singular Spectrum Analysis in Astrometry and Geodynamics . In: de Le \'o n M, de Diego DM, Ros RM (eds) Mathematics and Astronomy: A Joint Long Journey, AIP, American Institute of Physics Conference Series, vol 1283, pp 319--328, https://doi.org/10.1063/1.3506397

  40. [40]

    Vondr \'a k J (1988) Is Chandler Frequency Constant? In: Babcock AK, Wilkins GA (eds) The Earth's Rotation and Reference Frames for Geodesy and Geodynamics, IAU Symposium, vol 128, p 359

  41. [41]

    297:899--906

    Vondr\'ak J, Ron C, Pe s ek I, C epek A (1995) New global solution of Earth orientation parameters from optical astrometry in 1900-1990 . 297:899--906

  42. [42]

    Acta Geodyn Geomater 7(3):245–251

    Vondr\'ak J, Ron C, S tefka A (2010) Earth orientation parameters based on EOC-4 astrometric catalog . Acta Geodyn Geomater 7(3):245–251

  43. [43]

    https://doi.org/10.1186/s40623-023-01944-y

    Yamaguchi R, Furuya M (2024) Can we explain the post-2015 absence of the Chandler wobble? Earth, Planets and Space 76(1):1. https://doi.org/10.1186/s40623-023-01944-y

  44. [44]

    In: Dick S, McCarthy D, Luzum B (eds) IAU Colloq

    Yatskiv Y (2000) Chandler Motion Observations . In: Dick S, McCarthy D, Luzum B (eds) IAU Colloq. 178: Polar Motion: Historical and Scientific Problems, Astronomical Society of the Pacific, Astronomical Society of the Pacific Conference Series, vol 208, pp 383--395

  45. [45]

    Journal of Geophysical Research (Solid Earth) 126(5):e2020JB021227

    Yi S, Sneeuw N (2021) Filling the Data Gaps Within GRACE Missions Using Singular Spectrum Analysis . Journal of Geophysical Research (Solid Earth) 126(5):e2020JB021227. https://doi.org/10.1029/2020JB021227

  46. [46]

    A new formalism for the floating-mean and Keplerian periodograms

    Zechmeister M, K \"u rster M (2009) The generalised Lomb-Scargle periodogram. A new formalism for the floating-mean and Keplerian periodograms . 496(2):577--584. https://doi.org/10.1051/0004-6361:200811296

  47. [47]

    Journal of Geodynamics 62:30--34

    Zotov L, Bizouard C (2012) On modulations of the Chandler wobble excitation . Journal of Geodynamics 62:30--34. https://doi.org/10.1016/j.jog.2012.03.010

  48. [48]

    Artificial Satellites 45(2):95--106

    Zotov LV (2010) Dynamical Modeling and Excitation Reconstruction as Fundamental of Earth Rotation Prediction . Artificial Satellites 45(2):95--106. https://doi.org/10.2478/v10018-010-0010-y

  49. [49]

    Moscow University Physics Bulletin 77(3):555--563

    Zotov LV, Sidorenkov NS, Bizouard C (2022) Anomalies of the Chandler Wobble in 2010s . Moscow University Physics Bulletin 77(3):555--563. https://doi.org/10.3103/S0027134922030134

  50. [50]

    , " * write output.state after.block = add.period write newline

    ENTRY address archive author booktitle chapter doi edition editor eid eprint howpublished institution journal key month note number organization pages publisher school series title type url volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all ...

  51. [51]

    write newline

    " write newline "" before.all 'output.state := FUNCTION add.period duplicate empty 'skip "." * add.blank if FUNCTION if.digit duplicate "0" = swap duplicate "1" = swap duplicate "2" = swap duplicate "3" = swap duplicate "4" = swap duplicate "5" = swap duplicate "6" = swap duplicate "7" = swap duplicate "8" = swap "9" = or or or or or or or or or FUNCTION ...