A Comprehensive Study on the Mid-Infrared Variability of Blazars

We present a comprehensive investigation of mid-infrared (MIR) flux variability at 3.4 μm (W1 band) for a large sample of 3816 blazars, using Wide-field Infrared Survey Explorer (WISE) data through December 2022. The sample consists of 1740 flat-spectrum radio quasars (FSRQs), 1281 BL Lac objects (BL Lacs), and 795 blazars of uncertain type (BCUs). Considering Fermi Large Area Telescope detection, we classify 2331 as Fermi blazars and 1485 as non-Fermi blazars. Additionally, based on synchrotron peak frequency, the sample includes 2264 low-synchrotron peaked (LSP), 512 intermediate-synchrotron peaked (ISP), and 655 high-synchrotron peaked (HSP) sources. We conduct a comparative analysis of short- and long-term intrinsic variability amplitude (σm), duty cycle (DC), and ensemble structure function (ESF) across blazar subclasses. The median short-term σm values were 0.181−0.106+0.153, 0.104−0.054+0.101, 0.135−0.076+0.154, 0.173−0.097+0.158, 0.177−0.100+0.156, 0.096−0.050+0.109, and 0.106−0.058+0.100 mag for FSRQs, BL Lacs, Fermi blazars, non-Fermi blazars, LSPs, ISPs, and HSPs, respectively. The median DC values were 71.03−22.48+14.17, 64.02−22.86+16.97, 68.96−25.52+15.66, 69.40−22.17+14.42, 71.24−21.36+14.25, 63.03−33.19+16.93, and 64.63−24.26+15.88 percent for the same subclasses. The median long-term σm values were 0.137−0.105+0.408, 0.171−0.132+0.206, 0.282−0.184+0.332, 0.071−0.062+0.143, 0.218−0.174+0.386, 0.173−0.132+0.208, and 0.101−0.077+0.161 mag for the same subclasses, respectively. Our results reveal significant differences in 3.4 μm flux variability among these subclasses. FSRQs (LSPs) exhibit larger σm and DC values compared to BL Lacs (ISPs and HSPs). Fermi blazars display higher long-term σm but lower short-term σm relative to non-Fermi blazars, while DC distributions between the two groups are similar. ESF analysis further confirms the greater variability of FSRQs, LSPs, and Fermi blazars across a wide range of time scales compared to BL Lacs, ISPs/HSPs, and non-Fermi blazars. These findings highlight a close correlation between MIR variability and blazar properties, providing valuable insights into the underlying physical mechanisms responsible for their emission.