Our analysis
A plain sine wave does nearly as well
A 70-year sine with only an amplitude and a phase matches the paper’s illustrative Figure 1 curve (0.21 against 0.20 ms). With its own two settings tuned, the gravitational model gets to about 0.10–0.13 ms, but a sine with a free period, or even a cubic curve, does as well.
The reason: after 30-year smoothing, 56 years of data hold only about four independent values. The model’s strength and relaxation time act as an amplitude and a phase knob, though the seismic history limits how far the phase can move.
Our analysis
What still counts
The gravitational curve’s shape comes from independent seismic data, its fitted values are physically sensible, and the other two torques point the wrong way in all 400 core-flow versions, even in a test that ignores long-term averages. For them, any strength at all makes the fit worse than predicting no change.
Much of the record is a net shortening of the day by about 2.5 ms since 1964, which the gravitational prediction reproduces. The Methods say a temporal mean is subtracted from all three predictions [1 Methods]. The released code removes the 1964–2019 average from the boundary torques, as the supplement says, but measures gravity from its long-term equilibrium, so only gravity can produce the net shortening [8 p. 2][1 Methods][2]. That choice matters for the fit: with gravity’s 1964–2019 average removed too, the Figure 1 fit falls to r = 0.81, and in our tests no retuned setting, even with a resisting torque, gets under the 0.2 ms cut-off (the best are about 0.3–0.45 ms). It doesn’t decide the direction: compared directly with the torque the record needs, gravity still scores +0.89 and the boundary torques −0.67 and −0.68, and measured from their 1880–2019 average the boundary torques run even more opposite (r ≈ −0.97).