Comprehensive Analysis of Optical brightness and Color Variability of Blazars in the ZTF Survey DR22

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.