Relativistic effects dominate the emission of blazar jets complicating our understanding of their intrinsic properties. Although many methods have been proposed to account for them, the variability Doppler factor method has been shown to describe the blazar populations best. We use a Bayesian hierarchical code called Magnetron to model the light curves of 973 sources observed by the Owens Valley Radio Observatory's 40-m telescope as a series of flares with an exponential rise and decay, and estimate their variability brightness temperature. Our analysis allows us to place the most stringent constraints on the equipartition brightness temperature i.e., the maximum achieved intrinsic brightness temperature in beamed sources which we found to be =2.78x10^11^K+/-26%. Using our findings we estimated the variability Doppler factor for the largest sample of blazars increasing the number of available estimates in the literature by almost an order of magnitude. Our results clearly show that {gamma}-ray loud sources have faster and higher amplitude flares than {gamma}-ray quiet sources. As a consequence they show higher variability brightness temperatures and thus are more relativistically beamed, with all of the above suggesting a strong connection between the radio flaring properties of the jet and {gamma}-ray emission.

Cone search capability for table J/ApJ/866/137/table2 (Variability brightness temperatures and beaming properties for the sources in our sample (T_var_ and delta_var_ for bright blazars))

Cone search capability for table J/ApJ/866/137/table3 (Variability brightness temperatures and beaming properties for the sources in our sample without a redshift estimate (Tvar and deltavar for bright blazars))