The Rotational Motion of Vesta

Total Page:16

File Type:pdf, Size:1020Kb

The Rotational Motion of Vesta A&A 556, A151 (2013) Astronomy DOI: 10.1051/0004-6361/201220732 & c ESO 2013 Astrophysics The rotational motion of Vesta N. Rambaux1,2 1 Université Pierre et Marie Curie, UPMC – Paris 06, 4 place Jussieu, 75005 Paris, France e-mail: [email protected] 2 IMCCE, Observatoire de Paris, CNRS UMR 8028, 77 Avenue Denfert-Rochereau, 75014 Paris, France Received 13 November 2012 / Accepted 4 March 2013 ABSTRACT Context. Vesta is the second largest body of the main asteroid belt, and it has been studied recently in great detail by the Dawn mission. It was the first time that this asteroid, or protoplanet, has been explored by a space mission, and it revealed a differentiated body. The knowledge of the rotational motion and, especially, its precession-nutation and length-of-day variations may add precious information on its interior. Aims. The objective of this paper is to present the first rotational model of Vesta based on the data acquired by the Dawn mission. The Dawn mission determined the orientation of Vesta with 0.01 degree accuracy, as well as the second degree of the gravity coefficients with a few percent accuracy. Methods. We built a semi-analytical model of the rotational motion of Vesta based on the orbital perturbations and the large triaxial shape of Vesta. The rotational motion is then described through the polar motion, precession, nutation, and length-of-day variations. The sensitivity of the precession-nutation is given as a function of the polar moment of inertia, which has not yet been determined. Results. We find that the amplitude of the nutation is about 2000 milli-arcsec at the semi-annual period, whereas the amplitude of the length-of-day variation is about 3 mas/days, in the semi-diurnal period. Finally, we show that the signature of the polar moment of inertia, which is crucial for constraining geophysical model, is on the order of 200 milli-arcsec for the semi-annual nutations for two extrema of polar moment with inertia values of 0.38–0.42. Key words. minor planets, asteroids: individual: Vesta – celestial mechanics 1. Introduction models used in this paper. Then, we explain the rotational the- ory, and in Sect. 5, we discuss the range of solutions. Finally, we The protoplanet Vesta is the second largest asteroid of the main conclude. asteroid belt with a mean radius of 262.7 ± 0.1 km (Russell et al. 2012), and since July 2011, it has belonged to the short list of as- teroids visited by a spacecraft. The Dawn mission orbited Vesta 2. Input data for almost one year and investigated the main characteristics of 2.1. Orientation data this body in great detail. Notably, the observational data allow assessing its origin as a remnant protoplanet from the forma- The Dawn mission has determined the orientation of the spin tion of the early solar system (Russell et al. 2012). Therefore, axis of Vesta during its mission from July 2011 to August 2012. it is crucial to determine its interior properties. One way to in- The spin axis orientation is described in terms of the right ascen- vestigate the interior structure is to accurately determine the sion α and declination δ. That means that the spin axis is given rotational motion, which is the objective of this paper. in a reference frame defined by the ICRF (International Celestial The Dawn mission has mapped the surface of Vesta with an Reference Frame). Table 1 presents the values of the right ascen- unprecedented resolution of 0.1 km that leads to determining sion and declination for Vesta obtained by Russell et al. (2012), the orientation of the spin axis at 0.01 degree accuracy (Russell Li et al. (2011), and Thomas et al. (1997). The data from Thomas et al. 2012). This accuracy improves the precision of the spin et al. (1997) comes from images obtained during a campaign position upon the Hubble Space Telescope (HST) measurements with the HST in December 1996, whereas Li et al. (2011)used by Thomas et al. (1997) and from the reanalysis of the com- combined ground data and HST data from 1986 to 2003. bined HST/ground data by Li et al. (2011) collected from 1983 to Figure 1 shows the projection of Vesta’s spin axis on the 2006 by a factor of 500. In parallel to Vesta’s mapping, the radio- ICRF XY plane reported by the different authors. We plot the science experiment provides the gravity field of Vesta (Russell ellipses of uncertainties related to each set of observations, in- et al. 2012) and thereby leading to constraints on the mass distri- cluding the five degrees uncertainties obtained by Thomas et al. bution inside Vesta and proving that the interior is differentiated. (1997) and Li et al. (2011), whereas the uncertainty ellipse ob- The objective of this paper is to present the first rotational tained by Dawn is 0.01 degree. The ellipse for the Dawn data model of Vesta based on the recent data obtained by Dawn and to is therefore too small to be visible on the scale of the figure. predict the expected precession-nutation and length-of-day vari- The figure also shows the projection of the orbital axis in the ations from various interior models. In return, the knowledge of XY plane located at the coordinates 0.124, –0.375. The determi- the rotational motion can be used to infer the interior properties nation of the orbital axis comes from the Horizons ephemerides of this protoplanet. The article is divided into four parts. First, we (Acton et al. 1996). In this figure, the position of the spin axis present the input data coming from Dawn and from geophysical with respect to the orbital axis can be described by two spherical Article published by EDP Sciences A151, page 1 of 8 A&A 556, A151 (2013) Table 1. Estimation of the orientation of the spin axis of Vesta. α (deg) β (deg) Ω0 (deg/day) Russell et al. (2012) 309.03 ± 0.01 42.23 ± 0.01 1617.333119 ± 0.000003 Li et al. (2011) 305.8 ± 3.1 41.4 ± 1.5 – Thomas et al. (1997) 301 ± 541± 5 1617.332776 -0.35 where M and R are the mass and radius of Vesta, respec- / 2 -0.4 tively. Because Vesta has a non zero C22, the values of A MR 2 ρ and B/MR are different. In addition, C22 is large, around 0.4 -0.45 of the C20, meaning that Vesta has a large triaxial shape. To ob- λ -0.5 tain the set of the three moments of inertia, it is necessary to add an external constraint on the value of the polar moment of / 2 ICRF -0.55 inertia C MR . According to Russell et al. (2012), Vesta is dif- Y ferentiated into two layers, a mantle and a core. In addition, they -0.6 suggest that the core could still be in hydrostatic equilibrium. As -0.65 a consequene, we compute the expected core’s flattening by us- ing Clairaut equation for a two-layer Vesta. We obtain the core’s -0.7 flattening equal to (ac − cc)/ac = 0.122, where ac and cc are the -0.75 equatorial and polar axis of the oblate core. Then, the total polar 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 moment of inertia of Vesta is expressed as X ICRF π π 4 2 2 8 4 Fig. 1. XY C = abc a + b ρm + a cc(ρc − ρm), (2) Projection in the ICRF plane of the spin axis of Vesta from 15 15 c Thomas et al. (1997) (green curve), Li et al. (2011) (blue curve) and Russell et al. (2012) (red curve). The ellipses represent the measured where a, b,andc are the major, intermediate, and minor axes uncertainties and in the case of Dawn the ellipse is not visible on the of the ellipsoidal shape of Vesta measured by the Dawn mis- scale of the figure. The black point represents the orientation of the ρ ρ λ ρ sion, and c and m are the density of the core and of the man- orbital axis of Vesta. The angles and represent the longitude and tle. The mass of the ellipsoid is equal to M = 4π/3abcρ projection of the obliquity of Vesta into the ICRF from the orbital axis. mean with ρmean the mean density of the body. We compute the polar moment of inertia of Vesta for two sets of interior models: ρc = 3 3 7100 kg/m , rc = 113 km and ρc = 7800 kg/m , rc = 107 km angles, the projected obliquity ρ and longitude angle λ repre- (Russell et al. 2012). By conserving the mass, we obtain a den- / 3 / 3 sented in the figure. The values of ρ for each set of observa- sity for the mantle of 3148.2 kg m and 3148.4 kg m that is in tion are in very good agreement, whereas the longitude angles the range of mantle density determined in Zuber et al. (2011). / 2 = . present a larger error than the ellipse of uncertainties. The longi- The resulting polar moment of inertia is C MR 0 4168 / 2 = . tude angle uncertainty may be related to the angle of illumination and C MR 0 4161, which are very close. For the present at the surface of Vesta by the Sun. Thus, the uncertainty in the study, we fix the value of the polar moment of inertia equal / 2 = .
Recommended publications
  • Dynamical Adjustments in IAU 2000A Nutation Series Arising from IAU 2006 Precession A
    A&A 604, A92 (2017) Astronomy DOI: 10.1051/0004-6361/201730490 & c ESO 2017 Astrophysics Dynamical adjustments in IAU 2000A nutation series arising from IAU 2006 precession A. Escapa1; 2, J. Getino3, J. M. Ferrándiz2, and T. Baenas2 1 Dept. of Aerospace Engineering, University of León, 24071 León, Spain 2 Dept. of Applied Mathematics, University of Alicante, PO Box 99, 03080 Alicante, Spain e-mail: [email protected] 3 Dept. of Applied Mathematics, University of Valladolid, 47011 Valladolid, Spain Received 20 January 2017 / Accepted 23 May 2017 ABSTRACT The adoption of International Astronomical Union (IAU) 2006 precession model, IAU 2006 precession, requires IAU 2000A nutation to be adjusted to ensure compatibility between both theories. This consists of adding small terms to some nutation amplitudes relevant at the microarcsecond level. Those contributions were derived in previously published articles and are incorporated into current astronomical standards. They are due to the estimation process of nutation amplitudes by Very Long Baseline Interferometry (VLBI) and to the changes induced by the J2 rate present in the precession theory. We focus on the second kind of those adjustments, and develop a simple model of the Earth nutation capable of determining all the changes arising in the theoretical construction of the nutation series in a dynamical consistent way. This entails the consideration of three main classes of effects: the J2 rate, the orbital coefficients rate, and the variations induced by the update of some IAU 2006 precession quantities. With this aim, we construct a first order model for the nutations of the angular momentum axis of the non-rigid Earth.
    [Show full text]
  • The Spin, the Nutation and the Precession of the Earth's Axis Revisited
    The spin, the nutation and the precession of the Earth’s axis revisited The spin, the nutation and the precession of the Earth’s axis revisited: A (numerical) mechanics perspective W. H. M¨uller [email protected] Abstract Mechanical models describing the motion of the Earth’s axis, i.e., its spin, nutation and its precession, have been presented for more than 400 years. Newton himself treated the problem of the precession of the Earth, a.k.a. the precession of the equinoxes, in Liber III, Propositio XXXIX of his Principia [1]. He decomposes the duration of the full precession into a part due to the Sun and another part due to the Moon, to predict a total duration of 26,918 years. This agrees fairly well with the experimentally observed value. However, Newton does not really provide a concise rational derivation of his result. This task was left to Chandrasekhar in Chapter 26 of his annotations to Newton’s book [2] starting from Euler’s equations for the gyroscope and calculating the torques due to the Sun and to the Moon on a tilted spheroidal Earth. These differential equations can be solved approximately in an analytic fashion, yielding Newton’s result. However, they can also be treated numerically by using a Runge-Kutta approach allowing for a study of their general non-linear behavior. This paper will show how and explore the intricacies of the numerical solution. When solving the Euler equations for the aforementioned case numerically it shows that besides the precessional movement of the Earth’s axis there is also a nu- tation present.
    [Show full text]
  • Research on the Precession of the Equinoxes and on the Nutation of the Earth’S Axis∗
    Research on the Precession of the Equinoxes and on the Nutation of the Earth’s Axis∗ Leonhard Euler† Lemma 1 1. Supposing the earth AEBF (fig. 1) to be spherical and composed of a homogenous substance, if the mass of the earth is denoted by M and its radius CA = CE = a, the moment of inertia of the earth about an arbitrary 2 axis, which passes through its center, will be = 5 Maa. ∗Leonhard Euler, Recherches sur la pr´ecession des equinoxes et sur la nutation de l’axe de la terr,inOpera Omnia, vol. II.30, p. 92-123, originally in M´emoires de l’acad´emie des sciences de Berlin 5 (1749), 1751, p. 289-325. This article is numbered E171 in Enestr¨om’s index of Euler’s work. †Translated by Steven Jones, edited by Robert E. Bradley c 2004 1 Corollary 2. Although the earth may not be spherical, since its figure differs from that of a sphere ever so slightly, we readily understand that its moment of inertia 2 can be nonetheless expressed as 5 Maa. For this expression will not change significantly, whether we let a be its semi-axis or the radius of its equator. Remark 3. Here we should recall that the moment of inertia of an arbitrary body with respect to a given axis about which it revolves is that which results from multiplying each particle of the body by the square of its distance to the axis, and summing all these elementary products. Consequently this sum will give that which we are calling the moment of inertia of the body around this axis.
    [Show full text]
  • Positional Astronomy Coordinate Systems
    Positional Astronomy Observational Astronomy 2019 Part 2 Prof. S.C. Trager Coordinate systems We need to know where the astronomical objects we want to study are located in order to study them! We need a system (well, many systems!) to describe the positions of astronomical objects. The Celestial Sphere First we need the concept of the celestial sphere. It would be nice if we knew the distance to every object we’re interested in — but we don’t. And it’s actually unnecessary in order to observe them! The Celestial Sphere Instead, we assume that all astronomical sources are infinitely far away and live on the surface of a sphere at infinite distance. This is the celestial sphere. If we define a coordinate system on this sphere, we know where to point! Furthermore, stars (and galaxies) move with respect to each other. The motion normal to the line of sight — i.e., on the celestial sphere — is called proper motion (which we’ll return to shortly) Astronomical coordinate systems A bit of terminology: great circle: a circle on the surface of a sphere intercepting a plane that intersects the origin of the sphere i.e., any circle on the surface of a sphere that divides that sphere into two equal hemispheres Horizon coordinates A natural coordinate system for an Earth- bound observer is the “horizon” or “Alt-Az” coordinate system The great circle of the horizon projected on the celestial sphere is the equator of this system. Horizon coordinates Altitude (or elevation) is the angle from the horizon up to our object — the zenith, the point directly above the observer, is at +90º Horizon coordinates We need another coordinate: define a great circle perpendicular to the equator (horizon) passing through the zenith and, for convenience, due north This line of constant longitude is called a meridian Horizon coordinates The azimuth is the angle measured along the horizon from north towards east to the great circle that intercepts our object (star) and the zenith.
    [Show full text]
  • On the IAU 2000 Nutation Consistency with IAU 2006 Precession Draft Note
    On the IAU 2000 nutation consistency with IAU 2006 precession Draft note∗ Alberto Escapay and Nicole Capitainez Abstract The Earth precession-nutation model of the International Astronomical Union (IAU) is composed of the IAU 2006 precession and IAU 2000 nu- tation. The IAU 2006 precession, which is consistent with both dynamical theory and the IAU 2000 nutation, was adopted to replace the precession part of the IAU 2000 precession-nutation. In that process, it was noticed that very slight adjustments were required to the IAU 2000 nutation ampli- tudes in order to ensure consistency at the micro arcsecond level with the IAU 2006 precession. The formulae for these adjustments provided by Cap- itaine et al. (2005) were implanted in part of the standards and software for computing nutation (e.g., SOFA and IERS Conventions 2010). However, IAU has not adopted any resolution on that direction, so formally such an inconsistency remains. Moreover, it has been shown recently (Escapa et al. 2017) that a few additional terms should be added to the 2005 expressions. We examine the drawbacks that have arisen as a consequence of the lack of such a resolution and propose different options to address them within the framework of IAU 2006 precession and IAU 2000 nutation. They run from just supplementing current IAU resolutions to clarify the content and termi- nology of the IAU precession-nutation model, to adopting a potential new resolution that would also ensure dynamical consistency between precession and nutation. ∗This draft will be submitted to publication yDepartment of Aerospace Engineering, University of Leon,´ E-24071 Leon,´ Spain; [email protected] zSYRTE, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universites,´ F-75014 Paris, France; [email protected] 1 1 Introduction Nowadays, the International Astronomical Union (IAU) precession-nutation model is based on IAU 2000A nutation (Mathews et al.
    [Show full text]
  • SOFA Tools for Earth Attitude
    International Astronomical Union Standards Of Fundamental Astronomy SOFA Tools for Earth Attitude Software version 18 Document revision 1.64 Version for C programming language http://www.iausofa.org 2021 April 18 MEMBERS OF THE IAU SOFA BOARD (2021) John Bangert United States Naval Observatory (retired) Steven Bell Her Majesty’s Nautical Almanac Office Nicole Capitaine Paris Observatory Maria Davis United States Naval Observatory (IERS) Micka¨el Gastineau Paris Observatory, IMCCE Catherine Hohenkerk Her Majesty’s Nautical Almanac Office (chair, retired) Li Jinling Shanghai Astronomical Observatory Zinovy Malkin Pulkovo Observatory, St Petersburg Jeffrey Percival University of Wisconsin Wendy Puatua United States Naval Observatory Scott Ransom National Radio Astronomy Observatory Nick Stamatakos United States Naval Observatory Patrick Wallace RAL Space (retired) Toni Wilmot Her Majesty’s Nautical Almanac Office (trainee) Past Members Wim Brouw University of Groningen Mark Calabretta Australia Telescope National Facility William Folkner Jet Propulsion Laboratory Anne-Marie Gontier Paris Observatory George Hobbs Australia Telescope National Facility George Kaplan United States Naval Observatory Brian Luzum United States Naval Observatory Dennis McCarthy United States Naval Observatory Skip Newhall Jet Propulsion Laboratory Jin Wen-Jing Shanghai Observatory © Copyright 2013-20 International Astronomical Union. All Rights Reserved. Reproduction, adaptation, or translation without prior written permission is prohibited, except as al- lowed under the copyright laws. CONTENTS iii Contents 1 INTRODUCTION 1 1.1 The SOFA software ................................... 1 1.2 Quick start ....................................... 1 1.3 Abbreviations ...................................... 1 2 CELESTIAL COORDINATES 3 2.1 Stellar directions .................................... 3 2.2 Precession-nutation ................................... 3 2.3 Evolution of celestial reference systems ........................ 4 2.4 The IAU 2000 changes ................................
    [Show full text]
  • 3D Rigid Body Dynamics: Tops and Gyroscopes
    J. Peraire, S. Widnall 16.07 Dynamics Fall 2008 Version 2.0 Lecture L30 - 3D Rigid Body Dynamics: Tops and Gyroscopes 3D Rigid Body Dynamics: Euler Equations in Euler Angles In lecture 29, we introduced the Euler angles as a framework for formulating and solving the equations for conservation of angular momentum. We applied this framework to the free-body motion of a symmetrical body whose angular momentum vector was not aligned with a principal axis. The angular moment was however constant. We now apply Euler angles and Euler’s equations to a slightly more general case, a top or gyroscope in the presence of gravity. We consider a top rotating about a fixed point O on a flat plane in the presence of gravity. Unlike our previous example of free-body motion, the angular momentum vector is not aligned with the Z axis, but precesses about the Z axis due to the applied moment. Whether we take the origin at the center of mass G or the fixed point O, the applied moment about the x axis is Mx = MgzGsinθ, where zG is the distance to the center of mass.. Initially, we shall not assume steady motion, but will develop Euler’s equations in the Euler angle variables ψ (spin), φ (precession) and θ (nutation). 1 Referring to the figure showing the Euler angles, and referring to our study of free-body motion, we have the following relationships between the angular velocities along the x, y, z axes and the time rate of change of the Euler angles. The angular velocity vectors for θ˙, φ˙ and ψ˙ are shown in the figure.
    [Show full text]
  • The Precession and Nutation of Deformable Bodies
    The Precession and Nutation of Deformable Bodies Zdengk Kopal - -. ..I > / IPAQES) (NASA CR OR TUX OR AD NUMBER) MATH E MATI C S R E S E A R C H DECEMBER 1966 D1-82-0590 THE PRECESSION AND NUTATION OF DEFORMABLE BODIES by Zdengk Kopal I This research was supported in part by the National Aeronautics and Space Administration under Contract No. NASW-1470. Mathematical Note No. 496 Mathematics Research Laboratory E BOEING SCIENTIFIC RESEARCH LABORATORIES December 1966 ABSTRACT The aim of the present communication has been to set up the Eulerian system of equations which governs the motion of a self-gravitating de- formable body (regarded as a compressible fluid of arbitrarily high viscosity) about its own center of gravity in an arbitrary external field of force. If the latter were particularized to represent the tidal attraction of the Sun and the Moon, this motion would represent the luni-solar precession and nutation of a fluid Earth; if, on the other hand, the external field of force were governed by the Earth (or the Sun), the motion would define the physical librations of the Moon regarded as a deformable body. All these (and other) cases arising in the solar system will be treated in due course. The specific aim of this first of a series of reports in which these problems will be discussed will be to establish the explicit form of the system of differential equations which are basic to our problem. One specific aspect of their solution--namely, dynamical tides on deformable bodies and the consequent dissipation of energy--will be deferred to a second report of this series; while reports 111 and IV will be concerned with particular cases of the precession and librations of the Earth and the Moon.
    [Show full text]
  • Capitaine N., Soffel M.: on the Definition and Use of the Ecliptic In
    ON THE DEFINITION AND USE OF THE ECLIPTIC IN MODERN ASTRONOMY N. CAPITAINE1, M. SOFFEL2 1 SYRTE, Observatoire de Paris, CNRS, UPMC, 61, av. de l'Observatoire, 75014 Paris, France e-mail: [email protected] 2 Lohrmann Observatory, Dresden Technical University, 01062 Dresden, Germany e-mail: michael.so®[email protected] ABSTRACT. The ecliptic was a fundamental reference plane for astronomy from antiquity to the realization and use of the FK5 reference system. The situation has changed considerably with the adoption of the International Celestial Reference system (ICRS) by the IAU in 1998 and the IAU resolutions on reference systems that were adopted from 2000 to 2009. First, the ICRS has the property of being independent of epoch, ecliptic or equator. Second, the IAU 2000 resolutions, which speci¯ed the systems of space-time coordinates within the framework of General Relativity, for the solar system (the Barycentric Celestial Reference System, BCRS) and the Earth (the Geocentric Celestial Reference System, GCRS), did not refer to any ecliptic and did not provide a de¯nition of a GCRS ecliptic. These resolutions also provided the de¯nition of the pole of the nominal rotation axis (the Celestial intermediate pole, CIP) and of new origins on the equator (the Celestial and Terrestrial intermediate origins, CIO and TIO), which do not require the use of an ecliptic. Moreover, the models and standards adopted by the IAU 2006 and IAU 2009 resolutions are largely referred to the ICRS, BCRS, GCRS as well as to the new pole and origins. Therefore, the ecliptic has lost much of its importance.
    [Show full text]
  • Introduction to Astronomy ! AST0111-3 (Astronomía) ! ! ! ! ! ! ! ! ! ! ! ! Semester 2014B Prof
    Introduction to Astronomy ! AST0111-3 (Astronomía) ! ! ! ! ! ! ! ! ! ! ! ! Semester 2014B Prof. Thomas H. Puzia Theme Our Sky 1. Celestial Sphere 2. Diurnal Movement 3. Annual Movement 4. Lunar Movement 5. The Seasons 6. Eclipses Precession • The axis of the Earth (and of the celestial sphere) precesses like the axis of rotation of a top (trompo). The period of this movement is very long: P = 26000 yr. • The Earth is not perfectly spherical, it is wider at the Equator by ~43 km. The inclination of the rotational axis and the asymmetric gravitational pull of the Sun and the Moon produce a change (“torque”) in the axis direction. • Thus the Earth's equatorial plane changes position gradually. This precession of the equinoxes was discovered by Hipparchus in the second century BC. • Example: Today, the spring equinox is in Pisces, but in the year 2600 it will be in Aquarius. • Nutation: a small oscillatory motion superimposed on the Earth's axis precession. The Nutation of Equinoxes Nutation (discovered by James Bradley in 1748) is generally described as the sum of higher-order terms of Earth’s polar motion due to some time-variable nature of !tidal forces that act on Earth’s body (“Precise Geoid”). Nutation is generally: split into vector terms parallel and perpendicular to the direction of precession split into short- and long-period terms due to various effects, such as time-dependent distances of moon, sun, jupiter et al., variable tilt of orbits e.g. moon vs. earth orbit, ocean currents, location of Earth crust relative to her NiFe core, etc. largest nutation component (17x9 arcsec) has a period 6798 days or 18.6 years, while the second- largest (1.3x0.6 arcsec) has a period of 183 days.
    [Show full text]
  • The Specification of Nutation in the Iau System of Astronomical Constants
    THE SPECIFICATION OF NUTATION IN THE IAU SYSTEM OF ASTRONOMICAL CONSTANTS G. A. Wilkins Royal Greenwich Observatory SUMMARY The principal purpose of the IAU system of astronomical constants is to provide a self-consistent set of constants for use in the computation of the international ephemerides of the Sun9 Moon, planets and stars and in the reduction of observations of these bodies. At present nutation is computed from a theory of the rotation of the Earth as a rigid body and only the coefficient of the principal term in obliquity is specified in the system of constants. Such a simple specification will not be ade­ quate for use with the more precise observations that are becoming available, and it appears that it will be necessary to adopt a new model of the Earth and to develop a new theory of nutation which will take in­ to account the elastic properties of the Earth. The new model should be consistent with other constants of the IAU system, and with the model used in other branches of geophysics. The new specification of nutation should be formally adopted by the IAU in 1979 so that it can be used in the published ephemerides for 1984 onwards. INTRODUCTION This paper is intended to provide a general introduction to IAU Sympos­ ium No 78 on "Nutation and the Earth1s Rotation". One of the purposes of this symposium is to provide for the discussion of the problems in­ volved in the development and adoption for international use of a new theory of the nutation of the Earth1s axis of rotation under the action of perturbing forces.
    [Show full text]
  • SOFA Tools for Earth Attitude
    International Astronomical Union Standards Of Fundamental Astronomy SOFA Tools for Earth Attitude Software version 13 Document revision 1.6 Version for Fortran programming language http://www.iausofa.org 2017 July 13 MEMBERS OF THE IAU SOFA BOARD (2007) John Bangert United States Naval Observatory Wim Brouw University of Groningen Mark Calabretta Australia Telescope National Facility Anne-Marie Gontier Paris Observatory Catherine Hohenkerk Her Majesty's Nautical Almanac Office Wen-Jing Jin Shanghai Observatory Zinovy Malkin Inst. Appl. Astron., St Petersburg Dennis McCarthy United States Naval Observatory Jeffrey Percival University of Wisconsin Patrick Wallace Rutherford Appleton Laboratory ⃝c Copyright 2007-16 International Astronomical Union. All Rights Reserved. Reproduction, adaptation, or translation without prior written permission is prohibited, except as al- lowed under the copyright laws. CONTENTS iii Contents 1 INTRODUCTION 1 1.1 The SOFA software ................................... 1 1.2 Quick start ....................................... 1 1.3 Abbreviations ...................................... 1 2 CELESTIAL COORDINATES 3 2.1 Stellar directions .................................... 3 2.2 Precession-nutation ................................... 3 2.3 Evolution of celestial reference systems ........................ 4 2.4 The IAU 2000 changes ................................. 7 2.5 Frame bias ....................................... 7 2.6 CIO and TIO ...................................... 8 2.7 Equation of the origins ................................
    [Show full text]