ESA recently called for new "Science Ideas" to be investigated in terms of feasibility and technological developments -- for technologies not yet sufficiently mature. These ideas may in the future become candidates for M or L class missions within the ESA Science Program. With the launch of Gaia in December 2013, Europe entered a new era of space astrometry following in the footsteps of the very successful Hipparcos mission from the early 1990s. Gaia is the successor to Hipparcos, both of which operated in optical wavelengths, and Gaia is two orders of magnitude more accurate in the five astrometric parameters and is surveying four orders of magnitude more stars in a vast volume of the Milky Way. The combination of the Hipparcos/Tycho-2 catalogues with the first early Gaia data release will give improved proper motions over a long ~25 year baseline. The final Gaia solution will also establish a new optical reference frame by means of quasars, by linking the optical counterparts of radio (VLBI) sources defining the orientation of the reference frame, and by using the zero proper motion of quasars to determine a non-rotating frame. A weakness of Gaia is that it only operates at optical wavelengths. However, much of the Galactic centre and the spiral arm regions, important for certain studies, are obscured by interstellar extinction and this makes it difficult for Gaia to deeply probe. Traditionally, this problem is overcome by switching to the infra-red but this was not possible with Gaia's CCDs. Additionally, to scan the entire sky and make global absolute parallax measurements the spacecraft must have a constant rotation and this requires that the CCDs operate in TDI mode, increasing their complexity.

A new all-sky visible and Near-InfraRed (NIR) space astrometry mission with a wavelength cutoff in the K-band is not just focused on a single or small number of key science cases. Instead, it is extremely broad, answering key science questions in nearly every branch of astronomy while also providing a dense and accurate visible-NIR reference frame needed for future astronomy facilities. For almost 2 billion common stars the combination of Gaia and a new all-sky NIR astrometry mission would provide much improved proper motions, answering key science questions -- from the solar system and stellar systems, including exoplanet systems, to compact galaxies, quasars, neutron stars, binaries and dark matter substructures. The addition of NIR will result in up to 8 billion newly measured stars in some of the most obscured parts of our Galaxy, and crucially reveal the very heart of the Galactic bulge region. In this white paper we argue that rather than improving on the accuracy, a greater overall science return can be achieved by going deeper than Gaia and by expanding the wavelength range to the NIR.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Gaia Data Release 2 (Gaia DR2) contains results for 1693 million sources in the magnitude range 3 to 21 based on observations collected by the European Space Agency Gaia satellite during the first 22 months of its operational phase. We describe the input data, models, and processing used for the astrometric content of Gaia DR2, and the validation of these results performed within the astrometry task. Some 320 billion centroid positions from the pre-processed astrometric CCD observations were used to estimate the five astrometric parameters (positions, parallaxes, and proper motions) for 1332 million sources, and approximate positions at the reference epoch J2015.5 for an additional 361 million mostly faint sources. Special validation solutions were used to characterise the random and systematic errors in parallax and proper motion. For the sources with five-parameter astrometric solutions, the median uncertainty in parallax and position at the reference epoch J2015.5 is about 0.04 mas for bright (G<14 mag) sources, 0.1 mas at G=17 mag, and 0.7 mas at G=20 mag. In the proper motion components the corresponding uncertainties are 0.05, 0.2, and 1.2 mas/yr, respectively. The optical reference frame defined by Gaia DR2 is aligned with ICRS and is non-rotating with respect to the quasars to within 0.15 mas/yr. From the quasars and validation solutions we estimate that systematics in the parallaxes depending on position, magnitude, and colour are generally below 0.1 mas, but the parallaxes are on the whole too small by about 0.03 mas. Significant spatial correlations of up to 0.04 mas in parallax and 0.07 mas/yr in proper motion are seen on small (<1 deg) and intermediate (20 deg) angular scales. Important statistics and information for the users of the Gaia DR2 astrometry are given in the appendices.

The second release of Gaia data (Gaia DR2) contains the astrometric parameters for more than half a million quasars. This set defines a kinematically non-rotating reference frame in the optical domain referred to as the Gaia-CRF2. The Gaia-CRF2 is the first realisation of a non-rotating global optical reference frame that meets the ICRS prescriptions, meaning that it is built only on extragalactic sources. It consists of the positions of a sample of 556 869 sources in Gaia DR2, obtained from a positional cross-match with the ICRF3-prototype and AllWISE AGN catalogues. The sample constitutes a clean, dense, and homogeneous set of extragalactic point sources in the magnitude range G from 16 to 21 mag with accurately known optical positions. The median positional uncertainty is 0.12 mas for G < 18 mag and 0.5 mas at G = 20 mag. Large-scale systematics are estimated to be in the range 20 to 30 muas. The accuracy claims are supported by the parallaxes and proper motions of the quasars in Gaia DR2. The optical positions for a subset of 2820 sources in common with the ICRF3-prototype show very good overall agreement with the radio positions, but several tens of sources have significantly discrepant positions.

This paper discusses the effects of gravitational waves on high-accuracy astrometric observations such as those delivered by Gaia. Depending on the frequency of gravitational waves, two regimes for the influence of gravitational waves on astrometric data are identified: the regime when the effects of gravitational waves directly influence the derived proper motions of astrometric sources and the regime when those effects mostly appear in the residuals of the standard astrometric solution. The paper is focused on the second regime while the known results for the first regime are briefly summarized. The deflection of light due to a plane gravitational wave is then discussed. Starting from a model for the deflection we derive the corresponding partial derivatives and summarize some ideas for the search strategy of such signals in high-accuracy astrometric data. In order to reduce the dimensionality of the parameter space the use of vector spherical harmonics is suggested and explained. The explicit formulas for the VSH expansion of the astrometric signal of a plain gravitational wave are derived. Finally, potential sensitivity of Gaia astrometric data is discussed. Potential astrophysical sources of gravitational waves that can be interesting for astrometric detection are identified.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Determination of absolute parallaxes by means of a scanning astrometric satellite such as Hipparcos or Gaia relies on the short-term stability of the so-called basic angle between the two viewing directions. Uncalibrated variations of the basic angle may produce systematic errors in the computed parallaxes. 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. The changes in observables produced by small perturbations of the basic angle, attitude, and parallaxes are calculated analytically. We then look for a combination of perturbations that has no net effect on the observables. 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 such 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.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

As part of the data processing for Gaia Data Release~1 (Gaia DR1) a special astrometric solution was computed, the so-called auxiliary quasar solution. This gives positions for selected extragalactic objects, including radio sources in the second realisation of the International Celestial Reference Frame (ICRF2) that have optical counterparts bright enough to be observed with Gaia. A subset of these positions was used to align the positional reference frame of Gaia DR1 with the ICRF2. We describe the properties of the Gaia auxiliary quasar solution for a subset of sources matched to ICRF2, and compare their optical and radio positions at the sub-mas level. Their formal standard errors are better than 0.76~milliarcsec (mas) for 50% of the sources and better than 3.35~mas for 90%. Optical magnitudes are obtained in Gaia's unfiltered photometric G band. The comparison with the radio positions of the defining sources shows no systematic differences larger than a few tenths of a mas. The fraction of questionable solutions, not readily accounted for by the statistics, is less than 6%. Normalised differences have extended tails requiring case-by-case investigations for around 100 sources, but we have not seen any difference indisputably linked to an optical-radio offset in the sources.

Gaia Data Release 1 (Gaia DR1) contains astrometric results for more than 1 billion stars brighter than magnitude 20.7 based on observations collected by the Gaia satellite during the first 14 months of its operational phase. We give a brief overview of the astrometric content of the data release and of the model assumptions, data processing, and validation of the results. For stars in common with the Hipparcos and Tycho-2 catalogues, complete astrometric single-star solutions are obtained by incorporating positional information from the earlier catalogues. For other stars only their positions are obtained by neglecting their proper motions and parallaxes. The results are validated by an analysis of the residuals, through special validation runs, and by comparison with external data. Results. For about two million of the brighter stars (down to magnitude ~11.5) we obtain positions, parallaxes, and proper motions to Hipparcos-type precision or better. For these stars, systematic errors depending e.g. on position and colour are at a level of 0.3 milliarcsecond (mas). For the remaining stars we obtain positions at epoch J2015.0 accurate to ~10 mas. Positions and proper motions are given in a reference frame that is aligned with the International Celestial Reference Frame (ICRF) to better than 0.1 mas at epoch J2015.0, and non-rotating with respect to ICRF to within 0.03 mas/yr. The Hipparcos reference frame is found to rotate with respect to the Gaia DR1 frame at a rate of 0.24 mas/yr. Based on less than a quarter of the nominal mission length and on very provisional and incomplete calibrations, the quality and completeness of the astrometric data in Gaia DR1 are far from what is expected for the final mission products. The results nevertheless represent a huge improvement in the available fundamental stellar data and practical definition of the optical reference frame.

These are lecture notes of the lecture on elementary celestial mechanics read at Technische Universitaet Dresden, Dresden, Germany from 2000 till 2011.

The Gaia all-sky astrometric survey is challenged by several issues affecting the spacecraft stability. Amongst them, we find the focus evolution, straylight and basic angle variations Contrary to pre-launch expectations, the image quality is continuously evolving, during commissioning and the nominal mission. Payload decontaminations and wavefront sensor assisted refocuses have been carried out to recover optimum performance. Straylight and basic angle variations several orders of magnitude greater than foreseen were found and studied during commissioning by the Gaia scientists (payload experts). Building on their investigations, an ESA-Airbus DS working group was established during the early nominal mission and worked on a detailed root cause analysis. In parallel, Gaia scientists have also continued analysing the data, most notably comparing the BAM signal to global astrometric solutions, with remarkable agreement. In this contribution, a status review of these issues will be provided, with emphasis on the mitigation schemes and the lessons learned for future space missions where extreme stability is a key requirement.

Comparison of stellar catalogues with position and proper motion components using a decomposition on a set of orthogonal vector spherical harmonics. We show the theoretical and practical advantages of this technique as a result of invariance properties and the independence of the decomposition from a prior model. We describe the mathematical principles used to perform the spectral decomposition, evaluate the level of significance of the multipolar components and examine the transformation properties under space rotation. The principles are illustrated with a characterisation of the systematic effects in the FK5 catalogue compared to Hipparcos and with an application to the extraction of the rotation and dipole acceleration in the astrometric solution of QSOs expected from Gaia.

A highly precise model for the motion of a rigid Earth is indispensable to reveal the effects of non-rigidity in the rotation of the Earth from observations. To meet the accuracy goal of modern theories of Earth rotation of 1 microarcsecond (muas) it is clear, that for such a model also relativistic effects have to be taken into account. The largest of these effects is the so called geodetic precession. In this paper we will describe this effect and the standard procedure to deal with it in modeling Earth rotation up to now. With our relativistic model of Earth rotation Klioner et al. (2001) we are able to give a consistent post-Newtonian treatment of the rotational motion of a rigid Earth in the framework of General Relativity. Using this model we show that the currently applied standard treatment of geodetic precession is not correct. The inconsistency of the standard treatment leads to errors in all modern theories of Earth rotation with a magnitude of up to 200 muas for a time span of one century.

Although the formulas for the light deflection due to quadrupole gravitational field of deflecting bodies are well known, the formulas are rather complicated, so that massive computations of quadrupole light deflection (e.g., in the framework of astrometric survey missions like Gaia) are time-consuming. Considering an observer situated within a few million kilometers from the Earth (clearly the most practical case), we derive the simplest possible form of the relevant formulas still having numerical accuracy of 1 micro-arcsecond. This form leads to simple upper estimates for the quadrupole light deflection in various cases allowing one to relate the magnitude of the actual quadrupole deflection with the corresponding monopole deflection due to the same body. These upper estimates can be used to decide if, for a given configuration, the actual quadrupole deflection should be computed for a given accuracy goal.

Efficient computation of the quadrupole light deflection for quasars/quasars and solar system objects within the framework of the baseline Gaia relativity model (GREM) is discussed. Two refinements have been achieved with the goal to improve the performance of the model: First, the quadrupole deflection formulas for both cases are simplified as much as possible considering the Gaia nominal orbit (only approximate minimal distances between Gaia and the giant planets were used here), physical parameters of the giant planets and the envisaged accuracy of 1 microarcsecond for individual systematic effects. The recommended formulas are given by Eq.(40) for stars/quasars and by Eq.(81) for solar system objects. Second, simple expressions for the upper estimate of the quadrupole light deflection have been found allowing, with a few additional arithmetical operations, to judge a priori if the quadrupole light deflection should be computed or not for a given source and for a given requested accuracy. The recommended criteria are given by Eq.(45) for stars/quasars and by Eq.(92) for solar system objects. Additionally, the quadrupole Shapiro effect for solar system objects is reconsidered. A strict upper bound of quadrupole Shapiro effect for solar system objects is given in Eq.(109).

This paper summarizes the analysis of the consequences of the violation of the Local Lorentz Invariance (LLI) on astrometric observations. We demonstrate that from the point of view of the LLI astrometric observations represent an experiment of Michelson-Morley type. The future high-accuracy astrometric projects (e.g., Gaia) will be used to test the LLI.

Numerical integration of the differential equations of light propagation in the Schwarzschild metric shows that in some situations relevant for practical observations the well-known post-Newtonian solution for light propagation has an error up to 16 microarcsecond. The aim of this work is to demonstrate this fact, identify the reason for this error and to derive an analytical formula accurate at the level of 1 microarcsecond as needed for high-accuracy astrometric projects (e.g., Gaia). An analytical post-post-Newtonian solution for the light propagation for both Cauchy and boundary problems is given for the Schwarzschild metric augmented by the PPN and post-linear parameters $\beta$, $\gamma$ and $\epsilon$. Using analytical upper estimates of each term we investigate which post-post-Newtonian terms may play a role for an observer in the solar system at the level of 1 microarcsecond and conclude that only one post-post-Newtonian term remains important for this numerical accuracy. In this way, an analytical solution for the boundary problem for light propagation is derived. That solution contains terms of both post-Newtonian and post-post-Newtonian order, but is valid for the given numerical level of 1 microarcsecond. The derived analytical solution has been verified using the results of a high-accuracy numerical integration of differential equations of light propagation and found to be correct at the level well below 1 microarcsecond for arbitrary observer situated within the solar system. Furthermore, the origin of the post-post-Newtonian terms relevant for the microarcsecond accuracy is elucidated. We demonstrate that these terms result from an inadequate choice of the impact parameter in the standard post-Newtonian formulas.

Relativistic modelling of rotational motion of extended bodies represents one of the most complicated problems of Applied Relativity. The relativistic reference systems of IAU (2000) give a suitable theoretical framework for such a modelling. Recent developments in the post-Newtonian theory of Earth rotation in the limit of rigidly rotating multipoles are reported below. All components of the theory are summarized and the results are demonstrated. The experience with the relativistic Earth rotation theory can be directly applied to model the rotational motion of other celestial bodies. The high-precision theories of rotation of the Moon, Mars and Mercury can be expected to be of interest in the near future.

In [2] we have seen that among the post-post-Newtonian terms in the light deflection formulas there are "enhanced" ones that may become much larger than the other post-post-Newtonian terms. It is demonstrated here that these "enhanced" terms result from an inadequate choice of the impact parameter. Introducing another impact parameter, that can be considered as coordinate-independent, we demonstrate that all "enhanced" terms disappear from the formulas.

Numerical integration of the differential equations of light propagation in the Schwarzschild metric shows that in some extreme situations relevant for practical observations (e.g. for Gaia) the well-known standard post-Newtonian formula for the boundary problem has an error up to 16 \muas. The aim of this note is to identify the reason for this error and to derive an extended formula accurate at the level of 1 \muas as needed e.g. for Gaia. The analytical parametrized post-post-Newtonian solution for light propagation derived by \citet{report1} gives the solution for the boundary problem with all analytical terms of order $\OO4$ taken into account. Giving an analytical upper estimates of each term we investigate which post-post-Newtonian terms may play a role for an observer in the solar system at the level of 1 \muas. We conclude that only one post-post-Newtonian term remains important for this numerical accuracy and derive a simplified analytical solution for the boundary problem for light propagation containing all the terms that are indeed relevant at the level of 1 \muas. The derived analytical solution has been verified using the results of a high-accuracy numerical integration of differential equations of light propagation and found to be correct at the level well below 1 \muas for arbitrary observer situated within the solar system.

This note suggests nomenclature for dealing with the units of various astronomical quantities that are used with the relativistic time scales TT, TDB, TCB and TCG. It is suggested to avoid wordings like "TDB units" and "TT units" and avoid contrasting them to "SI units". The quantities intended for use with TCG, TCB, TT or TDB should be called "TCG-compatible", "TCB-compatible", "TT-compatible" or "TDB-compatible", respectively. The names of the units second and meter for numerical values of all these quantities should be used with out any adjectives. This suggestion comes from a special discussion forum created within IAU Commission 52 "Relativity in Fundamental Astronomy".

A rigorous analytical solution of light propagation in Schwarzschild metric in post-post Newtonian approximation has been presented in \cite{report1}. In \cite{report2} it has been asserted that the sum of all those terms which are of order ${{\cal O} (\frac{m^2}{d^2})}$ and ${{\cal O}(\frac{m^2}{d_\sigma^2})}$ is not greater than $15/4 \pi \frac{m^2}{d^2}}$ and $15/4 \pi \frac{m^2}{d_\sigma^2}}$, respectively. All these terms can be neglected on microarcsecond level of accuracy, leading to considerably simplified analytical transformations of light propagation. In this report, we give formal mathematical proofs for the inequalities used in the appendices of \cite{report2}.

An analytical solution for light propagation in the post-post-Newtonian approximation is given for the Schwarzschild metric in harmonic gauge augmented by PPN and post-linear parameters $\beta$, $\gamma$ and $\epsilon$. The solutions of both Cauchy and boundary problem are given. The Cauchy problem is posed using the initial position of the photon $\ve{x}_0 = \ve{x}(t_0)$ and its propagation direction \ve{\sigma} at minus infinity: $\ve{\sigma} = {1\over c} \lim\limits_{t \to -\infty}\dot{\ve{x}}(t)$. An analytical expression for the total light deflection is given. The solutions for $t - t_0$ and $\ve{\sigma}$ are given in terms of boundary conditions $\ve{x}_0 = \ve{x} (t_0)$ and $\ve{x} = \ve{x}(t)$.

ASTROD I is a planned interplanetary space mission with multiple goals. The primary aims are: to test General Relativity with an improvement in sensitivity of over 3 orders of magnitude, improving our understanding of gravity and aiding the development of a new quantum gravity theory; to measure key solar system parameters with increased accuracy, advancing solar physics and our knowledge of the solar system and to measure the time rate of change of the gravitational constant with an order of magnitude improvement and the anomalous Pioneer acceleration, thereby probing dark matter and dark energy gravitationally. It is an international project, with major contributions from Europe and China and is envisaged as the first in a series of ASTROD missions. ASTROD I will consist of one spacecraft carrying a telescope, four lasers, two event timers and a clock. Two-way, two-wavelength laser pulse ranging will be used between the spacecraft in a solar orbit and deep space laser stations on Earth, to achieve the ASTROD I goals. A second mission, ASTROD II is envisaged as a three-spacecraft mission which would test General Relativity to one part per billion, enable detection of solar g-modes, measure the solar Lense-Thirring effect to 10 parts per million, and probe gravitational waves at frequencies below the LISA bandwidth. In the third phase (ASTROD III or Super-ASTROD), larger orbits could be implemented to map the outer solar system and to probe primordial gravitational-waves at frequencies below the ASTROD II bandwidth.

This work gives an explicit exact expression for the Thomas precession arising in the framework of Special Theory of Relativity as the spatial rotation resulting from two subsequence Lorentz boosts. The final result for the orthogonal matrix of Thomas precession is given by Eqs. (21)--(25). A trivial calculation leads to the compact formula (26) for the angle of rotation due to Thomas precession. In the framework of Gaia the special-relativistic Thomas precession is an important step in the derivation of an aberrational formula with the Mansouri-Sexl parameters. The latter formula will be used to test the Local Lorentz Invariance with Gaia data as will be explained elsewhere.

The problem of determination of the orbital velocity of an astrometric satellite from its own observational data is studied. It is well known that data processing of microarcsecond-level astrometric observations imposes very stringent requirements on the accuracy of the orbital velocity of the satellite (a velocity correction of 1.45 mm/s implies an aberrational correction of 1 $\mu$as). Because of a number of degeneracies the orbital velocity cannot be fully restored from observations provided by the satellite. Seven constraints that must be applied on a velocity parameterization are discussed and formulated mathematically. It is shown what part of velocity can be recovered from astrometric data by a combined fit of both velocity parameters and astrometric parameters of the sources. Numerical simulations show that, with the seven constraints applied, the velocity and astrometric parameters can be reliably estimated from observational data. It is also argued that the idea to improve the velocity of an astrometric satellite from its own observational data is only useful if the a priori information on the orbital velocity justifies the applicability of the velocity constraints. The proposed model takes into account only translational motion of the satellite and ignores any satellite-specific parameters. Therefore, the results of this study are equally applicable to both scanning missions similar to Gaia, and pointing ones like SIM, provided that enough sources were observed sufficiently uniformly.

For relativistic modelling of high-accuracy astronomical data several time scales are used: barycentric and geocentric coordinate times, TCB and TCG, as well as two additional time scales, TDB and TT, that are defined as linear functions of TCB and TCG, respectively. The paper is devoted to a concise but still detailed explanation of the reasons and the implications of the relativistic scalings of astronomical quantities induced by the time scales TDB and TT. We consequently distinguish between quantities and their numerical values expressed in some units. It is argued that the scaled time scales, the scaled spatial coordinates and the scaled masses should be considered as distinct quantities which themselves can be expressed in any units, and not as numerical values of the same quantities expressed in some different, non-SI units (``TDB units'' and ``TT units''). Along the same lines of argumentation the system of astronomical units is discussed in the relativistic framework. The whole freedom in the definitions of the systems of astronomical units for TCB and TDB is demonstrated. A number of possible ways to freeze the freedom are shown and discussed. It is argued that in the future one should think about converting AU into a defined quantity by fixing its value in SI meters.

Context. High accuracy astrometric instruments like Gaia aiming at an accuracy of 1 microarcsecond cannot be considered as point-like observers in the framework of relativistic modelling of observable quantities. Aims. Special-relativistic effects on the imaging by a non-point-like arbitrarily moving optical instrument are discussed. Methods. A special-relativistic reflection law for a mirror of arbitrary shape and motion is derived in the limit of geometrical optics. The aberration patterns are computed with ray tracing using a full pecial-relativistic model for two simple rotating optical instrument. Results. It was found that the effect of special-relativistic reflection law on the photocenters of aberration patterns of an optical system rotating with a moderate angular velocity of $60 \arcsec/{\rm s}$ may be at the level of 1 microarcsecond if the system involves mirrors significantly inclined relative to the optical axis. Conclusions. Special-relativistic optical modelling of the future astrometric instrument is generally speaking indispensable if a level of a few microarcseconds is envisaged.

This paper represents a first attempt of embedding the Barycentric Celestial Reference System, the fundamental relativistic reference system to be used for the modeling of Gaia observations, into some cosmological background. The general Robertson-Walker metric is transformed into local coordinates where the cosmological effects are represented as tidal potentials. A version of a cosmological BCRS is then suggested to lowest order. The effects of cosmological background on the motion of the solar system are estimated and found to be completely negligible. The relation to the de Sitter and Schwarzschild-de Sitter solutions is discussed.

After a short review of experimental foundations of metric theories of gravity, the choice of general relativity as a theory to be used for the routine modeling of Gaia observations is justified. General principles of relativistic modeling of astronomical observations are then sketched and compared to the corresponding Newtonian principles. The fundamental reference system -- Barycentric Celestial Reference System, which has been chosen to be the relativistic reference system underlying the future Gaia reference frame is presented. Principal relativistic effects in each constituent of a relativistic model of astronomical observations are briefly elucidated. The structure of a relativistic model of positional observations which can be used as a standard relativistic model for Gaia is sketched. The physical meaning of the Gaia reference frame is discussed. It is discussed also how Gaia observations can be used to verify general relativity.

A relativistic definition of the physically adequate proper reference system of a test observer is suggested within the framework of the PPN formalism. According to the nomenclature accepted within the GAIA project this reference system is called Center-of-Mass Reference System (CoMRS). The interrelation between the suggested definition of the CoMRS and the Resolutions 2000 on relativity of the International Astronomical Union (IAU) are elucidated. The tetrad representation of the CoMRS at its origin is also explicated. It is demonstrated how to use that tetrad representation to calculate the relation between the observed direction of a light ray and the corresponding coordinate direction in the Barycentric Celestial Reference System of the IAU. It is argued that the kinematically non-rotating CoMRS is the natural choice of the reference system where the attitude of the observer (e.g. of the GAIA satellite) should be modeled. The relativistic equations of rotational motion of a satellite relative to its CoMRS are briefly discussed. A simple algorithm for the attitude description of the satellite is proposed.

One of the most subtle points in the modern relativistic models for microarcsecond astrometrical observations is the treatment of the influence of translational motion of gravitating bodies on the light propagation. This paper describes numerical simulations of the light propagation in the gravitational field of moving gravitating bodies as well as summarizes the underlying theory. The simulations include high-precision numerical integrations of both post-Newtonian and post-Minkowskian differential equations of light propagation and a detailed comparison of the results of the numerical integrations with various available approximate analytical formulas. The simulations has been performed both for hypothetical bodies with various parameters of trajectories as well as for all the major bodies of the solar system using the JPL ephemeris DE405/LE405 to calculate their motion. It is shown that for the accuracy of $\sim 0.2$ $\mu$as it is sufficient to use the well-known solution for the light propagation in the field of a motionless mass monopole and substitute in that solution the position of the body at the moment of closest approach between the actual trajectory of the body and the unperturbed light path (as it was first suggested by Hellings (1986)). For a higher accuracy one should use either the post-Newtonian solution for uniformly moving bodies (Klioner & Kopeikin, 1992) or the post-Minkowskian solution for arbitrarily moving bodies (Kopeikin & Schaefer, 1999). For astrometric observations performed from within the solar system these two solutions guarantee the accuracy of $\sim 0.002$ $\mu$as and are virtually indistinguishable from each other.

We show how to derive the equations of light propagation in the gravitational field of uniformly moving mass monopoles without formulating and integrating the differential equations of light propagation in that field. The well-known equations of light propagation in the gravitational field of a motionless mass monopole are combined with a suitable Lorentz transformation. The possibility to generalize this technique for the more complicated case of uniformly moving body of arbitrary multipole structure is discussed.

This paper discusses the IAU Resolutions B1.3, B1.4, B1.5 and B1.9 (2000) that were adopted during the 24th General Assembly in Manchester, 2000 and provides details and explanations for these Resolutions. It is explained why they present significant progress over the corresponding IAU 1991 Resolutions and why they are necessary in the light of present accuracies in astrometry, celestial mechanics and metrology. In fact most of these Resolutions are consistent with astronomical models and software already in use. The metric tensors and gravitational potentials of both the Barycentric Celestial Reference System and Geocentric Celestial Reference System are defined and discussed. The necessity and relevance of the two celestial reference systems are explained. The transformations of coordinates and gravitational potentials are discussed. Potential coefficients parameterizing the post-Newtonian gravitational potentials are expounded. Simplified versions of the time transformations suitable for modern clock accuracies are elucidated. Various approximations used in the Resolutions are explicated and justified. Some models (e.g. for higher spin moments) that serve the purpose for estimating orders of magnitude have actually never been published before.

A set of equations describing the rotational motion of the Earth relative to the GCRS is formulated in the approximation of rigidly rotating multipoles. The external bodies are supposed to be mass monopoles. The derived set of formulas is supposed to form the theoretical basis for a practical post-Newtonian theory of Earth precession and nutation.

This paper is devoted to a practical model for relativistic reduction of positional observations with an accuracy of 1 \muas which is expected to be attained in the future space astrometry missions. All relativistic effects which are caused by the gravitational field of the Solar system and which are of practical importance for this accuracy level are thoroughly calculated and discussed. The model includes relativistic modeling of the motion of observer, modeling of relativistic aberration and gravitational light deflection as well as a relativistic treatment of parallax and proper motion suitable for the accuracy of 1 \muas. The model is formulated both for remote sources (stars, quasars, etc.) and for Solar system objects (asteroids, etc.). The suggested model is formulated within the framework of Parametrized Post-Newtonian Formalism with parameters $\beta$ and $\gamma$. However, for general relativity ($\beta=\gamma=1$) the model is fully compatible with the IAU Resolutions (2000) on relativity in celestial mechanics, astrometry and metrology. The model is presented in a form suitable for implementation in a software system for data processing or simulation. The changes which should be applied to the model to attain the accuracy of 0.1 \muas are reviewed. Potentially important relativistic effects caused by additional gravitational fields which are generated outside of the Solar system are also briefly discussed.