Impact of basic angle variations on the parallax zero point for a scanning astrometric satellite
^{1} Pulkovo Observatory, Pulkovskoye shosse 65, 196140 SaintPetersburg, Russia
email: ag.butkevich@gmail.com
^{2} Lohrmann Observatory, Technische Universität Dresden, 01062 Dresden, Germany
email: sergei.klioner@tudresden.de
^{3} Lund Observatory, Department of Astronomy and Theoretical Physics, Lund University, Box 43, 22100 Lund, Sweden
email: lennart; david@astro.lu.se
^{4} Institute of Astronomy, Madingley Road, Cambridge CB3 OHA, UK
email: fvl@ast.cam.ac.uk
Received: 14 March 2017
Accepted: 16 April 2017
Context. Determination of absolute parallaxes by means of a scanning astrometric satellite such as Hipparcos or Gaia relies on the shortterm stability of the socalled basic angle between the two viewing directions. Uncalibrated variations of the basic angle may produce systematic errors in the computed parallaxes.
Aims. We examine the coupling between a global parallax shift and specific variations of the basic angle, namely those related to the satellite attitude with respect to the Sun.
Methods. The changes in observables produced by small perturbations of the basic angle, attitude, and parallaxes were calculated analytically. We then looked for a combination of perturbations that had no net effect on the observables.
Results. In the approximation of infinitely small fields of view, it is shown that certain perturbations of the basic angle are observationally indistinguishable from a global shift of the parallaxes. If these kinds of perturbations exist, they cannot be calibrated from the astrometric observations but will produce a global parallax bias. Numerical simulations of the astrometric solution, using both direct and iterative methods, confirm this theoretical result. For a given amplitude of the basic angle perturbation, the parallax bias is smaller for a larger basic angle and a larger solar aspect angle. In both these respects Gaia has a more favourable geometry than Hipparcos. In the case of Gaia, internal metrology is used to monitor basic angle variations. Additionally, Gaia has the advantage of detecting numerous quasars, which can be used to verify the parallax zero point.
Key words: methods: data analysis / methods: statistical / space vehicles: instruments / catalogs / astrometry / parallaxes
© ESO, 2017
1. Introduction
The European Space Agency’s space astrometry mission Gaia, aiming to determine astrometric parameters for at least one billion stars with accuracies reaching ten microarcseconds (de Bruijne 2012), was launched in December 2013 (Gaia Collaboration 2016b). Gaia is based on similar observation principles as the highly successful pioneering astrometric mission Hipparcos (ESA 1997). In particular, both satellites use an optical system providing two viewing directions separated by a wide angle, referred to as the basic angle (Perryman et al. 2001). The goal of the present paper is to show how certain timedependent variations of the basic angle can bias the parallax zero point of an astrometric solution derived from observations by a scanning satellite utilizing this sort of design.
The data processing for an astrometric satellite should, as far as possible, be based on the principle of selfcalibration (Lindegren & Bastian 2011). This means that the same observational data are used to determine both the scientifically interesting astrometric parameters and the socalled “nuisance parameters” that describe the instrument calibration, satellite attitude, and other relevant parts of the observation model (Lindegren et al. 2012). The selfcalibration is, however, of limited applicability when variations of different parameters do not produce fully independent effects in the observables. These kinds of situations can occur when the variation of certain parameters leads to changes in the observables that resemble the changes produced by the variation of some other parameters. The more similar the changes in the observables are, the stronger the correlation between the parameters. If the changes are identical, the problem of parameter estimation is degenerate, meaning that the same set of observables is equally well described by different sets of parameter values.
As long as the degeneracy involves only nuisance parameters, it has no effect on the astrometric solution and only leads to an arbitrary but unimportant shift of the respective nuisance parameters. However, if there is a degeneracy between the astrometric and nuisance parameters, the selfcalibration process, in general, leads to biased astrometry. The celestial reference frame is an example of a complete degeneracy between the astrometric and attitude parameters, which can only be lifted by means of external data, in this case the positions and proper motions of a number of extragalactic objects (Kovalevsky et al. 1997; Lindegren et al. 2012). Concerning the instrument calibration parameters, it is possible to formulate the calibration model in such a way that (near)degeneracies are avoided among its parameters, as well as between the calibration and attitude parameters.
However, it is still possible that the actual physical variations of the instrument contain components that are degenerate with the astrometric parameters. By definition, these kinds of variations cannot be detected internally by the astrometric solution, for example from an analysis of the residuals, but only through a comparison with external data, such as astrophysical information or independent metrology. In Gaia the latter is chosen, as is detailed in Sect. 4.1.
The basic angle is an important example of a quantity that can vary in a way that cannot be fully calibrated from observations. Already in the early years of the Hipparcos project it was realised that certain periodic variations of the basic angle, caused by a nonuniform heating of the satellite by the solar radiation, lead to a global shift of the parallaxes (Lindegren 1977). Subsequent analyses (Arenou et al. 1995; van Leeuwen 2005) concluded that the possible effect on the Hipparcos parallaxes was negligible, suggesting a very good shortterm stability of the basic angle in that satellite.
For Gaia the situation is different. The much higher accuracy targeted by this mission necessitates a very careful consideration of possible biases introduced by uncalibrated instrumental effects, including basic angle variations. This is even more evident in view of the very significant, ~1 mas amplitude, basic angle variations measured by the onboard metrology system of Gaia (Gaia Collaboration 2016b; Lindegren et al. 2016). In this context, the neardegeneracy between a global parallax zero point error and a possible basic angle variation induced by solar radiation is particularly relevant. The theoretical analysis of the problem presented here expands and clarifies earlier analytical results by Lindegren (1977, 2004) and van Leeuwen (2005).
An analytical treatment of the problem is given in Sect. 2. In Sect. 3 we present the results of numerical experiments that confirm the theoretical expectations. In Sect. 4 we consider the practical implications of results. Some concluding remarks are given in Sect. 5.
2. Theory
In this section we consider how small perturbations of various parameters change the observed quantities. We first demonstrate that, to first order in the small angles, arbitrary variations of observables are equivalent to certain variations of the basic angle and attitude (Sects. 2.1–2.4). Then we find the changes of observables due to a common shift of all parallaxes (Sect. 2.6). Combining these results, we derive in, Sects. 2.7–2.8, the specific variations of the basic angle and attitude that are observationally indistinguishable from a common shift of the parallaxes.
2.1. Reference system
To study the coupling between the instrument parameters and parallax, it is convenient to make use of the rotating reference system aligned with the fields of view. This system, known as the Scanning Reference System (SRS) in the Gaia nomenclature (Lindegren et al. 2012), is represented by the instrument axes x, y, z (Fig. 1), with z directed along the nominal spin axis of the satellite, x bisecting the two viewing directions separated by the basic angle Γ, and y = z × x. The direction towards an object is specified by the unit vector (1)with the instrument angles ϕ and h describing the position of the object with respect to the SRS (Fig. 1). For a star in the preceding field of view (PFoV) ϕ ≃ + Γ/2, while in the following field of view (FFoV) ϕ ≃ −Γ/2.
2.2. Field angles
An observation consists of a measurement of the coordinates of a stellar image in the focal plane at a particular time. In practice, the measurement is expressed in detector coordinates, such as pixels, but we consider here an idealised system providing a direct measurement of the two field angles g and h in the relevant field of view. Whereas the acrossscan field angle h coincides with the corresponding instrument angle, the alongscan field angle g is reckoned from the centre of the corresponding field of view in the direction of the satellite rotation (Fig. 1). Projected on the sky, the fieldofview centre defines two viewing directions separated by the basic angle Γ. Thus, (2)where subscripts p and f denote values for the preceding and following field of view, respectively. We assume that the instrument is ideal except for the basic angle Γ, which can deviate from its nominal (conventional) value Γ_{c} by a timedependent variation: (3)It is important to note that the alongscan field angle g, as defined here, is not the same as the alongscan field angle η normally used in the context of Gaia data processing (Lindegren et al. 2012). Whereas η is measured from a fixed, conventional origin at ϕ = ± Γ_{c}/ 2, our g is measured from the actual, variable field centre at ϕ = ± Γ(t)/2. This difference motivates the change in notation from η to g. For consistency, a corresponding change is made in the acrossscan direction, although our h is the same as the acrossscan field angle ζ used in the Gaia data processing.
Fig. 1 Definition of the instrument axes x, y, z of the Scanning Reference System (SRS), the basic angle Γ, and the field angles g and h specifying the observed direction to a star (u) in either field of view. ϕ is the alongscan instrument angle of the star. In the SRS the direction to the solar system barycentre, b, is specified by the angles ξ and Ω. 

Open with DEXTER 
2.3. Variations of the field angles due to a change in the basic angle
Any increase or decrease of the basic angle makes the fields of view move farther from each other or closer together. This, in turn, changes the observed field angle g for a given stellar image. However, because the attitude, that is the celestial pointing of the SRS axes, is unchanged, the value of ϕ for a given star is not affected by the basic angle. For example, in the preceding field of view, an increase of the basic angle causes the observed image to be shifted with respect to the centre of the field of view so that the observed alongscan field angle g_{p} is decreased. The opposite effect takes place in the following field of view. The acrossscan field angles h_{p} and h_{f} are obviously not affected. The variations of the field angles caused by the basicangle variation δΓ are therefore (4)This agrees with Eq. (2), taking into account that δϕ = 0.
2.4. Variations of the field angles due to a change in the attitude
A quaternion representation is used to parametrise the attitude of Gaia (Lindegren et al. 2012, Appendix A). Here it is more convenient to describe small changes in the attitude by means of three small angles δ_{x}, δ_{y}, and δ_{z} representing the rotations around the corresponding SRS axes. Because the direction u to the star is regarded here as fixed, the corresponding changes in the observed field angles are found to be (5)In these and the following equations, we neglect terms of the second and higher orders in δ_{x}, δ_{y}, δ_{z}, and δΓ. To this approximation, we can use Γ_{c} instead of Γ in the trigonometric factors.
2.5. Combined changes in the basic angle and attitude
Combining Eqs. (4) and (5) we see that a simultaneous change of the basic angle by δΓ and of the attitude by δ_{x}, δ_{y}, δ_{z} gives the following total change of the field angles: (6)An exact inversion of this system of equations gives (7)The first two equations in (7) show that arbitrary small changes in the acrossscan field angles δh_{p} and δh_{f} to first order can be represented as changes in the attitude by δ_{x} and δ_{y}. Similarly, the last two equations show that arbitrary changes in the alongscan field angles δg_{p} and δg_{f} can be represented as a combination of a change of the basic angle δΓ and a change in the attitude by δ_{z}.
In general, an arbitrary perturbation of the observed stellar positions, being a smooth function of time and stellar position, clearly results in a smooth, timedependent variation of δg_{p}, δg_{f}, δh_{p}, and δg_{f}. From Eq. (7) it follows that this kind of a perturbation is observationally indistinguishable from a certain timedependent variation of δ_{x}, δ_{y}, δ_{z}, and δΓ.
2.6. Variations of the field angles due to a change in the parallax
The position of the barycentre of the solar system with respect to the instrument can be specified by a distance R (in au) and two angular coordinates. We took the angular coordinates to be ξ and Ω defined as in Fig. 1. According to the scanning law, ξ is nearly constant while Ω is increasing with time as the satellite spins. The barycentric position of the satellite is (8)The observed direction u to a star is given by Eq. (4) of Lindegren et al. (2012) as a function of the astrometric parameters of the star. Linearisation yields the change in the direction caused by a small change of the parallax δϖ: (9)We now assume that the direction to a star is changed only from a change of its parallax, whereas the basic angle and attitude are kept constant. The fixed basic angle implies δϕ = δg. The fixed attitude means that x, y, and z are constant, so that Eq. (1) gives the change in direction (10)where (11)are unit vectors in the directions of increasing ϕ and h, respectively. They are evidently orthogonal to each other and to u. Equating δu from (9) and (10), and taking the dot product with v and w, gives (12)Substituting Eqs. (8) and (11) we find (13)Up to this point the derived formulae are valid throughout the field of view to first order in the small variations denoted with a δ. To proceed, we now consider a star observed at the centre of either field of view, so that g = h = 0 and ϕ = ± Γ_{c}/ 2. In this case the variations of the field angles caused by δϖ become (14)Considering only stars at the centre of either field of view effectively means that we neglected the finite size of the field of view. In both Hipparcos and Gaia the halfsize of the field of view is Φ < 10^{2} rad. Because  g ,  h  < Φ, neglected terms in Eq. (14) are of the order of Φ × δ, where δ represents any of the quantities δΓ, δ_{x}, etc. Equation (14) is therefore expected to be accurate to <1% at any point in the field of view. The implications of this approximation are further discussed in Sect. 4.4.
2.7. Relation between the changes in parallax, basic angle, and attitude
Substituting Eq. (14) into Eq. (7) we readily obtain a relation between the change in parallax and the corresponding changes in basic angle and attitude: (15)These equations should be interpreted as follows: a change in parallax by δϖ is observationally indistinguishable (to order Φ × δ) from a simultaneous change of the attitude by δ_{x}, δ_{y}, δ_{z} and of the basic angle by δΓ. The formulae were derived for one star, but if δϖ is the same for all stars, they hold for all observations of all stars. Equation (15) therefore defines the specific variations of the attitude and basic angle that mimic a global shift in parallax. Remarkably, the rotation around the x axis is not affected by the global parallax change, whereas the rotation around the y axis is independent of Ω and therefore, in practice, almost constant.
In a global astrometric solution, all the attitude and stellar parameters, including ϖ, are simultaneously fitted to the observations of g and h. A specific variation in the basic angle of the form δΓ(t) ∝ cosΩsinξR then leads to a global shift of the fitted parallaxes, together with some timedependent attitude errors δ_{y}, δ_{z}. Because the effects of this kind of a basic angle variation are fully absorbed by the attitude parameters and parallaxes, the variation is completely degenerate with the stellar and attitude model and cannot be detected from the residuals of the fit.
For a satellite in orbit around the Earth, such as Hipparcos, or near L_{2}, such as Gaia, R is approximately constant. To achieve a stable thermal regime of the instrument for a scanning astrometric satellite, one typically chooses a scanning law with a constant angle between the direction to the Sun and the spin axis, the socalled solar aspect angle. This means that angle ξ is nearly constant as well (see Sect. 4.5). In the next section we consider the case when R and ξ are exactly constant. Nevertheless, in reality both R and ξ are somewhat timedependent, and this case is discussed in Sect. 4.3.
Because ξ and R are nearly constant, the degenerate component of the basic angle variation is essentially of the form cosΩ, which is periodic with the satellite spin period relative to the solarsystem barycentre. This leads to a fundamental design requirement for a scanning astrometry satellite, namely that the basic angle should not have significant periodic variations with a period close to the period of rotation of the satellite, and especially not of the form cosΩ.
2.8. Harmonic representation of the variations
From Eq. (15) it is seen that a global parallax shift corresponds to variations of δΓ and δ_{z} proportional to cosΩ and sinΩ, respectively, whereas δ_{y} and δ_{x} are constant. These quantities correspond to terms of order k = 0 and 1 in the more general harmonic series (16)Specifically, if except for , we find the following relations between the amplitude of the basic angle variation, the global shift of the parallaxes, and the nonzero harmonics of the attitude errors: (17)The numerical values correspond to the mean parameters relevant for Gaia, that is Γ_{c} = 106.5°, ξ = 45°, and R = 1.01 au.
It is not the purpose of this paper to investigate the possible effects of other harmonics of the basic angle variation. However, it can be mentioned that is the only harmonic parameter that is degenerate with the attitude and stellar parameters. This means that , , and , for k> 1 can be determined as parameters in the calibration model.
3. Results of numerical simulations
In this section we present the results of numerical tests performed to check the above conclusions. To study different aspects of the problem, we made use of two different, though complimentary, solutions: a direct solution, where inversion of the normal matrix provides full covariance information; and an iterative solution, with separate updates for different groups of unknowns, similar to the method used in the actual processing of Gaia data (Lindegren et al. 2012). Although the direct solution can only handle a relatively small number of stars, and therefore is of limited practical use, it enables us to investigate important mathematical properties of the problem and to examine the correlations between all the parameters. By contrast, the iterative solution cannot provide this kind of information, but is more realistic in terms of the number of stars and has been successfully employed in the processing of real Gaia data.
3.1. Smallscale direct solutions
For the direct solutions a special simulation software was developed. It simulates the observations of a small number of stars and the reconstruction of their astrometric parameters based on conventional leastsquares fitting. The normal equations for the unknown parameters are accumulated and the normal matrix is inverted using singular value decomposition (Golub & van Loan 1996). This decomposition allowed us to study mathematical properties of the problem, especially details of its degeneracy. The simulations included 10^{4} stars uniformly distributed over the celestial sphere. Observations were generated with a harmonic perturbation of the basic angle as in Eq. (16d) with mas. To study the effects of the basic angle variations in pure form, no noise was added to the observations. The solutions always include five astrometric parameters per star: two components of the position, the parallax, and two components of the proper motion. Additional parameters representing the variations in attitude and basic angle were introduced as required by various types of solutions.
The first test was to check the theoretical predictions in Eq. (17). To this end, a solution was made including only the star (S) and attitude (A) parameters. The attitude part of the solution consisted of nine parameters: the harmonic amplitudes of δ_{x}, δ_{y}, and δ_{z} as given by the first three equations of (16) for k ≤ 1. The results of this solution, summarised in Table 1, are in very good agreement with the theoretical expectations. Small deviations from the predicted values may be caused, for example, by the limited number of stars, the finite size of the field of view, and numerical rounding errors.
Another test was the singular value analysis of the normal matrix. In addition to the solution described above (of type SA), we computed three other solutions with different sets of unknowns but always using the same observations. Solution S included only the star parameters as unknowns, solution SB included also the harmonic coefficients of δΓ in the last equation of (16) for k ≤ 1, and solution SBA included all three sets of unknowns.
The results, summarised in Table 2, again confirm the theoretical predictions. The problem is close to being degenerate only in case SBA where all three kinds of parameters, star, basic angle, and attitude, are fitted simultaneously. In this case there is one singular value, ~ 10^{5}, much smaller than the second smallest value, ~ 10^{2}. This indicates that the problem has a rank deficiency of one, which is clearly caused by the (near)degeneracy between the global parallax shift and the specific basic angle and attitude variations described by Eq. (17). In all other cases, SB, SA, and S, no isolated small singular values appear. The problem is then formally wellconditioned, although, as we saw in case SA, the solutions may be biased by the basic angle variations.
Parallax shift and attitude harmonics obtained in the smallscale direct solution (of type SA) with mas basicangle variation.
Singular values σ_{i} of the normal matrix for the different types of direct solutions.
The circumstance that the smallest singular value in case SBA is not zero is partly attributable to rounding errors in a solution involving ≥50 000 unknowns. However, even in exact arithmetics, the smallscale direct solution does not represent a fully degenerate problem because (i) the strict degeneracy only occurs if one neglects the finite size of the fields of view, and (ii) some parameters of the mission are slightly timedependent, but were assumed to be constant by the harmonic representations in Eq. (16). The question of the time dependence is further addressed in Sect. 4.3.
The attitude parameters used in the smallscale solutions are not representative of any practically useful attitude model, but were chosen solely to verify the expected degeneracy with the basic angle variation and parallax zero point. In particular, the harmonic model of δ_{x}, δ_{y}, δ_{z} in Eq. (16) cannot describe a solid rotation of the reference frame, which explains why Table 2 does not show the sixfold degeneracy between the attitude and stellar parameters normally expected from the unconstrained reference frame. This simplification is removed in the largescale simulations described below, which use a fully realistic attitude model.
3.2. Largescale iterative solutions
To test the effect of the basic angle variation in an iterative solution, we made use of the Gaia AGISLab simulation software (Holl et al. 2012, Appendix B). This tool allows us to simulate independent astrometric solutions in a reasonable time, based on the same principles as the astrometric global iterative solution (AGIS; Lindegren et al. 2012) used for Gaia but employing a smaller number of primary stars and several other timesaving simplifications.
To investigate the parallax zero point we did a set of tests for different values of the basic angle in the range from 30° to 150°, including the Hipparcos and Gaia values of 58° and 106.5°, respectively. The nominal Gaia scanning law and realistic geometry of the Gaia fields of view were assumed in these simulations, which include one million bright, G = 13, stars uniformly distributed on the sky and observed during five years with no dead time. The nominal alongscan observation noise of 95 μas per CCD was assumed, based on the estimated centroiding performance of Gaia for stars of G = 13 mag. According to de Bruijne (2012), this corresponds to an expected endofmission precision of around 10 μas for the parallaxes. Additionally, basic angle variations with an amplitude mas were included. The unknowns consisted of the five astrometric parameters per star and the attitude parameters based on a Bspline representation of the quaternion components (Lindegren et al. 2012) using a knot interval of 30 s. The AGISLab simulations therefore incorporate many detailed features from the observations and data reductions of an actual scanning astrometric mission, including a number of higherorder effects neglected in our analytical approach.
The parallaxes obtained in the iterative solutions are offset from their “true” values, assumed in the simulation, by random and systematic errors. Results for the mean offset δϖ are shown by the points in Fig. 2. The curve is the theoretical relation from the first equation in (17). The points deviate from the curve by at most 1 μas, or ≲0.1%. The random errors, measured by the sample standard deviation of the offsets of the individual parallaxes, were 7.3 μas in all the experiments, are practically independent of the basic angle and in rough agreement with the expected endofmission precision.
Fig. 2 Parallax zero point shift for a basic angle variation of the form cosΩ with amplitude 1 mas. The dots show results of the largescale iterative solution for several different basic angles, including the Hipparcos and Gaia values; the solid curve shows the theoretical relation from Eq. (17). 

Open with DEXTER 
4. Discussion
In the preceding sections we show that a global shift of the parallaxes is observationally indistinguishable from a certain time variation of the basic angle. The relevant relation, strictly valid at the centre of the field of view, is given by the last identity in Eq. (15). Here we proceed to discuss some practical implications of this result.
4.1. Physical relevance of the cos Ω dependence
Elementary design principles have led to the choice of a nearly constant solar aspect angle ξ (Fig. 1) for both Hipparcos (43°) and Gaia (45°). Moreover, for a satellite in orbit around the Earth or close to the second Lagrange (L_{2}) point of the SunEarthMoon system, the barycentric distance R is always close to 1 au. The form of basic angle variation that is degenerate with parallax is then essentially proportional to cosΩ. This result is highly significant in relation to expected thermal variations of the instrument. The oblique solar illumination of the rotating satellite may produce basic angle variations that are periodic with the spin period relative to the Sun, meaning that they are of the general harmonic form of the last line in Eq. (16). A nearly constant, nonzero coefficient could therefore be a very realistic physical consequence of the way the satellite is operated.
Knowledge of the close coupling between a possible thermal impact on the instrument and the parallax zero point has resulted in very strict engineering specifications for the acceptable amplitude of shortterm variations of the basic angle in both Hipparcos and Gaia. In the case of Gaia it was known already at an early design phase that the basic angle variations cannot be fully avoided passively and need to be measured. Therefore Gaia includes a dedicated laserinterferometric metrology system, the basic angle monitor (BAM; Mora et al. 2014), to measure the shortterm variations. According to BAM measurements, during the first year of the nominal operations of Gaia, the amplitude of the cosΩ term, referred to 1.01 au and epoch 2015.0, was about 0.848 mas (Lindegren et al. 2016). Uncorrected, this kind of a large variation would lead to a parallax bias of 0.741 mas according to Eq. (17). For Gaia Data Release 1 (Gaia Collaboration 2016a) the observations were corrected for the basic angle variations based on a harmonic model fitted to the BAM measurements.
4.2. Dependence on Γ_{c}
From Eq. (17) it is seen that the parallax shift is inversely proportional to sin(Γ_{c}/ 2). For a constant amplitude of the cosΩ term in the basic angle variation, the parallax shift therefore decreases with increasing basic angle (cf. Fig. 2). From this point of view the optimum basic angle is therefore 180°. This value would however be very bad for the overall conditioning and precision of the astrometric solution (Makarov 1997, 1998; Lindegren & Bastian 2011), which instead favours a value around 90°. Moreover, the sensitivity to the cosΩ variation is only a factor 1.4 smaller at 180° than at 90°. The Gaia value Γ_{c} = 106.5° is therefore a reasonable choice.
That the sensitivity to the basic angle variation increases with decreasing Γ_{c} can be understood from simple arguments. Consider the alongscan effects of the parallax shift δϖ. As long as the nominal basic angle Γ_{c} is large, the effects in the two fields of view are significantly different. However, the smaller the basic angle is, the more similar are the effects in the two fields of view. This can be seen from Eq. (14) but is also obvious without any formula. We emphasise that it is the basic angle variations that induce field angle perturbations and that the solution tries to find parallaxes and attitude parameters that fit the perturbed field angles. It is then obvious that a smaller basic angle will require a larger parallax shift to absorb a basic angle variation of given amplitude.
4.3. Time dependence of R and ξ
As noted above, the barycentric distance R and the angle ξ are not strictly constant but functions of time. In this case, Eq. (15) gives the particular time dependencies of δΓ, δ_{x}, δ_{y}, and δ_{z} that are degenerate with a global parallax shift. In particular, (18)where C is constant, is indistinguishable from a parallax shift of .
Any other form of the variation δΓ(t) is not completely degenerate with δϖ. It may therefore contain components that can be detected by analysis of the residuals and subsequently eliminated by means of additional calibration terms. However, an arbitrary variation δΓ(t) in general also contains a component of the form of Eq. (18), which will result in some parallax shift. This shift can be estimated by projecting the variation onto the function on the righthand side of Eq. (18) in the leastsquares sense: (19)where the angular brackets denote averaging over time. If R and ξ are constant, the factor (Rsinξ)^{1} can be taken out of the averages. If, in addition, δΓ(t) is strictly periodic in Ω, we recover the first equality in Eq. (17).
For Gaia, which operates in the vicinity of L_{2}, its barycentric distance R varies between approximately 0.99 and 1.03 au as a combination of the eccentric heliocentric orbit of the Earth, the Lissajour orbit around L_{2}, and the timedependent offset between the Sun and the solar system’s barycentre. As discussed in Sect. 4.5 below, ξ varies by about 1% from its nominal value of 45°.
4.4. Effect of the finite size of the fields of view
To derive Eqs. (14) and (15) we neglected the finite size of the field of view by considering only observations at the field centre (g = h = 0). This was necessary to obtain an exact relation between the parallax shift and the basic angle and attitude variations. In a finite field of view, additional terms appear due to the variation of the parallax factor across the field, which cannot be represented by a unique set of basic angle variations. These terms are of the order of Φ ≃ 10^{2} times smaller than the basic angle variations. Here Φ is half the size of the field of view. As a consequence, a basic angle variation of the form of Eq. (18) is not strictly degenerate with δϖ and the attitude angles when a finite field of view is considered. However, if the instrument has also periodic optical distortions separately in each field of view that need to be calibrated, the corresponding, more complex calibration model may contribute to the degeneracy and, in the worst case, restore complete degeneracy.
4.5. Dependence on heliotropic coordinates
The spherical coordinates R, ξ, Ω introduced in Sect. 2.6 define the position of the solar system barycentre in the SRS. This is what matters for calculating the parallax effect, which depends on the observer’s displacement from the barycentre. On the other hand, the physical relevance of the cosΩ modulation is connected with the changing illumination of the satellite by the Sun, which depends on the heliocentric distance R_{h} and the heliotropic angles ξ_{h} and Ω_{h} that represent the proper direction towards the centre of the Sun at the time of observation. The difference between the heliotropic and barytropic coordinates is at most about 0.01 au and 0.01 rad, respectively. This is small but should be taken into account for an accurate modelling of the basic angle variations. In this context it can be noted that the expected thermal impact on the satellite scales as , whereas the parallax factor scales as R. On the other hand, the scanning law is chosen to keep ξ_{h} as constant as possible, whereas ξ can vary at the level of 1%. If the basic angle varies as , it is no longer strictly of the form of Eq. (18). The resulting parallax bias can be estimated by means of Eq. (19).
4.6. Breaking the degeneracy?
If the basic angle varies as a consequence of the changing solar illumination of the rotating satellite, we expect to see a parallax bias according to Eq. (19). However, as discussed above, the degeneracy with the parallax zero point is not perfect, and in principle this opens up a possibility to calibrate the basic angle variations from the observations. At least three different effects that contribute to breaking the degeneracy could be used: the finite size of the field of view (Sect. 4.4), the time dependence of R due to the eccentricity of the Earth’s orbit (Sect. 4.3), and the difference between the barytropic and heliotropic angles (Sect. 4.5). Unfortunately all three effects only come in at a level of a few per cent of the variation, or less, which makes the result very sensitive to small errors in the calibration model. Moreover, the finite field of view is of little use if we have to calibrate complex periodic variations of optical distortions independently in each field of view. The best chance may be offered by the time variation of R, for which the ratio of the parallax effect to the illumination strength goes as R^{3} and consequently varies by ± 5% over the year. Thus the hope of breaking the degeneracy purely from the observations themselves, based on the selfcalibration principle, is somewhat limited.
5. Conclusions
We present an analysis of the effect of basic angle variations on the global shift of parallaxes derived from observations by a scanning astrometric satellite with two fields of view. The method of small perturbations was used to derive the changes in the four observables, the across and alongfield angles in both fields of view, resulting from perturbations of four instrument parameters, the basic angle and three components of the attitude. Conversely, any given perturbation of the four observables can equally be represented by a specific combination of the instrument parameters. Applying this technique to the perturbations induced by a change in the parallax, we derive the timedependent variations of the instrument parameters that exactly mimic a global shift of the parallaxes.
These relations confirm previous findings that an uncorrected variation of the basic angle of the form , with Ω being the barycentric spin phase, leads to a global shift of the parallax zero point of for the parameters of the Gaia design. Results of numerical simulations are in complete agreement with the analytical formulae.
In general, periodic variations of the basic angle can be expected from the thermal impact of solar radiation on the spinning satellite (Lindegren 1977, 2004). Those periodic variations are typically related to the heliotropic spin phase Ω_{h}, which is close to the barycentric spin phase Ω. If the thermallyinduced variations contain a significant component that is proportional to cosΩ_{h}, their effect on the observations is practically indistinguishable from a global shift of the parallaxes. Although the degeneracy is not perfect, it is difficult to break except by using of other kinds of data or external information. In the case of Gaia this includes, in particular, direct measurement of the basic angle variations by means of laser metrology (BAM). The use of astrophysical information such as parallaxes of pulsating stars (Windmark et al. 2011; Gould & Kollmeier 2017) and quasars is vital for verifying the successful determination of the parallax zero point.
Acknowledgments
The authors acknowledge useful discussions with many colleagues within the Gaia community, of which we especially wish to mention Ulrich Bastian, Anthony Brown, Jos de Bruijne, Uwe Lammers, François Mignard, and Timo Prusti. The authors warmly thank the anonymous referee for valuable comments and suggestions. The work at Technische Universität Dresden was partially supported by the BMWi grants 50 QG 0601, 50 QG 0901 and 50 QG 1402 awarded by the Deutsche Zentrum für Luft und Raumfahrt e.V. (DLR).
References
 Arenou, F., Lindegren, L., Froeschle, M., et al. 1995, A&A, 304, 52 [NASA ADS] [Google Scholar]
 de Bruijne, J. H. J. 2012, Ap&SS, 341, 31 [NASA ADS] [CrossRef] [Google Scholar]
 ESA, ed. 1997, The Hipparcos and TYCHO catalogues. Astrometric and photometric star catalogues derived from the ESA HIPPARCOS Space Astrometry Mission, ESA SP, 1200 [Google Scholar]
 Gaia Collaboration (Brown, A. G. A., et al.) 2016a, A&A, 595, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gaia Collaboration (Prusti, T., et al.) 2016b, A&A, 595, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Golub, G. H., & van Loan, C. F. 1996, Matrix computations (The John Hopkins Univ. Press) [Google Scholar]
 Gould, A., & Kollmeier, J. A. 2017, Journal of the Korean Astronomical Society, 50, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Holl, B., Lindegren, L., & Hobbs, D. 2012, A&A, 543, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kovalevsky, J., Lindegren, L., Perryman, M. A. C., et al. 1997, A&A, 323, 620 [NASA ADS] [Google Scholar]
 Lindegren, L. 1977, Thermal stability and the determination of parallaxes, Tech. rep., Lund Observatory [Google Scholar]
 Lindegren, L. 2004, Scientific requirements for basic angle stability monitoring, Tech. rep., Lund Observatory [Google Scholar]
 Lindegren, L., & Bastian, U. 2011, in EAS Publ. Ser., 45, 109 [Google Scholar]
 Lindegren, L., Lammers, U., Hobbs, D., et al. 2012, A&A, 538, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lindegren, L., Lammers, U., Bastian, U., et al. 2016, A&A, 595, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Makarov, V. V. 1997, in Hipparcos – Venice ’97, eds. R. M. Bonnet, E. Høg, P. L. Bernacca, et al., ESA SP, 402, 823 [NASA ADS] [Google Scholar]
 Makarov, V. V. 1998, A&A, 340, 309 [NASA ADS] [Google Scholar]
 Mora, A., Biermann, M., Brown, A. G. A., et al. 2014, in Space Telescopes and Instrumentation 2014: Optical, Infrared, and Millimeter Wave, Proc. SPIE, 9143, 91430X [CrossRef] [Google Scholar]
 Perryman, M. A. C., de Boer, K. S., Gilmore, G., et al. 2001, A&A, 369, 339 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 van Leeuwen, F. 2005, A&A, 439, 805 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Windmark, F., Lindegren, L., & Hobbs, D. 2011, A&A, 530, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Parallax shift and attitude harmonics obtained in the smallscale direct solution (of type SA) with mas basicangle variation.
Singular values σ_{i} of the normal matrix for the different types of direct solutions.
All Figures
Fig. 1 Definition of the instrument axes x, y, z of the Scanning Reference System (SRS), the basic angle Γ, and the field angles g and h specifying the observed direction to a star (u) in either field of view. ϕ is the alongscan instrument angle of the star. In the SRS the direction to the solar system barycentre, b, is specified by the angles ξ and Ω. 

Open with DEXTER  
In the text 
Fig. 2 Parallax zero point shift for a basic angle variation of the form cosΩ with amplitude 1 mas. The dots show results of the largescale iterative solution for several different basic angles, including the Hipparcos and Gaia values; the solid curve shows the theoretical relation from Eq. (17). 

Open with DEXTER  
In the text 