Glitch Size Research Articles

A multiband study of pulsar glitches with Fermi-LAT and Parkes

ABSTRACT Pulsar glitch is a phenomenon characterized by abrupt changes in the spin period over less than a minute. We present a comprehensive analysis of glitches in four gamma-ray pulsars by combining the timing observation data from Fermi-Large Area Telescope (Fermi-LAT) and Parkes 64 m radio telescope. The timing data of five pulsars, namely PSRs J1028–5819, J1420–6048, J1509–5850, J1709–4429 (B1706–44), and J1718–3825, are examined over 14 yr of observations for each. A total of 12 glitches are identified in four pulsars, including a previously unreported glitch. That is, a new small glitch is identified for PSR J1718–3825 in MJD $\sim$ 59121(8), with a fractional glitch size of $\Delta \nu /\nu \sim 1.9(2) \times 10^{-9}$. For PSR J1420–6048, our investigation confirms the presence of two linear recovery terms during the evolution of $\dot{\nu }$ following glitches 4, 6, and 8. Moreover, an exponential recovery process was identified after glitch 8, with a recovery fraction (Q) of $Q = 0.0131(5)$ and a corresponding time-scale of $\tau _{\rm d} = 100(6)$ d. Regarding the fourth glitch of PSR J1709–4429, our analysis reveals the presence of two exponential recovery terms with degree of recovery and decay time-scales Q1 = 0.0104(5), $\tau _{\rm d1}=72(4)$ d, and Q2 = 0.006(1), $\tau _{\rm d2}=4.2(6)$ d, respectively. For the remaining previously reported glitches, we also refine the glitch epochs and recovery process through precise fitting of the timing data. We discuss how multiband data of glitches can help better characterize the glitch recoveries and constrain the underlying physics of glitch events. Our findings demonstrate that the accumulation of observational data reveals the rich complexity of the glitch phenomenon, aiding in the search for a well-established interpretation.

Glitches and Glitching Clusters in Rotation-powered Pulsars

Abstract The study of pulsar glitch phenomena serves as a valuable probe into the dynamic properties of matter under extreme high-density conditions, offering insights into the physics within neutron stars. Providing theoretical explanations for the diverse manifestations observed in different pulsars has proven to be a formidable challenge. By analyzing the distribution of glitch sizes and waiting times, along with the evolution of cumulative glitch sizes over time, we have uncovered a long-term clustering phenomenon for pulsar glitches. This perspective allows us to approach the distinct glitch representations in various pulsars from a unified standpoint, connecting the same periodicity of observational data to the randomness. Without relying on specific physical models, we utilized the coefficient of variation to numerically determine optimal clustering numbers and clustering periods for sample pulsars. Our analysis involving 27 pulsars has revealed a clear linear relationship between the glitch cluster period and characteristic age. Of interest, the cumulative distribution of functions of cluster sizes and interval times have the same patterns, which can be synchronously fitted by Gaussian processes. These results may indicate novel understandings of glitches and the resulting processes.

Timing irregularities and glitches from the pulsar monitoring campaign at IAR

Pulsars have a very stable rotation overall. However, sudden increases in their rotation frequency, known as glitches, perturb their evolution. While many observatories commonly detect large glitches, small glitches are harder to detect because of the lack of daily cadence observations over long periods of time (years). We aim to explore and characterise the timing behaviour of young pulsars on daily timescales, looking for small glitches and other irregularities, in order to further our comprehension of the real distribution of glitch sizes. Our findings have consequences for the theoretical modelling of the glitch mechanism. We observed six pulsars with up to daily cadence between December 2019 and January 2024 with the two antennas of the Argentine Institute of Radio Astronomy (IAR). We used standard pulsar timing tools to obtain the times of arrival of the pulses and to characterise the pulsar's rotation. We developed an algorithm to look for small timing events in the data and calculate the changes in the frequency (nu ) and its derivative ($ at those epochs. We find that the rotation of all pulsars in this dataset is affected by small step changes in nu and $ Among them, we find three new glitches that have not been reported before: two glitches in PSR J1048$-$5832 with relative sizes of $ $ and $ $, and one glitch in the Vela pulsar with a size of $ $. We also report new decay terms on the 2021 Vela giant glitch, and on the 2022 giant glitches in PSR J0742$-$2822 and PSR J1740$-$3015, respectively. In addition, we find that the red noise contribution significantly diminished in PSR J0742$-$2822 after its giant glitch in 2022. Our results highlight the importance of high-cadence monitoring with an exhaustive analysis of the residuals to better characterise the distribution of glitch sizes and to deepen our understanding of the mechanisms behind glitches, red noise, and timing irregularities.

Measuring glitch recoveries and braking indices with Bayesian model selection

ABSTRACT For a selection of 35 pulsars with large spin-up glitches ($\Delta {\nu }/\nu \ge 10^{-6}$), which are monitored by the Jodrell Bank Observatory, we analyse 157 glitches and their recoveries. All parameters are measured consistently and we choose the best model to describe the post-glitch recovery based on Bayesian evidence. We present updated glitch epochs, sizes, changes of spin down rate, exponentially recovering components (amplitude and corresponding time-scale) when present, as well as pulsars’ second frequency derivatives and their glitch-associated changes if detected. We discuss the different observed styles of post-glitch recovery as well as some particularly interesting sources. Several correlations are revealed between glitch parameters and pulsar spin parameters, including a very strong correlation between a pulsar’s interglitch $|\ddot{\nu }|$ and $\dot{\nu }$, as well as between the glitch-induced spin-down rate change $\Delta \dot{\nu }_{\rm p}$ that does not relax exponentially and $\dot{\nu }$. We find that the ratio $\left|\Delta \dot{\nu }_{\mathrm{p}}/\ddot{\nu }\right|$ can be used as an estimate of glitch recurrence times, especially for those pulsars for which there are indications of a characteristic glitch size and interglitch waiting time. We calculate the interglitch braking index n and find that pulsars with large glitches typically have n greater than 3, suggesting that internal torques dominate the rotational evolution between glitches. The external torque, for example, from electromagnetic dipole radiation, could dominate the observed $\ddot{\nu }$ for the youngest pulsars ($\lesssim 10^{4}\,\,\mathrm{yr}$), which may be expected to display $n\sim 3$.

Pulsar glitch in a strangeon star model – III. The recovery

ABSTRACT Strangeon star model has passed various observational tests, such as the massive pulsars and the tidal deformability during binary mergers. Pulsar glitch, as a useful probe for studying the interior structure of pulsars, has also been studied in strangeon star model in our previous papers, including the recovery coefficient, the waiting time of glitches, and glitch activity. In this paper, the recovery process of a glitch is described in the strangeon star model, based on the starquake picture established before (in Paper I). After the starquake, the inner motion of the stellar matter would reduce the tangential pressure in the cracked places at the equatorial plane. The recovery (increase) of the tangential pressure would be achieved by a viscous flow towards the cracked places at equatorial plane, which leads to the exponential recovery of the spin frequency. A uniform viscous flow can reproduce the single exponential decay observed in some glitches, and the viscous time-scale τ and the depth h of the cracking place below the surface can be fitted by the recovery data. It is found that h increases with glitch size Δν/ν, which is expected in the glitch scenario of strangeon stars. The magnitude of the recovery predicted in this recovery model is also consistent with that derived from observations. The single exponential decay reproduced by a uniform viscous flow can be generalized to two or more exponentials by the multicomponent of viscous flows.

A new small glitch in Vela discovered with a hidden Markov model

ABSTRACT A striking feature of the Vela pulsar (PSR J0835−4510) is that it undergoes sudden increases in its spin frequency, known as glitches, with a fractional amplitude of the order of 10−6 approximately every 900 d. Glitches of smaller magnitudes are also known to occur in Vela. Their distribution in both time and amplitude is less well constrained but equally important for understanding the physical process underpinning these events. In order to better understand these small glitches in Vela, an analysis of high-cadence observations from the Mount Pleasant Observatory is presented. A hidden Markov model (HMM) is used to search for small, previously undetected glitches across 24 yr of observations covering MJD 44929 to MJD 53647. One previously unknown glitch is detected around MJD 48636 (1992 January 15), with fractional frequency jump Δf/f = (8.19 ± 0.04) × 10−10 and frequency derivative jump $\Delta \dot{f}/\dot{f} = (2.98 \pm 0.01) \times 10^{-4}$ . Two previously reported small glitches are also confidently redetected, and independent estimates of their parameters are reported. Excluding these events, 90 per cent confidence frequentist upper limits on the sizes of missed glitches are also set, with a median upper limit of $\Delta f^{90~{{\% }}}/f = 1.35 \times 10^{-9}$. Upper limits of this kind are enabled by the semi-automated and computationally efficient nature of the HMM, and are crucial to informing studies that are sensitive to the lower end of the glitch size distribution.

Stability of interlinked neutron vortex and proton flux-tube arrays in a neutron star – III. Proton feedback

ABSTRACT The coupled, time-dependent Gross–Pitaevskii and Ginzburg–Landau equations are solved simultaneously in three dimensions to investigate the equilibrium state and far-from-equilibrium, spin-down dynamics of an interpenetrating neutron superfluid and proton type-II superconductor, as an idealized description of the outer core of a neutron star. The simulations generalize previous calculations without the time-dependent Ginzburg–Landau equation, where proton feedback is absent. If the angle θ between the rotation and magnetic axes does not equal zero, the equilibrium state consists of geometrically complicated neutron vortex and proton flux-tube tangles, as the topological defects pin to one another locally but align with different axes globally. During spin-down, new types of motion are observed. For θ = 0, entire vortices pair rectilinearly with flux tubes and move together while pinned. For θ ≠ 0, vortex segments pair with segments from one or more flux tubes, and the paired segments move together while pinned. The degree to which proton feedback impedes the deceleration of the crust is evaluated as a function of θ and the pinning strength, η. Key geometric properties of vortex-flux-tube tangles, such as filament length, mean curvature, and polarity are analysed. It is found that proton feedback smooths the deceleration of the crust, reduces the rotational glitch sizes, and stabilizes the vortex tangle dynamics. The dimensionless control parameters in the simulations are mutually ordered to match what is expected in a real neutron star, but their central values and dynamic ranges differ from reality by many orders of magnitude due to computational limitations.

Glitches due to quasineutron-vortex scattering in the superfluid inner crust of a pulsar

We revisit the mechanism of vortex unpinning caused by the neutron-vortex scattering \cite{prad1} in the inner crust of a pulsar. The strain energy released by the crustquake is assumed to be absorbed in some part of the inner crust and causes pair-breaking quasi-neutron excitations from the existing free neutron superfluid in the bulk of the inner crust. The scattering of these quasi-neutrons with the vortex core normal neutrons unpins a large number of vortices from the thermally affected regions and results in pulsar glitches. We consider the geometry of a cylindrical shell of the affected pinning region to study the implications of the vortex unpinning in the context of pulsar glitches. We find that a pulsar can release about $\sim 10^{11} - 10^{13}$ vortices by this mechanism. These numbers are equivalent to the glitch size of orders $\sim 10^{-11} - 10^{-9}$ for Vela-like pulsars with the characteristic age $\tau \simeq 10^4$ years. We also suggest a possibility of a vortex avalanche triggered by the movement of the unpinned vortices. A rough estimate of the glitch size caused by an avalanche shows an encouraging result.

Prospects for detecting transient quasi-monochromatic gravitational waves from glitching pulsars with current and future detectors

ABSTRACT Pulsars are rotating neutron stars that emit periodic electromagnetic radiation. While pulsars generally slow down as they lose energy, some also experience glitches: spontaneous increases of their rotational frequency. According to several models, these glitches can also lead to the emission of long-duration transient gravitational waves (GWs). We present detection prospects for such signals by comparing indirect energy upper limits on GW strain for known glitches with the sensitivity of current and future ground-based GW detectors. We first consider the optimistic case of generic constraints based on the glitch size and find that realistic matched-filter searches in the fourth LIGO–Virgo–KAGRA observing run (O4) could make a detection, or set constraints below these indirect upper limits, for equivalents of 36 out of 726 previously observed glitches, and 74 in the O5 run. With the third-generation Einstein Telescope or Cosmic Explorer, 35–40 per cent of glitches would be accessible. When specializing to a scenario where transient mountains produce the post-glitch GW emission, following the Yim & Jones model, the indirect upper limits are stricter. Out of the smaller set of 119 glitches with measured healing parameter, as needed for predictions under that model, only 6 glitches would have been within reach for O4 and 14 for O5, with a similar percentage as before with third-generation detectors. We also discuss how this model matches the observed glitch population.

Prospects for detecting and localizing short-duration transient gravitational waves from glitching neutron stars without electromagnetic counterparts

Neutron stars are known to show accelerated spin-up of their rotational frequency called a glitch. Highly magnetized rotating neutron stars (pulsars) are frequently observed by radio telescopes (and in other frequencies), where the glitch is observed as irregular arrival times of pulses which are otherwise very regular. A glitch in an isolated neutron star can excite the fundamental (f)-mode oscillations which can lead to gravitational wave generation. Electromagnetic observations of pulsars (and hence pulsar glitches) require the pulsar to be oriented so that the jet is pointed toward the detector, but this is not a requirement for gravitational wave emission which is more isotropic and not jetlike. Hence, gravitational wave observations have the potential to uncover nearby neutron stars where the jet is not pointed towards the Earth. In this work, we study the prospects of finding glitching neutron stars using a generic all-sky search for short-duration gravitational wave transients. The analysis covers the high-frequency range from $1-4$ kHz of LIGO-Virgo detectors for signals up to a few seconds. We set upper limits for the third observing run of the LIGO-Virgo detectors and present the prospects for upcoming observing runs of LIGO, Virgo, KAGRA, and LIGO India. We find the detectable glitch size will be around $10^{-5}$ Hz for the fifth observing run for pulsars with spin frequencies and distances comparable to the Vela pulsar. We also present the prospects of localizing the direction in the sky of these sources with gravitational waves alone, which can facilitate electromagnetic follow-up. We find that for the five detector configuration, the localization capability for a glitch size of $10^{-5}$ Hz is around $132\,\mathrm{deg}^{2}$ at $1\sigma$ confidence for $50\%$ of events with distance and spin frequency as that of Vela.

Detection of 16 Small Glitches in Nine Pulsars

Pulsar timing measurements with a 26 m radio telescope at Nanshan between 2000 and 2014 were used to search for glitch events. The data span of nine pulsars ranges from 11.6 to 14.2 yr, and 16 new glitch events were identified in nine pulsars. Glitch parameters were determined through fitting the timing residuals data. All 16 glitches have a small fractional size. Six new glitches have been detected in PSR J1833−0827, making it another frequent glitching pulsar. Some of the 16 glitches may experience exponential or linear recovery, but it is unlikely for us to make further analyses with the large gap in the data set. All the glitch rates obtained from Nanshan are higher than that from Jodrell Bank Observatory. The small glitch size and high glitch rate could possibly attribute to the high observation cadence.

Pulse Profile Variations Associated with Two Glitches of PSR B1822–09

We reported two new glitches of PSR B1822−09 that were detected at the Shanghai Tian Ma Radio Telescope. The glitches occurred around MJD 58,025 and 58,474.5, respectively. The shapes of the integrated mean pulse profiles in both the radio-bright (B-mode) and the radio-quiet (Q-mode) modes changed after two glitches. Such changes are probably related to the trigger mechanisms of the glitches. According to the Gaussian fitting to the integrated mean pulse profiles, variations of the integrated mean pulse profiles can be attributed to variations of the Gaussian components. The fitting also indicates that there may be changes in only a part of the Gaussian components after glitches. We proposed an interpretation of the relation between the precursor and the interpulse of PSR B1822−09 according to variations of Gaussian components. We analyzed the cumulative distributions of the glitch sizes and the waiting times of all 14 glitches of PSR B1822−09. The cumulative distribution of the glitch sizes can be fitted well by a power law with an index α of 0.985 ± 0.005. The cumulative distribution of the waiting times follows a Poisson model with a mean waiting time λ of 466 ± 66 days. The correlation coefficient between the waiting times and the sizes of the corresponding preceding glitches is 0.906. In contrast, there is no apparent correlation between the waiting times and the sizes of the corresponding trailing glitches.

On the fractional glitch recoveries and pulsar spin properties

AbstractPulsar glitches are sudden jumps in the spin frequency of an otherwise steadily spinning down neutron star. A glitch, in many pulsars, is followed by a recovery phase in which the pulsars return fully or partially to the pre‐glitch state on a wide range of time scales. This paper statistically investigates the behavior of the glitch recovery parameter, Q, and its relationship with the pulsar rotational parameters. The Q for the different sub‐classes of the glitch sizes appear to arise from different distributions, suggesting a possible different path of recovery or interior dynamics for the recovery of glitches of different sizes. While the Q has moderate relation with characteristic age, rotational frequency, fractional moment of inertia, and coupling parameter, it shows a very weak correlation with spin‐down rate. The implications are discussed. Glitch and glitch recovery appear to have different underlying physics, but both point to some dynamical changes within the interior of the neutron star.

On the pulsar spin frequency derivatives and the glitch activity

The number of sudden spin-ups in radio pulsars known as pulsar glitches has increased over the years. Though a consensus has not been reached with regards to the actual cause of the phenomenon, the electromagnetic braking torque on the crust quantified via the magnitude of pulsar spin frequency first derivative, $ \dot{\nu} $ is a key factor in mechanisms put across toward the understanding of the underlying principles involved. The glitch size has been used to establish a quantity used to constrain the mean possible change in pulsar spin frequency $ (\nu) $ per year due to a glitch known as the `glitch activity'. Traditionally, the glitch activity parameter $ A_{g} $ is calculated from the cumulative glitch sizes in a pulsar at a certain observational time span. In this analysis, we test the possibility of of quantifying the $ A_{g} $ with the pulsars main spin frequency derivatives (i.e. $ \dot{\nu} $ and $\ddot{\nu} $). In this approach, the ratio of the frequency derivatives, i.e. $ |\ddot{\nu}|/\dot{\nu}^{2} $ is seen to constrains the glitch activity in radio pulsars. The glitch size is found to be independent of the magnitude of the ratio, however, based on the recorded glitch events, the lower end of $ |\ddot{\nu}|/\dot{\nu}^{2} $ distribution appear to have more glitches. The minimum inter-glitch time interval in the ensemble of pulsars scale with the ratio as $t_{g} \sim 3.35(|\ddot{\nu}|/\dot{\nu}^{2})^{0.23} $. The $ A_{g} $ quantified in this analysis supports the idea of neutron star inner-crust superfluid being the reservoir of momentum transferred during glitches. It suggests that the moment of inertia of the inner-crust to be at most 10 % of the entire neutron star moment of inertia.

Revisiting the Post-glitch Relaxation of the 2000 Vela Glitch with the Neutron Star Equation of States in the Brueckner and Relativistic Brueckner Theories

Abstract We revisit the short-term post-glitch relaxation of the Vela 2000 glitch in the simple two-component model of the pulsar glitch by making use of the latest realistic equations of states from the microscopic Brueckner and the relativistic Brueckner theories for neutron stars, which can reconcile with the available astrophysical constraints. We show that to fit both the glitch size and the post-glitch jumps in frequency derivatives approximately 1 minute after the glitch, the mass of the Vela pulsar is necessarily small, and there may be demands for a stiff equation of state (which results in a typical stellar radius larger than ∼12.5 km) and a strong suppression of the pairing gap in the nuclear medium. We discuss the implications of this result on the understanding of pulsar glitches.

Long-term Statistics of Pulsar Glitches Due to History-dependent Avalanches

Abstract Stress accumulation-relaxation meta-models of pulsar glitches make precise, microphysics-agnostic predictions of long-term glitch statistics, which can be falsified by existing and future timing data. Previous meta-models assume that glitches are triggered by an avalanche process, e.g., involving superfluid vortices, and that the probability density function (PDF) of the avalanche sizes is history independent and specified exogenously. Here, a recipe is proposed to generate the avalanche sizes endogenously in a history-dependent manner, by tracking the thresholds of occupied vortex pinning sites as a function of time. Vortices unpin spasmodically from sites with thresholds below a global, time-dependent stress and repin at sites with thresholds above the global stress, imbuing the system with long-term memory. The meta-model predicts PDFs, auto-, and cross-correlations for glitch sizes and waiting times, which are provisionally inconsistent with current observations, unlike some previous meta-models (e.g., state-dependent Poisson process), whose predictions are consistent. The theoretical implications are intriguing, albeit uncertain, because history-dependent avalanches embody faithfully the popular, idealized understanding in the literature of how vortex unpinning operates as a driven, stochastic process. The meta-model predicts aftershocks, which occur with larger than average sizes and longer than average waiting times after the largest, system-resetting glitches. This prediction will be tested, once more data are generated by the next generation of pulsar timing campaigns.

30 glitches in 18 radio pulsars

Pulsar timing is a classic technology of detecting irregularities in pulsar rotation. We carried out this method for 18 young radio pulsars, with long-term timing observations obtained between 2007 and 2015 using the Parkes 64-m radio telescope. As a result, 30 glitches were identified, ranging from 0.75 × 10−9 to 8.6 × 10−6 in the relative glitch sizes Δν/ν, where ν = 1/P is the pulse frequency. These glitches are composed of 26 new glitches and four published glitches with new exponential recoveries. All pulsars exhibit normal glitches, and six pulsars were observed to undergo a glitch event for the first time. We discuss the properties and implications for neutron-star physics of these glitches, and show that they are in agreement with previous work, except that the cumulative probability distributions of the mean waiting times for PSRs J0537–6910, J1341–6220 and J1740–3015 are not in consonance with the Poisson model.

Constraining mechanism associated with fast radio burst and glitch from SGR J1935

ABSTRACT The discovery of fast radio burst (FRB) 200428 from galactic SGR J1935+2154 makes it possible to measure rotational changes accompanied by FRBs and to test several FRB models which may be simultaneously associated with glitches. Inspired by this idea, we present order of magnitude calculations to the scenarios proposed. FRB models such as global starquakes, crust fractures, and collisions between pulsars and asteroids/comets are discussed. For each mechanism, the maximum glitch sizes are constrained by the isotropic energy release during the X-ray burst and/or the SGR J1935+2154-like radio burst rate. Brief calculations show that, the maximum glitch sizes for different mechanisms differ by order(s) of magnitude. If glitches are detected to be coincident with FRBs from galactic magnetars in the future, glitch behaviours (such as glitch size, rise time-scale, the recovery coefficient, and spin-down rate offset) are promising to serve as criterions to distinguish glitch mechanisms and in turn to constrain FRB models.

Effects of periodicity in observation scheduling on parameter estimation of pulsar glitches

ABSTRACT In certain pulsar timing experiments, where observations are scheduled approximately periodically (e.g. daily), timing models with significantly different frequencies (including but not limited to glitch models with different frequency increments) return near-equivalent timing residuals. The average scheduling aperiodicity divided by the phase error due to time-of-arrival uncertainties is a useful indicator when the degeneracy is important. Synthetic data are used to explore the effect of this degeneracy systematically. It is found that phase-coherent tempo2 or temponest-based approaches are biased sometimes towards reporting small glitch sizes regardless of the true glitch size. Local estimates of the spin frequency alleviate this bias. A hidden Markov model is free from bias towards small glitches and announces explicitly the existence of multiple glitch solutions but sometimes fails to recover the correct glitch size. Two glitches in the UTMOST public data release are reassessed, one in PSR J1709−4429 at MJD 58178 and the other in PSR J1452−6036 at MJD 58600. The estimated fractional frequency jump in PSR J1709−4429 is revised upward from Δf/f = (54.6 ± 1.0) × 10−9 to (2432.2 ± 0.1) × 10−9 with the aid of additional data from the Parkes radio telescope. We find that the available UTMOST data for PSR J1452−6036 are consistent with Δf/f = 270 × 10−9 + N/(fT) with N = 0, 1, and 2, where $T \approx 1\, \text{sidereal day}$ is the observation scheduling period. Data from the Parkes radio telescope can be included, and the N = 0 case is selected unambiguously with a combined data set.

Small glitches and other rotational irregularities of the Vela pulsar

Context. Glitches are sudden increases in the rotation rate ν of neutron stars, which are thought to be driven by the neutron superfluid inside the star. The Vela pulsar presents a comparatively high rate of glitches, with 21 events reported since observations began in 1968. These are amongst the largest known glitches (17 of them have sizes Δν/ν ≥ 10−6) and exhibit very similar characteristics. This similarity, combined with the regularity with which large glitches occur, has turned Vela into an archetype of this type of glitching behaviour. The properties of its smallest glitches, on the other hand, are not clearly established. Aims. We explore the population of small-amplitude, rapid rotational changes in the Vela pulsar and determine the rate of occurrence and sizes of its smallest glitches. This will help advance our understanding of the actual distribution of glitch sizes and inter-glitch waiting times in this pulsar, which has implications for theoretical models of the glitch mechanism. Methods. High-cadence observations of the Vela pulsar were taken between 1981 and 2005 at the Mount Pleasant Radio Observatory. An automated systematic search was carried out that investigated whether a significant change of spin frequency ν and/or the spin-down rate ν̇ takes place at any given time. Results. We find two glitches that have not been reported before, with respective sizes Δν/ν of (5.55 ± 0.03) × 10−9 and (38 ± 4) × 10−9. The latter is followed by an exponential-like recovery with a characteristic timescale of 31 d. In addition to these two glitch events, our study reveals numerous events of all possible signatures (i.e. combinations of Δν and Δν̇ signs), all of them small with |Δν|/ν < 10−9, which contribute to the Vela timing noise. Conclusions. The Vela pulsar presents an under-abundance of small glitches compared to many other glitching pulsars, which appears genuine and not a result of observational biases. In addition to typical glitches, the smooth spin-down of the pulsar is also affected by an almost continuous activity that can be partially characterised by small step-like changes in ν, ν̇ or both. Simulations indicate that a continuous wandering of the rotational phase, following a red spectrum, could mimic such step-like changes in the timing residuals.