New Metrics for Identifying Variables and Transients in Large Astronomical Surveys

Abstract

A key science goal of large sky surveys such as those conducted by the Vera C. Rubin Observatory and precursors to the Square Kilometre Array is the identification of variable and transient objects. One approach is the statistical analysis of the time series of the changing brightness of sources, that is, their light curves. However, finding adequate statistical representations of light curves is challenging because of data quality issues such as sparsity of observations, irregular sampling, and other nuisance factors inherent in astronomical data collection. The wide diversity of objects that a large-scale survey will observe also means that making parametric assumptions about the shape of light curves is problematic. We present a Gaussian process (GP) regression approach for characterising light curve variability that addresses these challenges. Our approach makes no assumptions about the shape of a light curve and, therefore, is general enough to detect a range of variable source types. In particular, we propose using the joint distribution of GP amplitude hyperparameters to distinguish variable and transient candidates from nominally stable ones and apply this approach to 6394 radio light curves from the ThunderKAT survey. We compare our results with two variability metrics commonly used in radio astronomy, namely ${\eta }_{\nu }$ and $V_{\nu }$, and show that our approach has better discriminatory power and interpretability. Finally, we conduct a rudimentary search for transient sources in the ThunderKAT dataset to demonstrate how our approach might be used as an initial screening tool. Computational notebooks in Python and R are available to help facilitate the deployment of this framework to other surveys.

26 pages, 13 figures

References in corpus (28)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Rank-normalization, folding, and localization: An improved $\hat{R}$ for assessing convergence of MCMC
• FINDCHIRP: an algorithm for detection of gravitational waves from inspiraling compact binaries
• Non-linear X-ray variability in X-ray binaries and active galaxies
• The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background
• Inferring probabilistic stellar rotation periods using Gaussian processes
• The Caltech-NRAO Stripe 82 Survey (CNSS) Paper I: The Pilot Radio Transient Survey In 50 deg${}^2$
• New advances in the Gaussian-process approach to pulsar-timing data analysis
• Transiting exoplanet candidates from K2 Campaigns 5 and 6
• SuperRAENN: A Semi-supervised Supernova Photometric Classification Pipeline Trained on Pan-STARRS1 Medium Deep Survey Supernovae
• A survey for transients and variables with the Murchison Widefield Array 32-tile prototype at 154 MHz
• Disentangling Time Series Spectra with Gaussian Processes: Applications to Radial Velocity Analysis
• Improving gravitational-wave parameter estimation using Gaussian process regression
• Looking at Blazar Light Curve Periodicities with Gaussian Processes
• 21 new long-term variables in the GX 339-4 field: two years of MeerKAT monitoring
• A Gaussian Process Regression Reveals No Evidence for Planets Orbiting Kapteyn's Star
• Serendipitous discovery of radio flaring behaviour from a nearby M dwarf with MeerKAT
• Density estimation with Gaussian processes for gravitational-wave posteriors
• Detecting the periodicity of highly irregularly sampled light-curves with Gaussian processes: the case of SDSSJ025214.67-002813.7
• The Murchison Widefield Array Transients Survey (MWATS). A search for low frequency variability in a bright Southern hemisphere sample
• A Matched Filter Technique For Slow Radio Transient Detection And First Demonstration With The Murchison Widefield Array
• Bursts from Space: MeerKAT - The first citizen science project dedicated to commensal radio transients
• What can Gaussian Processes really tell us about supernova lightcurves? Consequences for Type II(b) morphologies and genealogies
• Using Gaussian Processes to detect AGN flares
• Finding radio transients with anomaly detection and active learning based on volunteer classifications