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 (12)

• 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
• The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background
• Inferring probabilistic stellar rotation periods using Gaussian processes
• New advances in the Gaussian-process approach to pulsar-timing data analysis
• 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
• Serendipitous discovery of radio flaring behaviour from a nearby M dwarf with MeerKAT
• Detecting the periodicity of highly irregularly sampled light-curves with Gaussian processes: the case of SDSSJ025214.67-002813.7
• 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
• Using Gaussian Processes to detect AGN flares