Radio Variability Research Articles

Radio Variability and Broad-Band Spectra of Infrared Galaxies with and without OH Megamaser Emission

We study the radio variability of galaxies with and without sources of hydroxyl (OH) megamaser radiation based on the continuum radio measurements conducted in 2019-2022 with the radio telescope RATAN-600 at frequencies of 2.3, 4.7, 8.2, and 11.2 GHz. Presumably, radio continuum emission significantly affects the megamaser radiation brightness, therefore, such a characteristic as the variability of radio emission is important for determining the OHM galaxies parameters. With additional data from the literature, the parameters of radio variability on a time scale up to 30 years were estimated. The median values of the variability index for 48 OHM galaxies are in the range $V_{S}=0.08$-$0.17$, and for 30 galaxies without OH emission they are $V_{S}=0.08$-$0.28$. For some individual galaxies in both samples, flux density variations reach 30-50%. These sources either are commonly associated with AGNs or reveal active star formation. Generally, the variability of luminous infrared galaxies with and without OH megamaser emission is moderate and of the same order of magnitude on long time scales. From estimating the spectral energy distribution parameters in a broad frequency range (from MHz to THz), we determined the spectral index below 50 GHz and the color temperatures of dust components for megamaser and control sample galaxies. At a level of $\rho<0.05$, there are no statistically significant differences in the distribution of these parameters for the two samples, as well there are no statistically significant correlations between the dust color temperatures and the variability index or luminosity in the OH line.

Jet kinematics in the transversely stratified jet of 3C 84

3C 84 (NGC 1275) is one of the brightest radio sources in the millimetre radio bands, which led to a plethora of very-long-baseline interferometry (VLBI) observations at numerous frequencies over the years. They reveal a two-sided jet structure, with an expanding but not well-collimated parsec-scale jet, pointing southward. High-resolution millimetre-VLBI observations allow the study and imaging of the jet base on a sub-parsec scale. This could facilitate the investigation of the nature of the jet origin, also in view of the previously detected two-railed jet structure and east-west oriented core region seen withRadioAstronat 22 GHz. We produced VLBI images of this core and inner jet region, observed over the past twenty years at 15, 43, and 86 GHz. We determined the kinematics of the inner jet and ejected features at 43 and 86 GHz and compared their ejection times with radio andγ-ray variability. For the moving jet features, we find an average velocity ofβappavg = 0.055−0.22c (μavg = 0.04 − 0.18 mas yr−1). From the time-averaged VLBI images at the three frequencies, we measured the transverse jet width along the bulk flow. On the ≤1.5 parsec scale, we find a clear trend of the jet width being frequency dependent, with the jet being narrower at higher frequencies. This stratification is discussed in the context of a spine-sheath scenario, and we compare it to other possible interpretations. From quasi-simultaneous observations at 43 and 86 GHz, we obtain spectral index maps, revealing a time-variable orientation of the spectral index gradient due to structural variability of the inner jet.

Compact and Variable Radio Emission from an Active Galaxy with Supersoft X-Ray Emission

RX J1301.9+2747 is a unique active galaxy with a supersoft X-ray spectrum that lacks significant emission at energies above 2 keV. In addition, it is one of few galaxies displaying quasiperiodic X-ray eruptions that recur on a timescale of 13–20 ks. We present multiepoch radio observations of RX J1301.9+2747 using GMRT, Very Large Array (VLA), and Very Long Baseline Array (VLBA). The VLBA imaging at 1.6 GHz reveals a compact radio emission unresolved at a scale of <0.7 pc, with a brightness temperature of T b > 5 × 107 K. The radio emission is variable by more than a factor of 2.5 over a few days, based on the data taken from VLA monitoring campaigns. The short-term radio variability suggests that the radio emitting region has a size as small as 8 × 10−4 pc, resulting in an even higher brightness temperature of T b ∼ 1012 K. A similar limit on the source size can be obtained if the observed flux variability is not intrinsic and caused by the interstellar scintillation effect. The overall radio spectrum is steep with a time-averaged spectral index α = −0.78 ± 0.03 between 0.89 and 14 GHz. These observational properties rule out a thermal or star formation origin of the radio emission, and appear to be consistent with the scenario of episodic jet ejections driven by a magnetohydrodynamic process. Simultaneous radio and X-ray monitoring observations down to a cadence of hours are required to test whether the compact and variable radio emission is correlated with the quasiperiodic X-ray eruptions.

Quasi-simultaneous observations of radio and X-ray variability in three radio-quiet Seyfert galaxies

ABSTRACT Radio variability in some radio-quiet (RQ) active galactic nuclei suggests emission from regions close to the central engine, possibly the outer accretion disc corona. If the origins of the radio and the X-ray emission are physically related, their emission may be temporarily correlated, possibly with some time delays. We present the results of quasi-simultaneous radio and X-ray monitoring of three RQ Seyfert galaxies, Mrk 110, Mrk 766, and NGC 4593, carried out with the Very Large Array at 8.5 GHz over a period of about 300 d, and with the Rossi X-ray Timing Explorer at 2–10 keV over a period of about 2000 d. The radio core variability is likely detected in the highest resolution (A configuration) observations of Mrk 110 and NGC 4593, with a fractional variability amplitude of 6.3 per cent and 9.5 per cent, respectively. A cross-correlation analysis suggests an apparently strong (Pearson r = −0.89) and highly significant correlation (p = 1 × 10−6) in Mrk 110, with the radio lagging the X-ray by 56 d. However, a further analysis of the r values distribution for physically unrelated long time delays reveals that this correlation is not significant. This occurs since the Pearson correlation assumes white noise, while both the X-ray and the radio light curves follow red noise, which dramatically increases the chance, by a factor of ∼103, to get extremely high r values in uncorrelated data sets. A significantly longer radio monitoring with a higher sampling rate, preferably with a high-resolution fixed radio array, is required in order to reliably detect a delay.

Results of studying the radio and optical variability properties of MRK 501 active galaxy

MRK 501 is a giant radio galaxy with an active core in the constellation Hercules. It is a powerful gamma source. MRK 501 has gamma cyclic variability as well as high X-ray and optical activity. MRK 501 is a quite powerful variable radio source, with a flux density at 15 GHz about 1.2 Jy. Of particular interest is the presence of cyclic brightness variations in gamma and optical ranges, which, according to modern models, may be caused by the presence of a supermassive double black hole in the core of the galaxy MRK 501. In this work, optical observatories of three countries, Ukraine (Mayaki, Heavenly Owl), Slovakia (Vihorlat) and Latvia (Baldone), took part in a joint project to study the rapid variability of AGN in radio and optical bands. The analysis of observations obtained from different observatories in V , R, I bands showed the presence of a wave-like, presumably cyclic, variability in the light curves with an average quasi-period of 46.7 days. Taking into account that a quasi-period of about 23 days was detected in the gamma range, it is most likely that we have found a harmonic (doubled value) of this period in the optical range, which is a new result. In addition, an analysis of long-term variability of MRK 501, from historical light curves of the AAVSO catalog, showed the presence of two longer cycles of variability with quasiperiods of about 1 year and about 5–6 years. It is shown that these quasi-periods, as well as their amplitudes, are unevenly distributed over an observation time interval from 2003 to 2021 (18 years) and change significantly over time.

Constraining the radio properties of the z = 6.44 QSO VIK J2318−3113

The recent detection of the quasi-stellar object (QSO) VIKING J231818.3−311346 (hereafter VIK J2318−3113) at redshift z = 6.44 in the Rapid ASKAP Continuum Survey (RACS) uncovered its radio-loud nature, making it one of the most distant known to date in this class. By using data from several radio surveys of the Galaxy And Mass Assembly 23h field and from a dedicated follow-up, we were able to constrain the radio spectrum of VIK J2318−3113 in the observed range ∼0.1–10 GHz. At high frequencies (0.888–5.5 GHz in the observed frame) the QSO presents a steep spectrum (αr = 1.24, with Sν ∝ ν−αr), while at lower frequencies (0.4–0.888 GHz in the observed frame) it is nearly flat. The overall spectrum can be modelled by either a curved function with a rest-frame turnover around 5 GHz, or with a smoothly varying double power law that is flat below a rest-frame break frequency of about 20 GHz and above which it significantly steepens. Based on the model adopted, we estimated that the radio jets of VIK J2318−3113 must be a few hundred years old in the case of a turnover, or less than a few × 104 years in the case of a break in the spectrum. Having multiple observations at two frequencies (888 MHz and 5.5 GHz), we further investigated the radio variability previously reported for this source. We found that the marginally significant flux density variations are consistent with the expectations from refractive interstellar scintillation, even though relativistic effects related to the orientation of the source may still play a non-negligible role. Further radio and X-ray observations are required to conclusively discern the nature of this variation.

Radiative shock oscillation model for the long-term flares of Sgr A*

ABSTRACT We examine time-dependent 2D relativistic radiation magnetohydrodynamic (MHD) flows to develop the shock oscillation model for the long-term flares of Sagittarius A* (Sgr A*). Adopting modified flow parameters in addition to the previous studies, we confirm quasi-periodic flares with periods of ∼5 and 10 d that are compatible with observations by Chandra, Swift, and XMM–Newton monitoring of Sgr A*. Using a simplified two-temperature model of ions and electrons, we find that the flare due to synchrotron emission lags that of bremsstrahlung emission by 1–2 h that are qualitatively comparable to the time lags of 1–5 h reported in several simultaneous observations of radio and X-ray variability in Sgr A*. The synchrotron emission is confined in a core region of 3Rg size with the strong magnetic field, while the bremsstrahlung emission mainly originates in a distant region of 10–20Rg behind the oscillating shock, where Rg is the Schwarzschild radius. The time lag is estimated as the transit time of the acoustic wave between the above two regions. The time-averaged distribution of radiation shows a strong anisotropic nature along the rotational axis but isotropic distribution in the radial direction. A high-velocity jet with ∼0.6c along the rotational axis is intermittently found in a narrow funnel region with a collimation angle ∼15°. The shock oscillating model explains well the flaring rate and the time lag between radio and X-ray emissions for the long-term flares of Sgr A*.

21 new long-term variables in the GX 339−4 field: two years of MeerKAT monitoring

ABSTRACT We present 21 new long-term variable radio sources found commensally in 2 yr of weekly MeerKAT monitoring of the low-mass X-ray binary GX 339−4. The new sources are vary on time-scales of weeks to months and have a variety of light-curve shapes and spectral index properties. Three of the new variable sources are coincident with multiwavelength counterparts; and one of these is coincident with an optical source in deep MeerLICHT images. For most sources, we cannot eliminate refractive scintillation of active galactic nuclei as the cause of the variability. These new variable sources represent 2.2 ± 0.5 per cent of the unresolved sources in the field, which is consistent with the 1–2 per cent variability found in past radio variability surveys. However, we expect to find short-term variable sources in the field and these 21 new long-term variable sources. We present the radio light curves and spectral index variability of the new variable sources, as well as the absolute astrometry and matches to coincident sources at other wavelengths.

New Tests of Milli-lensing in the Blazar PKS 1413 + 135

Abstract Symmetric achromatic variability (SAV) is a rare form of radio variability in blazars that has been attributed to gravitational milli-lensing by a ∼102–105 M ⊙ mass condensate. Four SAVs have been identified between 1980 and 2020 in the long-term radio monitoring data of the blazar PKS 1413 + 135. We show that all four can be fitted with the same, unchanging, gravitational lens model. If SAV is due to gravitational milli-lensing, PKS 1413 + 135 provides a unique system for studying active galactic nuclei with unprecedented microarcsecond resolution, as well as for studying the nature of the milli-lens itself. We discuss two possible candidates for the putative milli-lens: a giant molecular cloud hosted in the intervening edge-on spiral galaxy, and an undetected dwarf galaxy with a massive black hole. We find a significant dependence of SAV crossing time on frequency, which could indicate a fast shock moving in a slower underlying flow. We also find tentative evidence for a 989 day periodicity in the SAVs, which, if real, makes possible the prediction of future SAVs: the next three windows for possible SAVs begin in 2022 August, 2025 May, and 2028 February.

FRAMEx. II. Simultaneous X-Ray and Radio Variability in Active Galactic Nuclei—The Case of NGC 2992

Abstract Using simultaneous Very Long Baseline Array and Neil Gehrels Swift Observatory X-ray Telescope observations of the active galactic nucleus (AGN) in NGC 2992 over a six-month observing campaign, we observed a large drop in core 5 cm radio luminosity, by a factor of &gt;3, in tandem with a factor of &gt;5 increase in 2–10 keV X-ray luminosity. While NGC 2992 has long been an important object for studies of X-ray variability, our study is the first simultaneous X-ray and radio variability campaign for this object. We observe that the X-ray spectral index does not change over the course of the flare, consistent with a change in the bulk amount of Comptonizing plasma, potentially due to a magnetic reconnection event in the accretion disk. The drop in apparent radio luminosity can be explained by a change in free–free absorption, which we calculate to correspond to an ionized region with a physical extent and electron density consistent with the broad-line region (BLR). Our results are consistent with magnetic reconnection events in the dynamic accretion disk creating outbursts of ionizing material, increasing Compton up-scattering of UV accretion disk photons and feeding material into the BLR. These findings present an important physical picture for the dynamical relationship between X-ray and radio emission in AGNs.

Radio and γ-Ray Activity in the Jet of the Blazar S5 0716+714

Abstract We explore the connection between the γ-ray and radio emission in the jet of the blazar 0716+714 by using 15, 37, and 230 GHz radio and 0.1–200 GeV γ-ray light curves spanning 10.5 yr (2008–2019). We find significant positive and negative correlations between radio and γ-ray fluxes in different time ranges. The time delays between radio and γ-ray emission suggest that the observed γ-ray flares originated from multiple regions upstream of the radio core, within a few parsecs from the central engine. Using time-resolved 43 GHz Very Long Baseline Array maps we identified 14 jet components moving downstream along the jet. Their apparent speeds range from 6c to 26c, and they show notable variations in their position angles upstream from the stationary component (∼0.53 mas from the core). The brightness temperature declines as a function of distance from the core according to a power law that becomes shallower at the location of the stationary component. We also find that the periods at which significant correlations between radio and γ-ray emission occur overlap with the times when the jet was oriented to the north. Our results indicate that the passage of a propagating disturbance (or shock) through the radio core and the orientation of the jet might be responsible for the observed correlation between the radio and γ-ray variability. We present a scenario that connects the positive correlation and the unusual anticorrelation by combining the production of a flare and a dip at γ-rays by a strong moving shock at different distances from the jet apex.

Testing the Potential for Radio Variability in Disks around T Tauri Stars with Observations and Chemical Modeling

Abstract A multiwavelength observing campaign of the T Tauri star (TTS) GM Aur was undertaken in 2019 December. This campaign obtained Swift X-ray and NUV fluxes, HST NUV spectra, LCOGT u ′ g ′ r ′ i ′ and TESS photometry, CHIRON Hα spectra, Atacama Large Millimeter/submillimeter Array (ALMA) 13CO and C18O line fluxes, and Very Large Array (VLA) 3 cm continuum fluxes taken contemporaneously over one month. The X-ray to optical observations were presented previously. Here we present the ALMA and VLA data and make comparisons to GM Aur’s accretion and X-ray properties. We report no variability in the observed millimeter CO emission. Using disk chemistry models, we show that the magnitude of the changes seen in the FUV luminosity of GM Aur could lead to a variation of up to ∼6% in CO line emission and changes in the X-ray luminosity could lead to larger changes of ∼25%. However, the FUV and X-ray luminosity increases must last at least 100 years in order to induce changes, which seems implausible in the TTS stage; also, these changes would be too small to be detectable by ALMA. We report no variability in the 3 cm emission observed by the VLA, showing that changes of less than a factor of ∼3 in the accretion rates of TTSs do not lead to detectable changes in the mass-loss rate traced by the jet at centimeter wavelengths. We conclude that changes typically seen in the FUV and X-ray luminosities of TTSs do not lead to observable changes in millimeter CO line emission or jet centimeter continuum emission.

Daily variability at milli-arcsecond scales in the radio-quiet NLSy1 Mrk 110

ABSTRACT The origin of radio emission in the majority of active galactic nuclei is still poorly understood. Various competing mechanisms are likely involved in the production of radio emission and precise diagnostic tools are needed to disentangle them, of which variability is among the most powerful. For the first time, we show evidence for significant radio variability at 5 GHz at milli-arcsecond scales on days to weeks time-scales in the highly accreting and extremely radio-quiet (RQ) Narrow Line Seyfert 1 (NLSy1) Mrk 110. The simultaneous Swift/XRT light curve indicates stronger soft than hard X-ray variability. The short-term radio variability suggests that the GHz emitting region has a size smaller than ∼180 Schwarzschild radii. The high brightness temperature and the radio and X-ray variability rule out a star-formation and a disc wind origin. Synchrotron emission from a low-power jet and/or an outflowing corona is then favoured.

CHILES VERDES: Radio Variability at an Unprecedented Depth and Cadence in the COSMOS Field

Abstract Although it is well established that some extragalactic radio sources are time-variable, the properties of this radio variability, and its connection with host galaxy properties, remain to be explored—particularly for faint sources. Here we present an analysis of radio variable sources from the CHILES Variable and Explosive Radio Dynamic Evolution Survey (CHILES VERDES)—a partner project of the 1.4 GHz COSMOS H i Large Extragalactic Survey. CHILES VERDES provides an unprecedented combination of survey depth, duration, and cadence, with 960 hr of 1–2 GHz continuum VLA data obtained over 209 epochs between 2013 and 2019 in a 0.44 deg2 section of the well-studied extragalactic deep field, COSMOS. We identified 18 moderate-variability sources (showing 10%–30% flux density variation) and 40 lower-variability sources (2%–10% flux density variation). They are mainly active galactic nuclei (AGNs) with radio luminosities in the range of 1022–1027 W Hz−1, based on cross-matching with COSMOS multiwavelength catalogs. The moderate-variability sources span redshifts z = 0.22–1.56, have mostly flat radio spectra (α &gt; −0.5), and vary on timescales ranging from days to years. The lower-variability sources have similar properties, but generally have higher radio luminosities than the moderate-variability sources, extending to z = 2.8, and have steeper radio spectra (α &lt; −0.5). No star-forming galaxy showed statistically significant variability in our analysis. The observed variability likely originates from scintillation on short (∼week) timescales, and Doppler-boosted intrinsic AGN variability on long (month–year) timescales.

Milliarcsecond Core Size Dependence of the Radio Variability of Blazars

AbstractStudying blazar radio variability on timescales ranging from months to years provides information on the sub-parsec-scale structures of the jets, and the physics of the central active galactic nuclei. In this study, we focus on the radio variability of 1158 blazars observed at 15 GHz through the Owens Valley Radio Observatory Blazar Monitoring Program, where these sources have been observed about twice a week for over a decade. We investigate the dependence of the variability amplitudes and timescales, derived using a simple model fit to the structure function, on the milliarcsecond radio core sizes measured by Very Long Baseline Interferometry. The most compact sources exhibit larger variability amplitudes and shorter variability timescales than the more extended sources; this could be explained by light travel-time effects.

Millimeter-sized Dust Grains Surviving the Water-sublimating Temperature in the Inner 10 au of the FU Ori Disk

Abstract Previous observations have shown that the ≲10 au, ≳400 K hot inner disk of the archetypal accretion outburst young stellar object, FU Ori, is dominated by viscous heating. To constrain dust properties in this region, we have performed radio observations toward this disk using the Karl G. Jansky Very Large Array in 2020 June–July, September, and November. We also performed complementary optical photometric monitoring observations. We found that the dust thermal emission from the hot inner disk mid-plane of FU Ori has been approximately stationary and the maximum dust grain size is ≳1.6 mm in this region. If the hot inner disk of FU Ori, which is inward of the 150–170 K water snowline, is turbulent (e.g., corresponding to a Sunyaev &amp; Shakura viscous α t ≳ 0.1), or if the actual maximum grain size is still larger than the lower limit we presently constrain, then as suggested by the recent analytical calculations and the laboratory measurements, water-ice-free dust grains may be stickier than water-ice-coated dust grains in protoplanetary disks. Additionally, we find that the free–free emission and the Johnson B- and V-band magnitudes of these binary stars were brightening in 2016–2020. The optical and radio variability might be related to the dynamically evolving protostellar- or disk-accretion activities. Our results highlight that the hot inner disks of outbursting objects are important laboratories for testing models of dust grain growth. Given the active nature of such systems, to robustly diagnose the maximum dust grain sizes, it is important to carry out coordinated multiwavelength radio observations.

The MAVERIC Survey: Variable Jet-accretion Coupling in Luminous Accreting Neutron Stars in Galactic Globular Clusters

Abstract Accreting neutron stars in low-mass X-ray binaries show outflows—and sometimes jets—in the general manner of accreting black holes. However, the quantitative link between the accretion flow (traced by X-rays) and outflows and/or jets (traced by radio emission) is much less well understood for neutron stars than for black holes, other than the general observation that neutron stars are fainter in the radio at a given X-ray luminosity. We use data from the deep MAVERIC radio continuum survey of Galactic globular clusters for a systematic radio and X-ray study of six luminous (L X &gt; 1034 erg s−1) persistent neutron star X-ray binaries in our survey, as well as two other transient systems also captured by our data. We find that these neutron star X-ray binaries show an even larger range in radio luminosity than previously observed. In particular, in quiescence at L X ∼ 3 × 1034 erg s−1, the confirmed neutron star binary GRS 1747–312 in Terzan 6 sits near the upper envelope of the black hole radio/X-ray correlation, and the persistently accreting neutron star systems AC 211 (in M15) and X1850–087 (in NGC 6712) show unusual radio variability and luminous radio emission. We interpret AC 211 as an obscured “Z source” that is accreting at close to the Eddington limit, while the properties of X1850–087 are difficult to explain, and motivate future coordinated radio and X-ray observations. Overall, our results show that neutron stars do not follow a single relation between inflow and outflow, and confirm that their accretion dynamics are more complex than for black holes.

The variability of brightest cluster galaxies at high radio frequencies

ABSTRACT Variability of a galaxy’s core radio source can be a significant consequence of active galactic nucleus accretion. However, this variability has not been well studied, particularly at high radio frequencies. As such, we report on a campaign monitoring the high radio frequency variability of 20 nearby, cool-core brightest cluster galaxies. From our representative sample, we show that most vary significantly on time-scales of approximately 1 yr and longer. Our highest cadence observations are at 15 GHz and are from the Owens Valley Radio Observatory. They have a median time interval of 7 d and mostly span between 8 and 13 yr. We apply a range of variability detection techniques to the sources’ light curves to analyse changes on week to decade long time-scales. Most notably, at least half of the sources show 20 per cent peak to trough variability on 3 yr time-scales, while at least a third vary by 60 per cent on 6 yr time-scales. Significant variability, which is important to studies of the Sunyaev–Zel’dovich Effect in the radio/sub-mm, is therefore a common feature of these sources. We also show how the variability relates to spectral properties at frequencies of up to 353 GHz using data from the Korean VLBI network, the NIKA2 instrument of the IRAM 30-m telescope, and the SCUBA-2 instrument of the James Clerk Maxwell Telescope.

What defines a compact symmetric object? A carefully vetted sample of compact symmetric objects

AbstractCompact symmetric objects (CSOs), young jetted‐active galactic nuclei of overall projected size &lt;1 kpc, are of great interest due to their youth and evolution. The classification was introduced to distinguish between ∼95% of powerful compact extragalactic radio sources in flux density limited samples that are dominated by asymmetric emission due to relativistic beaming from jets aligned close to the line of sight and ∼5% of objects that are not. The original classification criteria were (a) overall projected diameter smaller than 1 kpc, (b) identified center of activity, and (c) symmetric jet structure about the center. There is confusion and erosion of the value of the CSO classification due to misclassifications. Many jets contain compact bright features outside the core, resulting in a gigahertz peaked spectrum total spectrum and a “compact double” appearance, and some objects with jet axes aligned close to the line of sight appear symmetric because the approaching jet is projected on both sides of the core. To eliminate the confusion, we propose adding (d) slow radio variability and (e) low apparent velocity of bright features moving along the jets to the above CSO criteria. We are compiling a catalog of CSOs using these five criteria to eliminate the confusion of Doppler boosting.

Identifying changing jets through their radio variability

Context. Supermassive black holes can launch highly relativistic jets with velocities reaching Lorentz factors of as high as Γ &gt; 50. How the jets accelerate to such high velocities and where along the jet they reach terminal velocity are open questions that are tightly linked to their structure as well as their launching and dissipation mechanisms. Aims. Changes in the beaming factor along the jets could potentially reveal jet acceleration, deceleration, or bending. We aim to (1) quantify the relativistic effects in multiple radio frequencies and (2) study possible jet velocity–viewing angle variations at parsec scales. Methods. We used the state-of-the-art code Magnetron to model light curves from the University of Michigan Radio Observatory and the Metsähovi Radio Observatory’s monitoring programs in five frequencies covering about 25 years of observations in the 4.8 to 37 GHz range for 61 sources. We supplement our data set with high-frequency radio observations in the 100–340 GHz range from ALMA, CARMA, and SMA. For each frequency we estimate the Doppler factor which we use to quantify possible changes in the relativistic effects along the jets. Results. The majority of our sources do not show any statistically significant difference in their Doppler factor across frequencies. This is consistent with constant velocity in a conical jet structure, as expected at parsec scales. However, our analysis reveals 17 sources where relativistic beaming changes as a function of frequency. In the majority of cases, the Doppler factor increases towards lower frequencies. Only 1253–053 shows the opposite behavior. By exploring their jet properties we find that the jet of 0420–014 is likely bent across the 4.8–340 GHz range. For 0212+735, the jet is likely parabolic, and still accelerating in the 4.8–37 GHz range. We discuss possible interpretations for the trends found in the remaining sources.