Variability Of Blazars Research Articles

Long-term optical-infrared color variability of blazars

The long-term optical and infrared color variability of blazars has been investigated with monitoring data from the Small and Moderate Aperture Research Telescope System (SMARTS). The sample in this study consists of 49 flat spectrum radio quasars (FSRQs) and 22 BL Lacertae objects (BL Lacs). The fractional variability amplitudes of each source have been calculated in both optical R band and infrared J band. Overall, the variability amplitudes of FSRQs are larger than those of BL Lacs. The results also suggest that the variability amplitude of most FSRQs is larger at a lower energy band (J band) than at a higher one (R band), but the variability amplitude of BL Lacs is larger at the higher energy band. Both types of blazars display color variation along with variability in brightness. However, they show different variation behaviors in general. In the whole data set, 35 FSRQs exhibit redder-when-brighter trends, and 11 FSRQs exhibit opposite trends; 11 BL Lacs follow bluer-when-brighter trends, and seven BL Lacs follow opposite trends. A detailed examination indicates that there are 10 blazars showing redder-when-brighter trends in their low state, and bluer-when-brighter or stable-when-brighter trends in their high state. Some more complicated color behaviors have also been detected in several blazars. The non-thermal jet emission and the thermal emission from the accretion disk are employed to explain the observed color behaviors.

FOURIER ANALYSIS OF BLAZAR VARIABILITY: KLEIN–NISHINA EFFECTS AND THE JET SCATTERING ENVIRONMENT

The strong variability of blazars can be characterized by power spectral densities (PSDs) and Fourier frequency-dependent time lags. In previous work, we created a new theoretical formalism for describing the PSDs and time lags produced via a combination of stochastic particle injection and emission via the synchrotron, synchrotron self-Compton, and external Compton (EC) processes. This formalism used the Thomson cross section and simple $\delta$-function approximations to model the synchrotron and Compton emissivities. Here we expand upon this work, using the full Compton cross section and detailed and accurate emissivities. Our results indicate good agreement between the PSDs computed using the $\delta$-function approximations and those computed using the accurate expressions, provided the observed photons are produced primarily by electrons with energies exceeding the lower limit of the injected particle population. Breaks are found in the PSDs at frequencies corresponding to the cooling timescales of the electrons primarily responsible for the observed emission, and the associated time lags are related to the difference in electron cooling timescales between the two energy channels, as expected. If the electron cooling timescales can be determined from the observed time lags and/or the observed EC PSDs, then one could in principle use the method developed here to determine the energy of the external seed photon source for EC, which is an important unsolved problem in blazar physics.

KINEMATICS OF AND EMISSION FROM HELICALLY ORBITING BLOBS IN A RELATIVISTIC MAGNETIZED JET

We present a general relativistic (GR) model of jet variability in active galactic nuclei due to orbiting blobs in helical motion along a funnel or cone shaped magnetic surface anchored to the accretion disk near the black hole. Considering a radiation pressure driven flow in the inner region, we find that it stabilizes the flow, yielding Lorentz factors ranging between 1.1 and 7 at small radii for reasonable initial conditions. Assuming these as inputs, simulated light curves (LCs) for the funnel model include Doppler and gravitational shifts, aberration, light bending, and time delay. These LCs are studied for quasi-periodic oscillations (QPOs) and the power spectral density (PSD) shape and yield an increased amplitude ($\sim$ 12 %); a beamed portion and a systematic phase shift with respect to that from a previous special relativistic model. The results strongly justify implementing a realistic magnetic surface geometry in Schwarzschild geometry to describe effects on emission from orbital features in the jet close to the horizon radius. A power law shaped PSD with a typical slope of $-2$ and QPOs with timescales in the range of $(1.37 - 130.7)$ days consistent with optical variability in Blazars, emerges from the simulations for black hole masses $M_{\bullet} = (0.5 - 5) \times 10^8 M_{\odot}$ and initial Lorentz factors $\gamma_{jet,i} = 2 - 10$. The models presented here can be applied to explain radio, optical, and X-ray variability from a range of jetted sources including active galactic nuclei, X-ray binaries and neutron stars.

THROUGH THE RING OF FIRE:γ-RAY VARIABILITY IN BLAZARS BY A MOVING PLASMOID PASSING A LOCAL SOURCE OF SEED PHOTONS

Blazars exhibit flares across the electromagnetic spectrum. Many $\gamma$-ray flares are highly correlated with flares detected at optical wavelengths; however, a small subset appears to occur in isolation, with little or no variability detected at longer wavelengths. These "orphan" $\gamma$-ray flares challenge current models of blazar variability, most of which are unable to reproduce this type of behavior. We present numerical calculations of the time-variable emission of a blazar based on a proposal by Marscher et al. (2010) to explain such events. In this model, a plasmoid ("blob") propagates relativistically along the spine of a blazar jet and passes through a synchrotron-emitting ring of electrons representing a shocked portion of the jet sheath. This ring supplies a source of seed photons that are inverse-Compton scattered by the electrons in the moving blob. The model includes the effects of radiative cooling, a spatially varying magnetic field, and acceleration of the blob's bulk velocity. Synthetic light curves produced by our model are compared to the observed light curves from an orphan flare that was coincident with the passage of a superluminal knot through the inner jet of the blazar PKS 1510$-$089. In addition, we present Very Long Baseline Array polarimetric observations that point to the existence of a jet sheath in PKS 1510$-$089, thus providing further observational support for the plausibility of our model. An estimate of the bolometric luminosity of the sheath within PKS 1510$-$089 is made, yielding $L_{\rm sh} \sim 3 \times 10^{45} ~ \rm erg ~ \rm s^{-1}$. This indicates that the sheath within PKS 1510$-$089 is potentially a very important source of seed photons.

SDSS J14584479+3720215: A BENCHMARKJHKSBLAZAR LIGHT CURVE FROM THE 2MASS CALIBRATION SCANS

Active galactic nuclei (AGNs) are well-known to exhibit flux variability across a wide range of wavelength regimes, but the precise origin of the variability at different wavelengths remains unclear. To investigate the relatively unexplored near-IR variability of the most luminous AGNs, we conduct a search for variability using well sampled JHKs-band light curves from the 2MASS survey calibration fields. Our sample includes 27 known quasars with an average of 924 epochs of observation over three years, as well as one spectroscopically confirmed blazar (SDSSJ14584479+3720215) with 1972 epochs of data. This is the best-sampled NIR photometric blazar light curve to date, and it exhibits correlated, stochastic variability that we characterize with continuous auto-regressive moving average (CARMA) models. None of the other 26 known quasars had detectable variability in the 2MASS bands above the photometric uncertainty. A blind search of the 2MASS calibration field light curves for AGN candidates based on fitting CARMA(1,0) models (damped-random walk) uncovered only 7 candidates. All 7 were young stellar objects within the {\rho} Ophiuchus star forming region, five with previous X-ray detections. A significant {\gamma}-ray detection (5{\sigma}) for the known blazar using 4.5 years of Fermi photon data is also found. We suggest that strong NIR variability of blazars, such as seen for SDSSJ14584479+3720215, can be used as an efficient method of identifying previously-unidentified {\gamma}-ray blazars, with low contamination from other AGN.

ANTARES constrains a blazar origin of two IceCube PeV neutrino events

The source(s) of the neutrino excess reported by the IceCube Collaboration is unknown. The TANAMI Collaboration recently reported on the multiwavelength emission of six bright, variable blazars which are positionally coincident with two of the most energetic IceCube events. Objects like these are prime candidates to be the source of the highest-energy cosmic rays, and thus of associated neutrino emission. We present an analysis of neutrino emission from the six blazars using observations with the ANTARES neutrino telescope.The standard methods of the ANTARES candidate list search are applied to six years of data to search for an excess of muons --- and hence their neutrino progenitors --- from the directions of the six blazars described by the TANAMI Collaboration, and which are possibly associated with two IceCube events. Monte Carlo simulations of the detector response to both signal and background particle fluxes are used to estimate the sensitivity of this analysis for different possible source neutrino spectra. A maximum-likelihood approach, using the reconstructed energies and arrival directions of through-going muons, is used to identify events with properties consistent with a blazar origin.Both blazars predicted to be the most neutrino-bright in the TANAMI sample (1653$-$329 and 1714$-$336) have a signal flux fitted by the likelihood analysis corresponding to approximately one event. This observation is consistent with the blazar-origin hypothesis of the IceCube event IC14 for a broad range of blazar spectra, although an atmospheric origin cannot be excluded. No ANTARES events are observed from any of the other four blazars, including the three associated with IceCube event IC20. This excludes at a 90\% confidence level the possibility that this event was produced by these blazars unless the neutrino spectrum is flatter than $-2.4$.

A SEARCH FOR FAST X-RAY VARIABILITY FROM ACTIVE GALACTIC NUCLEI USINGSWIFT

Blazars are a class of active galactic nuclei (AGNs) known for their very rapid variabilty in the high energy regions of the electromagnetic spectrum. Despite this known fast variability, X-ray observations have generally not revealed variability in blazars with rate doubling or halving timescales less than approximately 15 min. Since its launch, the Swift X-ray Telescope has obtained 0.2-10 keV X-ray data on 143 AGNs, including blazars, through intense target of opportunity observations that can be analyzed in a multiwavelength context and used to model jet parameters, particularly during flare states. We have analyzed this broad Swift data set in a search for short timescale variability in blazars that could limit the size of the emission region in the blazar jet. While we do find several low-significance possible flares with potential indications of rapid variability, we find no strong evidence for rapid ($<$15 minutes) doubling or halving times in flares in the soft X-ray energy band for the AGNs analyzed.

An active state of the BL Lacertae object Markarian 421 detected by INTEGRAL in April 2013

Multiwavelength variability of blazars offers indirect insight into their powerful engines and on the mechanisms through which energy is propagated from the centre down the jet. The BL Lac object Mkn 421 is a TeV emitter, a bright blazar at all wavelengths, and therefore an excellent target for variability studies. Mkn 421 was observed by INTEGRAL and Fermi-LAT in an active state on 16-21 April 2013. Well sampled optical, soft, and hard X-ray light curves show the presence of two flares. The average flux in the 20-100 keV range is 9.1e-11 erg/s/cm2 (~4.5 mCrab) and the nuclear average apparent magnitude, corrected for Galactic extinction, is V ~12.2. In the time-resolved X-ray spectra (3.5-60 keV), which are described by broken power laws and, marginally better, by log-parabolic laws, we see a hardening that correlates with flux increase, as expected in refreshed energy injections in a population of electrons that later cool via synchrotron radiation. The hardness ratios between the JEM-X fluxes in two different bands and between the JEM-X and IBIS/ISGRI fluxes confirm this trend. During the observation, the variability level increases monotonically from the optical to the hard X-rays, while the large LAT errors do not allow a significant assessment of the MeV-GeV variability. The cross-correlation analysis during the onset of the most prominent flare suggests a monotonically increasing delay of the lower frequency emission with respect to that at higher frequency, with a maximum time-lag of about 70 minutes, that is however not well constrained. The spectral energy distributions from the optical to the TeV domain are satisfactorily described by homogeneous models of blazar emission based on synchrotron radiation and synchrotron self-Compton scattering, except in the state corresponding to the LAT softest spectrum and highest flux.

Time correlation between the radio and gamma-ray activity in blazars and the production site of the gamma-ray emission

In order to determine the location of the gamma-ray emission site in blazars, we investigate the time-domain relationship between their radio and gamma-ray emission. Light-curves for the brightest detected blazars from the first 3 years of the mission of the Fermi Gamma-ray Space Telescope are cross-correlated with 4 years of 15GHz observations from the OVRO 40-m monitoring program. The large sample and long light-curve duration enable us to carry out a statistically robust analysis of the significance of the cross-correlations, which is investigated using Monte Carlo simulations including the uneven sampling and noise properties of the light-curves. Modeling the light-curves as red noise processes with power-law power spectral densities, we find that only one of 41 sources with high quality data in both bands shows correlations with significance larger than 3-sigma (AO 0235+164), with only two more larger than even 2.25-sigma (PKS 1502+106 and B2 2308+34). Additionally, we find correlated variability in Mrk 421 when including a strong flare that occurred in July-September 2012. These results demonstrate very clearly the difficulty of measuring statistically robust multiwavelength correlations and the care needed when comparing light-curves even when many years of data are used. This should be a caution. In all four sources the radio variations lag the gamma-ray variations, suggesting that the gamma-ray emission originates upstream of the radio emission. Continuous simultaneous monitoring over a longer time period is required to obtain high significance levels in cross-correlations between gamma-ray and radio variability in most blazars.

Models for Very Rapid High-Energy γ-Ray Variability in Blazars

We present a family of models for rapid γ-ray variability in blazars based on a two-component jet. Fast variability occurs when relativistic electron–positron pairs interact with small-scale perturbations in the inner jet. These perturbations are produced by Kelvin–Helmholtz instabilities. We fit the minute-scale strong variability displayed by the blazar PKS 2155−304 and present synthetic light-curves of the kind that might be observed with forthcoming instruments such as the Cherenkov telescope array.

Symmetric and triangle-shaped variability of blazars

AbstractSymmetric and triangle-shaped flux variability in X-ray and gamma-ray light curves is observed from many blazars. We derived the X-ray spectrum changing in time by using a kinetic equation of high energy electrons. Giving linearly changing the injection of low energy electrons into accelerating and emitting region, we obtained the preliminary results that represent the characteristic X-ray variability of the linear flux increase with hardening in the rise phase and the linear decrease with softening in the decay phase.

Detection of significant cm to sub-mm band radio and -ray correlated variability in Fermi bright blazars

The exact location of the gamma-ray emitting region in blazars is still controversial. In order to attack this problem we present first results of a cross-correlation analysis between radio (11 cm to 0.8 mm wavelength, F-GAMMA program) and gamma-ray (0.1-300 GeV) ~ 3.5 year light curves of 54 Fermi-bright blazars. We perform a source stacking analysis and estimate significances and chance correlations using mixed source correlations. Our results reveal: (i) the first highly significant multi-band radio and gamma-ray correlations (radio lagging gamma rays) when averaging over the whole sample, (ii) average time delays (source frame: 76+/-23 to 7+/-9 days), systematically decreasing from cm to mm/sub-mm bands with a frequency dependence tau_r,gamma (nu) ~ nu^-1, in good agreement with jet opacity dominated by synchrotron self-absorption, (iii) a bulk gamma-ray production region typically located within/upstream of the 3 mm core region (tau_3mm,gamma=12+/-8 days), (iv) mean distances between the region of gamma-ray peak emission and the radio "tau=1 photosphere" decreasing from 9.8+/-3.0 pc (11 cm) to 0.9+/-1.1 pc (2 mm) and 1.4+/-0.8 pc (0.8 mm), (v) 3 mm/gamma-ray correlations in 9 individual sources at a significance level where one is expected by chance (probability: 4 times 10^-6), (vi) opacity and "time lag core shift" estimates for quasar 3C 454.3 providing a lower limit for the distance of the bulk gamma-ray production region from the supermassive black hole (SMBH) of ~ 0.8-1.6 pc, i.e. at the outer edge of the Broad Line Region (BLR) or beyond. A 3 mm tau=1 surface at ~ 2-3 pc from the jet-base (i.e. well outside the "canonical BLR") finally suggests that BLR material extends to several pc distances from the SMBH.

STOCHASTIC MODELING OF THEFERMI/LAT γ-RAY BLAZAR VARIABILITY

We study the gamma-ray variability of 13 blazars observed with the Fermi Large Area Telescope (LAT). These blazars have the most complete light curves collected during the first 4 years of the Fermi sky survey. We model them with the Ornstein-Uhlenbeck (OU) process or a mixture of the OU processes. The OU process has power spectral density (PSD) proportional to 1/f^alpha with alpha changing at a characteristic time scale, tau_0, from 0 (tau>>tau_0) to 2 (tau<<tau_0). The PSD of the mixed OU process has two characteristic time scales and an additional intermediate region with 0<alpha<2. We show that the OU model provides a good description of the Fermi/LAT light curves of three blazars in our sample. For the first time we constrain a characteristic gamma-ray time scale of variability in two BL Lac sources, 3C 66A and PKS 2155-304 (tau_0=25 day and tau_0=43 day, respectively, in the observer's frame), which are longer than the soft X-ray time scales detected in blazars and Seyfert galaxies. We find that the mixed OU process approximates the light curves of the remaining 10 blazars better than the OU process. We derive limits on their long and short characteristic time scales, and infer that their Fermi/LAT PSDs resemble a power-law function. We constrain the PSD slopes for all but one source in the sample. We find hints for sub-hour Fermi/LAT variability in four flat spectrum radio quasars. We discuss the implications of our results for theoretical models of blazar variability.

SEED PHOTON FIELDS OF BLAZARS IN THE INTERNAL SHOCK SCENARIO

We extend our approach of modeling spectral energy distribution (SED) and lightcurves of blazars to include external Compton (EC) emission due to inverse Compton scattering of an external anisotropic target radiation field. We describe the time-dependent impact of such seed photon fields on the evolution of multifrequency emission and spectral variability of blazars using a multi-zone time-dependent leptonic jet model, with radiation feedback, in the internal shock model scenario. We calculate accurate EC-scattered high-energy spectra produced by relativistic electrons throughout the Thomson and Klein-Nishina regimes. We explore the effects of varying the contribution of (1) a thermal Shakura-Sunyaev accretion disk, (2) a spherically symmetric shell of broad-line clouds, the broad line region (BLR), and (3) a hot infrared emitting dusty torus (DT), on the resultant seed photon fields. We let the system evolve to beyond the BLR and within the DT and study the manifestation of the varying target photon fields on the simulated SED and lightcurves of a typical blazar. The calculations of broadband spectra include effects of gamma-gamma absorption as gamma-rays propagate through the photon pool present inside the jet due to synchrotron and inverse Compton processes, but neglect gamma-gamma absorption by the BLR and DT photon fields outside the jet. Thus, our account of gamma-gamma absorption is a lower limit to this effect. Here, we focus on studying the impact of parameters relevant for EC processes on high-energy (HE) emission of blazars.

Connection between optical and γ-ray variability in blazars

We use optical data from the Palomar Transient Factory (PTF) and the Catalina Real-Time Transient Survey (CRTS) to study the variability of gamma-ray detected and non-detected objects in a large population of active galactic nuclei (AGN) selected from the Candidate Gamma-Ray Blazar Survey and Fermi Gamma-Ray Space Telescope catalogs. Our samples include 714 sources with PTF data and 1244 sources with CRTS data. We calculate the intrinsic modulation index to quantify the optical variability amplitude in these samples. We find the gamma-ray detected objects to be more variable than the non-detected ones. The flat spectrum radio quasars (FSRQs) are more variable than the BL Lac objects in our sample, but the significance of the difference depends on the sample used. When dividing the objects based on their synchrotron peak frequency, we find the low synchrotron peaked (LSP) objects to be significantly more variable than the high synchrotron peaked (HSP) ones, explaining the difference between the FSRQs and BL Lacs. This could be due to the LSPs being observed near their electron energy peak, while in the HSPs the emission is caused by lower energy electrons, which cool more slowly. We also find a significant correlation between the optical and gamma-ray fluxes that is stronger in the HSP BL Lacs than in the FSRQs. The FSRQs in our sample are also more Compton dominated than the HSP BL Lacs. These findings are consistent with models where the gamma-ray emission of HSP objects is produced by the synchrotron self-Compton mechanism, while the LSP objects need an additional external Compton component that increases the scatter in the flux-flux correlation.

A THREE-YEAR MULTI-WAVELENGTH STUDY OF THE VERY-HIGH-ENERGY γ-RAY BLAZAR 1ES 0229+200

The high-frequency-peaked BL Lacertae object 1ES 0229+200 is a relatively distant (z = 0.1396), hard-spectrum (Gamma ~ 2.5), very-high-energy-emitting (E > 100 GeV) gamma-ray blazar. Very-high-energy measurements of this active galactic nucleus have been used to place constraints on the intensity of the extragalactic background light and the intergalactic magnetic field. A multi-wavelength study of this object centered around very-high-energy observations by VERITAS is presented. This study obtained, over a period of three years, an 11.7 standard deviation detection and an average integral flux F(E>300 GeV) = (23.3 +- 2.8_stat +- 5.8_sys) x 10^-9 photons m^-2 s^-1, or 1.7% of the Crab Nebula's flux (assuming the Crab Nebula spectrum measured by H.E.S.S). Supporting observations from Swift and RXTE are analyzed. The Swift observations are combined with previously published Fermi observations and the very-high-energy measurements to produce an overall spectral energy distribution which is then modeled assuming one-zone synchrotron-self-Compton emission. The chi^2 probability of the TeV flux being constant is 1.6%. This, when considered in combination with measured variability in the X-ray band, and the demonstrated variability of many TeV blazars, suggests that the use of blazars such as 1ES 0229+200 for intergalactic magnetic field studies may not be straightforward and challenges models that attribute hard TeV spectra to secondary gamma-ray production along the line of sight.

The uncorrelated long-term γ-ray and X-ray variability of blazars and its implications on disc–jet coupling

We examine the long term (~10 years) gamma-ray variability of blazars observed by EGRET and Fermi and find that for six sources the average flux varied by more than an order of magnitude. For two of these sources (PKS 0208-512 and PKS 0528+134), there were extensive observations (at various observing periods) by EGRET. Hence these dramatic variations are not due to a single short time-scale flare, but reflect long term changes in the average flux. Over the last twenty years, these two sources were also the target of several X-ray observatories (e.g. ROSAT, ASCA, RXTE, BeppoSAX, Chandra, Suzaku, XMM-Newton and Swift). While the ratios of the average gamma-ray fluxes between EGRET and Fermi observations are 22.9 +/- 1.9 and 12.6 +/- 1.5, their estimated 2-10 keV X-ray flux do not show such dramatic variations. The X-ray emission from such flat spectrum radio quasars (FSRQs) are believed to be due to synchrotron self Compton, while gamma-rays originate from inverse Comptonization of external soft photons from an accretion disk and/or broad line region. We argue that in this scenario, the only explanation for the uncorrelated variability is that there was an order of magnitude decrease in the external soft photons, while the jet parameters remained more or less constant. This result indicates that perhaps the accretion and jet processes are not tightly coupled in these sources.

Numerical study of broadband spectra caused by internal shocks in magnetized relativistic jets of blazars

The internal-shocks scenario in relativistic jets has been used to explain the variability of blazars’ outflow emission. Recent simulations have shown that the magnetic field alters the dynamics of these shocks producing a whole zoo of spectral energy density patterns. However, the role played by magnetization in such high-energy emission is still not entirely understood. With the aid of Fermi’s second LAT AGN catalog, a comparison with observations in the γ-ray band was performed, in order to identify the effects of the magnetic field.

The minijets-in-a-jet statistical model and the rms-flux correlation

The flux variability of blazars at very high energies does not have a clear origin. Flux variations on time scales down to the minute suggest that variability originates in the jet, where a relativistic boost can shorten the observed time scale, while the linear relation between the flux and its RMS or the skewness of the flux distribution suggests that the variability stems from multiplicative processes, which are associated in some models with the accretion disk. We study the RMS-flux relation and emphasize its link to Pareto distributions, characterized by a power-law probability density function. Such distributions are naturally generated within a minijets- in-a-jet statistical model, in which boosted emitting regions are isotropically oriented within the bulk relativistic flow of a jet. We prove that, within this model, the flux of a single minijet is proportional to its RMS. This relation still holds when considering a large number of emitting regions, for which the distribution of the total flux is skewed and could be interpreted as being log-normal. The minijets-in-a-jet statistical model reconciles the fast variations and the statistical properties of the flux of blazars at very high energies.

CHARACTERIZING THE OPTICAL VARIABILITY OF BRIGHT BLAZARS: VARIABILITY-BASED SELECTION OFFERMIACTIVE GALACTIC NUCLEI

We investigate the use of optical photometric variability to select and identify blazars in large-scale time-domain surveys, in part to aid in the identification of blazar counterparts to the ~30% of gamma-ray sources in the Fermi 2FGL catalog still lacking reliable associations. Using data from the optical LINEAR asteroid survey, we characterize the optical variability of blazars by fitting a damped random walk model to individual light curves with two main model parameters, the characteristic timescales of variability (tau), and driving amplitudes on short timescales (sigma). Imposing cuts on minimum tau and sigma allows for blazar selection with high efficiency E and completeness C. To test the efficacy of this approach, we apply this method to optically variable LINEAR objects that fall within the several-arcminute error ellipses of gamma-ray sources in the Fermi 2FGL catalog. Despite the extreme stellar contamination at the shallow depth of the LINEAR survey, we are able to recover previously-associated optical counterparts to Fermi AGN with E > 88% and C = 88% in Fermi 95% confidence error ellipses having semimajor axis r < 8'. We find that the suggested radio counterpart to Fermi source 2FGL J1649.6+5238 has optical variability consistent with other gamma-ray blazars, and is likely to be the gamma-ray source. Our results suggest that the variability of the non-thermal jet emission in blazars is stochastic in nature, with unique variability properties due to the effects of relativistic beaming. After correcting for beaming, we estimate that the characteristic timescale of blazar variability is ~3 years in the rest-frame of the jet, in contrast with the ~320 day disk flux timescale observed in quasars. The variability-based selection method presented will be useful for blazar identification in time-domain optical surveys, and is also a probe of jet physics.