Gaia astrometry for stars with too few observations - a Bayesian approach

TL;DR

This paper introduces a Bayesian method to estimate star positions and uncertainties in Gaia data when observations are too few for full five-parameter solutions, improving error estimates for under-observed stars.

Contribution

It develops a Bayesian approach incorporating prior information to improve astrometric solutions for stars with limited observations in Gaia data.

Findings

Bayesian method provides realistic error estimates for under-observed stars

Optimal priors depend on magnitude and galactic coordinates

Method demonstrated successfully on simulated Gaia data

Abstract

Gaia's astrometric solution aims to determine at least five parameters for each star, together with appropriate estimates of their uncertainties and correlations. This requires at least five distinct observations per star. In the early data reductions the number of observations may be insufficient for a five-parameter solution, and even after the full mission many stars will remain under-observed, including faint stars at the detection limit and transient objects. In such cases it is reasonable to determine only the two position parameters. Their formal uncertainties would however grossly underestimate the actual errors, due to the neglected parallax and proper motion. We aim to develop a recipe to calculate sensible formal uncertainties that can be used in all cases of under-observed stars. Prior information about the typical ranges of stellar parallaxes and proper motions is incorporated in the astrometric solution by means of Bayes' rule. Numerical simulations based on the Gaia Universe Model Snapshot (GUMS) are used to investigate how the prior influences the actual errors and formal uncertainties when different amounts of Gaia observations are available. We develop a criterion for the optimum choice of priors, apply it to a wide range of cases, and derive a global approximation of the optimum prior as a function of magnitude and galactic coordinates. The feasibility of the Bayesian approach is demonstrated through global astrometric solutions of simulated Gaia observations. With an appropriate prior it is possible to derive sensible positions with realistic error estimates for any number of available observations. Even though this recipe works also for well-observed stars it should not be used where a good five-parameter astrometric solution can be obtained without a prior. Parallaxes and proper motions from a solution using priors are always biased and should not be used.