Blazar 3C Research Articles

3C 66A: Variability in 2007–2015

The results of photometric (BV RIJHK) and polarimetric (R)monitoring of the blazar 3C 66A performed at the St. Petersburg State University and the Central AstronomicalObservatory of the Russian Academy of Sciences in 2007–2015, radio observations performed by the Boston University team with the Very Long Baseline Array at 43 GHz, and a gamma-ray light curve based on observations with the Fermi SpaceObservatory are presented. Color variations of the object are studied. Changes in the optical spectral energy distribution are observed at some times, indicating the appearance and disappearance of individual variable sources. A variable source with a degree of polarization of 36% is identified, which is responsible for the polarization variations observed during one episode. The correlations between the variations in the different spectral ranges indicate that the optical and gamma-ray radiation originates near the radio core detected at 43 GHz. The presence of five superluminal components emerging from the core is detected.

Clues on high-energy emission mechanism from blazar 3C 454.3 during 2015 August flare

We perform a detailed spectral study of a recent flaring activity from the Flat Spectrum Radio Quasar (FSRQ), 3C\,454.3, observed simultaneously in optical, UV, X-ray and $\gamma$-ray energies during 16 to 28 August, 2015. The source reached its peak $\gamma$-ray flux of $(1.9\pm0.2)\times\,10^{-05} \; {\rm ph\,cm^{-2}\,s^{-1}}$ on 22 August. The time averaged broadband spectral energy distribution (SED) is obtained for three time periods, namely "flaring state"; covering the peak $\gamma$-ray flux, "post flaring state"; immediately following the peak flare and "quiescent state"; separated from the flaring event and following the post flaring state. The SED corresponding to the flaring state is investigated using different emission models involving synchrotron, synchrotron self Compton (SSC) and external Compton (EC) mechanisms. Our study suggests that the X-ray and $\gamma$-ray emission from 3C\,454.3 cannot be attributed to a single emission mechanism and instead, one needs to consider both SSC and EC mechanisms. Moreover, the target photon energy responsible for the EC process corresponds to an equivalent temperature of 564 K, suggesting that the flare location lies beyond the broad line emitting region of the FSRQ. SED fitting of the other two flux states further supports these inferences.

A connection between γ-ray and parsec-scale radio flares in the blazar 3C 273

We present a comprehensive 5-43 GHz VLBA study of the blazar 3C 273 initiated after an onset of a strong $\gamma$-ray flare in this source. We have analyzed the kinematics of new-born components, light curves, and position of the apparent core to pinpoint the location of the $\gamma$-ray emission. Estimated location of the $\gamma$-ray emission zone is close to the jet apex, 2 pc to 7 pc upstream from the observed 7 mm core. This is supported by ejection of a new component. The apparent core position was found to be inversely proportional to frequency. The brightness temperature in the 7 mm core reached values up to at least $10^{13}$ K during the flare. This supports the dominance of particle energy density over that of magnetic field in the 7 mm core. Particle density increased during the radio flare at the apparent jet base, affecting synchrotron opacity. This manifested itself as an apparent core shuttle along the jet during the 7 mm flare. It is also shown that a region where optical depth decreases from $\tau\sim1$ to $\tau<<1$ spans over several parsecs along the jet. The jet bulk flow speed estimated at the level of 12c on the basis of time lags between 7 mm light curves of stationary jet features is 1.5 times higher than that derived from VLBI apparent kinematics analysis.

OUP accepted manuscript

We present results based on the observations of the blazar 3C 66A from 2005 November 06 to 2016 February 14 in the BVR and I broadbands using 1.2m telescope of the Mt. Abu InfraRed Observatory (MIRO). The source was observed on 160 nights out of which on 89 nights it was monitored for more than 1 hr to check for the presence of any intra-day variability (IDV). The blazar 3C 66A exhibited significant variations in the optical flux on short and long term time scales. However, unlike highly variable S5 0716+71, it showed IDV duty cycle of about 8% only. Our statistical studies suggest the IDV time scales ranging from $\sim$ 37 min to about 3.12 hours and, in one case, a possibility of the quasi-periodic variations with characteristic timescale of $\sim$ 1.4 hrs. The IDV amplitudes in R$-$band were found to vary from 0.03 mag to as much as 0.6 mag, with larger amplitude of variation when the source was relatively fainter. The typical rate of the flux variation was estimated to be $\sim$0.07 mag hr$^{-1}$ in both, the rising and the falling phases. However, the rates of the brightness variation as high as 1.38 mag/hr were also detected. The shortest timescale of the variation, 37 min, sets an upper limit of $6.92 \times 10^{14}$ cm on the size of the emission region and about $3.7 \times 10^8\ \mathrm{M}_{\odot}$ as an estimate of the mass of the black hole, assuming the origin of the rapid optical variability is in close vicinity of the central SMBH. The long-term study suggests a mild bluer-when-brighter behavior, typical for BL Lacs.

Multifrequency Study of the Blazar 3C 454.3

This work is devoted to multi-frequency studying of the blazar 3C 454.3. The study includes spectroscopic and photometric observations in the optical, IR, and gamma-ray bands. We investigate whether a correlation exists in the light curves at different wavelengths. We have carried out observations of the optical spectrum (from 4000 to 7000 Angstroms) between 2007 and 2009, and identified MgII [2800 Angstroms] and FeII emission line features. We have obtained optical magnitudes and color indices of the quasar and performed a correlation between the optical, IR, and gamma-ray light curves. We have found statistically significant correlations between the light curves at different wavelengths.

Physical Conditions and Variability Processes in AGN Jets through Multi-Frequency Linear and Circular Radio Polarization Monitoring

Radio polarimetry is an invaluable tool to investigate the physical conditions and variability processes in active galactic nuclei (AGN) jets. However, detecting their linear and circular polarization properties is a challenging endeavor due to their low levels and possible depolarization effects. We have developed an end-to-end data analysis methodology to recover the polarization properties of unresolved sources with high accuracy. It has been applied to recover the linear and circular polarization of 87 AGNs measured by the F-GAMMA program from July 2010 to January 2015 with a mean cadence of 1.3 months. Their linear polarization was recovered at four frequencies between 2.64 and 10.45 GHz and the circular polarization at 4.85 and 8.35 GHz. The physical conditions required to reproduce the observed polarization properties and the processes which induce their variability were investigated with a full-Stokes radiative transfer code which emulates the synchrotron emission of modeled jets. The model was used to investigate the conditions needed to reproduce the observed polarization behavior for the blazar 3C 454.3, assuming that the observed variability is attributed to evolving internal shocks propagating downstream.

HIGH-ENERGY NEUTRINOS FROM RECENT BLAZAR FLARES

ABSTRACT The energy density of cosmic neutrinos measured by IceCube matches the one observed by Fermi in extragalactic photons that predominantly originate in blazars. This has inspired attempts to match Fermi sources with IceCube neutrinos. A spatial association combined with a coincidence in time with a flaring source may represent a smoking gun for the origin of the IceCube flux. In 2015 June, the Fermi Large Area Telescope observed an intense flare from blazar 3C 279 that exceeded the steady flux of the source by a factor of 40 for the duration of a day. We show that IceCube is likely to observe neutrinos, if indeed hadronic in origin, in data that are still blinded at this time. We also discuss other opportunities for coincident observations that include a recent flare from blazar 1ES 1959+650 that previously produced an intriguing coincidence with AMANDA observations.

Multiwavelength temporal variability of the blazar 3C 454.3 during 2014 activity phase

We present a multi-wavelength temporal analysis of the blazar 3C 454.3 during the high $\gamma$-ray active period from May-December, 2014. Except for X-rays, the period is well sampled at near-infrared (NIR)-optical by the \emph{SMARTS} facility and the source is detected continuously on daily timescale in the \emph{Fermi}-LAT $\gamma$-ray band. The source exhibits diverse levels of variability with many flaring/active states in the continuously sampled $\gamma$-ray light curve which are also reflected in the NIR-optical light curves and the sparsely sampled X-ray light curve by the \emph{Swift}-XRT. Multi-band correlation analysis of this continuous segment during different activity periods shows a change of state from no lags between IR and $\gamma$-ray, optical and $\gamma$-ray, and IR and optical to a state where $\gamma$-ray lags the IR/optical by $\sim$3 days. The results are consistent with the previous studies of the same during various $\gamma$-ray flaring and active episodes of the source. This consistency, in turn, suggests an extended localized emission region with almost similar conditions during various $\gamma$-ray activity states. On the other hand, the delay of $\gamma$-ray with respect to IR/optical and a trend similar to IR/optical in X-rays along with strong broadband correlations favor magnetic field related origin with X-ray and $\gamma$-ray being inverse Comptonized of IR/optical photons and external radiation field, respectively.

Exploring the nature of the broadband variability in the flat spectrum radio quasar 3C 273

Detailed investigation of broadband flux variability in the blazar 3C 273 allows us to probe the location and size of emission regions and their physical conditions. We report the results on correlation studies of the flaring activity observed between 2008 and 2012. The observed broadband variations were investigated using the structure function and the discrete correlation function, and power spectral density analysis (PSD) methods. The PSD analysis showed that the optical/IR light curve slopes are consistent with the slope of white noise processes, while, the PSD slopes at radio, X-ray and gamma-ray energies are consistent with red-noise processes. The flux variations at gamma-ray and mm-radio bands are found to be significantly correlated. Using the estimated time lag of (110\pm27) days between gamma-ray and radio light curves, we constrained the location of the gamma-ray emission region at a de-projected distance of 1.2\pm0.9 pc from the jet apex. Flux variations at X-ray bands were found to have a significant correlation with variations at both radio and \gamma-rays energies. The correlation between X-rays and gamma-rays light curves suggests presence of two components responsible for the X-ray emission. A negative time lag of -(50\pm20) days, where the X-rays are leading the emission, suggests X-rays are emitted closer to the jet apex from a compact region at a distance of ~(0.5\pm0.4) pc from the jet apex. A positive time lag of (110\pm20) days suggests jet-base origin of the other X-ray component at ~(4--5)~pc from the jet apex. The flux variations at radio frequencies were found to be well correlated with each other such that the variations at higher frequencies are leading the lower frequencies, which could be expected in the standard shock-in-jet model.

COLLISION-INDUCED MAGNETIC RECONNECTION AND A UNIFIED INTERPRETATION OF POLARIZATION PROPERTIES OF GRBs AND BLAZARS

ABSTRACT The jet composition and energy dissipation mechanism of gamma-ray bursts (GRBs) and blazars are fundamental questions that remain not fully understood. One plausible model is to interpret the γ-ray emission of GRBs and optical emission of blazars as synchrotron radiation of electrons accelerated from the collision-induced magnetic dissipation regions in Poynting-flux-dominated jets. The polarization observation is an important and independent information to test this model. Based on our recent 3D relativistic MHD simulations of collision-induced magnetic dissipation of magnetically dominated blobs, here we perform calculations of the polarization properties of the emission in the dissipation region and apply the results to model the polarization observational data of GRB prompt emission and blazar optical emission. We show that the same numerical model with different input parameters can reproduce well the observational data of both GRBs and blazars, especially the 90° polarization angle (PA) change in GRB 100826A and the 180° PA swing in blazar 3C279. This supports a unified model for GRB and blazar jets, suggesting that collision-induced magnetic reconnection is a common physical mechanism to power the relativistic jet emission from events with very different black hole masses.

A HARD GAMMA-RAY FLARE FROM 3C 279 IN 2013 DECEMBER

ABSTRACT The blazar 3C 279 exhibited twin γ-ray flares of similar intensity in 2013 December and 2014 April. In this work, we present a detailed multi-wavelength analysis of the 2013 December flaring event. Multi-frequency observations reveal the uncorrelated variability patterns with X-ray and optical–UV fluxes peaking after the γ-ray maximum. The broadband spectral energy distribution (SED) at the peak of the γ-ray activity shows a rising γ-ray spectrum but a declining optical–UV flux. This observation along with the detection of uncorrelated variability behavior rules out the one-zone leptonic emission scenario. We, therefore, adopt two independent methodologies to explain the SED: a time-dependent lepto-hadronic modeling and a two-zone leptonic radiative modeling approach. In the lepto-hadronic modeling, a distribution of electrons and protons subjected to a randomly orientated magnetic field produces synchrotron radiation. Electron synchrotron is used to explain the IR to UV emission while proton synchrotron emission is used to explain the high-energy γ-ray emission. A combination of both electron synchrotron self-Compton emission and proton synchrotron emission is used to explain the X-ray spectral break seen during the later stage of the flare. In the two-zone modeling, we assume a large emission region emitting primarily in IR to X-rays and γ-rays to come primarily from a fast-moving compact emission region. We conclude by noting that within a span of four months, 3C 279 has shown the dominance of a variety of radiative processes over each other and this reflects the complexity involved in understanding the physical properties of blazar jets in general.

MONITORING OF GAMMA-RAY BRIGHT AGN: THE MULTI-FREQUENCY POLARIZATION OF THE FLARING BLAZAR 3C 279

We present results of long-term multi-wavelength polarization observations of the powerful blazar 3C~279 after its $\gamma$-ray flare on 2013~December 20. We followed up this flare with single-dish polarization observations using two 21-m telescopes of the Korean VLBI Network. Observations carried out weekly from 2013~December~25 to 2015~January~11, at 22~GHz, 43~GHz, 86~GHz simultaneously, as part of the Monitoring Of GAmma-ray Bright AGN (MOGABA) program. We measured 3C~279 total flux densities of 22--34~Jy at 22~GHz, 15--28~Jy (43~GHz), and 10--21~Jy (86~GHz), showing mild variability of $\leq 50\,\%$ over the period of our observations. The spectral index between 22~GHz and 86~GHz ranged from $-0.13$ to $-0.36$. Linear polarization angles were 27$^{\circ}$--38$^{\circ}$, 30$^{\circ}$--42$^{\circ}$, and 33$^{\circ}$--50$^{\circ}$ at 22~GHz, 43~GHz, and 86~GHz, respectively. The degree of linear polarization was in the range of 6--12\,\%, and slightly decreased with time at all frequencies. We investigated Faraday rotation and depolarization of the polarized emission at 22--86~GHz, and found Faraday rotation measures (RM) of $-300$ to $-1200$~rad~m$^{-2}$ between 22~GHz and 43~GHz, and $-800$ to $-5100$~rad~m$^{-2}$ between 43~GHz and 86~GHz. The RM values follow a power law with a mean power law index $a$ of $2.2$, implying that the polarized emission at these frequencies travels through a Faraday screen in or near the jet. We conclude that the regions emitting polarized radio emission may be different from the region responsible for the 2013 December $\gamma$-ray flare and are maintained by the dominant magnetic field perpendicular to the direction of the radio jet at milliarcsecond scales.

The nature of the γ-ray flare associated with blazar 3C 454.3

We have studied the simultaneous spectral energy distributions (SEDs) of the 2009 December flare and those of the quiescent state of blazar 3C 454.3 by constructing a multi-component model. We find that all six SEDs can be explained by a one-zone leptonic model involving synchrotron self-Compton (SSC) plus external Compton emission from an accretion disk (ECD) and that from a broad-line region (ECC). X-ray emission is dominated by the SSC mechanism, and the γ-ray spectrum is well represented by a combination of ECD and ECC processes. Our results indicate that the energy density of the magnetic field and electrons decrease with distance from the central engine, and the Doppler factor increases with the blob moving outward in the development of the 2009 December flare. The increase in the observed flux density is possibly due to the increase in the Doppler factor of the blob. The relation between the Doppler factor δb and the distance from the central black hole suggests the magnetically driven jets span a sub-pc scale, and the relation between the magnetic field B' and the dimension of the emission region R'b is in good agreement with what is required by conservation of magnetic flux. The weak “harder-when-brighter” behavior of the γ-ray spectrum could be a result of the increase in Doppler factor during the outward motion of the blob. The parameters during the quiescent state obviously deviate from those during the flare state. We propose that the flare was likely caused by the ejection of a new blob. The gamma-ray emissions in different states are associated with the evolution of the blob.

First 230 GHz VLBI fringes on 3C 279 using the APEX Telescope

We report about a 230 GHz very long baseline interferometry (VLBI) fringe finder observation of blazar 3C 279 with the APEX telescope in Chile, the phased submillimeter array (SMA), and the SMT of the Arizona Radio Observatory (ARO). We installed VLBI equipment and measured the APEX station position to 1 cm accuracy (1 sigma). We then observed 3C 279 on 2012 May 7 in a 5 hour 230 GHz VLBI track with baseline lengths of 2800 M$\lambda$ to 7200 M$\lambda$ and a finest fringe spacing of 28.6 micro-arcseconds. Fringes were detected on all baselines with SNRs of 12 to 55 in 420 s. The correlated flux density on the longest baseline was ~0.3 Jy/beam, out of a total flux density of 19.8 Jy. Visibility data suggest an emission region <38 uas in size, and at least two components, possibly polarized. We find a lower limit of the brightness temperature of the inner jet region of about 10^10 K. Lastly, we find an upper limit of 20% on the linear polarization fraction at a fringe spacing of ~38 uas. With APEX the angular resolution of 230 GHz VLBI improves to 28.6 uas. This allows one to resolve the last-photon ring around the Galactic Center black hole event horizon, expected to be 40 uas in diameter, and probe radio jet launching at unprecedented resolution, down to a few gravitational radii in galaxies like M 87. To probe the structure in the inner parsecs of 3C 279 in detail, follow-up observations with APEX and five other mm-VLBI stations have been conducted (March 2013) and are being analyzed.

Frequency-dependent core shifts and parameter estimation for the blazar 3C 454.3

We study the core shift effect in the parsec scale jet of the blazar 3C 454.3 using the 4.8 GHz - 36.8 GHz radio light curves obtained from three decades of continuous monitoring. From a piecewise Gaussian fit to each flare, time lags $\Delta t$ between the observation frequencies $\nu$ and spectral indices $\alpha$ based on peak amplitudes $A$ are determined. From the fit $\Delta t \propto \nu^{1/k_r}$, $k_r = 1.10 \pm 0.18$ indicating equipartition between the magnetic field energy density and the particle energy density. From the fit $A \propto \nu^\alpha$, $\alpha$ is in the range $-0.24$ to $1.52$. A mean magnetic field strength at 1 pc, $B_1 = 0.5 \pm 0.2$ G, and at the core, $B_{\rm core} = 46 \pm 16$ mG, are inferred, consistent with previous estimates. The measure of core position offset is $\Omega_{r\nu} = 6.4 \pm 2.8$ pc GHz$^{1/k_r}$ when averaged over all frequency pairs. Based on the statistical trend shown by the measured core radius $r_{\rm core}$ as a function of $\nu$, we infer that the synchrotron opacity model may not be valid for all cases. A Fourier periodogram analysis yields power law slopes in the range $-1.6$ to $-3.5$ describing the power spectral density shape and gives bend timescales in the range $0.52 - 0.66~$yr. This result, and both positive and negative $\alpha$, indicate that the flares originate from multiple shocks in a small region. Important objectives met in our study include: the demonstration of the computational efficiency and statistical basis of the piecewise Gaussian fit; consistency with previously reported results; evidence for the core shift dependence on observation frequency and its utility in jet diagnostics in the region close to the resolving limit of very long baseline interferometry observations.

RAPID VARIABILITY OF BLAZAR 3C 279 DURING FLARING STATES IN 2013−2014 WITH JOINTFERMI-LAT,NuSTAR,SWIFT, AND GROUND-BASED MULTI-WAVELENGTH OBSERVATIONS

We report the results of a multi-band observing campaign on the famous blazar 3C 279 conducted during a phase of increased activity from 2013 December to 2014 April, including first observations of it with NuSTAR. The $\gamma$-ray emission of the source measured by Fermi-LAT showed multiple distinct flares reaching the highest flux level measured in this object since the beginning of the Fermi mission, with $F(E > 100\,{\rm MeV})$ of $10^{-5}$ photons cm$^{-2}$ s$^{-1}$, and with a flux doubling time scale as short as 2 hours. The $\gamma$-ray spectrum during one of the flares was very hard, with an index of $\Gamma_\gamma = 1.7 \pm 0.1$, which is rarely seen in flat spectrum radio quasars. The lack of concurrent optical variability implies a very high Compton dominance parameter $L_\gamma/L_{\rm syn} > 300$. Two 1-day NuSTAR observations with accompanying Swift pointings were separated by 2 weeks, probing different levels of source activity. While the 0.5$-$70 keV X-ray spectrum obtained during the first pointing, and fitted jointly with Swift-XRT is well-described by a simple power law, the second joint observation showed an unusual spectral structure: the spectrum softens by $\Delta\Gamma_{\rm X} \simeq 0.4$ at $\sim$4 keV. Modeling the broad-band SED during this flare with the standard synchrotron plus inverse Compton model requires: (1) the location of the $\gamma$-ray emitting region is comparable with the broad line region radius, (2) a very hard electron energy distribution index $p \simeq 1$, (3) total jet power significantly exceeding the accretion disk luminosity $L_{\rm j}/L_{\rm d} \gtrsim 10$, and (4) extremely low jet magnetization with $L_{\rm B}/L_{\rm j} \lesssim 10^{-4}$. We also find that single-zone models that match the observed $\gamma$-ray and optical spectra cannot satisfactorily explain the production of X-ray emission.

Multiband variability in the blazar 3C 273 with XMM–Newton

We have undertaken a nearly simultaneous optical/UV and X-ray variability study of the flat spectrum radio quasar, 3C 273 using data available from the XMM–Newton satellite mission from 2000 June to 2012 July. Here, we focus on the multiwavelength flux variability on both intra-day and long time-scales of this very well known radio-loud source. We found high flux variability over long time-scales in all bands for which observations were made. The optical/UV variability amplitude was more than twice than that in the X-ray bands. There is some frequency dependence of the variability in optical/UV bands in the sense that the variability amplitude increases with increasing frequency; however, the X-ray emissions disagreewiththistrendasthevariabilityamplitudedecreasesfromsofttohardX-raybands.On intra-day time-scales 3C 273 showed small amplitude variability in X-ray bands. A hardness ratio analysis in the X-rayregime indicates that theparticle acceleration mechanism dominates the cooling mechanism during most of the ∼12 yr span of these observations.

POLARIZATION SWINGS REVEAL MAGNETIC ENERGY DISSIPATION IN BLAZARS

The polarization signatures of the blazar emissions are known to be highly variable. In addition to small fluctuations of the polarization angle around a mean value, sometimes large (> 180^o) polarization angle swings are observed. We suggest that such p henomena can be interpreted as arising from light-travel-time effects within an underlying axisymmetric emission region. We present the first simultaneous fitting of the multi-wavelength spectrum, variability and time-dependent polarization features of a correlated optical and gamma-ray flaring event of the prominent blazar 3C279, which was accompanied by a drastic change of its polarization signatures. This unprecedented combination of spectral, variability, and polarization information in a coherent physical model allows us to place stringent constraints on the particle acceleration and magnetic-field topology in the relativistic jet of a blazar, strongly favoring a scenario in which magnetic energy dissipation is the primary driver of the flare event.

MULTI-WAVELENGTH OBSERVATIONS OF 3C 279 DURING THE EXTREMELY BRIGHT GAMMA-RAY FLARE IN 2014 MARCH–APRIL

The well studied blazar 3C 279 underwent a giant $\gamma$-ray outburst in 2014 March-April. The measured $\gamma$-ray flux (1.21 $\pm$ 0.10 $\times$ 10$^{-5}$ ph cm$^{-2}$ s$^{-1}$ in 0.1-300 GeV energy range) is the highest detected from 3C 279 by Fermi Large Area Telescope. Hour scale $\gamma$-ray flux variability are observed, with a flux doubling time as short as 1.19 $\pm$ 0.36 hours detected during one flare. The $\gamma$-ray spectrum is found to be curved at peak of the flare suggesting low probability of detecting very high energy (VHE; E $>$ 100 GeV) emission, which is further confirmed by the Very Energetic Radiation Imaging Telescope Array System observations. The $\gamma$-ray flux increased by more than an order in comparison to low activity state and the flare consists of multiple sub-structures having fast rise and slow decay profile. The flux enhancement is seen in all the wavebands though at a lesser extent compared to $\gamma$-rays. During the flare, a considerable amount of the kinetic jet power gets converted to $\gamma$-rays and the jet becomes radiatively efficient. A one zone leptonic emission model is used to reproduce the flare and we find increase in the bulk Lorentz factor as a major cause of the outburst. From the observed fast variability, lack of VHE detection, and the curved $\gamma$-ray spectrum, we conclude that the location of the emission region cannot be far out from the broad line region (BLR) and contributions from both BLR and torus photons are required to explain the observed $\gamma$-ray spectrum.

TIME DEPENDENT HADRONIC MODELING OF FLAT SPECTRUM RADIO QUASARS

We introduce a new time-dependent lepto-hadronic model for blazar emission that takes into account the radiation emitted by secondary particles, such as pions and muons, from photo hadronic interactions. Starting from a baseline parameter set guided by a fit to the spectral energy distribution of the blazar 3C 279, we perform a parameter study to investigate the effects of perturbations of the input parameters to mimic different flaring events to study the resulting lightcurves in the optical, X-ray, high energy (HE: E > 100 MeV) and very-high-energy (VHE: E > 100 GeV) gamma-rays as well as the neutrino emission associated with charged-pion and muon decay. We find that flaring events from an increase in the efficiency of Fermi II acceleration will produce a positive correlation between all bandpasses and a marked plateau in the HE gamma-ray lightcurve. We also predict a distinctive dip in the HE lightcurve for perturbations caused by a change in the proton injection spectral index. These plateaus / dips could be a tell tale signature of hadronic models for perturbations that lead to more efficient acceleration of high energy protons in parameter regimes where pion and muon synchrotron emission is non-negligible.