100 random realisations (blue) and medians for each epoch of the time-domain reconstructions of the DM variations delays at 1.4 GHz radio frequency modelled as a Gaussian process for PSRs J0613–0200 (left) and J1909–3744 (right), with the “DR2full+”. The black dashed lines display the last epochs of DR2full (that is, EPTA only). The inclusion of InPTA data allows to measure sharp changes at after MJD 59000. The delays are obtained from marginalising the timing model parameters which comprise the DM constant and the first two time derivatives of the dispersion measure (DM1, DM2).