Comprehensive Analysis of Optical brightness and Color Variability of Blazars in the ZTF Survey DR22
Pith reviewed 2026-06-28 09:17 UTC · model grok-4.3
The pith
BL Lacs show a bluer-when-brighter trend while FSRQs show redder-when-brighter, with brighter states favoring bluer trends in both.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The analysis establishes that BL Lacs display a BWB trend in 14.7 percent of cases versus 2.3 percent RWB, while FSRQs display BWB in 8.8 percent versus RWB in 14.1 percent; brighter states raise the chance of BWB behavior in both classes, and sources with a BWB trend in BL Lacs vary more than those with RWB while the opposite holds for FSRQs.
What carries the argument
The joint statistical dependence of color sign on blazar subclass, instantaneous brightness state, and measured variability amplitude, extracted from quasi-simultaneous g-r light curves.
If this is right
- BL Lacs that follow BWB vary more strongly than those that follow RWB.
- FSRQs that follow RWB vary more strongly than those that follow BWB.
- Both subclasses become more likely to show BWB behavior during their brighter intervals.
- Color evolution therefore cannot be described by subclass alone; brightness state and variability amplitude must also be included.
Where Pith is reading between the lines
- Jet emission models may need separate prescriptions for the particle populations or magnetic-field geometry that dominate at high versus low states.
- The same analysis applied to multi-band data from other surveys could test whether the type-dependent trends persist at radio or X-ray wavelengths.
- If the trends survive host subtraction, they supply an observational prior for simulations that track how synchrotron or inverse-Compton components shift with flux.
Load-bearing premise
The measured color trends arise from the blazar emission itself rather than residual host-galaxy light or from the particular way variability amplitude and color slope are calculated.
What would settle it
Repeating the identical source selection and trend classification after subtracting host-galaxy contributions or after switching to a different variability metric reverses the reported BWB and RWB fractions for either subclass.
Figures
read the original abstract
This study conducts a comprehensive analysis of brightness and color variability in blazars, utilizing over six years of quasi-simultaneous g-band and r-band data from 1149 sources in the ZTF Data Release 22 (DR22), including 589 BL Lacs and 560 FSRQs. We quantify the amplitude of variability and the fractional root mean square (rms) variability for each source and statistically assess the overall and short-term color behaviors across different subclasses; examine the distribution of brightness variability characteristics across different blazar types and investigate how the extent of variability correlates with color trends. We found BL Lacs tend to exhibit a BWB (bluer when brighter) trend, while FSRQs display a RWB (redder when brighter) trend; BL Lacs with negligible host-galaxy contamination exhibit a BWB trend fraction of 14.7% (68/462) compared to 2.3% (11/462) for RWB trend, while FSRQs show 8.8% (49/560) BWB trend versus 14.1% (79/560) RWB trend. By statistically investigating how color behavior depends on brightness state across different timescales, we find that brighter states in both BL Lacs and FSRQs are more likely to exhibit BWB trend. Our results also show that BL Lacs with a BWB trend exhibit higher variability than those with a RWB trend, whereas FSRQs with a RWB trend display significantly greater variability than those with a BWB trend. These results suggest that blazar color variability depends jointly on source type, brightness state, and variability amplitude, highlighting the complexity of color evolution in blazars.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes brightness and color variability in 1149 blazars (589 BL Lacs and 560 FSRQs) from over six years of ZTF DR22 g- and r-band data. It quantifies variability amplitudes and fractional rms variability, assesses color trends, and reports that BL Lacs predominantly exhibit BWB trends (14.7% vs 2.3% RWB in the negligible host-contamination subset of 462 sources) while FSRQs exhibit RWB trends (8.8% BWB vs 14.1% RWB); brighter states favor BWB trends in both classes, and variability amplitude correlates differently with trends by source type. The central claim is that blazar color variability depends jointly on source type, brightness state, and variability amplitude.
Significance. The large sample of over 1100 sources provides substantial statistical power for measuring color-magnitude trends across blazar subclasses and brightness states. If the host-contamination controls and trend definitions prove robust, these empirical fractions and correlations would offer useful constraints on jet emission models versus host or disk contributions.
major comments (3)
- [Abstract] Abstract: the BWB/RWB fractions are reported only after restricting BL Lacs to the 462/589 sources with 'negligible host-galaxy contamination,' with no equivalent control or quantitative test described for the remaining BL Lacs or for any FSRQs; because host light can preferentially redden fainter states, this is load-bearing for the joint-dependence claim on source type.
- [Results on color trends] Results section on color trends: the thresholds used to classify BWB versus RWB trends (slope sign, significance cut, or fractional rms variability threshold) receive no sensitivity analysis; different choices could alter the reported 14.7 % vs 2.3 % and 8.8 % vs 14.1 % numbers and thereby the claimed joint dependence on type, state, and amplitude.
- [Brightness state analysis] Section on brightness-state dependence: the definitions of 'brighter states,' the separation of short-term versus long-term trends, and the exact statistical test for the BWB preference in brighter states are not specified in sufficient detail to evaluate whether the reported trend is robust to binning or metric choices.
minor comments (1)
- [Abstract] Abstract: the phrase 'statistically assess the overall and short-term color behaviors across different subclasses' is used without indicating the timescale binning or metric employed.
Simulated Author's Rebuttal
We thank the referee for the constructive comments on our manuscript. We address each major comment below and indicate planned revisions to improve clarity and robustness.
read point-by-point responses
-
Referee: [Abstract] Abstract: the BWB/RWB fractions are reported only after restricting BL Lacs to the 462/589 sources with 'negligible host-galaxy contamination,' with no equivalent control or quantitative test described for the remaining BL Lacs or for any FSRQs; because host light can preferentially redden fainter states, this is load-bearing for the joint-dependence claim on source type.
Authors: We acknowledge the referee's concern regarding the host-contamination control. Section 3.2 of the manuscript describes the quantitative criterion (host contribution fraction <5% based on SED fitting) used to select the 462 BL Lacs. FSRQs are reported for the full sample because their nuclear emission dominates across the observed flux range, rendering host effects minimal; this is supported by literature on FSRQ SEDs. For the remaining 127 BL Lacs, we will add a supplementary table and brief discussion in the revised Results section comparing color-trend fractions with and without the contamination cut to test robustness. The abstract will be updated to explicitly note the sample restriction for BL Lacs. revision: partial
-
Referee: [Results on color trends] Results section on color trends: the thresholds used to classify BWB versus RWB trends (slope sign, significance cut, or fractional rms variability threshold) receive no sensitivity analysis; different choices could alter the reported 14.7 % vs 2.3 % and 8.8 % vs 14.1 % numbers and thereby the claimed joint dependence on type, state, and amplitude.
Authors: We agree that sensitivity tests would strengthen the classification. The current thresholds are a linear slope sign with >2σ significance on the color-magnitude fit and a minimum fractional rms variability of 0.1. In the revised manuscript we will add a dedicated subsection (or appendix) performing sensitivity analysis by varying the significance cut (1σ–3σ) and rms threshold (0.05–0.2), reporting the resulting BWB/RWB fractions for both classes. This will directly address potential dependence on threshold choice. revision: yes
-
Referee: [Brightness state analysis] Section on brightness-state dependence: the definitions of 'brighter states,' the separation of short-term versus long-term trends, and the exact statistical test for the BWB preference in brighter states are not specified in sufficient detail to evaluate whether the reported trend is robust to binning or metric choices.
Authors: We apologize for insufficient detail in the original text. Brighter states are defined as individual epochs where the g-band flux exceeds the source's median flux by >20%; short-term trends are computed within single ZTF observing seasons while long-term trends use the full six-year baseline. The preference for BWB in brighter states is assessed via a two-proportion z-test on the fraction of sources showing positive color-magnitude slopes in the brighter versus fainter subsamples. The revised Methods and Results sections will include explicit equations, the exact z-test implementation, and additional robustness checks using alternative definitions (e.g., 10% or 30% flux thresholds and different seasonal binning). revision: yes
Circularity Check
No circularity; purely empirical observational analysis
full rationale
The paper reports direct statistical measurements of variability amplitudes (fractional rms) and color-magnitude trends (BWB/RWB fractions) computed from ZTF g/r light curves for 1149 sources. No equations, fitted models, or derivations are presented whose outputs reduce by construction to the inputs; all reported fractions and correlations are counts and comparisons performed on the observed data after applying selection cuts. No self-citation load-bearing steps, ansatzes, or uniqueness theorems appear in the derivation chain. The analysis is self-contained against the external ZTF dataset.
Axiom & Free-Parameter Ledger
free parameters (1)
- Thresholds for defining BWB and RWB trends
axioms (2)
- domain assumption The g and r band data are sufficiently simultaneous to measure color without lag-induced bias.
- domain assumption Host galaxy contamination is negligible or correctly accounted for in the BL Lac subsample.
Reference graph
Works this paper leans on
-
[1]
Astropy Collaboration, Price-Whelan, A. M., Sip˝ ocz, B. M., et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f
-
[2]
Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214, doi: 10.1086/170270
-
[3]
Bellm, E. C., Kulkarni, S. R., Graham, M. J., et al. 2019, PASP, 131, 018002, doi: 10.1088/1538-3873/aaecbe
-
[4]
Bonning, E., Urry, C. M., Bailyn, C., et al. 2012, ApJ, 756, 13, doi: 10.1088/0004-637X/756/1/13 B¨ ottcher, M. 2019, Galaxies, 7, 20, doi: 10.3390/galaxies7010020
-
[5]
Britzen, S., Fendt, C., Witzel, G., et al. 2018, MNRAS, 478, 3199, doi: 10.1093/mnras/sty1026
-
[6]
Chandra, S., Baliyan, K. S., Ganesh, S., & Joshi, U. C. 2011, ApJ, 731, 118, doi: 10.1088/0004-637X/731/2/118
-
[7]
Frank, J., King, A., & Raine, D. J. 2002, Accretion Power in Astrophysics: Third Edition
2002
-
[8]
Gaur, H., Gupta, A. C., Strigachev, A., et al. 2012, MNRAS, 425, 3002, doi: 10.1111/j.1365-2966.2012.21583.x
-
[9]
Gu, M.-F., & Ai, Y. L. 2011, A&A, 528, A95, doi: 10.1051/0004-6361/201016280
-
[10]
F., Lee, C.-U., Pak, S., Yim, H
Gu, M. F., Lee, C.-U., Pak, S., Yim, H. S., & Fletcher, A. B. 2006, A&A, 450, 39, doi: 10.1051/0004-6361:20054271
-
[11]
2004, A&A, 421, 877, doi: 10.1051/0004-6361:20040235
Guetta, D., Ghisellini, G., Lazzati, D., & Celotti, A. 2004, A&A, 421, 877, doi: 10.1051/0004-6361:20040235
-
[12]
Heidt, J., & Wagner, S. J. 1996, A&A, 305, 42, doi: 10.48550/arXiv.astro-ph/9506032
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.astro-ph/9506032 1996
-
[13]
2017, AJ, 154, 42, doi: 10.3847/1538-3881/aa799a
Hong, S., Xiong, D., & Bai, J. 2017, AJ, 154, 42, doi: 10.3847/1538-3881/aa799a
-
[14]
Hovatta, T., Valtaoja, E., Tornikoski, M., & L¨ ahteenm¨ aki, A. 2009, A&A, 494, 527, doi: 10.1051/0004-6361:200811150
-
[15]
Hu, S. M., Chen, X., Guo, D. F., Jiang, Y. G., & Li, K. 2014, MNRAS, 443, 2940, doi: 10.1093/mnras/stu1373
-
[16]
2011, PASJ, 63, 639, doi: 10.1093/pasj/63.3.327
Ikejiri, Y., Uemura, M., Sasada, M., et al. 2011, PASJ, 63, 639, doi: 10.1093/pasj/63.3.327
-
[17]
Isler, J. C., Urry, C. M., Coppi, P., et al. 2017, ApJ, 844, 107, doi: 10.3847/1538-4357/aa79fc
-
[18]
Jiang, Y.-F., & Blaes, O. 2020, ApJ, 900, 25, doi: 10.3847/1538-4357/aba4b7
-
[19]
2015, MNRAS, 449, 94, doi: 10.1093/mnras/stv241
Kokubo, M. 2015, MNRAS, 449, 94, doi: 10.1093/mnras/stv241
-
[20]
Larionov, V. M., Villata, M., & Raiteri, C. M. 2010, A&A, 510, A93, doi: 10.1051/0004-6361/200913536
-
[21]
, archivePrefix = "arXiv", eprint =
Li, S.-L., & Cao, X. 2008, MNRAS, 387, L41, doi: 10.1111/j.1745-3933.2008.00480.x
-
[22]
2016, Ap&SS, 361, 345, doi: 10.1007/s10509-016-2934-6
Mao, L., & Zhang, X. 2016, Ap&SS, 361, 345, doi: 10.1007/s10509-016-2934-6
-
[23]
Marscher, A. P., & Gear, W. K. 1985, ApJ, 298, 114, doi: 10.1086/163592
-
[24]
2015, Ap&SS, 357, 75, doi: 10.1007/s10509-015-2254-2 13
Massaro, E., Maselli, A., Leto, C., et al. 2015, Ap&SS, 357, 75, doi: 10.1007/s10509-015-2254-2 13
-
[25]
Massaro, E., Nesci, R., Maesano, M., Montagni, F., & D’Alessio, F. 1998, MNRAS, 299, 47, doi: 10.1046/j.1365-8711.1998.01696.x
-
[26]
Negi, V., Joshi, R., Chand, K., et al. 2022, MNRAS, 510, 1791, doi: 10.1093/mnras/stab3591
-
[27]
Papadakis, I. E., Villata, M., & Raiteri, C. M. 2007, A&A, 470, 857, doi: 10.1051/0004-6361:20077516
-
[28]
Pereyra, N. A., Vanden Berk, D. E., Turnshek, D. A., et al. 2006, ApJ, 642, 87, doi: 10.1086/500919
-
[29]
2005, MNRAS, 364, 1105, doi: 10.1111/j.1365-2966.2005.09655.x
Pian, E., Falomo, R., & Treves, A. 2005, MNRAS, 361, 919, doi: 10.1111/j.1365-2966.2005.09216.x
-
[30]
Poon, H., Fan, J. H., & Fu, J. N. 2009, ApJS, 185, 511, doi: 10.1088/0067-0049/185/2/511
-
[31]
Potter, W. J., & Cotter, G. 2015, MNRAS, 453, 4070, doi: 10.1093/mnras/stv1657
-
[32]
M., Villata M., Acosta-Pulido J
Raiteri, C. M., Villata, M., Acosta-Pulido, J. A., et al. 2017, Nature, 552, 374, doi: 10.1038/nature24623
-
[33]
Safna, P. Z., Stalin, C. S., Rakshit, S., & Mathew, B. 2020, MNRAS, 498, 3578, doi: 10.1093/mnras/staa2622
-
[34]
Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103, doi: 10.1088/0004-637X/737/2/103
work page internal anchor Pith review doi:10.1088/0004-637x/737/2/103 2011
-
[35]
I., & Sunyaev, R
Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337
1973
-
[36]
2001, MNRAS, 322, 247, doi: 10.1046/j.1365-8711.2001.04072.x 14E
Spada, M., Ghisellini, G., Lazzati, D., & Celotti, A. 2001, MNRAS, 325, 1559, doi: 10.1046/j.1365-8711.2001.04557.x
-
[37]
M., Scarpa, R., O’Dowd, M., et al
Urry, C. M., Scarpa, R., O’Dowd, M., et al. 2000, ApJ, 532, 816, doi: 10.1086/308616
-
[38]
2003, ApJ, 590, 123, doi: 10.1086/374889
Vagnetti, F., Trevese, D., & Nesci, R. 2003, ApJ, 590, 123, doi: 10.1086/374889
-
[39]
Vaughan, S., Edelson, R., Warwick, R. S., & Uttley, P. 2003, MNRAS, 345, 1271, doi: 10.1046/j.1365-2966.2003.07042.x
-
[40]
Villata, M., Raiteri, C. M., Kurtanidze, O. M., et al. 2002, A&A, 390, 407, doi: 10.1051/0004-6361:20020662
-
[41]
Villata, M., Raiteri, C. M., Kurtanidze, O. M., et al. 2004, A&A, 421, 103, doi: 10.1051/0004-6361:20035895
-
[42]
2007, AJ, 133, 1599, doi: 10.1086/511773
Wu, J., Zhou, X., Ma, J., et al. 2007, AJ, 133, 1599, doi: 10.1086/511773
-
[43]
2014, Research in Astronomy and Astrophysics, 14, 933, doi: 10.1088/1674-4527/14/8/004
Zhang, B.-K., Zhao, X.-Y., Wang, C.-X., & Dai, B.-Z. 2014, Research in Astronomy and Astrophysics, 14, 933, doi: 10.1088/1674-4527/14/8/004
-
[44]
2024, MNRAS, 529, 3699, doi: 10.1093/mnras/stae765
Zhang, X., Xiong, D.-r., Gao, Q.-g., et al. 2024, MNRAS, 529, 3699, doi: 10.1093/mnras/stae765
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.