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
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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.
- [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)
- [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.
- [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
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
-
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
-
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
-
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
- Comparison against independent historical proxies for 1858.9-1860.9, as no such proxies are known to exist.
Circularity Check
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
free parameters (1)
- linear amplitude trends
axioms (2)
- domain assumption Polar motion consists of bias plus Chandler and annual wobbles whose amplitudes may change linearly
- domain assumption Singular spectrum analysis extracts the dominant components of polar motion from the observed series
Reference graph
Works this paper leans on
-
[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]
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]
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]
Chandler SC (1894 b ) Variation of latitude from the Greenwich muralcircle Observations, 1836-51 . 14:57--60. https://doi.org/10.1086/102062
-
[5]
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
work page 2000
- [6]
-
[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)
work page 1972
-
[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
work page 2000
-
[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]
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]
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]
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]
Golyandina N, Nekrutkin V, Zhigljavsky A (2001) Analysis of Time Series Structure: SSA and related techniques . Chapman and Hall/CRC
work page 2001
-
[14]
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]
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]
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]
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]
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]
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
work page 1996
-
[20]
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
work page 1999
-
[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]
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
work page 1982
-
[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)
work page 1974
-
[24]
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]
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]
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]
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
work page 2019
-
[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)
work page 1961
-
[29]
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
work page 2010
-
[30]
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)
work page 1970
-
[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)
work page 1970
-
[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
work page 1990
-
[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]
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]
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
work page 1975
-
[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
work page doi:10.5194/npg-22-371-201510.5194/npgd-1-1947-2014 2015
-
[37]
Thackeray WG (1893) Latitude variation and Greenwich observations, 1851-91 . 53:120. https://doi.org/10.1093/mnras/53.3.120
-
[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]
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]
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
work page 1988
-
[41]
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
work page 1995
-
[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
work page 2010
-
[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]
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
work page 2000
-
[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]
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]
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]
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]
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]
, " * 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]
" 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 ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.