JBE phase residuals R(t) with respect to local fits ϕ_j(t) calculated for the 10 chosen episodes (see Fig. 3). Gray vertical lines are plotted at the dates of reported glitches (Espinoza et al. 2011). The braking index during episodes is shown in orange. Sinusoidal fits to residuals, discussed below, are shown in blue and cyan to distinguish episodes. Gray points at the bottom show ten times the difference between R(t) and sinusoidal fits R_j(t) at data points; the horizontal dashed lines bracket these final residuals with their standard deviation of 0.057 turns; the distribution of final residuals has wider wings than a Gaussian.