Accurately predicting Earth’s rotation rate, as represented by Length of Day (LOD) variations, is essential for applications such as satellite navigation, climate studies, geophysical research, and disaster prevention. However, predicting LOD is challenging due to its sensitivity to various geophysical and meteorological factors. Current methods, including statistical approaches, often struggle with short-term forecasting accuracy. In this study, we use Monte Carlo Singular Spectrum Analysis (MCSSA) to distinguish between deterministic and non-deterministic components within the LOD time series. The deterministic components are extended using the SSA prediction algorithm. To enhance robustness, we refine Allen and Smith’s methodology (testing significance of eigenmodes against an autoregressive (AR) (1) noise null hypothesis) by integrating an autoregressive moving average (ARMA) model to account for noise, providing valuable insights into the non-deterministic behaviors present in the series. We comprehensively evaluate our methodology through a comparative analysis. For long-term prediction (365 days), we compare our method against the combined LS and autoregressive (AR) method. For short-term prediction (next 10 days), we compare it against the results of the second Earth Orientation Parameters Prediction Comparison Campaign (second EOP PCC). Using the IERS 20 C04 time series, our hybrid model demonstrates a superior long-term prediction accuracy with a mean absolute error (MAE) of 0.201 ms/day on the 365th day. Additionally, the short-term prediction performance is comparable to the second EOP PCC results. These results illustrate that the proposed method efficiently predicts LOD, showing significant improvement in long-term accuracy and robustness in short-term forecasting.

Earth orientation parameters (EOP) are critical for applications in orbit determination, astronomy, space geodesy, and geophysics. Accurate predictions of EOP rely on the identification of both deterministic periodic patterns and noise characteristics. This study addresses these requirements by analyzing the Polar Motion (PM) and Length Of Day (LOD) time series to determine its stochastic model structure, estimated using least squares variance component estimation (LS-VCE). With this model, deterministic periodic patterns were extracted through least squares harmonic estimation (LS-HE) and validated against colored noise components to identify significant signals. Using the IERS 14 C04 data from January 1, 2000, to December 31, 2019, the study identified the noise as power-law with a spectral index of −1.5, suggesting non-stationary fractional Brownian characteristics. LS-HE detected dominant frequencies in PM–notably the Chandler and annual signals–and in LOD, with annual, semiannual, 14-day, and 9-day signals. Building on these findings, short-, mid-, and long-term prediction model were developed for the PM and LOD time series from September 2021 to December 2022. The predictive model combines the LS-HE-extracted signals and the noise model to generate forecasts. These predictions were compared with other models from the Second Earth Orientation Parameter Prediction Comparison Campaign, demonstrating competitive accuracy, particularly for the initial forecast days. The results validate that combining LS-HE with a realistic noise model provides an effective approach for short-term of the PM and LOD forecasting, meeting the accuracy goals of geodetic and geophysical applications.