The discovery of the accelerating expansion of the Universe, thought to be driven by a mysterious form of `dark energy' constituting most of the Universe, has further revived the interest in testing Einstein's theory of General Relativity. At the very foundation of Einstein's theory is the geodesic motion of a small, structureless test-particle. Depending on the physical context, a star, planet or satellite can behave very nearly like a test-particle, so geodesic motion is used to calculate the advance of the perihelion of a planet's orbit, the dynamics of a binary pulsar system and of an Earth orbiting satellite. Verifying geodesic motion is then a test of paramount importance to General Relativity and other theories of fundamental physics. On the basis of the first few months of observations of the recently launched satellite LARES, its orbit shows the best agreement of any satellite with the test-particle motion predicted by General Relativity. That is, after modelling its known non-gravitational perturbations, the LARES orbit shows the smallest deviations from geodesic motion of any artificial satellite. LARES-type satellites can thus be used for accurate measurements and for tests of gravitational and fundamental physics. Already with only a few months of observation, LARES provides smaller scatter in the determination of several low-degree geopotential coefficients (Earth gravitational deviations from sphericity) than available from observations of any other satellite or combination of satellites.

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.