  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.
